[deal.ii] Mesh refinement + Periodic boundary conditions

[Jorge R.G.]

Dear all,

I'm trying to solve a diffusion-like problem with periodic boundary conditions. 
I know how to enforce the periodic boundary conditions by adding constraints to 
a ConstraintMatrix object. I also know how to do an adaptive mesh refinement 
following the tutorials. However, I have problems when I try to use both 
periodic conditions and mesh refinement. Sometimes, the resulting mesh is not 
smooth across the boundary (the refinement level of two adjacent cells across 
the periodic boundary is not the same), and I have difficulties building the 
ConstraintMatrix that represents the periodic conditions.

How can I enforce a smooth refinement of the mesh accross the periodic 
boundary? How can I build hanging_node_constraints across this boundary?

I'm a little lost here. I guess I can extend the code in the functions

bool Triangulation< dim >::prepare_coarsening_and_refinement 

and

DoFTools::make_hanging_node_constraints

Has anybody done this before? 

Thanks,

Jorge



      ¡Todo sobre la Liga Mexicana de fútbol! Estadisticas, resultados, 
calendario, fotos y más:&lt;
http://espanol.sports.yahoo.com/
_______________________________________________