[Getdp] Solving 3D magnetostatic problem

Stepan Roucka stepan at roucka.eu
Sun Jun 28 20:11:35 CEST 2009


Just to let you know, that I wrote a conversion script
from msh vector data to vtk unstructured xml file,
so no advice is needed anymore.

Stepan Roucka

2009/6/24 Stepan Roucka <stepan at roucka.eu>:
> Well, I rewrote the scripts and now it works. I still don't know
> where the problem was. I'll probably add this to the wiki,
> because 3D magnetostatic computations are quite useful IMO.
>
> However, now I have another problem: Is it possible to
> export the resulting vector field to some format readable
> by ParaView (vtk?)? If not, are there some conversion scripts?
> What format is the best for writing my own conversion utility?
>
> Thanks
> Stepan Roucka
>
> 2009/6/19 Stepan Roucka <stepan at roucka.eu>:
>> 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
>>
>