Hi! I'm not using a mesh.transformation(M), this is a complete test program that reproduces it (copied from the documentation):
************************ 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()); const int numPoints = mymesh.points_index().size(); std::cout << mymesh.points()[numPoints-1] << std::endl; ************************* At the end it prints the coordinates of the last point, which I expected to be [1, 1, 1] since those are the limits set in 'vect', but I get: [10, 10, 10] I compile it with g++ 4.4.4 (no compilation warnings): g++ Test.cpp -o Test.elf -Wall -Wextra -O3 -lgetfem and use GetFEM++ 4.1 in the current Debian Sid OS, and a 64bit computer with a Intel T4400 dual core processor. - Torquil On 24/09/10 09:14, Yves Renard wrote: > > 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) ? > > 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 >> [email protected] >> https://mail.gna.org/listinfo/getfem-users > > _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
