Re: [Getfem-users] Fwd: Re: combining 2D and 1D elements

[Yves Renard]

On mardi 27 avril 2010, Till Heinemann wrote:
> Notes:
>
> // First associate FEM mf_bubble (among others) with my 2D mesh
> // MESH_FILE='structured:GT="GT_LINEAR_QK(2)";SIZES=[1,1];NOISED=0';
> // N=2; NX=1;
> // BUBBLE_FEM_TYPE = 'FEM_PK(1, 1)'; % for 1D elements
> // FEM_TYPE = 'FEM_QK(2,1)'; % for 2D elements
> elastostatic_problem(void) : mim(mesh),mf_u(mesh),
>              mf_bubble(mesh),mim_bubble(mesh) {}
>
> // ... then create convexes from faces of convex 0
> mesh.add_faces_of_convex(0);
> // ... ! cannot call mf_u.nb_dof() if updating is enabled
> getfem::pfem pf_disable_update=0;
> mf_u.set_finite_element(mf_u.linked_mesh().convex_index(N),pf_u);
> mf_u.set_auto_add(pf_disable_update);
> mf_bubble.set_finite_element(mf_bubble.linked_mesh().convex_index(N-1),pf_b
>ubble); mf_bubble.set_auto_add(pf_disable_update);
> mf_u.nb_dof(); // works (unless I update it)
> mf_bubble.nb_dof(); // DOES NOT WORK.
>
> console output:
> terminate called after throwing an instance of 'gmm::gmm_error'
>   what():  Error in ../src/getfem/getfem_mesh.h, line 222
> bgeot::pgeometric_trans getfem::mesh::trans_of_convex(bgeot::size_type)
> const:
> No geometric transformation or nonexisting element
>
> What am I doing wrong?

The function  mesh.add_faces_of_convex(0); should not be used in this context. 
This is a low level function of the object bgeot::mesh which should not be 
used for a getfem::mesh because it does not define the geometric 
transformation. May be this function could be re-written for getfem::mesh but 
specifying the geometric transformation.

Yves.



-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  Pole de Mathematiques, INSA-Lyon             fax : (33) 04.72.43.85.29
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------