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

[Jed Brown]

On Mon, Nov 12, 2012 at 8:43 AM, Alexander Grayver
<agrayver at gfz-potsdam.de>wrote:

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

No, this was missed in the rename. It's fixed here.

https://bitbucket.org/petsc/petsc-dev/changeset/cfa58f9da86656f6ddf509089493a99a1d503d6f


This should probably be back-ported to 3.3, though it will break an
existing Fortran user that forgot to fix that call (or noticed the bug and
didn't report it).

The petsc.cs.iit.edu server is down right now and although we are migrating
to bitbucket, we haven't decided how to manage the release repository there
yet.

https://bitbucket.org/petsc/petsc-dev/issue/22/backport-bug-fixes-to-33-once-release


>
> *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);
> }*
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121112/ff7393df/attachment-0001.html>