126bd1501SBarry Smith #if !defined(PETSCSYSTYPES_H) 226bd1501SBarry Smith #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 111cdc6ee08SBarry Smith (except on very rare BLAS/LAPACK implementations that support 64 bit integers see the notes 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: 117cdc6ee08SBarry Smith ./configure automatically determines the size of the integers used by BLAS/LAPACK except when --with-batch is used 118cdc6ee08SBarry Smith in that situation one must know (by some other means) if the integers used by BLAS/LAPACK are 64 bit and if so pass the flag --known-64-bit-blas-indice 119cdc6ee08SBarry Smith 120cdc6ee08SBarry Smith MATLAB ships with BLAS and LAPACK that use 64 bit integers, for example if you run ./configure with, the option 121df4397b0SStefano Zampini --with-blaslapack-lib=[/Applications/MATLAB_R2010b.app/bin/maci64/libmwblas.dylib,/Applications/MATLAB_R2010b.app/bin/maci64/libmwlapack.dylib] 122df4397b0SStefano Zampini 123cdc6ee08SBarry Smith MKL ships with both 32 and 64 bit integer versions of the BLAS and LAPACK. If you pass the flag -with-64-bit-blas-indices PETSc will link 124cdc6ee08SBarry Smith against the 64 bit version, otherwise it use the 32 bit version 125df4397b0SStefano Zampini 126cdc6ee08SBarry Smith OpenBLAS can be built to use 64 bit integers. The ./configure options --download-openblas -with-64-bit-blas-indices will build a 64 bit integer version 127df4397b0SStefano Zampini 128df4397b0SStefano Zampini External packages such as hypre, ML, SuperLU etc do not provide any support for passing 64 bit integers to BLAS/LAPACK so cannot 129cdc6ee08SBarry Smith be used with PETSc when PETSc links against 64 bit integer BLAS/LAPACK. ./configure will generate an error if you attempt to link PETSc against any of 130cdc6ee08SBarry Smith these external libraries while using 64 bit integer BLAS/LAPACK. 131df4397b0SStefano Zampini 132df4397b0SStefano Zampini Level: intermediate 133df4397b0SStefano Zampini 134df4397b0SStefano Zampini .seealso: PetscMPIInt, PetscInt, PetscBLASIntCast() 135df4397b0SStefano Zampini 136df4397b0SStefano Zampini M*/ 137df4397b0SStefano Zampini #if defined(PETSC_HAVE_64BIT_BLAS_INDICES) 138df4397b0SStefano Zampini typedef PetscInt64 PetscBLASInt; 139df4397b0SStefano Zampini #else 140df4397b0SStefano Zampini typedef int PetscBLASInt; 141df4397b0SStefano Zampini #endif 142df4397b0SStefano Zampini 143df4397b0SStefano Zampini /*E 144df4397b0SStefano Zampini PetscBool - Logical variable. Actually an int in C and a logical in Fortran. 145df4397b0SStefano Zampini 146df4397b0SStefano Zampini Level: beginner 147df4397b0SStefano Zampini 148df4397b0SStefano Zampini Developer Note: 149df4397b0SStefano 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 150df4397b0SStefano Zampini boolean values. It is not easy to have a simple macro that that will work properly in all circumstances with all three mechanisms. 151df4397b0SStefano Zampini 152df4397b0SStefano Zampini .seealso: PETSC_TRUE, PETSC_FALSE, PetscNot() 153df4397b0SStefano Zampini E*/ 154df4397b0SStefano Zampini typedef enum { PETSC_FALSE,PETSC_TRUE } PetscBool; 155df4397b0SStefano Zampini 156df4397b0SStefano Zampini /*MC 157df4397b0SStefano Zampini PetscReal - PETSc type that represents a real number version of PetscScalar 158df4397b0SStefano Zampini 159df4397b0SStefano Zampini 160df4397b0SStefano Zampini Notes: 161df4397b0SStefano Zampini For MPI calls that require datatypes, use MPIU_REAL as the datatype for PetscScalar and MPIU_SUM, MPIU_MAX, etc. for operations. 162df4397b0SStefano Zampini They will automatically work correctly regardless of the size of PetscReal. 163df4397b0SStefano Zampini 164df4397b0SStefano Zampini See PetscScalar for details on how to ./configure the size of PetscReal. 165df4397b0SStefano Zampini 166df4397b0SStefano Zampini Level: beginner 167df4397b0SStefano Zampini 168df4397b0SStefano Zampini .seealso: PetscScalar, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT 169df4397b0SStefano Zampini M*/ 170df4397b0SStefano Zampini 171df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 172df4397b0SStefano Zampini typedef float PetscReal; 173df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL_DOUBLE) 174df4397b0SStefano Zampini typedef double PetscReal; 175df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 176df4397b0SStefano Zampini # if defined(__cplusplus) 177df4397b0SStefano Zampini extern "C" { 178df4397b0SStefano Zampini # endif 179df4397b0SStefano Zampini # include <quadmath.h> 180df4397b0SStefano Zampini # if defined(__cplusplus) 181df4397b0SStefano Zampini } 182df4397b0SStefano Zampini # endif 183df4397b0SStefano Zampini typedef __float128 PetscReal; 184df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FP16) 185df4397b0SStefano Zampini typedef __fp16 PetscReal; 186df4397b0SStefano Zampini #endif /* PETSC_USE_REAL_* */ 187df4397b0SStefano Zampini 188df4397b0SStefano Zampini /*MC 189df4397b0SStefano Zampini PetscComplex - PETSc type that represents a complex number with precision matching that of PetscReal. 190df4397b0SStefano Zampini 191df4397b0SStefano Zampini Synopsis: 192df4397b0SStefano Zampini #include <petscsys.h> 193df4397b0SStefano Zampini PetscComplex number = 1. + 2.*PETSC_i; 194df4397b0SStefano Zampini 195df4397b0SStefano Zampini Notes: 196df4397b0SStefano Zampini For MPI calls that require datatypes, use MPIU_COMPLEX as the datatype for PetscComplex and MPIU_SUM etc for operations. 197df4397b0SStefano Zampini They will automatically work correctly regardless of the size of PetscComplex. 198df4397b0SStefano Zampini 199df4397b0SStefano Zampini See PetscScalar for details on how to ./configure the size of PetscReal 200df4397b0SStefano Zampini 201df4397b0SStefano Zampini Complex numbers are automatically available if PETSc was able to find a working complex implementation 202df4397b0SStefano Zampini 203df4397b0SStefano Zampini Level: beginner 204df4397b0SStefano Zampini 205df4397b0SStefano Zampini .seealso: PetscReal, PetscScalar, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT, PETSC_i 206df4397b0SStefano Zampini M*/ 207df4397b0SStefano Zampini 208df4397b0SStefano Zampini #if defined(__cplusplus) && defined(PETSC_HAVE_CXX_COMPLEX) && !defined(PETSC_USE_REAL___FLOAT128) 209df4397b0SStefano Zampini # if !defined(PETSC_SKIP_COMPLEX) 210df4397b0SStefano Zampini /* C++ support of complex number */ 211df4397b0SStefano Zampini # define PETSC_HAVE_COMPLEX 1 212df4397b0SStefano Zampini # if defined(PETSC_HAVE_CUDA) && __CUDACC_VER_MAJOR__ > 6 213df4397b0SStefano Zampini /* complex headers in thrust only available in CUDA 7.0 and above */ 214df4397b0SStefano Zampini # define petsccomplexlib thrust 215df4397b0SStefano Zampini # include <thrust/complex.h> 216df4397b0SStefano Zampini # else 217df4397b0SStefano Zampini # define petsccomplexlib std 218df4397b0SStefano Zampini # include <complex> 219df4397b0SStefano Zampini # endif 220df4397b0SStefano Zampini # if defined(PETSC_USE_REAL_SINGLE) 221df4397b0SStefano Zampini typedef petsccomplexlib::complex<float> PetscComplex; 222df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL_DOUBLE) 223df4397b0SStefano Zampini typedef petsccomplexlib::complex<double> PetscComplex; 224df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FLOAT128) 225df4397b0SStefano Zampini typedef petsccomplexlib::complex<__float128> PetscComplex; /* Notstandard and not expected to work, use __complex128 */ 226df4397b0SStefano Zampini # endif /* PETSC_USE_REAL_ */ 227df4397b0SStefano Zampini # endif /* ! PETSC_SKIP_COMPLEX */ 22839829747SLisandro Dalcin # if !defined(PETSC_SKIP_CXX_COMPLEX_FIX) 22939829747SLisandro Dalcin # include <petsccxxcomplexfix.h> 23039829747SLisandro Dalcin # endif /* ! PETSC_SKIP_CXX_COMPLEX_FIX */ 231df4397b0SStefano Zampini #elif defined(PETSC_HAVE_C99_COMPLEX) && !defined(PETSC_USE_REAL___FP16) 232df4397b0SStefano Zampini # if !defined(PETSC_SKIP_COMPLEX) 233df4397b0SStefano Zampini # define PETSC_HAVE_COMPLEX 1 234df4397b0SStefano Zampini # include <complex.h> 235df4397b0SStefano Zampini # if defined(PETSC_USE_REAL_SINGLE) || defined(PETSC_USE_REAL___FP16) 236df4397b0SStefano Zampini typedef float _Complex PetscComplex; 237df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL_DOUBLE) 238df4397b0SStefano Zampini typedef double _Complex PetscComplex; 239df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FLOAT128) 240df4397b0SStefano Zampini typedef __complex128 PetscComplex; 241df4397b0SStefano Zampini # endif /* PETSC_USE_REAL_* */ 242df4397b0SStefano Zampini # endif /* !PETSC_SKIP_COMPLEX */ 243df4397b0SStefano Zampini #elif (defined(PETSC_USE_COMPLEX) && !defined(PETSC_SKIP_COMPLEX)) 244df4397b0SStefano Zampini # error "PETSc was configured --with-scalar-type=complex, but a language-appropriate complex library is not available" 245df4397b0SStefano Zampini #endif /* !PETSC_SKIP_COMPLEX */ 246df4397b0SStefano Zampini 247df4397b0SStefano Zampini /*MC 248df4397b0SStefano Zampini PetscScalar - PETSc type that represents either a double precision real number, a double precision 249df4397b0SStefano Zampini complex number, a single precision real number, a __float128 real or complex or a __fp16 real - if the code is configured 250df4397b0SStefano Zampini with --with-scalar-type=real,complex --with-precision=single,double,__float128,__fp16 251df4397b0SStefano Zampini 252df4397b0SStefano Zampini Notes: 253df4397b0SStefano 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. 254df4397b0SStefano Zampini 255df4397b0SStefano Zampini Level: beginner 256df4397b0SStefano Zampini 257df4397b0SStefano Zampini .seealso: PetscReal, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT, PetscRealPart(), PetscImaginaryPart() 258df4397b0SStefano Zampini M*/ 259df4397b0SStefano Zampini 260df4397b0SStefano Zampini #if (defined(PETSC_USE_COMPLEX) && !defined(PETSC_SKIP_COMPLEX)) 261df4397b0SStefano Zampini typedef PetscComplex PetscScalar; 262df4397b0SStefano Zampini #else /* PETSC_USE_COMPLEX */ 263df4397b0SStefano Zampini typedef PetscReal PetscScalar; 264df4397b0SStefano Zampini #endif /* PETSC_USE_COMPLEX */ 265df4397b0SStefano Zampini 266df4397b0SStefano Zampini /*E 267df4397b0SStefano Zampini PetscCopyMode - Determines how an array passed to certain functions is copied or retained 268df4397b0SStefano Zampini 269df4397b0SStefano Zampini Level: beginner 270df4397b0SStefano Zampini 271df4397b0SStefano Zampini $ PETSC_COPY_VALUES - the array values are copied into new space, the user is free to reuse or delete the passed in array 272df4397b0SStefano 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 273df4397b0SStefano Zampini $ delete the array. The array MUST have been obtained with PetscMalloc(). Hence this mode cannot be used in Fortran. 274df4397b0SStefano 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 275df4397b0SStefano Zampini the array but the user must delete the array after the object is destroyed. 276df4397b0SStefano Zampini 277df4397b0SStefano Zampini E*/ 278df4397b0SStefano Zampini typedef enum {PETSC_COPY_VALUES, PETSC_OWN_POINTER, PETSC_USE_POINTER} PetscCopyMode; 279df4397b0SStefano Zampini 280df4397b0SStefano Zampini /*MC 281df4397b0SStefano Zampini PETSC_FALSE - False value of PetscBool 282df4397b0SStefano Zampini 283df4397b0SStefano Zampini Level: beginner 284df4397b0SStefano Zampini 285df4397b0SStefano Zampini Note: 286df4397b0SStefano Zampini Zero integer 287df4397b0SStefano Zampini 288df4397b0SStefano Zampini .seealso: PetscBool, PETSC_TRUE 289df4397b0SStefano Zampini M*/ 290df4397b0SStefano Zampini 291df4397b0SStefano Zampini /*MC 292df4397b0SStefano Zampini PETSC_TRUE - True value of PetscBool 293df4397b0SStefano Zampini 294df4397b0SStefano Zampini Level: beginner 295df4397b0SStefano Zampini 296df4397b0SStefano Zampini Note: 297df4397b0SStefano Zampini Nonzero integer 298df4397b0SStefano Zampini 299df4397b0SStefano Zampini .seealso: PetscBool, PETSC_FALSE 300df4397b0SStefano Zampini M*/ 301df4397b0SStefano Zampini 302df4397b0SStefano Zampini /*MC 303df4397b0SStefano Zampini PetscLogDouble - Used for logging times 304df4397b0SStefano Zampini 305df4397b0SStefano Zampini Notes: 306df4397b0SStefano Zampini Contains double precision numbers that are not used in the numerical computations, but rather in logging, timing etc. 307df4397b0SStefano Zampini 3086b3bf505SPatrick Sanan Level: developer 3096b3bf505SPatrick Sanan 310df4397b0SStefano Zampini M*/ 311df4397b0SStefano Zampini typedef double PetscLogDouble; 312df4397b0SStefano Zampini 313df4397b0SStefano Zampini /*E 314df4397b0SStefano Zampini PetscDataType - Used for handling different basic data types. 315df4397b0SStefano Zampini 316df4397b0SStefano Zampini Level: beginner 317df4397b0SStefano Zampini 318df4397b0SStefano Zampini Notes: 319df4397b0SStefano Zampini Use of this should be avoided if one can directly use MPI_Datatype instead. 320df4397b0SStefano Zampini 321df4397b0SStefano Zampini Developer comment: 322df4397b0SStefano Zampini It would be nice if we could always just use MPI Datatypes, why can we not? 323df4397b0SStefano Zampini 32491e29162SBarry Smith If you change any values in PetscDatatype make sure you update their usage in 32591e29162SBarry Smith share/petsc/matlab/PetscBagRead.m 32691e29162SBarry Smith 327df4397b0SStefano Zampini .seealso: PetscBinaryRead(), PetscBinaryWrite(), PetscDataTypeToMPIDataType(), 328df4397b0SStefano Zampini PetscDataTypeGetSize() 329df4397b0SStefano Zampini 330df4397b0SStefano Zampini E*/ 3312ebc710fSVaclav Hapla typedef enum {PETSC_DATATYPE_UNKNOWN = 0, 3322ebc710fSVaclav Hapla PETSC_DOUBLE = 1, PETSC_COMPLEX = 2, PETSC_LONG = 3, PETSC_SHORT = 4, PETSC_FLOAT = 5, 33349de050aSVaclav Hapla PETSC_CHAR = 6, PETSC_BIT_LOGICAL = 7, PETSC_ENUM = 8, PETSC_BOOL = 9, PETSC___FLOAT128 = 10, 334ea127267SVaclav Hapla PETSC_OBJECT = 11, PETSC_FUNCTION = 12, PETSC_STRING = 13, PETSC___FP16 = 14, PETSC_STRUCT = 15, 3359e3e4c22SLisandro Dalcin PETSC_INT = 16, PETSC_INT64 = 17} PetscDataType; 336df4397b0SStefano Zampini 337df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 338df4397b0SStefano Zampini # define PETSC_REAL PETSC_FLOAT 3395117d392SLisandro Dalcin #elif defined(PETSC_USE_REAL_DOUBLE) 3405117d392SLisandro Dalcin # define PETSC_REAL PETSC_DOUBLE 341df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 342df4397b0SStefano Zampini # define PETSC_REAL PETSC___FLOAT128 343df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FP16) 344df4397b0SStefano Zampini # define PETSC_REAL PETSC___FP16 345df4397b0SStefano Zampini #else 346df4397b0SStefano Zampini # define PETSC_REAL PETSC_DOUBLE 347df4397b0SStefano Zampini #endif 3485117d392SLisandro Dalcin 3495117d392SLisandro Dalcin #if defined(PETSC_USE_COMPLEX) 3505117d392SLisandro Dalcin # define PETSC_SCALAR PETSC_COMPLEX 3515117d392SLisandro Dalcin #else 3525117d392SLisandro Dalcin # define PETSC_SCALAR PETSC_REAL 3535117d392SLisandro Dalcin #endif 3545117d392SLisandro Dalcin 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 .seealso: PetscRandomCreate(), PetscRandomGetValue(), PetscRandomType 475df4397b0SStefano Zampini S*/ 476df4397b0SStefano Zampini typedef struct _p_PetscRandom* PetscRandom; 477df4397b0SStefano Zampini 478df4397b0SStefano Zampini /* 479df4397b0SStefano Zampini In binary files variables are stored using the following lengths, 480df4397b0SStefano Zampini regardless of how they are stored in memory on any one particular 481df4397b0SStefano Zampini machine. Use these rather then sizeof() in computing sizes for 482df4397b0SStefano Zampini PetscBinarySeek(). 483df4397b0SStefano Zampini */ 484df4397b0SStefano Zampini #define PETSC_BINARY_INT_SIZE (32/8) 485df4397b0SStefano Zampini #define PETSC_BINARY_FLOAT_SIZE (32/8) 486df4397b0SStefano Zampini #define PETSC_BINARY_CHAR_SIZE (8/8) 487df4397b0SStefano Zampini #define PETSC_BINARY_SHORT_SIZE (16/8) 488df4397b0SStefano Zampini #define PETSC_BINARY_DOUBLE_SIZE (64/8) 489df4397b0SStefano Zampini #define PETSC_BINARY_SCALAR_SIZE sizeof(PetscScalar) 490df4397b0SStefano Zampini 491df4397b0SStefano Zampini /*E 492df4397b0SStefano Zampini PetscBinarySeekType - argument to PetscBinarySeek() 493df4397b0SStefano Zampini 494df4397b0SStefano Zampini Level: advanced 495df4397b0SStefano Zampini 496df4397b0SStefano Zampini .seealso: PetscBinarySeek(), PetscBinarySynchronizedSeek() 497df4397b0SStefano Zampini E*/ 498df4397b0SStefano Zampini typedef enum {PETSC_BINARY_SEEK_SET = 0,PETSC_BINARY_SEEK_CUR = 1,PETSC_BINARY_SEEK_END = 2} PetscBinarySeekType; 499df4397b0SStefano Zampini 500df4397b0SStefano Zampini /*E 501df4397b0SStefano Zampini PetscBuildTwoSidedType - algorithm for setting up two-sided communication 502df4397b0SStefano Zampini 503df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_ALLREDUCE - classical algorithm using an MPI_Allreduce with 504df4397b0SStefano Zampini $ a buffer of length equal to the communicator size. Not memory-scalable due to 505df4397b0SStefano Zampini $ the large reduction size. Requires only MPI-1. 506df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_IBARRIER - nonblocking algorithm based on MPI_Issend and MPI_Ibarrier. 507df4397b0SStefano Zampini $ Proved communication-optimal in Hoefler, Siebert, and Lumsdaine (2010). Requires MPI-3. 508df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_REDSCATTER - similar to above, but use more optimized function 509df4397b0SStefano Zampini $ that only communicates the part of the reduction that is necessary. Requires MPI-2. 510df4397b0SStefano Zampini 511df4397b0SStefano Zampini Level: developer 512df4397b0SStefano Zampini 513df4397b0SStefano Zampini .seealso: PetscCommBuildTwoSided(), PetscCommBuildTwoSidedSetType(), PetscCommBuildTwoSidedGetType() 514df4397b0SStefano Zampini E*/ 515df4397b0SStefano Zampini typedef enum { 516df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_NOTSET = -1, 517df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_ALLREDUCE = 0, 518df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_IBARRIER = 1, 519df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_REDSCATTER = 2 520df4397b0SStefano Zampini /* Updates here must be accompanied by updates in finclude/petscsys.h and the string array in mpits.c */ 521df4397b0SStefano Zampini } PetscBuildTwoSidedType; 522df4397b0SStefano Zampini 523df4397b0SStefano Zampini /*E 524df4397b0SStefano Zampini InsertMode - Whether entries are inserted or added into vectors or matrices 525df4397b0SStefano Zampini 526df4397b0SStefano Zampini Level: beginner 527df4397b0SStefano Zampini 528df4397b0SStefano Zampini .seealso: VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 529df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), 530df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd() 531df4397b0SStefano Zampini E*/ 5320c91829fSJunchao Zhang typedef enum {NOT_SET_VALUES, INSERT_VALUES, ADD_VALUES, MAX_VALUES, MIN_VALUES, INSERT_ALL_VALUES, ADD_ALL_VALUES, INSERT_BC_VALUES, ADD_BC_VALUES} InsertMode; 533df4397b0SStefano Zampini 534df4397b0SStefano Zampini /*MC 535df4397b0SStefano Zampini INSERT_VALUES - Put a value into a vector or matrix, overwrites any previous value 536df4397b0SStefano Zampini 537df4397b0SStefano Zampini Level: beginner 538df4397b0SStefano Zampini 539df4397b0SStefano Zampini .seealso: InsertMode, VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 540df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), ADD_VALUES, 541df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd(), MAX_VALUES 542df4397b0SStefano Zampini 543df4397b0SStefano Zampini M*/ 544df4397b0SStefano Zampini 545df4397b0SStefano Zampini /*MC 546df4397b0SStefano Zampini ADD_VALUES - Adds a value into a vector or matrix, if there previously was no value, just puts the 547df4397b0SStefano Zampini value into that location 548df4397b0SStefano Zampini 549df4397b0SStefano Zampini Level: beginner 550df4397b0SStefano Zampini 551df4397b0SStefano Zampini .seealso: InsertMode, VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 552df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), INSERT_VALUES, 553df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd(), MAX_VALUES 554df4397b0SStefano Zampini 555df4397b0SStefano Zampini M*/ 556df4397b0SStefano Zampini 557df4397b0SStefano Zampini /*MC 558df4397b0SStefano Zampini MAX_VALUES - Puts the maximum of the scattered/gathered value and the current value into each location 559df4397b0SStefano Zampini 560df4397b0SStefano Zampini Level: beginner 561df4397b0SStefano Zampini 562df4397b0SStefano Zampini .seealso: InsertMode, VecScatterBegin(), VecScatterEnd(), ADD_VALUES, INSERT_VALUES 563df4397b0SStefano Zampini 564df4397b0SStefano Zampini M*/ 565df4397b0SStefano Zampini 566*421aa1e7SJunchao Zhang /*MC 567*421aa1e7SJunchao Zhang MIN_VALUES - Puts the minimal of the scattered/gathered value and the current value into each location 568*421aa1e7SJunchao Zhang 569*421aa1e7SJunchao Zhang Level: beginner 570*421aa1e7SJunchao Zhang 571*421aa1e7SJunchao Zhang .seealso: InsertMode, VecScatterBegin(), VecScatterEnd(), ADD_VALUES, INSERT_VALUES 572*421aa1e7SJunchao Zhang 573*421aa1e7SJunchao Zhang M*/ 574*421aa1e7SJunchao Zhang 575*421aa1e7SJunchao Zhang 576df4397b0SStefano Zampini /*S 577df4397b0SStefano Zampini PetscSubcomm - A decomposition of an MPI communicator into subcommunicators 578df4397b0SStefano Zampini 579df4397b0SStefano Zampini Notes: 580df4397b0SStefano Zampini After a call to PetscSubcommSetType(), PetscSubcommSetTypeGeneral(), or PetscSubcommSetFromOptions() one may call 581df4397b0SStefano Zampini $ PetscSubcommChild() returns the associated subcommunicator on this process 582df4397b0SStefano Zampini $ PetscSubcommContiguousParent() returns a parent communitor but with all child of the same subcommunicator having contiguous rank 583df4397b0SStefano Zampini 584df4397b0SStefano Zampini Sample Usage: 585df4397b0SStefano Zampini PetscSubcommCreate() 586df4397b0SStefano Zampini PetscSubcommSetNumber() 587df4397b0SStefano Zampini PetscSubcommSetType(PETSC_SUBCOMM_INTERLACED); 588df4397b0SStefano Zampini ccomm = PetscSubcommChild() 589df4397b0SStefano Zampini PetscSubcommDestroy() 590df4397b0SStefano Zampini 591df4397b0SStefano Zampini Level: advanced 592df4397b0SStefano Zampini 593df4397b0SStefano Zampini Notes: 594df4397b0SStefano 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 595df4397b0SStefano Zampini $ PETSC_SUBCOMM_CONTIGUOUS - each new communicator contains a set of process with contiguous ranks in the original MPI communicator 596df4397b0SStefano Zampini $ PETSC_SUBCOMM_INTERLACED - each new communictor contains a set of processes equally far apart in rank from the others in that new communicator 597df4397b0SStefano Zampini 598df4397b0SStefano Zampini Example: Consider a communicator with six processes split into 3 subcommunicators. 599df4397b0SStefano 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 600df4397b0SStefano Zampini $ PETSC_SUBCOMM_INTERLACED - the first communicator contains rank 0,3, the second 1,4 and the third 2,5 601df4397b0SStefano Zampini 602df4397b0SStefano Zampini Developer Notes: 603df4397b0SStefano Zampini This is used in objects such as PCREDUNDANT to manage the subcommunicators on which the redundant computations 604df4397b0SStefano Zampini are performed. 605df4397b0SStefano Zampini 606df4397b0SStefano Zampini 607df4397b0SStefano Zampini .seealso: PetscSubcommCreate(), PetscSubcommSetNumber(), PetscSubcommSetType(), PetscSubcommView(), PetscSubcommSetFromOptions() 608df4397b0SStefano Zampini 609df4397b0SStefano Zampini S*/ 610df4397b0SStefano Zampini typedef struct _n_PetscSubcomm* PetscSubcomm; 611df4397b0SStefano Zampini typedef enum {PETSC_SUBCOMM_GENERAL=0,PETSC_SUBCOMM_CONTIGUOUS=1,PETSC_SUBCOMM_INTERLACED=2} PetscSubcommType; 612df4397b0SStefano Zampini 613df4397b0SStefano Zampini /*S 614df4397b0SStefano Zampini PetscHeap - A simple class for managing heaps 615df4397b0SStefano Zampini 616df4397b0SStefano Zampini Level: intermediate 617df4397b0SStefano Zampini 618df4397b0SStefano Zampini .seealso: PetscHeapCreate(), PetscHeapAdd(), PetscHeapPop(), PetscHeapPeek(), PetscHeapStash(), PetscHeapUnstash(), PetscHeapView(), PetscHeapDestroy() 619df4397b0SStefano Zampini S*/ 620df4397b0SStefano Zampini typedef struct _PetscHeap *PetscHeap; 621df4397b0SStefano Zampini 622df4397b0SStefano Zampini typedef struct _n_PetscShmComm* PetscShmComm; 623df4397b0SStefano Zampini typedef struct _n_PetscOmpCtrl* PetscOmpCtrl; 624df4397b0SStefano Zampini 625df4397b0SStefano Zampini /*S 626df4397b0SStefano Zampini PetscSegBuffer - a segmented extendable buffer 627df4397b0SStefano Zampini 628df4397b0SStefano Zampini Level: developer 629df4397b0SStefano Zampini 630df4397b0SStefano Zampini .seealso: PetscSegBufferCreate(), PetscSegBufferGet(), PetscSegBufferExtract(), PetscSegBufferDestroy() 631df4397b0SStefano Zampini S*/ 632df4397b0SStefano Zampini typedef struct _n_PetscSegBuffer *PetscSegBuffer; 633df4397b0SStefano Zampini 634df4397b0SStefano Zampini typedef struct _n_PetscOptionsHelpPrinted *PetscOptionsHelpPrinted; 635df4397b0SStefano Zampini 636df4397b0SStefano Zampini #endif 637