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
>
>
> First of all, Thank you very much for your reply
>
> This can only work if you know that the maximum is attained at one of the
> support_points
> of your FiniteElement. In general, you can use
> VectorTools::integrate_difference to find a good approximation
> to the maximum value of
Kyusik,
//by Han(To get psi0)
> Vector::iterator max;
> max=std::max_element(solution.begin(),solution.end());
>
This can only work if you know that the maximum is attained at one of the
support_points
of your FiniteElement. In general, you can use
VectorTools::integrate_difference to
