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 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 */ 228df4397b0SStefano Zampini #elif defined(PETSC_HAVE_C99_COMPLEX) && !defined(PETSC_USE_REAL___FP16) 229df4397b0SStefano Zampini # if !defined(PETSC_SKIP_COMPLEX) 230df4397b0SStefano Zampini # define PETSC_HAVE_COMPLEX 1 231df4397b0SStefano Zampini # include <complex.h> 232df4397b0SStefano Zampini # if defined(PETSC_USE_REAL_SINGLE) || defined(PETSC_USE_REAL___FP16) 233df4397b0SStefano Zampini typedef float _Complex PetscComplex; 234df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL_DOUBLE) 235df4397b0SStefano Zampini typedef double _Complex PetscComplex; 236df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FLOAT128) 237df4397b0SStefano Zampini typedef __complex128 PetscComplex; 238df4397b0SStefano Zampini # endif /* PETSC_USE_REAL_* */ 239df4397b0SStefano Zampini # endif /* !PETSC_SKIP_COMPLEX */ 240df4397b0SStefano Zampini #elif (defined(PETSC_USE_COMPLEX) && !defined(PETSC_SKIP_COMPLEX)) 241df4397b0SStefano Zampini # error "PETSc was configured --with-scalar-type=complex, but a language-appropriate complex library is not available" 242df4397b0SStefano Zampini #endif /* !PETSC_SKIP_COMPLEX */ 243df4397b0SStefano Zampini 244df4397b0SStefano Zampini /*MC 245df4397b0SStefano Zampini PetscScalar - PETSc type that represents either a double precision real number, a double precision 246df4397b0SStefano Zampini complex number, a single precision real number, a __float128 real or complex or a __fp16 real - if the code is configured 247df4397b0SStefano Zampini with --with-scalar-type=real,complex --with-precision=single,double,__float128,__fp16 248df4397b0SStefano Zampini 249df4397b0SStefano Zampini Notes: 250df4397b0SStefano 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. 251df4397b0SStefano Zampini 252df4397b0SStefano Zampini Level: beginner 253df4397b0SStefano Zampini 254df4397b0SStefano Zampini .seealso: PetscReal, PetscComplex, PetscInt, MPIU_REAL, MPIU_SCALAR, MPIU_COMPLEX, MPIU_INT, PetscRealPart(), PetscImaginaryPart() 255df4397b0SStefano Zampini M*/ 256df4397b0SStefano Zampini 257df4397b0SStefano Zampini #if (defined(PETSC_USE_COMPLEX) && !defined(PETSC_SKIP_COMPLEX)) 258df4397b0SStefano Zampini typedef PetscComplex PetscScalar; 259df4397b0SStefano Zampini #else /* PETSC_USE_COMPLEX */ 260df4397b0SStefano Zampini typedef PetscReal PetscScalar; 261df4397b0SStefano Zampini #endif /* PETSC_USE_COMPLEX */ 262df4397b0SStefano Zampini 263df4397b0SStefano Zampini /*E 264df4397b0SStefano Zampini PetscCopyMode - Determines how an array passed to certain functions is copied or retained 265df4397b0SStefano Zampini 266df4397b0SStefano Zampini Level: beginner 267df4397b0SStefano Zampini 268df4397b0SStefano Zampini $ PETSC_COPY_VALUES - the array values are copied into new space, the user is free to reuse or delete the passed in array 269df4397b0SStefano 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 270df4397b0SStefano Zampini $ delete the array. The array MUST have been obtained with PetscMalloc(). Hence this mode cannot be used in Fortran. 271df4397b0SStefano 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 272df4397b0SStefano Zampini the array but the user must delete the array after the object is destroyed. 273df4397b0SStefano Zampini 274df4397b0SStefano Zampini E*/ 275df4397b0SStefano Zampini typedef enum {PETSC_COPY_VALUES, PETSC_OWN_POINTER, PETSC_USE_POINTER} PetscCopyMode; 276df4397b0SStefano Zampini 277df4397b0SStefano Zampini /*MC 278df4397b0SStefano Zampini PETSC_FALSE - False value of PetscBool 279df4397b0SStefano Zampini 280df4397b0SStefano Zampini Level: beginner 281df4397b0SStefano Zampini 282df4397b0SStefano Zampini Note: 283df4397b0SStefano Zampini Zero integer 284df4397b0SStefano Zampini 285df4397b0SStefano Zampini .seealso: PetscBool, PETSC_TRUE 286df4397b0SStefano Zampini M*/ 287df4397b0SStefano Zampini 288df4397b0SStefano Zampini /*MC 289df4397b0SStefano Zampini PETSC_TRUE - True value of PetscBool 290df4397b0SStefano Zampini 291df4397b0SStefano Zampini Level: beginner 292df4397b0SStefano Zampini 293df4397b0SStefano Zampini Note: 294df4397b0SStefano Zampini Nonzero integer 295df4397b0SStefano Zampini 296df4397b0SStefano Zampini .seealso: PetscBool, PETSC_FALSE 297df4397b0SStefano Zampini M*/ 298df4397b0SStefano Zampini 299df4397b0SStefano Zampini /*MC 300df4397b0SStefano Zampini PetscLogDouble - Used for logging times 301df4397b0SStefano Zampini 302df4397b0SStefano Zampini Notes: 303df4397b0SStefano Zampini Contains double precision numbers that are not used in the numerical computations, but rather in logging, timing etc. 304df4397b0SStefano Zampini 3056b3bf505SPatrick Sanan Level: developer 3066b3bf505SPatrick Sanan 307df4397b0SStefano Zampini M*/ 308df4397b0SStefano Zampini typedef double PetscLogDouble; 309df4397b0SStefano Zampini 310df4397b0SStefano Zampini /*E 311df4397b0SStefano Zampini PetscDataType - Used for handling different basic data types. 312df4397b0SStefano Zampini 313df4397b0SStefano Zampini Level: beginner 314df4397b0SStefano Zampini 315df4397b0SStefano Zampini Notes: 316df4397b0SStefano Zampini Use of this should be avoided if one can directly use MPI_Datatype instead. 317df4397b0SStefano Zampini 318df4397b0SStefano Zampini Developer comment: 319df4397b0SStefano Zampini It would be nice if we could always just use MPI Datatypes, why can we not? 320df4397b0SStefano Zampini 321*91e29162SBarry Smith If you change any values in PetscDatatype make sure you update their usage in 322*91e29162SBarry Smith share/petsc/matlab/PetscBagRead.m 323*91e29162SBarry Smith 324df4397b0SStefano Zampini .seealso: PetscBinaryRead(), PetscBinaryWrite(), PetscDataTypeToMPIDataType(), 325df4397b0SStefano Zampini PetscDataTypeGetSize() 326df4397b0SStefano Zampini 327df4397b0SStefano Zampini E*/ 3282ebc710fSVaclav Hapla typedef enum {PETSC_DATATYPE_UNKNOWN = 0, 3292ebc710fSVaclav Hapla PETSC_DOUBLE = 1, PETSC_COMPLEX = 2, PETSC_LONG = 3, PETSC_SHORT = 4, PETSC_FLOAT = 5, 33049de050aSVaclav Hapla PETSC_CHAR = 6, PETSC_BIT_LOGICAL = 7, PETSC_ENUM = 8, PETSC_BOOL = 9, PETSC___FLOAT128 = 10, 331ea127267SVaclav Hapla PETSC_OBJECT = 11, PETSC_FUNCTION = 12, PETSC_STRING = 13, PETSC___FP16 = 14, PETSC_STRUCT = 15, 3329e3e4c22SLisandro Dalcin PETSC_INT = 16, PETSC_INT64 = 17} PetscDataType; 333df4397b0SStefano Zampini 334df4397b0SStefano Zampini #if defined(PETSC_USE_COMPLEX) 335df4397b0SStefano Zampini # define PETSC_SCALAR PETSC_COMPLEX 336df4397b0SStefano Zampini #else 337df4397b0SStefano Zampini # if defined(PETSC_USE_REAL_SINGLE) 338df4397b0SStefano Zampini # define PETSC_SCALAR PETSC_FLOAT 339df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FLOAT128) 340df4397b0SStefano Zampini # define PETSC_SCALAR PETSC___FLOAT128 341df4397b0SStefano Zampini # elif defined(PETSC_USE_REAL___FP16) 342df4397b0SStefano Zampini # define PETSC_SCALAR PETSC___FP16 343df4397b0SStefano Zampini # else 344df4397b0SStefano Zampini # define PETSC_SCALAR PETSC_DOUBLE 345df4397b0SStefano Zampini # endif 346df4397b0SStefano Zampini #endif 347df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 348df4397b0SStefano Zampini # define PETSC_REAL PETSC_FLOAT 349df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 350df4397b0SStefano Zampini # define PETSC_REAL PETSC___FLOAT128 351df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FP16) 352df4397b0SStefano Zampini # define PETSC_REAL PETSC___FP16 353df4397b0SStefano Zampini #else 354df4397b0SStefano Zampini # define PETSC_REAL PETSC_DOUBLE 355df4397b0SStefano Zampini #endif 356df4397b0SStefano Zampini #define PETSC_FORTRANADDR PETSC_LONG 357df4397b0SStefano Zampini 358df4397b0SStefano Zampini /*S 359df4397b0SStefano Zampini PetscToken - 'Token' used for managing tokenizing strings 360df4397b0SStefano Zampini 361df4397b0SStefano Zampini Level: intermediate 362df4397b0SStefano Zampini 363df4397b0SStefano Zampini .seealso: PetscTokenCreate(), PetscTokenFind(), PetscTokenDestroy() 364df4397b0SStefano Zampini S*/ 365df4397b0SStefano Zampini typedef struct _p_PetscToken* PetscToken; 366df4397b0SStefano Zampini 367df4397b0SStefano Zampini /*S 368df4397b0SStefano Zampini PetscObject - any PETSc object, PetscViewer, Mat, Vec, KSP etc 369df4397b0SStefano Zampini 370df4397b0SStefano Zampini Level: beginner 371df4397b0SStefano Zampini 372df4397b0SStefano Zampini Note: 373df4397b0SStefano Zampini This is the base class from which all PETSc objects are derived from. 374df4397b0SStefano Zampini 375df4397b0SStefano Zampini .seealso: PetscObjectDestroy(), PetscObjectView(), PetscObjectGetName(), PetscObjectSetName(), PetscObjectReference(), PetscObjectDereference() 376df4397b0SStefano Zampini S*/ 377df4397b0SStefano Zampini typedef struct _p_PetscObject* PetscObject; 378df4397b0SStefano Zampini 379df4397b0SStefano Zampini /*MC 380df4397b0SStefano Zampini PetscObjectId - unique integer Id for a PetscObject 381df4397b0SStefano Zampini 382df4397b0SStefano Zampini Level: developer 383df4397b0SStefano Zampini 384df4397b0SStefano Zampini Notes: 385df4397b0SStefano Zampini Unlike pointer values, object ids are never reused. 386df4397b0SStefano Zampini 387df4397b0SStefano Zampini .seealso: PetscObjectState, PetscObjectGetId() 388df4397b0SStefano Zampini M*/ 389df4397b0SStefano Zampini #if defined(PETSC_USING_F90) && !defined(PETSC_USE_FORTRANKIND) /* compaq F90 */ 390df4397b0SStefano Zampini typedef int PetscObjectId; 391df4397b0SStefano Zampini #else 392df4397b0SStefano Zampini typedef PetscInt64 PetscObjectId; 393df4397b0SStefano Zampini #endif 394df4397b0SStefano Zampini 395df4397b0SStefano Zampini /*MC 396df4397b0SStefano Zampini PetscObjectState - integer state for a PetscObject 397df4397b0SStefano Zampini 398df4397b0SStefano Zampini Level: developer 399df4397b0SStefano Zampini 400df4397b0SStefano Zampini Notes: 401df4397b0SStefano Zampini Object state is always-increasing and (for objects that track state) can be used to determine if an object has 402df4397b0SStefano 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. 403df4397b0SStefano Zampini 404df4397b0SStefano Zampini .seealso: PetscObjectId, PetscObjectStateGet(), PetscObjectStateIncrease(), PetscObjectStateSet() 405df4397b0SStefano Zampini M*/ 406df4397b0SStefano Zampini #if defined(PETSC_USING_F90) && !defined(PETSC_USE_FORTRANKIND) /* compaq F90 */ 407df4397b0SStefano Zampini typedef int PetscObjectState; 408df4397b0SStefano Zampini #else 409df4397b0SStefano Zampini typedef PetscInt64 PetscObjectState; 410df4397b0SStefano Zampini #endif 411df4397b0SStefano Zampini 412df4397b0SStefano Zampini /*S 413df4397b0SStefano Zampini PetscFunctionList - Linked list of functions, possibly stored in dynamic libraries, accessed 414df4397b0SStefano Zampini by string name 415df4397b0SStefano Zampini 416df4397b0SStefano Zampini Level: advanced 417df4397b0SStefano Zampini 418df4397b0SStefano Zampini .seealso: PetscFunctionListAdd(), PetscFunctionListDestroy(), PetscOpFlist 419df4397b0SStefano Zampini S*/ 420df4397b0SStefano Zampini typedef struct _n_PetscFunctionList *PetscFunctionList; 421df4397b0SStefano Zampini 422df4397b0SStefano Zampini /*E 423df4397b0SStefano Zampini PetscFileMode - Access mode for a file. 424df4397b0SStefano Zampini 425df4397b0SStefano Zampini Level: beginner 426df4397b0SStefano Zampini 427df4397b0SStefano Zampini $ FILE_MODE_READ - open a file at its beginning for reading 428df4397b0SStefano Zampini $ FILE_MODE_WRITE - open a file at its beginning for writing (will create if the file does not exist) 429df4397b0SStefano Zampini $ FILE_MODE_APPEND - open a file at end for writing 430df4397b0SStefano Zampini $ FILE_MODE_UPDATE - open a file for updating, meaning for reading and writing 431df4397b0SStefano Zampini $ FILE_MODE_APPEND_UPDATE - open a file for updating, meaning for reading and writing, at the end 432df4397b0SStefano Zampini 433df4397b0SStefano Zampini .seealso: PetscViewerFileSetMode() 434df4397b0SStefano Zampini E*/ 435df4397b0SStefano Zampini typedef enum {FILE_MODE_READ, FILE_MODE_WRITE, FILE_MODE_APPEND, FILE_MODE_UPDATE, FILE_MODE_APPEND_UPDATE} PetscFileMode; 436df4397b0SStefano Zampini 437df4397b0SStefano Zampini typedef void* PetscDLHandle; 438df4397b0SStefano Zampini typedef enum {PETSC_DL_DECIDE=0,PETSC_DL_NOW=1,PETSC_DL_LOCAL=2} PetscDLMode; 439df4397b0SStefano Zampini 440df4397b0SStefano Zampini /*S 441df4397b0SStefano Zampini PetscObjectList - Linked list of PETSc objects, each accessable by string name 442df4397b0SStefano Zampini 443df4397b0SStefano Zampini Level: developer 444df4397b0SStefano Zampini 445df4397b0SStefano Zampini Notes: 446df4397b0SStefano Zampini Used by PetscObjectCompose() and PetscObjectQuery() 447df4397b0SStefano Zampini 448df4397b0SStefano Zampini .seealso: PetscObjectListAdd(), PetscObjectListDestroy(), PetscObjectListFind(), PetscObjectCompose(), PetscObjectQuery(), PetscFunctionList 449df4397b0SStefano Zampini S*/ 450df4397b0SStefano Zampini typedef struct _n_PetscObjectList *PetscObjectList; 451df4397b0SStefano Zampini 452df4397b0SStefano Zampini /*S 453df4397b0SStefano Zampini PetscDLLibrary - Linked list of dynamics libraries to search for functions 454df4397b0SStefano Zampini 455df4397b0SStefano Zampini Level: advanced 456df4397b0SStefano Zampini 457df4397b0SStefano Zampini .seealso: PetscDLLibraryOpen() 458df4397b0SStefano Zampini S*/ 459df4397b0SStefano Zampini typedef struct _n_PetscDLLibrary *PetscDLLibrary; 460df4397b0SStefano Zampini 461df4397b0SStefano Zampini /*S 462df4397b0SStefano Zampini PetscContainer - Simple PETSc object that contains a pointer to any required data 463df4397b0SStefano Zampini 464df4397b0SStefano Zampini Level: advanced 465df4397b0SStefano Zampini 466df4397b0SStefano Zampini .seealso: PetscObject, PetscContainerCreate() 467df4397b0SStefano Zampini S*/ 468df4397b0SStefano Zampini typedef struct _p_PetscContainer* PetscContainer; 469df4397b0SStefano Zampini 470df4397b0SStefano Zampini /*S 471df4397b0SStefano Zampini PetscRandom - Abstract PETSc object that manages generating random numbers 472df4397b0SStefano Zampini 473df4397b0SStefano Zampini Level: intermediate 474df4397b0SStefano Zampini 475df4397b0SStefano Zampini Concepts: random numbers 476df4397b0SStefano Zampini 477df4397b0SStefano Zampini .seealso: PetscRandomCreate(), PetscRandomGetValue(), PetscRandomType 478df4397b0SStefano Zampini S*/ 479df4397b0SStefano Zampini typedef struct _p_PetscRandom* PetscRandom; 480df4397b0SStefano Zampini 481df4397b0SStefano Zampini /* 482df4397b0SStefano Zampini In binary files variables are stored using the following lengths, 483df4397b0SStefano Zampini regardless of how they are stored in memory on any one particular 484df4397b0SStefano Zampini machine. Use these rather then sizeof() in computing sizes for 485df4397b0SStefano Zampini PetscBinarySeek(). 486df4397b0SStefano Zampini */ 487df4397b0SStefano Zampini #define PETSC_BINARY_INT_SIZE (32/8) 488df4397b0SStefano Zampini #define PETSC_BINARY_FLOAT_SIZE (32/8) 489df4397b0SStefano Zampini #define PETSC_BINARY_CHAR_SIZE (8/8) 490df4397b0SStefano Zampini #define PETSC_BINARY_SHORT_SIZE (16/8) 491df4397b0SStefano Zampini #define PETSC_BINARY_DOUBLE_SIZE (64/8) 492df4397b0SStefano Zampini #define PETSC_BINARY_SCALAR_SIZE sizeof(PetscScalar) 493df4397b0SStefano Zampini 494df4397b0SStefano Zampini /*E 495df4397b0SStefano Zampini PetscBinarySeekType - argument to PetscBinarySeek() 496df4397b0SStefano Zampini 497df4397b0SStefano Zampini Level: advanced 498df4397b0SStefano Zampini 499df4397b0SStefano Zampini .seealso: PetscBinarySeek(), PetscBinarySynchronizedSeek() 500df4397b0SStefano Zampini E*/ 501df4397b0SStefano Zampini typedef enum {PETSC_BINARY_SEEK_SET = 0,PETSC_BINARY_SEEK_CUR = 1,PETSC_BINARY_SEEK_END = 2} PetscBinarySeekType; 502df4397b0SStefano Zampini 503df4397b0SStefano Zampini /*E 504df4397b0SStefano Zampini PetscBuildTwoSidedType - algorithm for setting up two-sided communication 505df4397b0SStefano Zampini 506df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_ALLREDUCE - classical algorithm using an MPI_Allreduce with 507df4397b0SStefano Zampini $ a buffer of length equal to the communicator size. Not memory-scalable due to 508df4397b0SStefano Zampini $ the large reduction size. Requires only MPI-1. 509df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_IBARRIER - nonblocking algorithm based on MPI_Issend and MPI_Ibarrier. 510df4397b0SStefano Zampini $ Proved communication-optimal in Hoefler, Siebert, and Lumsdaine (2010). Requires MPI-3. 511df4397b0SStefano Zampini $ PETSC_BUILDTWOSIDED_REDSCATTER - similar to above, but use more optimized function 512df4397b0SStefano Zampini $ that only communicates the part of the reduction that is necessary. Requires MPI-2. 513df4397b0SStefano Zampini 514df4397b0SStefano Zampini Level: developer 515df4397b0SStefano Zampini 516df4397b0SStefano Zampini .seealso: PetscCommBuildTwoSided(), PetscCommBuildTwoSidedSetType(), PetscCommBuildTwoSidedGetType() 517df4397b0SStefano Zampini E*/ 518df4397b0SStefano Zampini typedef enum { 519df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_NOTSET = -1, 520df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_ALLREDUCE = 0, 521df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_IBARRIER = 1, 522df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_REDSCATTER = 2 523df4397b0SStefano Zampini /* Updates here must be accompanied by updates in finclude/petscsys.h and the string array in mpits.c */ 524df4397b0SStefano Zampini } PetscBuildTwoSidedType; 525df4397b0SStefano Zampini 526df4397b0SStefano Zampini /*E 527df4397b0SStefano Zampini InsertMode - Whether entries are inserted or added into vectors or matrices 528df4397b0SStefano Zampini 529df4397b0SStefano Zampini Level: beginner 530df4397b0SStefano Zampini 531df4397b0SStefano Zampini .seealso: VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 532df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), 533df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd() 534df4397b0SStefano Zampini E*/ 535df4397b0SStefano 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; 536df4397b0SStefano Zampini 537df4397b0SStefano Zampini /*MC 538df4397b0SStefano Zampini INSERT_VALUES - Put a value into a vector or matrix, overwrites any previous value 539df4397b0SStefano Zampini 540df4397b0SStefano Zampini Level: beginner 541df4397b0SStefano Zampini 542df4397b0SStefano Zampini .seealso: InsertMode, VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 543df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), ADD_VALUES, 544df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd(), MAX_VALUES 545df4397b0SStefano Zampini 546df4397b0SStefano Zampini M*/ 547df4397b0SStefano Zampini 548df4397b0SStefano Zampini /*MC 549df4397b0SStefano Zampini ADD_VALUES - Adds a value into a vector or matrix, if there previously was no value, just puts the 550df4397b0SStefano Zampini value into that location 551df4397b0SStefano Zampini 552df4397b0SStefano Zampini Level: beginner 553df4397b0SStefano Zampini 554df4397b0SStefano Zampini .seealso: InsertMode, VecSetValues(), MatSetValues(), VecSetValue(), VecSetValuesBlocked(), 555df4397b0SStefano Zampini VecSetValuesLocal(), VecSetValuesBlockedLocal(), MatSetValuesBlocked(), INSERT_VALUES, 556df4397b0SStefano Zampini MatSetValuesBlockedLocal(), MatSetValuesLocal(), VecScatterBegin(), VecScatterEnd(), MAX_VALUES 557df4397b0SStefano Zampini 558df4397b0SStefano Zampini M*/ 559df4397b0SStefano Zampini 560df4397b0SStefano Zampini /*MC 561df4397b0SStefano Zampini MAX_VALUES - Puts the maximum of the scattered/gathered value and the current value into each location 562df4397b0SStefano Zampini 563df4397b0SStefano Zampini Level: beginner 564df4397b0SStefano Zampini 565df4397b0SStefano Zampini .seealso: InsertMode, VecScatterBegin(), VecScatterEnd(), ADD_VALUES, INSERT_VALUES 566df4397b0SStefano Zampini 567df4397b0SStefano Zampini M*/ 568df4397b0SStefano Zampini 569df4397b0SStefano Zampini /*S 570df4397b0SStefano Zampini PetscSubcomm - A decomposition of an MPI communicator into subcommunicators 571df4397b0SStefano Zampini 572df4397b0SStefano Zampini Notes: 573df4397b0SStefano Zampini After a call to PetscSubcommSetType(), PetscSubcommSetTypeGeneral(), or PetscSubcommSetFromOptions() one may call 574df4397b0SStefano Zampini $ PetscSubcommChild() returns the associated subcommunicator on this process 575df4397b0SStefano Zampini $ PetscSubcommContiguousParent() returns a parent communitor but with all child of the same subcommunicator having contiguous rank 576df4397b0SStefano Zampini 577df4397b0SStefano Zampini Sample Usage: 578df4397b0SStefano Zampini PetscSubcommCreate() 579df4397b0SStefano Zampini PetscSubcommSetNumber() 580df4397b0SStefano Zampini PetscSubcommSetType(PETSC_SUBCOMM_INTERLACED); 581df4397b0SStefano Zampini ccomm = PetscSubcommChild() 582df4397b0SStefano Zampini PetscSubcommDestroy() 583df4397b0SStefano Zampini 584df4397b0SStefano Zampini Level: advanced 585df4397b0SStefano Zampini 586df4397b0SStefano Zampini Concepts: communicator, create 587df4397b0SStefano Zampini 588df4397b0SStefano Zampini Notes: 589df4397b0SStefano 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 590df4397b0SStefano Zampini $ PETSC_SUBCOMM_CONTIGUOUS - each new communicator contains a set of process with contiguous ranks in the original MPI communicator 591df4397b0SStefano Zampini $ PETSC_SUBCOMM_INTERLACED - each new communictor contains a set of processes equally far apart in rank from the others in that new communicator 592df4397b0SStefano Zampini 593df4397b0SStefano Zampini Example: Consider a communicator with six processes split into 3 subcommunicators. 594df4397b0SStefano 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 595df4397b0SStefano Zampini $ PETSC_SUBCOMM_INTERLACED - the first communicator contains rank 0,3, the second 1,4 and the third 2,5 596df4397b0SStefano Zampini 597df4397b0SStefano Zampini Developer Notes: 598df4397b0SStefano Zampini This is used in objects such as PCREDUNDANT to manage the subcommunicators on which the redundant computations 599df4397b0SStefano Zampini are performed. 600df4397b0SStefano Zampini 601df4397b0SStefano Zampini 602df4397b0SStefano Zampini .seealso: PetscSubcommCreate(), PetscSubcommSetNumber(), PetscSubcommSetType(), PetscSubcommView(), PetscSubcommSetFromOptions() 603df4397b0SStefano Zampini 604df4397b0SStefano Zampini S*/ 605df4397b0SStefano Zampini typedef struct _n_PetscSubcomm* PetscSubcomm; 606df4397b0SStefano Zampini typedef enum {PETSC_SUBCOMM_GENERAL=0,PETSC_SUBCOMM_CONTIGUOUS=1,PETSC_SUBCOMM_INTERLACED=2} PetscSubcommType; 607df4397b0SStefano Zampini 608df4397b0SStefano Zampini /*S 609df4397b0SStefano Zampini PetscHeap - A simple class for managing heaps 610df4397b0SStefano Zampini 611df4397b0SStefano Zampini Level: intermediate 612df4397b0SStefano Zampini 613df4397b0SStefano Zampini Concepts: random numbers 614df4397b0SStefano Zampini 615df4397b0SStefano Zampini .seealso: PetscHeapCreate(), PetscHeapAdd(), PetscHeapPop(), PetscHeapPeek(), PetscHeapStash(), PetscHeapUnstash(), PetscHeapView(), PetscHeapDestroy() 616df4397b0SStefano Zampini S*/ 617df4397b0SStefano Zampini typedef struct _PetscHeap *PetscHeap; 618df4397b0SStefano Zampini 619df4397b0SStefano Zampini typedef struct _n_PetscShmComm* PetscShmComm; 620df4397b0SStefano Zampini typedef struct _n_PetscOmpCtrl* PetscOmpCtrl; 621df4397b0SStefano Zampini 622df4397b0SStefano Zampini /*S 623df4397b0SStefano Zampini PetscSegBuffer - a segmented extendable buffer 624df4397b0SStefano Zampini 625df4397b0SStefano Zampini Level: developer 626df4397b0SStefano Zampini 627df4397b0SStefano Zampini .seealso: PetscSegBufferCreate(), PetscSegBufferGet(), PetscSegBufferExtract(), PetscSegBufferDestroy() 628df4397b0SStefano Zampini S*/ 629df4397b0SStefano Zampini typedef struct _n_PetscSegBuffer *PetscSegBuffer; 630df4397b0SStefano Zampini 631df4397b0SStefano Zampini typedef struct _n_PetscOptionsHelpPrinted *PetscOptionsHelpPrinted; 632df4397b0SStefano Zampini 633df4397b0SStefano Zampini #endif 634