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

[Tetsuo Koyama]

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.