[petsc-users] Does PETSc take advantage of symmetry in iterative solvers

[Alexander Grayver]

On 12.11.2012 13:40, Jed Brown wrote:
> The option does not change the format, but it does change some 
> defaults (e.g., it makes ASM symmetric instead of using RASM). If you 
> want to use a symmetric _format_, use MATSBAIJ.

Jed,

Thanks.
Does it make sense that in fortran one has to call MatCreateMPISBAIJ?

/EXTERN_C_BEGIN
void PETSC_STDCALL matcreatempisbaij_(MPI_Comm *comm,PetscInt 
*bs,PetscInt *m,PetscInt *n,PetscInt *M,PetscInt *N,
          PetscInt *d_nz,PetscInt *d_nnz,PetscInt *o_nz,PetscInt 
*o_nnz,Mat *newmat,PetscErrorCode *ierr)
{
   CHKFORTRANNULLINTEGER(d_nnz);
   CHKFORTRANNULLINTEGER(o_nnz);
   *ierr = MatCreateSBAIJ(MPI_Comm_f2c(*(MPI_Fint 
*)&*comm),*bs,*m,*n,*M,*N,*d_nz,d_nnz,*o_nz,o_nnz,newmat);
}
/

-- 
Regards,
Alexander

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121112/1418a2fa/attachment.html>