The following code clamps each cell in a nalgebra square DMatrix. Is there a way to do this in a more “vectorized” fashion, say, using SIMD instructions?

When I run the code with any values for a to e and select integration or differentiation my code panics saying that there is a dimension mismatch and I’m not sure why. For reference here is the code I am using