[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