1df4397b0SStefano Zampini #if !defined(__PETSCSYSTYPES_H) 2df4397b0SStefano Zampini #define __PETSCSYSTYPES_H 3df4397b0SStefano Zampini 4df4397b0SStefano Zampini #include <petscconf.h> 5df4397b0SStefano Zampini #include <petscfix.h> 6df4397b0SStefano Zampini 7df4397b0SStefano Zampini /*MC 8df4397b0SStefano Zampini PetscErrorCode - datatype used for return error code from almost all PETSc functions 9df4397b0SStefano Zampini 10df4397b0SStefano Zampini Level: beginner 11df4397b0SStefano Zampini 12df4397b0SStefano Zampini .seealso: CHKERRQ, SETERRQ 13df4397b0SStefano Zampini M*/ 14df4397b0SStefano Zampini typedef int PetscErrorCode; 15df4397b0SStefano Zampini 16df4397b0SStefano Zampini /*MC 17df4397b0SStefano Zampini 18df4397b0SStefano Zampini PetscClassId - A unique id used to identify each PETSc class. 19df4397b0SStefano Zampini 20df4397b0SStefano Zampini Notes: 21df4397b0SStefano Zampini Use PetscClassIdRegister() to obtain a new value for a new class being created. Usually 22df4397b0SStefano Zampini XXXInitializePackage() calls it for each class it defines. 23df4397b0SStefano Zampini 24df4397b0SStefano Zampini Developer Notes: 25df4397b0SStefano Zampini Internal integer stored in the _p_PetscObject data structure. 26df4397b0SStefano Zampini These are all computed by an offset from the lowest one, PETSC_SMALLEST_CLASSID. 27df4397b0SStefano Zampini 28df4397b0SStefano Zampini Level: developer 29df4397b0SStefano Zampini 30df4397b0SStefano Zampini .seealso: PetscClassIdRegister(), PetscLogEventRegister(), PetscHeaderCreate() 31df4397b0SStefano Zampini M*/ 32df4397b0SStefano Zampini typedef int PetscClassId; 33df4397b0SStefano Zampini 34df4397b0SStefano Zampini /*MC 35df4397b0SStefano Zampini PetscMPIInt - datatype used to represent 'int' parameters to MPI functions. 36df4397b0SStefano Zampini 37df4397b0SStefano Zampini Level: intermediate 38df4397b0SStefano Zampini 39df4397b0SStefano Zampini Notes: 40df4397b0SStefano Zampini usually this is the same as PetscInt, but if PETSc was built with --with-64-bit-indices but 41df4397b0SStefano Zampini standard C/Fortran integers are 32 bit then this is NOT the same as PetscInt; it remains 32 bit. 42df4397b0SStefano Zampini 43df4397b0SStefano Zampini PetscMPIIntCast(a,&b) checks if the given PetscInt a will fit in a PetscMPIInt, if not it 44df4397b0SStefano Zampini generates a PETSC_ERR_ARG_OUTOFRANGE error. 45df4397b0SStefano Zampini 46df4397b0SStefano Zampini .seealso: PetscBLASInt, PetscInt, PetscMPIIntCast() 47df4397b0SStefano Zampini 48df4397b0SStefano Zampini M*/ 49df4397b0SStefano Zampini typedef int PetscMPIInt; 50df4397b0SStefano Zampini 51df4397b0SStefano Zampini /*MC 52df4397b0SStefano Zampini PetscEnum - datatype used to pass enum types within PETSc functions. 53df4397b0SStefano Zampini 54df4397b0SStefano Zampini Level: intermediate 55df4397b0SStefano Zampini 56df4397b0SStefano Zampini .seealso: PetscOptionsGetEnum(), PetscOptionsEnum(), PetscBagRegisterEnum() 57df4397b0SStefano Zampini M*/ 58df4397b0SStefano Zampini typedef enum { ENUM_DUMMY } PetscEnum; 59df4397b0SStefano Zampini 60df4397b0SStefano Zampini typedef short PetscShort; 61df4397b0SStefano Zampini typedef char PetscChar; 62df4397b0SStefano Zampini typedef float PetscFloat; 63df4397b0SStefano Zampini 64df4397b0SStefano Zampini /*MC 65df4397b0SStefano Zampini PetscInt - PETSc type that represents an integer, used primarily to 66df4397b0SStefano Zampini represent size of arrays and indexing into arrays. Its size can be configured with the option --with-64-bit-indices to be either 32-bit (default) or 64-bit. 67df4397b0SStefano Zampini 68df4397b0SStefano Zampini Notes: 69df4397b0SStefano Zampini For MPI calls that require datatypes, use MPIU_INT as the datatype for PetscInt. It will automatically work correctly regardless of the size of PetscInt. 70df4397b0SStefano Zampini 71df4397b0SStefano Zampini Level: beginner 72df4397b0SStefano Zampini 73df4397b0SStefano Zampini .seealso: PetscBLASInt, PetscMPIInt, PetscReal, PetscScalar, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT 74df4397b0SStefano Zampini M*/ 75df4397b0SStefano Zampini 76df4397b0SStefano Zampini #if defined(PETSC_HAVE_STDINT_H) 77df4397b0SStefano Zampini # include <stdint.h> 78df4397b0SStefano Zampini #endif 79df4397b0SStefano Zampini #if defined (PETSC_HAVE_INTTYPES_H) 80df4397b0SStefano Zampini # if !defined(__STDC_FORMAT_MACROS) 81df4397b0SStefano Zampini # define __STDC_FORMAT_MACROS /* required for using PRId64 from c++ */ 82df4397b0SStefano Zampini # endif 83df4397b0SStefano Zampini # include <inttypes.h> 84df4397b0SStefano Zampini # if !defined(PRId64) 85df4397b0SStefano Zampini # define PRId64 "ld" 86df4397b0SStefano Zampini # endif 87df4397b0SStefano Zampini #endif 88df4397b0SStefano Zampini 89df4397b0SStefano Zampini #if defined(PETSC_HAVE_STDINT_H) && defined(PETSC_HAVE_INTTYPES_H) && defined(PETSC_HAVE_MPI_INT64_T) /* MPI_INT64_T is not guaranteed to be a macro */ 90df4397b0SStefano Zampini typedef int64_t PetscInt64; 91df4397b0SStefano Zampini #elif (PETSC_SIZEOF_LONG_LONG == 8) 92df4397b0SStefano Zampini typedef long long PetscInt64; 93df4397b0SStefano Zampini #elif defined(PETSC_HAVE___INT64) 94df4397b0SStefano Zampini typedef __int64 PetscInt64; 95df4397b0SStefano Zampini #else 96df4397b0SStefano Zampini # error "cannot determine PetscInt64 type" 97df4397b0SStefano Zampini #endif 98df4397b0SStefano Zampini 99df4397b0SStefano Zampini #if defined(PETSC_USE_64BIT_INDICES) 100df4397b0SStefano Zampini typedef PetscInt64 PetscInt; 101df4397b0SStefano Zampini #else 102df4397b0SStefano Zampini typedef int PetscInt; 103df4397b0SStefano Zampini #endif 104df4397b0SStefano Zampini 105df4397b0SStefano Zampini /*MC 106df4397b0SStefano Zampini PetscBLASInt - datatype used to represent 'int' parameters to BLAS/LAPACK functions. 107df4397b0SStefano Zampini 108df4397b0SStefano Zampini Notes: 109df4397b0SStefano Zampini Usually this is the same as PetscInt, but if PETSc was built with --with-64-bit-indices but 110df4397b0SStefano Zampini standard C/Fortran integers are 32 bit then this is NOT the same as PetscInt it remains 32 bit 111df4397b0SStefano Zampini (except on very rare BLAS/LAPACK implementations that support 64 bit integers see the note below). 112df4397b0SStefano Zampini 113df4397b0SStefano Zampini PetscErrorCode PetscBLASIntCast(a,&b) checks if the given PetscInt a will fit in a PetscBLASInt, if not it 114df4397b0SStefano Zampini generates a PETSC_ERR_ARG_OUTOFRANGE error 115df4397b0SStefano Zampini 116df4397b0SStefano Zampini Installation Notes: 117df4397b0SStefano Zampini The 64bit versions of MATLAB ship with BLAS and LAPACK that use 64 bit integers for sizes etc, 118df4397b0SStefano Zampini if you run ./configure with the option 119df4397b0SStefano Zampini --with-blaslapack-lib=[/Applications/MATLAB_R2010b.app/bin/maci64/libmwblas.dylib,/Applications/MATLAB_R2010b.app/bin/maci64/libmwlapack.dylib] 120df4397b0SStefano Zampini but you need to also use --known-64-bit-blas-indices. 121df4397b0SStefano Zampini 122df4397b0SStefano Zampini MKL also ships with 64 bit integer versions of the BLAS and LAPACK, if you select those you must also ./configure with 123df4397b0SStefano Zampini --known-64-bit-blas-indices 124df4397b0SStefano Zampini 125df4397b0SStefano Zampini OpenBLAS can also be built to use 64 bit integers. The ./configure options --download-openblas -download-openblas-64-bit-blas-indices 126df4397b0SStefano Zampini will build a 64 bit integer version 127df4397b0SStefano Zampini 128df4397b0SStefano Zampini Developer Notes: 129df4397b0SStefano Zampini Eventually ./configure should automatically determine the size of the integers used by BLAS/LAPACK. 130df4397b0SStefano Zampini 131df4397b0SStefano Zampini External packages such as hypre, ML, SuperLU etc do not provide any support for passing 64 bit integers to BLAS/LAPACK so cannot 132df4397b0SStefano Zampini be used with PETSc if you have set PetscBLASInt to long int. 133df4397b0SStefano Zampini 134df4397b0SStefano Zampini Level: intermediate 135df4397b0SStefano Zampini 136df4397b0SStefano Zampini .seealso: PetscMPIInt, PetscInt, PetscBLASIntCast() 137df4397b0SStefano Zampini 138df4397b0SStefano Zampini M*/ 139df4397b0SStefano Zampini #if defined(PETSC_HAVE_64BIT_BLAS_INDICES) 140df4397b0SStefano Zampini typedef PetscInt64 PetscBLASInt; 141df4397b0SStefano Zampini #else 142df4397b0SStefano Zampini typedef int PetscBLASInt; 143df4397b0SStefano Zampini #endif 144df4397b0SStefano Zampini 145df4397b0SStefano Zampini /*E 146df4397b0SStefano Zampini PetscBool - Logical variable. Actually an int in C and a logical in Fortran. 147df4397b0SStefano Zampini 148df4397b0SStefano Zampini Level: beginner 149df4397b0SStefano Zampini 150df4397b0SStefano Zampini Developer Note: 151df4397b0SStefano Zampini Why have PetscBool , why not use bool in C? The problem is that K and R C, C99 and C++ all have different mechanisms for 152df4397b0SStefano Zampini boolean values. It is not easy to have a simple macro that that will work properly in all circumstances with all three mechanisms. 153df4397b0SStefano Zampini 154df4397b0SStefano Zampini .seealso: PETSC_TRUE, PETSC_FALSE, PetscNot() 155df4397b0SStefano Zampini E*/ 156df4397b0SStefano Zampini typedef enum { PETSC_FALSE,PETSC_TRUE } PetscBool; 157df4397b0SStefano Zampini 158df4397b0SStefano Zampini /*MC 159df4397b0SStefano Zampini PetscReal - PETSc type that represents a real number version of PetscScalar 160df4397b0SStefano Zampini 161df4397b0SStefano Zampini 162df4397b0SStefano Zampini Notes: 163df4397b0SStefano Zampini For MPI calls that require datatypes, use MPIU_REAL as the datatype for PetscScalar and MPIU_SUM, MPIU_MAX, etc. for operations. 164df4397b0SStefano Zampini They will automatically work correctly regardless of the size of PetscReal. 165df4397b0SStefano Zampini 166df4397b0SStefano Zampini See PetscScalar for details on how to ./configure the size of PetscReal. 167df4397b0SStefano Zampini 168df4397b0SStefano Zampini Level: beginner 169df4397b0SStefano Zampini 170df4397b0SStefano Zampini .seealso: PetscScalar, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT 171df4397b0SStefano Zampini M*/ 172df4397b0SStefano Zampini 173df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 174df4397b0SStefano Zampini typedef float PetscReal; 175df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL_DOUBLE) 176df4397b0SStefano Zampini typedef double PetscReal; 177df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 178df4397b0SStefano Zampini # if defined(__cplusplus) 179df4397b0SStefano Zampini extern "C" { 180df4397b0SStefano Zampini # endif 181df4397b0SStefano Zampini # include <quadmath.h> 182df4397b0SStefano Zampini # if defined(__cplusplus) 183df4397b0SStefano Zampini } 184df4397b0SStefano Zampini # endif 185df4397b0SStefano Zampini typedef __float128 PetscReal; 186df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FP16) 187df4397b0SStefano Zampini typedef __fp16 PetscReal; 188df4397b0SStefano Zampini #endif /* PETSC_USE_REAL_* */ 189df4397b0SStefano Zampini 190df4397b0SStefano Zampini /*MC 191df4397b0SStefano Zampini PetscComplex - PETSc type that represents a complex number with precision matching that of PetscReal. 192df4397b0SStefano Zampini 193df4397b0SStefano Zampini Synopsis: 194df4397b0SStefano Zampini #include <petscsys.h> 195df4397b0SStefano Zampini PetscComplex number = 1. + 2.*PETSC_i; 196df4397b0SStefano Zampini 197df4397b0SStefano Zampini Notes: 198df4397b0SStefano Zampini For MPI calls that require datatypes, use MPIU_COMPLEX as the datatype for PetscComplex and MPIU_SUM etc for operations. 199df4397b0SStefano Zampini They will automatically work correctly regardless of the size of PetscComplex. 200df4397b0SStefano Zampini 201df4397b0SStefano Zampini See PetscScalar for details on how to ./configure the size of PetscReal 202df4397b0SStefano Zampini 203df4397b0SStefano Zampini Complex numbers are automatically available if PETSc was able to find a working complex implementation 204df4397b0SStefano Zampini 205df4397b0SStefano Zampini Level: beginner 206df4397b0SStefano Zampini 207df4397b0SStefano Zampini .seealso: PetscReal, PetscScalar, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT, PETSC_i 208df4397b0SStefano Zampini M*/ 209df4397b0SStefano Zampini 210df4397b0SStefano Zampini #if defined(__cplusplus) && defined(PETSC_HAVE_CXX_COMPLEX) && !defined(PETSC_USE_REAL___FLOAT128) 211df4397b0SStefano Zampini # if !defined(PETSC_SKIP_COMPLEX) 212df4397b0SStefano Zampini /* C++ support of complex number */ 213df4397b0SStefano Zampini # define PETSC_HAVE_COMPLEX 1 214df4397b0SStefano Zampini # if defined(PETSC_HAVE_CUDA) && __CUDACC_VER_MAJOR__ > 6 215df4397b0SStefano Zampini /* complex headers in thrust only available in CUDA 7.0 and above */ 216df4397b0SStefano Zampini # define petsccomplexlib thrust 217df4397b0SStefano Zampini # include <thrust/complex.h> 218df4397b0SStefano Zampini # else 219df4397b0SStefano Zampini # define petsccomplexlib std 220df4397b0SStefano Zampini # include <complex> 221df4397b0SStefano Zampini # endif 222df4397b0SStefano Zampini # if defined(PETSC_USE_REAL_SINGLE) 223df4397b0SStefano Zampini typedef petsccomplexlib::complex<float> PetscComplex; 224df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL_DOUBLE) 225df4397b0SStefano Zampini typedef petsccomplexlib::complex<double> PetscComplex; 226df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FLOAT128) 227df4397b0SStefano Zampini typedef petsccomplexlib::complex<__float128> PetscComplex; /* Notstandard and not expected to work, use __complex128 */ 228df4397b0SStefano Zampini # endif /* PETSC_USE_REAL_ */ 229df4397b0SStefano Zampini # endif /* ! PETSC_SKIP_COMPLEX */ 230df4397b0SStefano Zampini #elif defined(PETSC_HAVE_C99_COMPLEX) && !defined(PETSC_USE_REAL___FP16) 231df4397b0SStefano Zampini # if !defined(PETSC_SKIP_COMPLEX) 232df4397b0SStefano Zampini # define PETSC_HAVE_COMPLEX 1 233df4397b0SStefano Zampini # include <complex.h> 234df4397b0SStefano Zampini # if defined(PETSC_USE_REAL_SINGLE) || defined(PETSC_USE_REAL___FP16) 235df4397b0SStefano Zampini typedef float _Complex PetscComplex; 236df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL_DOUBLE) 237df4397b0SStefano Zampini typedef double _Complex PetscComplex; 238df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FLOAT128) 239df4397b0SStefano Zampini typedef __complex128 PetscComplex; 240df4397b0SStefano Zampini # endif /* PETSC_USE_REAL_* */ 241df4397b0SStefano Zampini # endif /* !PETSC_SKIP_COMPLEX */ 242df4397b0SStefano Zampini #elif (defined(PETSC_USE_COMPLEX) && !defined(PETSC_SKIP_COMPLEX)) 243df4397b0SStefano Zampini # error "PETSc was configured --with-scalar-type=complex, but a language-appropriate complex library is not available" 244df4397b0SStefano Zampini #endif /* !PETSC_SKIP_COMPLEX */ 245df4397b0SStefano Zampini 246df4397b0SStefano Zampini /*MC 247df4397b0SStefano Zampini PetscScalar - PETSc type that represents either a double precision real number, a double precision 248df4397b0SStefano Zampini complex number, a single precision real number, a __float128 real or complex or a __fp16 real - if the code is configured 249df4397b0SStefano Zampini with --with-scalar-type=real,complex --with-precision=single,double,__float128,__fp16 250df4397b0SStefano Zampini 251df4397b0SStefano Zampini Notes: 252df4397b0SStefano Zampini For MPI calls that require datatypes, use MPIU_SCALAR as the datatype for PetscScalar and MPIU_SUM, MPIU_MAX etc for operations. They will automatically work correctly regardless of the size of PetscScalar. 253df4397b0SStefano Zampini 254df4397b0SStefano Zampini Level: beginner 255df4397b0SStefano Zampini 256df4397b0SStefano Zampini .seealso: PetscReal, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT, PetscRealPart(), PetscImaginaryPart() 257df4397b0SStefano Zampini M*/ 258df4397b0SStefano Zampini 259df4397b0SStefano Zampini #if (defined(PETSC_USE_COMPLEX) && !defined(PETSC_SKIP_COMPLEX)) 260df4397b0SStefano Zampini typedef PetscComplex PetscScalar; 261df4397b0SStefano Zampini #else /* PETSC_USE_COMPLEX */ 262df4397b0SStefano Zampini typedef PetscReal PetscScalar; 263df4397b0SStefano Zampini #endif /* PETSC_USE_COMPLEX */ 264df4397b0SStefano Zampini 265df4397b0SStefano Zampini /*E 266df4397b0SStefano Zampini PetscCopyMode - Determines how an array passed to certain functions is copied or retained 267df4397b0SStefano Zampini 268df4397b0SStefano Zampini Level: beginner 269df4397b0SStefano Zampini 270df4397b0SStefano Zampini $ PETSC_COPY_VALUES - the array values are copied into new space, the user is free to reuse or delete the passed in array 271df4397b0SStefano Zampini $ PETSC_OWN_POINTER - the array values are NOT copied, the object takes ownership of the array and will free it later, the user cannot change or 272df4397b0SStefano Zampini $ delete the array. The array MUST have been obtained with PetscMalloc(). Hence this mode cannot be used in Fortran. 273df4397b0SStefano Zampini $ PETSC_USE_POINTER - the array values are NOT copied, the object uses the array but does NOT take ownership of the array. The user cannot use 274df4397b0SStefano Zampini the array but the user must delete the array after the object is destroyed. 275df4397b0SStefano Zampini 276df4397b0SStefano Zampini E*/ 277df4397b0SStefano Zampini typedef enum {PETSC_COPY_VALUES, PETSC_OWN_POINTER, PETSC_USE_POINTER} PetscCopyMode; 278df4397b0SStefano Zampini 279df4397b0SStefano Zampini /*MC 280df4397b0SStefano Zampini PETSC_FALSE - False value of PetscBool 281df4397b0SStefano Zampini 282df4397b0SStefano Zampini Level: beginner 283df4397b0SStefano Zampini 284df4397b0SStefano Zampini Note: 285df4397b0SStefano Zampini Zero integer 286df4397b0SStefano Zampini 287df4397b0SStefano Zampini .seealso: PetscBool, PETSC_TRUE 288df4397b0SStefano Zampini M*/ 289df4397b0SStefano Zampini 290df4397b0SStefano Zampini /*MC 291df4397b0SStefano Zampini PETSC_TRUE - True value of PetscBool 292df4397b0SStefano Zampini 293df4397b0SStefano Zampini Level: beginner 294df4397b0SStefano Zampini 295df4397b0SStefano Zampini Note: 296df4397b0SStefano Zampini Nonzero integer 297df4397b0SStefano Zampini 298df4397b0SStefano Zampini .seealso: PetscBool, PETSC_FALSE 299df4397b0SStefano Zampini M*/ 300df4397b0SStefano Zampini 301df4397b0SStefano Zampini /*MC 302df4397b0SStefano Zampini PetscLogDouble - Used for logging times 303df4397b0SStefano Zampini 304df4397b0SStefano Zampini Notes: 305df4397b0SStefano Zampini Contains double precision numbers that are not used in the numerical computations, but rather in logging, timing etc. 306df4397b0SStefano Zampini 307*6b3bf505SPatrick Sanan Level: developer 308*6b3bf505SPatrick Sanan 309df4397b0SStefano Zampini M*/ 310df4397b0SStefano Zampini typedef double PetscLogDouble; 311df4397b0SStefano Zampini 312df4397b0SStefano Zampini /*E 313df4397b0SStefano Zampini PetscDataType - Used for handling different basic data types. 314df4397b0SStefano Zampini 315df4397b0SStefano Zampini Level: beginner 316df4397b0SStefano Zampini 317df4397b0SStefano Zampini Notes: 318df4397b0SStefano Zampini Use of this should be avoided if one can directly use MPI_Datatype instead. 319df4397b0SStefano Zampini 320df4397b0SStefano Zampini Developer comment: 321df4397b0SStefano Zampini It would be nice if we could always just use MPI Datatypes, why can we not? 322df4397b0SStefano Zampini 323df4397b0SStefano Zampini .seealso: PetscBinaryRead(), PetscBinaryWrite(), PetscDataTypeToMPIDataType(), 324df4397b0SStefano Zampini PetscDataTypeGetSize() 325df4397b0SStefano Zampini 326df4397b0SStefano Zampini E*/ 3272ebc710fSVaclav Hapla typedef enum {PETSC_DATATYPE_UNKNOWN = 0, 3282ebc710fSVaclav Hapla PETSC_DOUBLE = 1, PETSC_COMPLEX = 2, PETSC_LONG = 3, PETSC_SHORT = 4, PETSC_FLOAT = 5, 32949de050aSVaclav Hapla PETSC_CHAR = 6, PETSC_BIT_LOGICAL = 7, PETSC_ENUM = 8, PETSC_BOOL = 9, PETSC___FLOAT128 = 10, 330ea127267SVaclav Hapla PETSC_OBJECT = 11, PETSC_FUNCTION = 12, PETSC_STRING = 13, PETSC___FP16 = 14, PETSC_STRUCT = 15, 3312ebc710fSVaclav Hapla PETSC_INT = 16} PetscDataType; 332df4397b0SStefano Zampini 333df4397b0SStefano Zampini #if defined(PETSC_USE_COMPLEX) 334df4397b0SStefano Zampini # define PETSC_SCALAR PETSC_COMPLEX 335df4397b0SStefano Zampini #else 336df4397b0SStefano Zampini # if defined(PETSC_USE_REAL_SINGLE) 337df4397b0SStefano Zampini # define PETSC_SCALAR PETSC_FLOAT 338df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FLOAT128) 339df4397b0SStefano Zampini # define PETSC_SCALAR PETSC___FLOAT128 340df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FP16) 341df4397b0SStefano Zampini # define PETSC_SCALAR PETSC___FP16 342df4397b0SStefano Zampini # else 343df4397b0SStefano Zampini # define PETSC_SCALAR PETSC_DOUBLE 344df4397b0SStefano Zampini # endif 345df4397b0SStefano Zampini #endif 346df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 347df4397b0SStefano Zampini # define PETSC_REAL PETSC_FLOAT 348df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 349df4397b0SStefano Zampini # define PETSC_REAL PETSC___FLOAT128 350df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FP16) 351df4397b0SStefano Zampini # define PETSC_REAL PETSC___FP16 352df4397b0SStefano Zampini #else 353df4397b0SStefano Zampini # define PETSC_REAL PETSC_DOUBLE 354df4397b0SStefano Zampini #endif 355df4397b0SStefano Zampini #define PETSC_FORTRANADDR PETSC_LONG 356df4397b0SStefano Zampini 357df4397b0SStefano Zampini /*S 358df4397b0SStefano Zampini PetscToken - 'Token' used for managing tokenizing strings 359df4397b0SStefano Zampini 360df4397b0SStefano Zampini Level: intermediate 361df4397b0SStefano Zampini 362df4397b0SStefano Zampini .seealso: PetscTokenCreate(), PetscTokenFind(), PetscTokenDestroy() 363df4397b0SStefano Zampini S*/ 364df4397b0SStefano Zampini typedef struct _p_PetscToken* PetscToken; 365df4397b0SStefano Zampini 366df4397b0SStefano Zampini /*S 367df4397b0SStefano Zampini PetscObject - any PETSc object, PetscViewer, Mat, Vec, KSP etc 368df4397b0SStefano Zampini 369df4397b0SStefano Zampini Level: beginner 370df4397b0SStefano Zampini 371df4397b0SStefano Zampini Note: 372df4397b0SStefano Zampini This is the base class from which all PETSc objects are derived from. 373df4397b0SStefano Zampini 374df4397b0SStefano Zampini .seealso: PetscObjectDestroy(), PetscObjectView(), PetscObjectGetName(), PetscObjectSetName(), PetscObjectReference(), PetscObjectDereference() 375df4397b0SStefano Zampini S*/ 376df4397b0SStefano Zampini typedef struct _p_PetscObject* PetscObject; 377df4397b0SStefano Zampini 378df4397b0SStefano Zampini /*MC 379df4397b0SStefano Zampini PetscObjectId - unique integer Id for a PetscObject 380df4397b0SStefano Zampini 381df4397b0SStefano Zampini Level: developer 382df4397b0SStefano Zampini 383df4397b0SStefano Zampini Notes: 384df4397b0SStefano Zampini Unlike pointer values, object ids are never reused. 385df4397b0SStefano Zampini 386df4397b0SStefano Zampini .seealso: PetscObjectState, PetscObjectGetId() 387df4397b0SStefano Zampini M*/ 388df4397b0SStefano Zampini #if defined(PETSC_USING_F90) && !defined(PETSC_USE_FORTRANKIND) /* compaq F90 */ 389df4397b0SStefano Zampini typedef int PetscObjectId; 390df4397b0SStefano Zampini #else 391df4397b0SStefano Zampini typedef PetscInt64 PetscObjectId; 392df4397b0SStefano Zampini #endif 393df4397b0SStefano Zampini 394df4397b0SStefano Zampini /*MC 395df4397b0SStefano Zampini PetscObjectState - integer state for a PetscObject 396df4397b0SStefano Zampini 397df4397b0SStefano Zampini Level: developer 398df4397b0SStefano Zampini 399df4397b0SStefano Zampini Notes: 400df4397b0SStefano Zampini Object state is always-increasing and (for objects that track state) can be used to determine if an object has 401df4397b0SStefano Zampini changed since the last time you interacted with it. It is 64-bit so that it will not overflow for a very long time. 402df4397b0SStefano Zampini 403df4397b0SStefano Zampini .seealso: PetscObjectId, PetscObjectStateGet(), PetscObjectStateIncrease(), PetscObjectStateSet() 404df4397b0SStefano Zampini M*/ 405df4397b0SStefano Zampini #if defined(PETSC_USING_F90) && !defined(PETSC_USE_FORTRANKIND) /* compaq F90 */ 406df4397b0SStefano Zampini typedef int PetscObjectState; 407df4397b0SStefano Zampini #else 408df4397b0SStefano Zampini typedef PetscInt64 PetscObjectState; 409df4397b0SStefano Zampini #endif 410df4397b0SStefano Zampini 411df4397b0SStefano Zampini /*S 412df4397b0SStefano Zampini PetscFunctionList - Linked list of functions, possibly stored in dynamic libraries, accessed 413df4397b0SStefano Zampini by string name 414df4397b0SStefano Zampini 415df4397b0SStefano Zampini Level: advanced 416df4397b0SStefano Zampini 417df4397b0SStefano Zampini .seealso: PetscFunctionListAdd(), PetscFunctionListDestroy(), PetscOpFlist 418df4397b0SStefano Zampini S*/ 419df4397b0SStefano Zampini typedef struct _n_PetscFunctionList *PetscFunctionList; 420df4397b0SStefano Zampini 421df4397b0SStefano Zampini /*E 422df4397b0SStefano Zampini PetscFileMode - Access mode for a file. 423df4397b0SStefano Zampini 424df4397b0SStefano Zampini Level: beginner 425df4397b0SStefano Zampini 426df4397b0SStefano Zampini $ FILE_MODE_READ - open a file at its beginning for reading 427df4397b0SStefano Zampini $ FILE_MODE_WRITE - open a file at its beginning for writing (will create if the file does not exist) 428df4397b0SStefano Zampini $ FILE_MODE_APPEND - open a file at end for writing 429df4397b0SStefano Zampini $ FILE_MODE_UPDATE - open a file for updating, meaning for reading and writing 430df4397b0SStefano Zampini $ FILE_MODE_APPEND_UPDATE - open a file for updating, meaning for reading and writing, at the end 431df4397b0SStefano Zampini 432df4397b0SStefano Zampini .seealso: PetscViewerFileSetMode() 433df4397b0SStefano Zampini E*/ 434df4397b0SStefano Zampini typedef enum {FILE_MODE_READ, FILE_MODE_WRITE, FILE_MODE_APPEND, FILE_MODE_UPDATE, FILE_MODE_APPEND_UPDATE} PetscFileMode; 435df4397b0SStefano Zampini 436df4397b0SStefano Zampini typedef void* PetscDLHandle; 437df4397b0SStefano Zampini typedef enum {PETSC_DL_DECIDE=0,PETSC_DL_NOW=1,PETSC_DL_LOCAL=2} PetscDLMode; 438df4397b0SStefano Zampini 439df4397b0SStefano Zampini /*S 440df4397b0SStefano Zampini PetscObjectList - Linked list of PETSc objects, each accessable by string name 441df4397b0SStefano Zampini 442df4397b0SStefano Zampini Level: developer 443df4397b0SStefano Zampini 444df4397b0SStefano Zampini Notes: 445df4397b0SStefano Zampini Used by PetscObjectCompose() and PetscObjectQuery() 446df4397b0SStefano Zampini 447df4397b0SStefano Zampini .seealso: PetscObjectListAdd(), PetscObjectListDestroy(), PetscObjectListFind(), PetscObjectCompose(), PetscObjectQuery(), PetscFunctionList 448df4397b0SStefano Zampini S*/ 449df4397b0SStefano Zampini typedef struct _n_PetscObjectList *PetscObjectList; 450df4397b0SStefano Zampini 451df4397b0SStefano Zampini /*S 452df4397b0SStefano Zampini PetscDLLibrary - Linked list of dynamics libraries to search for functions 453df4397b0SStefano Zampini 454df4397b0SStefano Zampini Level: advanced 455df4397b0SStefano Zampini 456df4397b0SStefano Zampini .seealso: PetscDLLibraryOpen() 457df4397b0SStefano Zampini S*/ 458df4397b0SStefano Zampini typedef struct _n_PetscDLLibrary *PetscDLLibrary; 459df4397b0SStefano Zampini 460df4397b0SStefano Zampini /*S 461df4397b0SStefano Zampini PetscContainer - Simple PETSc object that contains a pointer to any required data 462df4397b0SStefano Zampini 463df4397b0SStefano Zampini Level: advanced 464df4397b0SStefano Zampini 465df4397b0SStefano Zampini .seealso: PetscObject, PetscContainerCreate() 466df4397b0SStefano Zampini S*/ 467df4397b0SStefano Zampini typedef struct _p_PetscContainer* PetscContainer; 468df4397b0SStefano Zampini 469df4397b0SStefano Zampini /*S 470df4397b0SStefano Zampini PetscRandom - Abstract PETSc object that manages generating random numbers 471df4397b0SStefano Zampini 472df4397b0SStefano Zampini Level: intermediate 473df4397b0SStefano Zampini 474df4397b0SStefano Zampini Concepts: random numbers 475df4397b0SStefano Zampini 476df4397b0SStefano Zampini .seealso: PetscRandomCreate(), PetscRandomGetValue(), PetscRandomType 477df4397b0SStefano Zampini S*/ 478df4397b0SStefano Zampini typedef struct _p_PetscRandom* PetscRandom; 479df4397b0SStefano Zampini 480df4397b0SStefano Zampini /* 481df4397b0SStefano Zampini In binary files variables are stored using the following lengths, 482df4397b0SStefano Zampini regardless of how they are stored in memory on any one particular 483df4397b0SStefano Zampini machine. Use these rather then sizeof() in computing sizes for 484df4397b0SStefano Zampini PetscBinarySeek(). 485df4397b0SStefano Zampini */ 486df4397b0SStefano Zampini #define PETSC_BINARY_INT_SIZE (32/8) 487df4397b0SStefano Zampini #define PETSC_BINARY_FLOAT_SIZE (32/8) 488df4397b0SStefano Zampini #define PETSC_BINARY_CHAR_SIZE (8/8) 489df4397b0SStefano Zampini #define PETSC_BINARY_SHORT_SIZE (16/8) 490df4397b0SStefano Zampini #define PETSC_BINARY_DOUBLE_SIZE (64/8) 491df4397b0SStefano Zampini #define PETSC_BINARY_SCALAR_SIZE sizeof(PetscScalar) 492df4397b0SStefano Zampini 493df4397b0SStefano Zampini /*E 494df4397b0SStefano Zampini PetscBinarySeekType - argument to PetscBinarySeek() 495df4397b0SStefano Zampini 496df4397b0SStefano Zampini Level: advanced 497df4397b0SStefano Zampini 498df4397b0SStefano Zampini .seealso: PetscBinarySeek(), PetscBinarySynchronizedSeek() 499df4397b0SStefano Zampini E*/ 500df4397b0SStefano Zampini typedef enum {PETSC_BINARY_SEEK_SET = 0,PETSC_BINARY_SEEK_CUR = 1,PETSC_BINARY_SEEK_END = 2} PetscBinarySeekType; 501df4397b0SStefano Zampini 502df4397b0SStefano Zampini /*E 503df4397b0SStefano Zampini PetscBuildTwoSidedType - algorithm for setting up two-sided communication 504df4397b0SStefano Zampini 505df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_ALLREDUCE - classical algorithm using an MPI_Allreduce with 506df4397b0SStefano Zampini $ a buffer of length equal to the communicator size. Not memory-scalable due to 507df4397b0SStefano Zampini $ the large reduction size. Requires only MPI-1. 508df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_IBARRIER - nonblocking algorithm based on MPI_Issend and MPI_Ibarrier. 509df4397b0SStefano Zampini $ Proved communication-optimal in Hoefler, Siebert, and Lumsdaine (2010). Requires MPI-3. 510df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_REDSCATTER - similar to above, but use more optimized function 511df4397b0SStefano Zampini $ that only communicates the part of the reduction that is necessary. Requires MPI-2. 512df4397b0SStefano Zampini 513df4397b0SStefano Zampini Level: developer 514df4397b0SStefano Zampini 515df4397b0SStefano Zampini .seealso: PetscCommBuildTwoSided(), PetscCommBuildTwoSidedSetType(), PetscCommBuildTwoSidedGetType() 516df4397b0SStefano Zampini E*/ 517df4397b0SStefano Zampini typedef enum { 518df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_NOTSET = -1, 519df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_ALLREDUCE = 0, 520df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_IBARRIER = 1, 521df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_REDSCATTER = 2 522df4397b0SStefano Zampini /* Updates here must be accompanied by updates in finclude/petscsys.h and the string array in mpits.c */ 523df4397b0SStefano Zampini } PetscBuildTwoSidedType; 524df4397b0SStefano Zampini 525df4397b0SStefano Zampini /*E 526df4397b0SStefano Zampini InsertMode - Whether entries are inserted or added into vectors or matrices 527df4397b0SStefano Zampini 528df4397b0SStefano Zampini Level: beginner 529df4397b0SStefano Zampini 530df4397b0SStefano Zampini .seealso: VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 531df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), 532df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd() 533df4397b0SStefano Zampini E*/ 534df4397b0SStefano Zampini typedef enum {NOT_SET_VALUES, INSERT_VALUES, ADD_VALUES, MAX_VALUES, INSERT_ALL_VALUES, ADD_ALL_VALUES, INSERT_BC_VALUES, ADD_BC_VALUES} InsertMode; 535df4397b0SStefano Zampini 536df4397b0SStefano Zampini /*MC 537df4397b0SStefano Zampini INSERT_VALUES - Put a value into a vector or matrix, overwrites any previous value 538df4397b0SStefano Zampini 539df4397b0SStefano Zampini Level: beginner 540df4397b0SStefano Zampini 541df4397b0SStefano Zampini .seealso: InsertMode, VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 542df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), ADD_VALUES, 543df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd(), MAX_VALUES 544df4397b0SStefano Zampini 545df4397b0SStefano Zampini M*/ 546df4397b0SStefano Zampini 547df4397b0SStefano Zampini /*MC 548df4397b0SStefano Zampini ADD_VALUES - Adds a value into a vector or matrix, if there previously was no value, just puts the 549df4397b0SStefano Zampini value into that location 550df4397b0SStefano Zampini 551df4397b0SStefano Zampini Level: beginner 552df4397b0SStefano Zampini 553df4397b0SStefano Zampini .seealso: InsertMode, VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 554df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), INSERT_VALUES, 555df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd(), MAX_VALUES 556df4397b0SStefano Zampini 557df4397b0SStefano Zampini M*/ 558df4397b0SStefano Zampini 559df4397b0SStefano Zampini /*MC 560df4397b0SStefano Zampini MAX_VALUES - Puts the maximum of the scattered/gathered value and the current value into each location 561df4397b0SStefano Zampini 562df4397b0SStefano Zampini Level: beginner 563df4397b0SStefano Zampini 564df4397b0SStefano Zampini .seealso: InsertMode, VecScatterBegin(), VecScatterEnd(), ADD_VALUES, INSERT_VALUES 565df4397b0SStefano Zampini 566df4397b0SStefano Zampini M*/ 567df4397b0SStefano Zampini 568df4397b0SStefano Zampini /*S 569df4397b0SStefano Zampini PetscSubcomm - A decomposition of an MPI communicator into subcommunicators 570df4397b0SStefano Zampini 571df4397b0SStefano Zampini Notes: 572df4397b0SStefano Zampini After a call to PetscSubcommSetType(), PetscSubcommSetTypeGeneral(), or PetscSubcommSetFromOptions() one may call 573df4397b0SStefano Zampini $ PetscSubcommChild() returns the associated subcommunicator on this process 574df4397b0SStefano Zampini $ PetscSubcommContiguousParent() returns a parent communitor but with all child of the same subcommunicator having contiguous rank 575df4397b0SStefano Zampini 576df4397b0SStefano Zampini Sample Usage: 577df4397b0SStefano Zampini PetscSubcommCreate() 578df4397b0SStefano Zampini PetscSubcommSetNumber() 579df4397b0SStefano Zampini PetscSubcommSetType(PETSC_SUBCOMM_INTERLACED); 580df4397b0SStefano Zampini ccomm = PetscSubcommChild() 581df4397b0SStefano Zampini PetscSubcommDestroy() 582df4397b0SStefano Zampini 583df4397b0SStefano Zampini Level: advanced 584df4397b0SStefano Zampini 585df4397b0SStefano Zampini Concepts: communicator, create 586df4397b0SStefano Zampini 587df4397b0SStefano Zampini Notes: 588df4397b0SStefano Zampini $ PETSC_SUBCOMM_GENERAL - similar to MPI_Comm_split() each process sets the new communicator (color) they will belong to and the order within that communicator 589df4397b0SStefano Zampini $ PETSC_SUBCOMM_CONTIGUOUS - each new communicator contains a set of process with contiguous ranks in the original MPI communicator 590df4397b0SStefano Zampini $ PETSC_SUBCOMM_INTERLACED - each new communictor contains a set of processes equally far apart in rank from the others in that new communicator 591df4397b0SStefano Zampini 592df4397b0SStefano Zampini Example: Consider a communicator with six processes split into 3 subcommunicators. 593df4397b0SStefano Zampini $ PETSC_SUBCOMM_CONTIGUOUS - the first communicator contains rank 0,1 the second rank 2,3 and the third rank 4,5 in the original ordering of the original communicator 594df4397b0SStefano Zampini $ PETSC_SUBCOMM_INTERLACED - the first communicator contains rank 0,3, the second 1,4 and the third 2,5 595df4397b0SStefano Zampini 596df4397b0SStefano Zampini Developer Notes: 597df4397b0SStefano Zampini This is used in objects such as PCREDUNDANT to manage the subcommunicators on which the redundant computations 598df4397b0SStefano Zampini are performed. 599df4397b0SStefano Zampini 600df4397b0SStefano Zampini 601df4397b0SStefano Zampini .seealso: PetscSubcommCreate(), PetscSubcommSetNumber(), PetscSubcommSetType(), PetscSubcommView(), PetscSubcommSetFromOptions() 602df4397b0SStefano Zampini 603df4397b0SStefano Zampini S*/ 604df4397b0SStefano Zampini typedef struct _n_PetscSubcomm* PetscSubcomm; 605df4397b0SStefano Zampini typedef enum {PETSC_SUBCOMM_GENERAL=0,PETSC_SUBCOMM_CONTIGUOUS=1,PETSC_SUBCOMM_INTERLACED=2} PetscSubcommType; 606df4397b0SStefano Zampini 607df4397b0SStefano Zampini /*S 608df4397b0SStefano Zampini PetscHeap - A simple class for managing heaps 609df4397b0SStefano Zampini 610df4397b0SStefano Zampini Level: intermediate 611df4397b0SStefano Zampini 612df4397b0SStefano Zampini Concepts: random numbers 613df4397b0SStefano Zampini 614df4397b0SStefano Zampini .seealso: PetscHeapCreate(), PetscHeapAdd(), PetscHeapPop(), PetscHeapPeek(), PetscHeapStash(), PetscHeapUnstash(), PetscHeapView(), PetscHeapDestroy() 615df4397b0SStefano Zampini S*/ 616df4397b0SStefano Zampini typedef struct _PetscHeap *PetscHeap; 617df4397b0SStefano Zampini 618df4397b0SStefano Zampini typedef struct _n_PetscShmComm* PetscShmComm; 619df4397b0SStefano Zampini typedef struct _n_PetscOmpCtrl* PetscOmpCtrl; 620df4397b0SStefano Zampini 621df4397b0SStefano Zampini /*S 622df4397b0SStefano Zampini PetscSegBuffer - a segmented extendable buffer 623df4397b0SStefano Zampini 624df4397b0SStefano Zampini Level: developer 625df4397b0SStefano Zampini 626df4397b0SStefano Zampini .seealso: PetscSegBufferCreate(), PetscSegBufferGet(), PetscSegBufferExtract(), PetscSegBufferDestroy() 627df4397b0SStefano Zampini S*/ 628df4397b0SStefano Zampini typedef struct _n_PetscSegBuffer *PetscSegBuffer; 629df4397b0SStefano Zampini 630df4397b0SStefano Zampini typedef struct _n_PetscOptionsHelpPrinted *PetscOptionsHelpPrinted; 631df4397b0SStefano Zampini 632df4397b0SStefano Zampini #endif 633