[Getdp] Solving 3D magnetostatic problem

[Stepan Roucka]

Hello everybody,
I am trying to solve for magnetic field of two coils
and two ferromagnetic rings. I adopted config files from
http://www.geuz.org/pipermail/getdp/2009/001193.html
to my geometry, but solver gives obviously wrong results with
message for GMRES:
GetDP   : Solver    : Generalized Minimum RESidual (GMRES)
GetDP   :     1  3.1602108e+03  1.0000000e+00
GetDP   :     2  3.1602108e+03  1.0000000e+00
GetDP   : Warning   : Iterative solver is facing a break-down

for CG:
GetDP   :   998  1.7430613e+16  5.5156487e+12
GetDP   :   999  1.8376004e+16  5.8148032e+12
GetDP   : Warning   : Iterative solver has iterated too many times

there is probably something wrong with my config file, but I can't
find what. Could somebody please take a look at it? I just
started learning principles of FEM, so any advice is welcome.

Thanks in advance
Stepan Roucka
-------------- next part --------------
A non-text attachment was scrubbed...
Name: attachment.geo
Type: application/octet-stream
Size: 6747 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20090619/7347f42d/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: attachment.pro
Type: application/octet-stream
Size: 2357 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20090619/7347f42d/attachment.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Magsta_a_3D.pro
Type: application/octet-stream
Size: 2638 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20090619/7347f42d/attachment-0001.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: solver.par
Type: application/octet-stream
Size: 4891 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20090619/7347f42d/attachment.par>