[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>