Ok, I got it. I had to assign the whole tetrahedra to the region. So there is no bug.
For those who are interested in details: http://dl.dropbox.com/u/2465545/mesh_refine_tetrahedra2.m and the diff: http://dl.dropbox.com/u/2465545/mesh_refine_tetrahedra2.m.diff Thanks, Mirko 2010/8/27 Renard Yves <[email protected]>: > > > > Dear Mirko, > > The algorithm of Bank refinement implemented in getfemç_mesh.cc is normally > designed to keep > the region number when refining elements (see the function > mesh::handle_region_refinement). > So, there is probably a bug somewhere ! > I you have a small example (matlab or c++) exhibiting the error, I could > have a look to > this function to see what is wrong. > > Yves. > > > > Mirko Windhoff <[email protected]> a écrit : > >> Dear GetFEM devs and users, >> >> I have a tetrahedral mesh with several region numbers, which I want to >> refine using the matlab interface function: >> gf_mesh_set(mesh M, 'refine'[, ivec CVIDs]) >> The refining adds new tetrahedra to the mesh, with regionnumber 0. >> Since the new tetrahedra are spread throughout the whole mesh, I can't >> determine to which region they belong (see screenshot). Then I >> selected only tetrahedra from one region. But unfortunately that >> did'nt help, because new tetrahedra are found also outside of the >> region. >> Is there any possibility to restrict refining to a certain region or >> to assign the correct region number to the refined tetrahedra? >> >> I think it would be possible to add fem data first based on the region >> number and refine afterwards, but I would like to refine the mesh >> first and the use the refined mesh to assign values based on the >> regions and a lookup volume. >> >> My refining function for my custom mesh structs: >> http://dl.dropbox.com/u/2465545/mesh_refine_tetrahedra.m >> Screenshot: >> http://dl.dropbox.com/u/2465545/refined_mesh.png >> >> Thank you for any hints. >> Mirko >> >> _______________________________________________ >> 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
