Hello Sebastian,
in what norm do you want to make it unit?
Best Regards,
Markus
Am 23.09.10 18:44, schrieb M. Sebastian Pauletti:
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
