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 > > [email protected] > > https://mail.gna.org/listinfo/getfem-users -- Yves Renard ([email protected]) 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 --------- _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
