[Getfem-users] estimate diameter of elements

[s.brambilla93]

Hi all,

 

i need to compute the Peclèt number for a transport problem, and I have some issues computing h, the maximum diameter between all the the elements of the mesh.

 

This iwhat i’ve done so far:

 

scalar_type h=0;

scalar_type temp=0;

for (getfem::mr_visitor i(mymesh.region( j )); !i.finished(); i++){

                temp=2*mymesh.convex_radius_estimate(i.cv);

                if(temp>h) h=temp;  

}

 

Mi principal issue is that i need another cycle, over the j-th  region, but i don’t know how to do it: i have created only regions for the faces of the boundary, i have not a region for the inner elements.

Is it possible to use an iterator for all the convexes in a mesh (and not for all the convexes in a region)?

 

Thank you for your help

Best regards

Stefano Brambilla