Dear Tetsuo, I think that the problem come from the fact that an initialisation bgeot::base_node(0.0) is understand as bgeot::base_node(0) i.e. converting the float in integer and gives a zero size point ... Best regards, Yves Le 28/01/2019 à 11:16, Tetsuo Koyama a écrit :
Dear GetFEM++ Project Thank you for useful mesh build library in GetFEM++. I'm wondering error occured when using add_point for 1D mesh. I checked newest git repository version and version 5.3. I want to confirm that it is bug or not. Please check if it is my mistake of use or specification of library. If it is bug, I'll fix it myself for my study. Code is... getfem::mesh mymesh; std::vector<bgeot::size_type> ind(3); ind[0] = mymesh.add_point(bgeot::base_node(0.0)); ind[1] = mymesh.add_point(bgeot::base_node(1.0)); ind[2] = mymesh.add_point(bgeot::base_node(2.0)); Error is... terminate called after throwing an instance of 'gmm::gmm_error' what(): Error in bgeot_node_tab.cc, line 106 gmm::size_type bgeot::node_tab::add_node(const base_node&, bgeot::scalar_type, bool): Nodes should have the same dimension Aborted (core dumped) Best Regard Tetsuo.
-- Yves Renard ([email protected]) tel : (33) 04.72.43.87.08 INSA-Lyon 20, rue Albert Einstein 69621 Villeurbanne Cedex, FRANCE http://math.univ-lyon1.fr/~renard ---------
