<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi all,<br>

<br>

suppose you have a P1 Lagrange vector-valued function, defined over a certain grid in R^3.

<br>

Given a solution of my problem represented by such a function, I need to know the value taken by the solution (which is naturally a 3-component vector) in each vertex and then calculate the lenght of such a vector.<br>

I need specifically this vertex-wise information, for experimental convergence purposes.<br>

<br>

How can I extract such a 3-component vector?<br>

<br>

<div>

<div style="font-family:Tahoma; font-size:13px">The instance of the solution is the following:<br>

<br>

DiscretePositionFunctionType discrete_position_solution ( "discrete_position_solution", dfpSpace );<br>

<br>

where the type is:<br>

<br>

typedef Dune::AdaptiveDiscreteFunction < DiscretePositionFunctionSpaceType > DiscretePositionFunctionType;<br>

<br>

and dfpSpace is an instance of the corresponding space type:<br>

<br>

typedef Dune::LagrangeDiscreteFunctionSpace< VectorPositionFunctionSpaceType, GridPartType, vectorPolynomialOrder > DiscretePositionFunctionSpaceType;<br>

<br>

<br>

Thanks a lot !!!<br>

<br>

Andrea<br>

__________________________________________________________<br>

<br>

Andrea Sacconi<br>

PhD student, Applied Mathematics<br>

AMMP Section, Department of Mathematics, Imperial College London,<br>

London SW7 2AZ, UK<br>

a.sacconi11@imperial.ac.uk<br>

</div>

</div>

</div>

</body>

</html>