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, [email protected] 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
> [email protected]
> https://mail.gna.org/listinfo/getfem-users
--
Yves Renard ([email protected]) 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
---------
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
