Changes: 2.0.22#

General:

  • The Fortran90-specific routines such as VecGetArrayF90() now work with the IBM Fortran compiler, the SGI Fortran 90 compiler, and the Cray T3E Fortran compiler.

  • The third argument for the Fortran versions of the routines VecGetArray()/VecRestoreArray(), ISGetIndices()/ISRestoreIndices(),             MatGetArray()/MatRestoreArray() and DAGetGlobalIndices() has been changed from integer to PetscOffset. See src/snes/examples/tutorials/ex5f.F for example usage.

  • Changed the location of the manual pages from docs/www/www.html and docs/www to docs/manualpages/manualpages.html and docs/manualpages.

  • The location of the mpiuni include files (used in the base.site file when not compiling with an MPI) has moved from src/mpiuni to src/sys/src/mpiuni

  • Added the option -log_summary_exclude <vec,mat,sles,snes> to limit the information printed in the summary table

  • Added the option -log_info_exclude <vec,mat,sles,snes>

  • Changed the calling sequences of PetscGetTime(),             PetscGetCPUTime(), PetscGetFlops() to return error codes like all other PETSc functions.

  • Changed the internal structure of the PETSc objects. This will not effect most users but makes PETSc more extensible for the future.

  • Removed the PETSc object child support and replaced it with

    These allow the user to attach objects and functions to any PETSc object.

  • Added a feeble example that demonstrates how to manage partitioning a grid in parallel in src/mat/examples/tutorials/ex2.c

  • The build command ‘make all’ in ${PETSC_DIR} will build the PETSc libraries, including the fortran interface.

AO (Application Orderings):

TS (Timestepping Solvers):

  • The calling sequence of TSSetType() has changed; see the manual page.

  • TSType is now a string, rather than enum type

SNES (Nonlinear Solvers):

  • The calling sequence of SNESSetType() has changed; see the manual page.

  • SNESType is now a string, rather than enum type

SLES (Linear Solvers):

  • See PC and KSP

KSP (Krylov Subspace Methods):

  • The calling sequence of KSPSetType() has changed; see the manual page.

  • KSPType is now a string, rather than enum type.

PC (Preconditioners):

  • Added two new PC classes, PCSLES and PCCOMPOSITE that allow use of any SLES object as a preconditioner and enable combining several preconditioners.

  • The calling sequence of PCSetType() has changed; see the manual page.

  • PCType is now a string, rather than enum type.

MAT (Matrices):

  • Added support for partitioning using the ParMETIS parallel partitioning package, see the manual page for PartitioningCreate()

  • Added MatGetColumnVector()

  • Added argument column size to MatGetSubMatrix()

  • Changed the memory allocation algorithm in MatLUFactorSymbolic() and MatILUFactorSymbolic() based on a suggestion by David Hysom of ODU.

  • The third argument for the Fortran version of the routines MatGetArray()/MatRestoreArray() has been changed from integer to PetscOffset.

DA (Distributed Arrays):

  • Changed DAGetDistributedVector() and DAGetLocalVector() to DACreateGlobalVector() and DACreateLocalVector(). You should destroy these vectors with a usual call to VecDestroy().

  • The third argument for the Fortran version of the routine DAGetGlobalIndices() has been changed from integer to PetscOffset.

VEC (Vectors):

  • Added VecSetBlockSize(), VecSetValuesBlocked(),             VecSetValuesBlockedLocal(),             VecSetLocalToGlobalMappingBlocked()

  • AddedVecCreatedShared(), which creates a parallel vector in shared memory on the SGI machines

  • Changed VEC_IGNORE_OFF_PROCESSOR_ENTRIES to VEC_IGNORE_OFF_PROC_ENTRIES

  • The third argument for the Fortran versions of the routines VecGetArray()/VecRestoreArray(), has been changed from integer to PetscOffset. See src/snes/examples/tutorials/ex5f.F for example usage.

  • VecGetArray(), VecRestoreArray() now work from Fortran on all machines, finally after 3 long years!

  • Changed the calling sequence of VecCreate(), now takes both local and global size, see the manual page.

IS (Index Sets):

  • The third argument for the Fortran versions of the routines ISGetIndices()/ISRestoreIndices() has been changed from integer to PetscOffset.

Draw (Graphics):

Viewers:

System:

Error Handling:

Event Logging:

Fortran Interface:

  • PETSC_NULL has been replaced with PETSC_NULL_INTEGER, PETSC_NULL_SCALAR, PETSC_NULL_DOUBLE or PETSC_NULL_CHARACTER depending on the context for usage. If the function expects an integer, a scalar (double or complex depending if you compile with BOPT=*_complex), a double precision array, or a string.

  • The Fortran include files are now located in the directory “include/finclude/.”, rather than “include/FINCLUDE/.h”. The include files that used to be in “include/finclude/.h” (not recommended) are now in “include/foldinclude/.h”

  • If you have stored PETSc objects in Fortran integer arrays, such as user(3) = vec you have to declare those arrays to be of type PetscFortranAddr instead of as integer. Note that PetscFortranAddr is simply integer on 32 bit machines and integer*8 on 64-bit machines.

  • The Fortran90-specific routines such as VecGetArrayF90() now work with the IBM Fortran compiler, the SGI Fortran 90 compiler, and the Cray T3E Fortran compiler.

  • PETSc objects can now be passed between C and Fortran directly. You do not have to use the commands PetscCObjectToFortranObject() and PetscFortranObjectToCObject() to translate the objects between languages.

  • Changed Fortran object argument from type int to type PetscFortranAddr in the routines PetscCObjectToFortranObject() and PetscFortranObjectToCObject().

  • Added the routine PetscInitializeFortran() to assist in mixed-language use of PETSc. See the manual page for details.