[Getdp] iterative method

Christophe Geuzaine Christophe.Geuzaine at ulg.ac.be
Mon Nov 5 10:36:31 CET 2001


I did not get the attachments. 

Try to precondition your system a bit more (i.e. increase the NbFill
parameter in the solver.par file).

Christophe

Lin Ji wrote:
> 
> Dear Christophe,
>     I am sorry to bother you again. I still can not make the Helmholtz
> equation solved correctly. It seems there is something wrong with the
> iterative method used to solve the linear system. I got the book 'Iterative
> Methods for Sparse Linear Systems'. But, it seems it requires quite a bit
> learning curve in order to understand the methods and the preconditioning
> technique. Can you take a look of my problem again? I would really appreciate
> it.
>     I attached 'hhltz.pro', 'hhltz.geo' and 'solver.par'. The following is
> the output of running getdp on hhltz.pro:
> 
> /personal/research/proj_rep/hhltz2D (66) ../../getdp-0.84/getdp hhltz -pre
> -cal -pos
> Loading   : Problem definition 'hhltz.pro'
> Info      : Available Resolutions
> (1) hhltz
> Choice: 1
> Info      : Available PostOperations
> (1) u
> Choice: 1
> Info      : Selected Resolution 'hhltz'
> Loading   : Geometric data 'hhltz.msh'
> Info      : System 'A' : Real
> P r e - P r o c e s s i n g . . .
> Operation : Treatment Formulation 'hhltz'
> Resources : cpu 11.670000 s / mem 0 kb
> 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[A]
> Solver    : Loading parameter file 'solver.par'
> Resources : cpu 47.360000 s / mem 0 kb
> Operation : Solve[A]
> Solver    : No scaling of system of equations
> Solver    : RCMK algebraic renumbering
> Solver    : N: 186723, NZ: 2140569, BW max/avg: 28/11, SW max: 186222
> Resources : cpu 52.100000 s / mem 0 kb
> Solver    : ILUTP (Float, fill-in = 20)
> Solver    : N: 186723, NZ: 6165203, BW max/avg: 40/33, SW max: 172178
> Resources : cpu 82.380000 s / mem 0 kb
> Solver    : Generalized Minimum RESidual (GMRES)
>    1  1.3449314e-01  1.0000000e+00
> Warning   : Iterative solver is facing a break-down
>    2  0.0000000e+00  0.0000000e+00
> Solver    : 2 Iterations / Residual: 0.134493
> Resources : cpu 83.890000 s / mem 0 kb
> Operation : SaveSolution[A]
> Resources : cpu 84.250000 s / mem 0 kb
> E n d   P r o c e s s i n g
> P o s t - P r o c e s s i n g . . .
> Info      : Selected PostProcessing 'hhltz'
> Info      : Selected Mesh 'hhltz.msh'
> Operation : PostOperation 1/1
>           > 'uh.pos'
> Info      : Initializing rapid search grid 20x20x1
> Resources : cpu 189.980000 s / mem 0 kb
> E n d   P o s t - P r o c e s s i n g
> E n d
> 
> Best regards,
> 
> Lin
> 
> --
> Dept. of Math., RPI
> 518-276,6904 (work)
> _______________________________________________
> GetDP mailing list
> GetDP at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp

-- 
Christophe Geuzaine

Tel: 32 (0) 4 366 37 10    http://geuz.org
Fax: 32 (0) 4 366 29 10    mailto:Christophe.Geuzaine at ulg.ac.be