You are right about partial_mesh_fem. It didn't work because of some implementation mistake.
I now have another problem: an integration method is mandatory when we build a new brick. But my integration method changes as the mesh changes. So after one iteration, the integration method I have given to the brick is not up-to-date. Then the system cannot be solved and I get the following error message 'error: the mesh_fem and mesh_im live on different meshes!'. Can I update the integration method of a brick? Thank you for your help. Kind regards, Anne-Lise Guilmin. -----Message d'origine----- De : Yves Renard [mailto:[email protected]] Envoyé : vendredi 30 avril 2010 15:31 À : [email protected] Cc : GUILMIN Anne-Lise Objet : Re: [Getfem-users] Open mesh On vendredi 30 avril 2010, GUILMIN Anne-Lise wrote: > Hello, > > Ok then I will go on with my own strategy. > I still have a problem with the partial mesh fem (used for lagrangian > multipliers). I update it the same way I do with mesh fem, but this leads > to a context error in the model. So I am thinking about creating new > variables for multipliers at each change of mesh. Is there any way to > destroy a variable of a model? Or just disable it? Normally, the partial mesh_fem should adapt itself automatically to the new situation. Perhaps, you can send me a piece of code exhibiting the problem. Yves. > > Thank you for your help. > Kind regards, > > Anne-Lise Guilmin. > > -----Message d'origine----- > De : Yves Renard [mailto:[email protected]] > Envoyé : jeudi 29 avril 2010 11:14 > À : [email protected] > Cc : GUILMIN Anne-Lise > Objet : Re: [Getfem-users] Open mesh > > > > Daer Anne-Lise, > > You are right, there is no specific tool for the moment in Getfem++ to deal > with such a situation. This is not a trivial task to produce a generic > algorithm for the automatic re-interpolation (or projection) of a field > after refinement or change of the mesh. So, for the moment, depending on > the situation, a specific strategy have to be developped for each > situation. > > Yves. > > On mercredi 28 avril 2010, GUILMIN Anne-Lise wrote: > > Dear all, > > > > I am working with an open mesh: I add convexes to the mesh before each > > new increment. > > I would like to keep the values of my model variables from one increment > > to another. > > As the mesh changes, the numbering of the dof may change too. Is there > > any predefined function to interpolate the old values on the new > > mesh_fem? > > You may have already met this problem when refining meshes. > > > > For now, I construct a new mesh_fem, then new vectors of variables from > > an interpolation of old ones on the new mesh_fem, and finally I > > overwrite the old meshfem and variables in the model with the new ones. > > This works well with my main variable but not with lagrangian > > multipliers: the computation stops due to an error of context of the > > model after updating the partial_mesh_fem on which are defined the > > multipliers. > > How could I update a partial_mesh_fem and still get a correct context > > for my model? > > > > Thank you for your help. > > > > Anne-Lise Guilmin > > __________________________ > > > > Ce message (et toutes ses pièces jointes éventuelles) est confidentiel et > > établi à l'intention exclusive de ses destinataires. Toute utilisation de > > ce message non conforme à sa destination, toute diffusion ou toute > > publication, totale ou partielle, est interdite, sauf autorisation > > expresse. L'IFP décline toute responsabilité au titre de ce message. > > > > This message and any attachments (the message) are confidential and > > intended solely for the addressees. Any unauthorised use or dissemination > > is prohibited. IFP should not be liable for this message. > > > > Visitez notre site Web / Visit our web site : http://www.ifp.fr > > __________________________ -- 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 --------- __________________________ Ce message (et toutes ses pièces jointes éventuelles) est confidentiel et établi à l'intention exclusive de ses destinataires. Toute utilisation de ce message non conforme à sa destination, toute diffusion ou toute publication, totale ou partielle, est interdite, sauf autorisation expresse. L'IFP décline toute responsabilité au titre de ce message. This message and any attachments (the message) are confidential and intended solely for the addressees. Any unauthorised use or dissemination is prohibited. IFP should not be liable for this message. Visitez notre site Web / Visit our web site : http://www.ifp.fr __________________________ _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
