[Getdp] help on tutorial

David Colignon David.Colignon at ulg.ac.be
Fri Aug 25 16:29:30 CEST 2006


Hi,

your mesh must be wrong. You get a 4x4 matrix to solve ... ( Solver    : N: 4, NZ: 4, BW max/avg: 1/1, SW max: 1)
Visualize the mesh with "gmsh Core.msh" to check.

Cheers,

Dave



Faisal Ahmed wrote:
> Hello Everybody
> 
> i am quite new in the field of computational electromagnetics & as well in
> getdp. i was trying to do simulation with one of the tutorials files named
> as 'Mag_av_2D.pro', for which the geometry is defined in file Core.geo.
> I suppose everybody knows which tutorial I'm talking about. I converted
> the Core.geo file to .msh format in GMSH. When I try to run the simulation
> with
> getdp CoreMassive -msh Core.msh -solve MagDyn_av_2D
> 
> I get error message as below:
> 
> Loading   : Problem definition 'CoreMassive.pro'
> Loading   : Problem definition 'Jacobian_Lib.pro'
> Loading   : Problem definition 'Integration_Lib.pro'
> Loading   : Problem definition 'MagDyn_av_2D.pro'
> Info      : Selected Resolution 'MagDyn_av_2D'
> Loading   : Geometric data 'Core.msh'
> Info      : System 'Sys_Mag' : Complex, Frequency = 1 Hz
> P r e - P r o c e s s i n g . . .
> Operation : Treatment Formulation 'Magnetodynamics_av_2D'
> Resources : cpu 0.008 s
> E n d   P r e - P r o c e s s i n g
> P r o c e s s i n g . . .
> Operation : Generate[Sys_Mag]
> Solver    : Loading parameter file 'solver.par'
> Info      : Setting System {A,b} to zero
> Resources : cpu 0.008 s
> Operation : Solve[Sys_Mag]
> Solver    : Scaling system of equations
> Warning   : Scaling aborted because of zero row (3)
> Solver    : RCMK algebraic renumbering
> Solver    : N: 4, NZ: 4, BW max/avg: 1/1, SW max: 1
> Resources : cpu 0.008 s
> Solver    : ILUT (Float, fill-in = 150)
> Error     : Zero row encountered in ILU
> 
> However the solver.par looks as follows:
> 
>             Matrix_Format            1
>           Matrix_Printing            0
>            Matrix_Storage            0
>                   Scaling            4
>     Renumbering_Technique            1
>            Preconditioner            1
>   Preconditioner_Position            2
>                   Nb_Fill           150
>     Permutation_Tolerance         0.001
>        Dropping_Tolerance            0
>     Diagonal_Compensation            0
>                Re_Use_ILU            1
>                 Algorithm           8
>               Krylov_Size           40
>           IC_Acceleration            1
>                 Re_Use_LU            0
>     Iterative_Improvement            0
>               Nb_Iter_Max         10000
>             Stopping_Test        1e-8
> 
> 
> It will be very nice to have some opininion for the solution of this
> problem. I hope I am not bothering everybody with a very small problem.
> 
> Regards
> Rookie
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp