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
[email protected]
https://mail.gna.org/listinfo/getfem-users
