Re: [petsc-users] Controlling MUMPS parameters inside PETSC

Evan: When PCLU is requested, matrix is treated as unsymmetric by mumps and petsc, thus the setting of id.sym=0 is actually used. If user requests PCCholesky, then mumps needs to be set as sym.id=1 (spd) or 2 (symmetric). In mumps interface, we do #if defined(PETSC_USE_COMPLEX) mumps->sym =

Re: [petsc-users] Controlling MUMPS parameters inside PETSC

Hi Hong, Thanks for your reply. When I use standalone MUMPS for a symmetric matrix, I pass a upper/lower triangular part to MUMPS. To do this, I use id.sym=2. Under PETSC environment, KSPCreate(PETSC_COMM_WORLD, ); KSPSetOperators(ksp, A, A); KSPSetType (ksp, KSPPREONLY); KSPGetPC(ksp, );

Re: [petsc-users] Controlling MUMPS parameters inside PETSC

Evan: > Dear PETSC-users, > > I use parallel direct solver MUMPS inside PETSC and need to control some > MUMPS parameters inside PETSC. For example, I want to set up MUMPS > parameters as shown below. > > ZMUMPS_STRUC_C id; > id.job=-1; /* Initialize mumps instance*/ > This is set in

[petsc-users] Controlling MUMPS parameters inside PETSC

Dear PETSC-users, I use parallel direct solver MUMPS inside PETSC and need to control some MUMPS parameters inside PETSC. For example, I want to set up MUMPS parameters as shown below. ZMUMPS_STRUC_C id; id.job=-1; /* Initialize mumps instance*/ id.par=1; /* 0: host is not involved in

Re: [petsc-users] Good recommendation on meshing package?

On Fri, Sep 29, 2017 at 11:06 AM, Zou, Ling wrote: > Hi all, > > I know this is a bit off topic on PETSc email list. > I would like to try some finite volume type of CFD algorithm with PETSc, > but I found it quite troublesome to manage mesh by myself. > I wonder if there is