EvalStrainOnGaussPts
double u(?)
double varespilon(?,3) -> output