Daniel, > If you know that your FiniteElement is interpolating (which I assume given > the way you compute the maximum) and that J_tor= x*solution+solution_max/y > is in your ansatz space (probably not) or an interpolation is sufficient, > you can use something like > > J_tor(local_dof_indices[i]) = > right_hand_side.value(fe_values.quadrature_point(q_index), > sol_tmp[q_index],*max) > > and hence just loop once over all the cells instead of solving a linear > system. Of course, you have to decide if the projection is that expensive > that you want to go with an interpolation instead. > > please excuse my ignorance, but I can't understand this part...
But, your conclusion is I don't need to use VectorTools::project because it does similar thing I did.. is it right..? Thank you. Kyusik. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
