Real InitialConditionTempl<RealVectorValue>::dotHelper(const GradientType & op1, const GradientType & op2) { return op1.contract(op2); }