[Getdp] example of fixed point iteration in GetDP

Frederic Trillaud ftrillaudp at gmail.com
Tue Jul 15 10:54:03 CEST 2014


Dear Ruth,

Thank you for your quick reply. I have the oneLab bundle Gmsh-Getdp
installed on my linux box. However, it does not come with examples. Is
the one mentioned can be found on the Getdp website or I should install
Getdp to access it?


Best,


Frederic


On mar, 2014-07-15 at 07:23 +0000, Ruth Vazquez Sabariego wrote:
> Hi Frederic,
> 
> 
> You can inspire yourself from the magnetostatic demo you get when
> downloading GetDP.
> 
> 
> There is a complete Newton-Raphson example, i.e. a term for completing
> the Jacobian matrix is included (JacNL).
> You can just remove or comment it for getting a fixed point iteration.
> The resolution does not change.
> 
> 
> HTH,
> Ruth
> 
> 
> 
> 
> 
> 
> 
> 
> 
> On 14 Jul 2014, at 18:29, Frederic Trillaud Pighi
> <ftrillaudp at pumas.iingen.unam.mx> wrote:
> 
> > Dear all,
> > 
> > 
> > 
> > I have a non-linear magnetodynamics problem that involves a fixed
> > point iteration. Would someone have an example to share or some
> > recommendation how to implement the algorithm?
> > 
> > 
> > 
> > The problem is the following (A formulation, sig: electrical
> > conductivity, Je: applied current):
> > 
> > 
> > 
> > Curl[nu*Curl(A)]+sig*dA/dt = Je
> > 
> > 
> > 
> > with,
> > 
> > 
> > 
> > sig = f(-dA/dt+Je/sig])
> > 
> > 
> > 
> > 
> > Best,
> > 
> > 
> > 
> > 
> > Frederic
> > 
> > _______________________________________________
> > getdp mailing list
> > getdp at geuz.org
> > http://www.geuz.org/mailman/listinfo/getdp
> 
>> Prof. Ruth V. Sabariego
> KU Leuven  
> Dept. Electrical Engineering ESAT/Electa
> http://www.esat.kuleuven.be/electa
> 
> 
> 
> 
> 
> 
> 
> 
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp