Hello,

I am trying to obtain the vertex global dofs using the DOFAccessor. I 
employ the following code or the same: 

const auto vertex_end = tria.end_vertex();
for (auto vertex = tria.begin_active_vertex(); vertex != vertex_end;
++vertex)
{
DoFAccessor<0, 1, 3, false> vertex_dofs(&tria,
vertex->level(),
vertex->index(),
&dof_handler);
std::cout << vertex_dofs.dof_index(0) << ", " << vertex_dofs.dof_index(1)
<< std::endl;
}

The error that  i obtain is as follows:

    dealii::TriaAccessor<0, 1, spacedim>::TriaAccessor(const 
dealii::Triangulation<1, spacedim>*, int, int, const AccessorData*) [with 
int spacedim = 3; dealii::TriaAccessor<0, 1, spacedim>::AccessorData = void]
The violated condition was: 
    (level == -2) && (index == -2)
Additional information: 
    This exception -- which is used in many places in the library --
    usually indicates that some condition which the author of the code
    thought must be satisfied at a certain point in an algorithm, is not
    fulfilled. An example would be that the first part of an algorithm
    sorts elements of an array in ascending order, and a second part of
    the algorithm later encounters an element that is not larger than the
    previous one.
    
    There is usually not very much you can do if you encounter such an
    exception since it indicates an error in deal.II, not in your own
    program. Try to come up with the smallest possible program that still
    demonstrates the error and contact the deal.II mailing lists with it
    to obtain help.

Can someone help me with this? The main objective is to obtain the global 
dof indices using the - dealii::TriaActiveIterator<dealii::TriaAccessor<0, 1, 
3>>

Thanks
Vinayak

-- 
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].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/d13673eb-ac1e-477a-8c63-9186daaa8d30n%40googlegroups.com.

Reply via email to