[Getdp] example of fixed point iteration in GetDP

Ruth Vazquez Sabariego Ruth.Sabariego at esat.kuleuven.be
Tue Jul 15 09:23:11 CEST 2014


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<mailto: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<mailto: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







-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140715/2ce6cf88/attachment.html>