FemSourceScalar
double f(?)
double rhs(?) -> output