Re: [Getfem-users] Open mesh

[GUILMIN Anne-Lise]

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?

Thank you for your help.
Kind regards,

Anne-Lise Guilmin. 

-----Message d'origine-----
De : Yves Renard [mailto:address@hidden
Envoyé : jeudi 29 avril 2010 11:14
À : address@hidden
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 (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

---------
__________________________

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
__________________________