[deal.II] How to clear data of triangulation completely

[Felix Lorenz]

Hello everybody,
I have the following problem:

Class C{
  ...

  Triangulation <1,2>boundary_triangulation; 

  FE_Q<1,2>  boundary_fe;

 DoFHandler <1,2>   boundary_dof_handler;
  MappingQGeneric <1,2>  boundary_map;
  .../
  void run ();
}

void C::run(){

 for (unsigned int cycle=0; cycle< ncycles; ++cycle)  //works for 
ncycles=1, error for ncycles > 1

{

   boundary_triangulation.create_triangulation (vertices, cells, 
SubCellData());// 
for every cycle I have here other vertices and cells to create 

// a new triangulation

  ... //here I do some finite element calculation, it works 
   ...

   boundary_dof_handler.clear();
   boundary_triangulation.clear();// Reset this triangulation into a virgin 
state by deleting all data. This operation is only allowed if no 
  //subscriptions to this object exist any 
more, such as DoFHandler objects using it
}
}

After the call of boundary_triangulation.clear() lots of vertices and cells 
are still stored by boundary_triangulation. Resulting to this I get the 
error in the function create_triangulation:

"You have a vertex in your triangulation " "at which more than two cells 
come together. "
 
How can I delete the data of the triangulation completely? Why does the 
clear() function not delete all data?

Best regards 
Felix



-- 
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.

Re: [deal.II] Compute jacobian at point (no quadrature point)

[Felix Lorenz]

Thank you very much. This is exactly what I wanted.

Best regards,
Felix

-- 
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.

[deal.II] Compute jacobian at point (no quadrature point)

[Felix Lorenz]

Hi everyone,

I want to compute the transformation of a gradient at a point (no 
quadrature point) on the reference cell to the real space. Mathematically 
formulated I want to compute:





How do I get the inverse Jacobian at this point?


Best regards,

Felix

-- 
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.

[deal.II] Compute jacobian at a point (no quadrature point)

[Felix Lorenz]

Hi everyone,

I want to compute the transformation of the gradient of the shape function 
at a point  (no quadrature point) on the reference cell to the real space. 
Mathematically written I want to compute:





   

: point in the reference cell

: jacobian

: gradient of the i-th shape function 


How can I get this  ?


Best regards,

Felix

-- 
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.