Sebastian,
vector fields created by polynomial interpolation cannot be made
uniformly of length one by adjusting the values in the support points,
unless they are constant.
What you need to do is: in the function, that uses the vector field, for
instance in quadrature points, correct it in each point.
Best,
Guido
On 9/23/2010 11:44 AM, M. Sebastian Pauletti wrote:
Dear all,
I have the following issue:
Suppose I have a finite element approximate solution to a vector
problem that is stored in
Vector<double> solution;
This solution is supposed to be a unit vector field, because of
approximation issues it is not exactly unit. As I want to feed this as
an input to another process I want to make it unit (up to machine
precision). So I would like to do something like
for ( it=solution.begin() ; it< solution.end(); it++ )
make_unit(*it);
Off course this does NOT work as there is no information in solution
as to what components in solution are the components of the vector
field for each dof.
Question: Is there a simple way to make a vector field solution a unit
vector field?
Thank you.
Best
-Sebastian
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
