KSP/PC choice

2007-07-27 Thread Tim Kröger
Dear Matt On Wed, 25 Jul 2007, Matthew Knepley wrote: -ksp_type preonly -pc_type lu -mat_type aijmumps -ksp_monitor It appears that you never call MatSetFromOptions() on the matrix you create, which would prevent the option from changing the type. Thank you very much for your help.

KSP/PC choice

2007-07-27 Thread gaate...@irisa.fr
Dear Tims, MUMPS is direct solver so you do not need Krylov subspace which are generally iterative solver. May be you are using MUMPS as solver in you problem not as preconditionner. Regards Dear Matt On Wed, 25 Jul 2007, Matthew Knepley wrote: -ksp_type preonly -pc_type lu -mat_type

KSP/PC choice

2007-07-27 Thread Matthew Knepley
You could, of course, embed MUMPS as a preconditioner in a Krylov iteration like GMRES and getthe effect that you want. Mat On 7/27/07, Stephan Kramer stephan.kramer at imperial.ac.uk wrote: Hi Tim, As Guy said MUMPS is a direct solver, i.e. it constructs the inverse of the matrix (in some

KSP/PC choice

2007-07-25 Thread Ben Tay
Hi, I am also trying to solve the NS eqns in 2d, finite volume mtd, using fractional step formulation. I am currently using KSPBCGS and PCILU to solve the momentum eqn matrix and HYPRE's AMG and the default ksptype to solve the poisson eqn. Do you ppl think this is a good choice? Or is there

KSP/PC choice

2007-07-25 Thread Tim Kröger
Dear Satish On Tue, 24 Jul 2007, Satish Balay wrote: Its best to send installation issues involving configure.log to petsc-maint at mcs.anl.gov - and not the list. Okay; sorry for that. --with-cc=gcc --with-fc=gfortran --with-shared=0 --download-mumps=1

KSP/PC choice

2007-07-25 Thread Matthew Knepley
On 7/25/07, Tim Kr?ger tim at cevis.uni-bremen.de wrote: Dear Satish But now, the next problem is that I am not able to use MUMPS. As described in the manual, I start my application using -ksp_type preonly -pc_type lu -mat_type aijmumps -ksp_monitor and I get It appears that you never

KSP/PC choice

2007-07-25 Thread Lisandro Dalcin
On 7/24/07, Ben Tay zonexo at gmail.com wrote: I am also trying to solve the NS eqns in 2d, finite volume mtd, using fractional step formulation. I am currently using KSPBCGS and PCILU to solve the momentum eqn matrix and HYPRE's AMG and the default ksptype to solve the poisson eqn. Do you

KSP/PC choice

2007-07-24 Thread Tim Kröger
Dear Lisandro, On Mon, 23 Jul 2007, Lisandro Dalcin wrote: On 7/23/07, Tim Kr?ger tim at cevis.uni-bremen.de wrote: What kind of stabilization for advection and pressure are you using? Streamline diffusion (for advection). I am not aware of the requirement to stabilize the pressure as

KSP/PC choice

2007-07-24 Thread gaate...@irisa.fr
Dear Tim, I affraid that you have same memory problem while using MUMPS. I have observed same memory problem while trying to solve some 3D Stokes equations with MUMPS. I think the better way to solve this type of problem is KSP with appropriate PC. Regards, Guy Dear Lisandro, On Mon, 23 Jul

KSP/PC choice

2007-07-24 Thread Satish Balay
On Tue, 24 Jul 2007, owner-petsc-users at mcs.anl.gov wrote: Date: Tue, 24 Jul 2007 09:47:37 +0200 (CEST) From: =?iso-8859-15?Q?Tim_Kr=F6ger?= tim at cevis.uni-bremen.de X-X-Sender: tim at elektrik To: petsc-users at mcs.anl.gov Subject: BOUNCE petsc-users at mcs.anl.gov: Message too

KSP/PC choice

2007-07-23 Thread Tim Kröger
Dear all, the more KSP and PC methods are available, the more difficult it is to choose a suitable one. I want to compute Navier-Stokes (using Finite Elements) for medium-sized Reynolds numbers (still laminar, though). I am currently using -ksp_type gmres -ksp_gmres_restart 30 -pc_type

KSP/PC choice

2007-07-23 Thread Lisandro Dalcin
On 7/23/07, Tim Kr?ger tim at cevis.uni-bremen.de wrote: Dear all, the more KSP and PC methods are available, the more difficult it is to choose a suitable one. I want to compute Navier-Stokes (using Finite Elements) for medium-sized Reynolds numbers (still laminar, though). I am currently

KSP/PC choice

2007-07-23 Thread Tim Kröger
Dear Lissandro, On Mon, 23 Jul 2007, Lisandro Dalcin wrote: On 7/23/07, Tim Kr?ger tim at cevis.uni-bremen.de wrote: the more KSP and PC methods are available, the more difficult it is to choose a suitable one. I want to compute Navier-Stokes (using Finite Elements) for medium-sized

KSP/PC choice

2007-07-23 Thread Matthew Knepley
1) Until you run out of memory, I would use sparse direct like MUMPS 2) After that, as long as you have the memory I would increase the GMRES vectors, say to 50 or 100. 3) After that I would try LGMRES which generally converges better on these problems. Matt On 7/23/07, Tim Kr?ger tim

KSP/PC choice

2007-07-23 Thread gaate...@irisa.fr
Dear Tims, The problem is how to solve more accuretly the local problem in ASM. You can also use KSP to solve local problem. I think you can look how to do it with PCASMGetSubKSP(). Best Regards, Guy Dear gaatenek On Mon, 23 Jul 2007, gaatenek at irisa.fr wrote: Do you monitoring the

KSP/PC choice

2007-07-23 Thread Lisandro Dalcin
On 7/23/07, Tim Kr?ger tim at cevis.uni-bremen.de wrote: Monolithic. (I am not quite sure about this word but I assume it means that I don't try to decouple the two equations.) That's the meaning for me (not sure if the word is completelly correct) What kind of stabilization for advection