Dear Danilo,
The definition of the region 10 is in fact inapropriate. Your definition m.set_region(10,m.faces_from_cvid(m.cvid())) means that you include faces of elements (edges in 2D), not elements themselves. The result is that the elliptic term is indeed integrated on edges, not on elements. And of course it makes a difference. You can define your region as follows domain = m.cvid() m.set_region(10, [domain, domain*0-1]) Indeed, a -1 is needed (0 for the matlab interface) to specify that the element itself is added. This is not very satisfactory. I think, I will modify the set_region method in order to be able to write m.set_region(10, m.cvid()). Thank you for your report, Yves. On vendredi 23 septembre 2011, Danilo Sergi wrote: > Dear Prof. Renard, > > > > I am grateful for your consideration of my problem. I have attached a > program for the heat equation as implemented in Matlab. By considering > only conduction (Laplace term), it is already seen a difference between > the case where the whole mesh is specified explicitly (lines 48 and > 108) and the default setting. > > > > I hope that the remark makes sense and can be useful. > > > > Sincerely yours, > > > > Danilo Sergi -- 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
