[Getdp] Solving 3D magnetostatic problem

Stepan Roucka stepan at roucka.eu
Fri Jun 19 01:04:45 CEST 2009


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: 6746 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: 2356 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: 2637 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: 4890 bytes
Desc: not available
Url : http://www.geuz.org/pipermail/getdp/attachments/20090619/7347f42d/attachment.par 


More information about the getdp mailing list