In attempting to modify Step 39 to do DG for a different type of equation (with a different error estimator), I have been computing the DG L^2 projection of a function I have already created along the lines of:
VectorTools::project (mapping, dof_handler, constraints, QGauss<dim>(dof_handler.get_fe().tensor_degree()+1), functionb<dim>(), b_h); under Step12<dim>::estimate(). This works fine but the problem arises when I try to somehow transfer this vector to Estimator<dim>::cell. I tried using a similar framework to that used to transferring the solution vector but it kept throwing up errors. Being new at C++ I'm somewhat unfamiliar with how to do it (despite trying lots of different ways) and any help on what code should be used at both ends is appreciated. Thanks!
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
