Linear algebra performance issue

[baggepinnen]

You do not need to update sigma in place when you use static arrays, simply return the new one

function stress(epsilon::SMatrix{3,3,Float64})
  (epsilon[1,1]+epsilon[2,2]+epsilon[3,3]) * UniformScaling(1.) 
end