Re: [Getfem-users] interface condition question

[Yves Renard]

Hello, 

For such a simple situation (matching meshes), you can also use the 
set_dof_partition method of the mesh_fem object. See line 678 of 
tests/crack.cc to define the partition on a domain and line 833 to obtain the 
constraint matrix (you can preferably use the gmm::range_basis algorithm to 
select the constraints in the whole mass matrix instead of simply looking at 
the diagonal terms).

For non-matching meshes, you have to use the (imperfect) interpolated fem 
object.

Best regards,

Yves.




On mercredi 1 juillet 2009, Ronan Perrussel wrote:
> Hello,
>
> it was not fully straightforward but you can find an example in:
> getfem++/contrib/inter_element_test,
> if you have matching meshes on the interface. If not it is certainly
> possible but I am not the expert to explain how to do.
>
> I hope it helps,
> Best regards,
> Ronan
>
> Jehanzeb Hameed a écrit :
> > Hello,
> >
> > I am investigating the possibility of using getfem for implementing a
> > 2 domain problem, with an interface between the domains. I want to use
> > continuous lagrange elements on each domain, and a jump condition at
> > the interface. The jump condition will take the form (u_1 - u_2, v_1 -
> > v_2), where u_1 and v_1 and test and trial functions for domain 1, and
> > u_2 and v_2 are test and trial functions for domain 2.  The integral
> > (u_1 - u_2, v_1 - v_2) is only over the interface between the domains.
> >
> > Can this be done relatively easily in getfem?
> >
> > Thanks,
> > -Jehanzeb
> >
> > _______________________________________________
> > Getfem-users mailing list
> > address@hidden
> > https://mail.gna.org/listinfo/getfem-users



-- 

  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

---------