Re: [Getfem-users] Build a mesh with quadrature points

[julien pommier]

Hi Giovani,
If you want to retrieve the coordinates of the quadrature points of a 
convex cv, you can do:

   /* get the approximate integration method for the convex */
   getfem::papprox_integration pai
     = mim.int_method_of_element(cv)->approx_method();
   /* prepare the geotrans_interpolation_context to apply the geometric 
transformation */
   bgeot::base_matrix G;
   bgeot::vectors_to_base_matrix(G, m.points_of_convex(cv));
   bgeot::geotrans_interpolation_context c(m.trans_of_convex(cv),
                       pai->point(0), G);
   /* apply the geotrans to each point of the integration method */
   for (size_type j = 0; j < pai->nb_points_on_convex(); ++j) {
     c.set_xref(pai->point(j));
     std::cout << " integration node " << j << " for convex " << cv << 
" is " << c.xreal() << "\n";
   }

Best regards,
Julien

Giovani wrote:
> Good day  gentlemen!
>
> Suppose I have a getfem::mesh object that has a getfem::mesh_fem 
> linked to it, and also a getfem::mesh_im that's also linked.
>
> I would like to build a new getfem::mesh, whose nodes are located at 
> the quadrature points of the previous mesh. Is there any simple way of 
> retrieving a list with the global coordinates of those quadrature points?
>
> Best Regards,
>
> Giovani M. Faccin
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users


--
Julien Pommier, bureau 111
GMM, INSA Toulouse, tél:05 61 55 93 42