On 02/07/2012 02:12 PM, Mohammad Mohsenie wrote:
hi there
in my code vector gradient is defined and calculated as :
Vector<double> gradient ;
for (unsigned int icell = 0; cell!=endc; ++cell, ++icell)
{ .
gradient(icell) = -(0.5*gr_1)- gr_2 ;
What are gr_1 and gr_2?
so now gradient is calculated in any cell . now i want to know how i can
have gradient in every point . i have an idea to at first specify the
cell which point belong to it so we will have gradient there but i dont
know can such procedure be done in dealii or not if yes how ؟
Using FEValues::get_function_gradients allows you to obtain the gradient
of a function at arbitrary quadrature points. This is in analogy to how
step-9 uses the FEValues::get_function_values function.
Best
W.
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
