[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] interface condition question
From: |
Yves Renard |
Subject: |
Re: [Getfem-users] interface condition question |
Date: |
Wed, 1 Jul 2009 10:53:57 +0200 |
User-agent: |
KMail/1.9.9 |
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
---------