Re: [Getfem-users] Structured mesh coordinate values

[Yves Renard]

Dear Torquil,

No, this really generates a mesh on the cube [0,1]^3. Are you sure you do not 
have an instruction

mesh.transformation(M);

after the mesh generation  (like in tests/laplacian.cc) ?

Another simple way to generate a structured mesh is to use getfem::import_mesh 
like this :


getfem::import_mesh("structured:GT_QK(3,1)";SIZES=[1,1,1];NSUBDIV=[10,10,10]", 
mesh);

This allows to use a the regular mesh from a parameter file instead of a 
filename.

Yves.


On jeudi 23 septembre 2010, Torquil Macdonald Sørensen wrote:
> Hi!
> 
> In the "Short User Documentation" it is written that to generate a 10x10x10
> mesh on the unit cube [0,1]^3, the following code would work:
> 
> getfem::mesh mymesh;
> bgeot::base_node org(0.0, 0.0, 0.0);
> std::vector<bgeot::base_small_vector> vect(3);
> vect[0] = bgeot::base_small_vector(1.0, 0.0, 0.0);
> vect[1] = bgeot::base_small_vector(0.0, 1.0, 0.0);
> vect[2] = bgeot::base_small_vector(0.0, 0.0, 1.0);
> std::vector<int> ref(3);
> ref[0] = ref[1] = ref[2] = 10;
> getfem::parallelepiped_regular_simplex_mesh(mymesh, 3, org, vect.begin(),
> ref.begin());
> 
> But when I print the coordinates of the points in this mesh, i.e.
> mymesh.points(), I get component values not within [0,1], but in [0,10].
> 
> Is this a mistake in the documentation, or do these values not represent
> the coordinates of the points (but instead "scaled coordinates")?
> 
> Thanks!
> Torquil
> 
> _______________________________________________
> 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

---------