// ---------------------------------------------------------------------------- double AdaptiveLinearVariationalSolver:: evaluate_goal(Form& M, std::shared_ptr<const Function> u) const { dolfin_assert(M.num_coefficients() > 0); M.set_coefficient(M.num_coefficients() - 1, u); return assemble(M); }