When using the Spooles Cholesky solver I am getting the error:
0:[0]PETSC ERROR: MatSetValues_SeqSBAIJ() line 792 in
src/mat/impls/sbaij/seq/sbaij.c 0:[0]PETSC ERROR: ! 0:[0]PETSC ERROR: Lower
triangular value cannot be set for sbaij format. Ignoring these values, run
with -mat_ignore_lower_triangular or call
MatSetOption(mat,MAT_IGNORE_LOWER_TRIANGULAR)!
0:[0]PETSC ERROR: MatSetValues() line 702 in src/mat/interface/matrix.c
When I attempt to call MatSetOption(mat,MAT_IGNORE_LOWER_TRIANGULAR,ierr)
via the Fortran interface, I get the compilation error that entity
mat_ignore_lower_triangular has undefined type. Can you tell me if this
option is supported via the Fortran interface?
This is certainly there in the development version, and we plan to release in a week
or so. However, a very easy fix is to update the list of enums at the top of
include/finclude/petscmat.h
with the enum from
include/petscmat.h
Thanks,
Matt
Note that using runtime option -mat_ignore_lower_triangular works, but I
would prefer to set this option in my code using the MatSetOption() routine.
Regards
Stephen
--
_______________________________________________________________________________
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