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