[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Using the Spooles Cholesky solver via PETSc in serial


   Stephen,

    This was an error in our last release and did not properly fixed
in a patch (why?). It will be in the next release.

   You can edit include/finclude/petscmat.h and add

      PetscEnum MAT_USE_COMPRESSEDROW,MAT_DO_NOT_USE_COMPRESSEDROW
      PetscEnum MAT_IGNORE_LOWER_TRIANGULAR,MAT_ERROR_LOWER_TRIANGULAR
      PetscEnum MAT_GETROW_UPPERTRIANGULAR

      parameter (MAT_USE_COMPRESSEDROW=89)
      parameter (MAT_DO_NOT_USE_COMPRESSEDROW=90)
      parameter (MAT_IGNORE_LOWER_TRIANGULAR=91)
      parameter (MAT_ERROR_LOWER_TRIANGULAR=92)
      parameter (MAT_GETROW_UPPERTRIANGULAR=93)

then it should work properly.

    Barry


On Fri, 2 Jun 2006, Stephen R Ball wrote:


Hi

I am using PETSc v2.3.1 (patch 9) via the Fortran interface. When I use the
Spooles CHOLESKY solver in serial I get the error:

	MatSetValues_MPIBAIJ_MatScaler() line 259 in
src/mat/impls/sbaij/mpi/mpisbaij.c
	Lower triangular value cannot be set for sbaij format. Ignoring
these values, run with -mat_ignore_lower_triagular or call
	MatSetOption(mat,MAT_IGNORE_LOWER_TRIANGULAR)
	MatSetValues() line 827 in src/mat/interface/matrix

The command line option works fine but when I try to use MatSetOption() with
the MAT_IGNORE_LOWER_TRIANGULAR flag I get the error:

	This name does not have a type, and must have an explicit type.
[MAT_IGNORE_LOWER_TRIANGULAR]

Can you tell me if use of the MAT_IGNORE_LOWER_TRIANGULAR flag is supported
in the Fortran interface for MatSetOption()?

Regards

Stephen R. Ball
--
_______________________________________________________________________________

The information in this email and in any attachment(s) is commercial in confidence. If you are not the named addressee(s) or if you receive this email in error then any distribution, copying or use of this communication or the information in it is strictly prohibited.  Please notify us immediately by email at admin.internet(at)awe.co.uk, and then delete this message from your computer.  While attachments are virus checked, AWE plc does not accept any liability in respect of any virus which is not detected.

AWE Plc
Registered in England and Wales
Registration No 02763902
AWE, Aldermaston, Reading, RG7 4PR