[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