Hello Everyone:

Is there a preferred way to obtain the actual coordinates of the point at
which the solution (and its derivatives and normals) are provided in the
DataPostprocessor class?  For example, if one has a material with varying
elastic moduli, computing the stress may involve evaluating a function to
get the moduli at that specific point.  If it isn't implemented already, how
difficult would it be to take this:

compute_derived_quantities_scalar (const std::vector< double
>                 & /*uh*/,
                                const std::vector< Tensor< 1, dim > >
& duh,
                                const std::vector< Tensor< 2, dim > >
& /* dduh */,
                                const std::vector< Point< dim > >
& /* normals */,
                                std::vector< Vector< double > >
& computed_quantities) const


and make it into something like this:


compute_derived_quantities_scalar (const std::vector< double
>                 & /*uh*/,
                                const std::vector< Tensor< 1, dim > >
& duh,
                                const std::vector< Tensor< 2, dim > >
& /* dduh */,
                                const std::vector< Point< dim > >
& /* normals */,
                                const std::vector< Point<dim> >
& support_points     <-------CHANGE HERE
                                std::vector< Vector< double > >
& computed_quantities) const


Thanks as always,

Jon Pitt
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to