Re: [Getfem-users] Error occuer when using add_point in 1D mesh

[Tetsuo Koyama]

Dear Yves

Thank you.

Best Regard Tetsuo.

> 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

2019年1月29日(火) 0:46 Tetsuo Koyama <address@hidden>:
>
> Dear GetFEM++ Project
>
> Hello. Sorry I could solved the problem by myself.
>
> getfem::mesh mymesh;
> std::vector<bgeot::size_type> ind(3);
> bgeot::base_node P1{1};
> bgeot::base_node P2{1};
> bgeot::base_node P3{1};
>
> P1[0] = 0.0;
> P2[0] = 1.0;
> P3[0] = 2.0;
>
> ind[0] = mymesh.add_point(P1);
> ind[1] = mymesh.add_point(P2);
> ind[2] = mymesh.add_point(P3);
>
> Thank you.
> Best Regard Tetsuo.
>
> 2019年1月29日(火) 0:25 Tetsuo Koyama <address@hidden>:
> >
> > Dear GetFEM++ Project
> >
> > Hello.
> > I understanded that bgeot::base_node(N) makes the N-dim base_node.
> > But I'm still thinking about how to make 1D base_node.
> >
> > Best Regard Tetsuo.
> >
> > 2019年1月28日(月) 19:16 Tetsuo Koyama <address@hidden>:
> > >
> > > 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.