Re: [Getfem-users] estimate diameter of elements

[Yves Renard]

Dear Stefano,

Yes, of course, it is possible to iterate on the the whole set of element of a mesh. It reads

for (getfem::dal_bv_visitor i(mymesh.convex_index()); !i.finished(); i++){

    temp=2*mymesh.convex_radius_estimate(i);

    if(temp>h) h=temp;

}

Best regards,

Yves.

Le 15/03/2017 à 15:43, address@hidden a écrit :

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

_______________________________________________
Getfem-users mailing list
address@hidden
https://mail.gna.org/listinfo/getfem-users

-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  Pole de Mathematiques, INSA-Lyon             fax : (33) 04.72.43.85.29
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------