Re: Fwd: Re: [Getfem-users] new_interpolated_fem problem

[Andriy Andreykiv]

Dear Yves,

        I have finaly managed to implement the application or Dirichlet 
boundary 
conditions in a weak sence. What you've suggested below works fine. 

Thank you very much. 

        In the attachment you can see an example picture of a problem where 
Dirichlet 
BC's weere applied along the cirle mesh (made of segments).

Best regards,
                        Andriy

>On Thursday 01 March 2007 18:38, Yves Renard wrote:

> Yes, due to renumeration of the dof, the correspondance is not trivial. It
>  can be obtained by the following code:
>
>      // the previous code which compute the mass matrix
>
>      getfem::mesh_fem mq_interpolate(ml.linked_mesh());
>      getfem::pfem ifem=getfem::new_interpolated_fem(mq,mim);
>      mq_interpolate.set_finite_element(ifem);
>      sparse_matrix M(6,2);
>      getfem::asm_mass_matrix(M,mim,mq_interpolate,ml);
>
>      // the code which compute the correspondance into org_index between
>      // dof of mq_interpolate and those of mq.
>      std::vector<size_type> org_index(mq_interpolate.nb_dof());
>      for (dal::bv_visitor icv(mq_interpolate.convex_index());
>         !icv.finished(); ++icv) {
>        for (size_type j = 0; j < mq_interpolate.nb_dof_of_element(icv);
> ++j) { org_index[mq_interpolate.ind_dof_of_element(icv)[j]]
>          = ifem->index_of_global_dof(icv, j);
>        }
>      }
>
> Can you try this and tell me if it works ?
>
> Yves.
>
> -------------------------------------------------------------------------
>   Yves Renard (address@hidden)       tel : (33) 04.72.43.80.11
>   Pole de Mathematiques, INSA de Lyon          fax : (33) 04.72.43.85.29
>   Institut Camille Jordan - CNRS UMR 5208
>   20, rue Albert Einstein
>   69621 Villeurbanne Cedex, FRANCE
>   http://math.univ-lyon1.fr/~renard
> -------------------------------------------------------------------------
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users

-- 
__________________________________________
Andriy Andreykiv
Delft University of Technology
Faculty of Mechanical Engineering 
Fundamentals of Microsystems
Mekelweg 2
2628 CD  Delft
The Netherlands

E-mail  : address@hidden
Tel.    :   (31) 15 2786818
Fax.    :   (31) 15 2789475
www     :   http://www-tm.wbmt.tudelft.nl/~andrico
private:  (31) 6 47376804

FictitiousDomainMethod.jpg
Description: JPEG image