FemToGaussPointsMfem
double u(?)
double out(?) -> output