PETSc

Documentation: Changes: 2.2.1

• Home
• Download
• Features
• Documentation
  • Manual pages and Users Manual
  • Referencing PETSc
  • Tutorials
  • Installation
  • SAWs
  • Changes
    • dev
    • 3.15
    • 3.14
    • 3.13
    • 3.12
    • 3.11
    • 3.10
    • 3.9
    • 3.8
    • 3.7
    • 3.6
    • 3.5
    • 3.4
    • 3.3
    • 3.2
    • 3.1
    • 3.0.0
    • 2.3.3
    • 2.3.2
    • 2.3.1
    • 2.3.0
    • 2.2.1
    • 2.2.0
    • 2.1.6
    • 2.1.5
    • 2.1.3
    • 2.1.2
    • 2.1.1
    • 2.1.0
    • 2.0.29
    • 2.0.28
    • 2.0.24
    • 2.0.22
    • 2.0.18-21
    • 2.0.17
    • 2.0.16
    • 2.0.15
  • Bug Reporting
  • Code Management
  • FAQ
  • License
  • Linear Solver Table
  • Nonlinear Solver Table
  • TAO Solver Table
• Applications/Publications
• Miscellaneous
• External Software
• Developers Site

General:

• Introduced 4 new PETSc data types: PetscInt, PetscErrorCode, PetscMPIInt and PetscBLASInt. For 99% of users these are just int or integer*4 and you do not need to change your code.
• For users with more than roughly 2 billion unknowns you can run configure with --with-64-bit-ints and then PetscInt will represent 64 bit integers, long long int in C and integer*8 in Fortran. But the the other 3 types remain 32 bit (i.e. int in C and integer*4 in Fortran). Now you can index into vectors and matrices with virtually unlimited sizes.

Vec:

• VecSetStashInitialSize() -> VecStashSetInitialSize()
• VecConvertMPIToAll() and VecConvertMPIToZero() have been replaced with VecScatterCreateToAll() and VecScatterCreateToZero()

IS:

• ISSum() now has a throughput argument instead of one in and one output.

Mat:

• MatSetStashInitialSize() -> MatStashSetInitialSize()
• added MatFactorInfoInitialize()
• added MatStashGetInfo()
• added MatGetVecs()
• added MatIsSymmetric(), MatIsTranspose() now take a tolerance flag for equality (use 0.0 for identical)
• added MatMPIAIJSetPreallocationCSR() and MatMPIBAIJSetPreallocationCSR()

PC:

• ICC now uses Manteuffel shift by default
• If matrix is symmetric, PC defaults to ICC when possible, if you told with MatSetOption(mat,MAT_SYMMETRIC)
• PCSetVector() is gone (no longer needed)
• PCNullSpaceAttach() is gone, use KSPSetNullSpace()
• MatGetSubmatrices() now preserves symmetric flag of original matrix
• hypre BoomerAMG (-pc_type hypre -pc_hypre_type boomeramg) now uses one 1 cycle of multigrid as preconditioner (Previously it used either a decrease in residual of 10^-7 or a maximum of 20 iterations

KSP:

• Note that a relative residual tolerance now means that ||r_k|| < tol*||b||
• Removed KSPSetRhs() and KSPSetSolution(), pass
• Added KSPSetNullSpace(), KSPGetNullSpace()

DA:

• Added DAGetGhostedCoordinates() and DAGetCoordinateDA() (see src/dm/da/examples/tutorials/ex3.c)
• Added DACoor2d and DACoor3d (see src/dm/da/examples/tutorials/ex3.c)
• DAGetInterpolation() now provides the correct interpolation for nonuniform grids using the coordinate information provided with DASetCoordinates(). (see src/dm/da/examples/tutorials/ex3.c)

config/configure.py:

• --with-c/f-blas-lapack -> --download-c/f-blas-lapack=yes
• --with-c/f-blas-lapack-if-needed -> --download-c/f-blas-lapack=ifneeded
• added --download-mpich=yes,no,ifneeded
• added --download-mpich-device=
• added --download-mpich-machines=[machines]
• added configure including download support for hypre, parmetis
• added configure support (without download) to Mumps, UMF, Superlu,SuperLU_dist,DSCPACK
• Made configure much more robust

SNES:

• Added -snes_converged_reason
• Added SNESSetRhs(snes,rhs) to allow for solving F(x) = rhs

TS:

DMMG:

• Added DMMGSetNullSpace()

SYS:

Fortran:

ExternalPackages: