Hi All,
I have successfully computed the secondary electric field for 3D case as
described in the paper by Grayver&Buerg, 2014. These field is declared as
BlockVector<double> solution, where the block correspond to real and imaginary
parts.
Now I would like to compute total electric E (secondary plus primary) and
magnetic field (coef* curlE)
1) I can use PostProcessor class to do this
template <int dim>
class Postprocessor : public DataPostprocessor<dim>
{
const SmartPointer<EPrimary<dim> > eprimary;
public:
Postprocessor (const SmartPointer<EPrimary<dim> > eprimary) :
eprimary (eprimary)
{
}
virtual
void
evaluate_vector_field( ...
...
}
But then I need to extract values of the total fields at set of receiver
locations.
If I had total fields in the form similar to solution, I could use
for (unsigned int i = 0; i < receiver_locations.size (); ++i)
{
VectorTools::point_value (dof_handler, solution,
receiver_locations[i],
solution_at_receiver); ...
}
But since total fields are computed with help of postprocessor, I do not know
how to extract values at receivers.
2) The other option is to follow the advice by Ross Kynch, and have loop over
cells, quadrature points and dofs to get local contributions to total electric
and magnetic field.
In this case I am not sure how then I get a global vector.
Please let me know what is the best way.
I would very much appreciate you reply,
Anna
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
