Re: [Libmesh-users] Where to Store the solution gradients?

[Roy S. <roy...@ic...>]

On Mon, 4 Jun 2007, John Peterson wrote:

> Before initializing the EquationSystems object, you'll probably want
> to add an additional vector via the System::add_vector function.
>
> -John
>
>  NumericVector<Number> & add_vector (const std::string& vec_name,
> 				      const bool projections=true);

This will let you store a copy of the solution field, but it gets
stored as degree of freedom coefficients.  LibMesh has no facility for
saving values at quadrature points; if you really have to store them
rather than recomputing them on the fly, you might try a
std::vector<std::vector<Real> >, indexed by element id and quadrature
point index, but that indexing will break after any mesh refinement.
---
Roy

> Bin writes:
> > Hi,
> >
> > After I use libmesh to solve an elasticity problem, I got displacement on
> > each nodes. However, the solution gradients (stress and strain) are also
> > important to me. My question is: after I compute these solution gradients on
> > gauss points in each element, where can I store them? What object in libmesh
> > can I use to store these gradients? I want to use them in the next step.
> >
> > Thank you.
> >
> > --
> > Best regards,
> >
> > Bin
>
> -------------------------------------------------------------------------
> This SF.net email is sponsored by DB2 Express
> Download DB2 Express C - the FREE version of DB2 express and take
> control of your XML. No limits. Just data. Click to get it now.
> http://sourceforge.net/powerbar/db2/
> _______________________________________________
> Libmesh-users mailing list
> Lib...@li...
> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>