Author Topic: Default linear solver  (Read 9215 times)

aleximos

  • Jr. Member
  • **
  • Posts: 34
Default linear solver
« on: September 11, 2012, 06:36:08 AM »
I would like to ask which is the default linear solver for FEAP by issuing this command?

TANG !form tangent only
LOOP iter 10 !perform 10 modified Newton iterations
   FORM !form residual
   SOLVe !solve linearized equations
NEXT iter !proceed to next iteration

FEAP_Admin

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 993
Re: Default linear solver
« Reply #1 on: September 11, 2012, 07:58:04 AM »
It is a multi-frontal profile solver.

The solver is basically described in the Z&T book vol 1 (appendix) and p344-350 Eng. Comput., 1985, Vol. 2.
The actual solver is a multi-column version of what you will find in these references.



aleximos

  • Jr. Member
  • **
  • Posts: 34
Re: Default linear solver
« Reply #2 on: September 11, 2012, 08:10:40 AM »
Can this solver be used in Parfeap?
I guess that it can, because the multifrontal solver is supposed to work on multiprocessors so it can be used for parallel computing, but as I read in the manual the methods that can be used (with the help of PETSC) are the CG, MINRES and GMRES with preconditioners.

FEAP_Admin

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 993
Re: Default linear solver
« Reply #3 on: September 11, 2012, 08:43:53 AM »
No.

parfeap requires you to use the PETSc solvers -- iterative or direct will work.