Dear all,
I am now using FE_Nedelec to solve the Maxwell's equation. I want to use
the AMS (Auxiliary-space Maxwell Solver provided in HYPRE) to precondition
the system. The solver needs the discrete gradient matrix as well as the
coordinates of the vertices in the mesh. For lowest order nedelec element,
the discrete gradient matrix is represented by the edges of the mesh in
terms of its vertices. So I need the indices of the two vertices of the
edge on which the DoF is defined.
For now, I would use the following method to get those information.
for (auto cell = dh.begin_active(); cell != dh.end(); ++cell) {
for (size_t i = 0; i < dofs_per_cell; ++i) {
auto line = cell->line(i);
size_t v0 = line->vertex_index(0);
size_t v1 = line->vertex_index(1);
}
}
But this method depends on the assumption that the order of DoFs in a cell
are same as the order of lines in a cell. I want to know whether this
assumption is valid. And is there a more general method?
Thank you in advance.
Best regards,
Ce Qin
--
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
