126bd1501SBarry Smith #if !defined(PETSCSYSTYPES_H) 226bd1501SBarry Smith #define PETSCSYSTYPES_H 3df4397b0SStefano Zampini 4df4397b0SStefano Zampini #include <petscconf.h> 5e1bf4ed2SJacob Faibussowitsch #include <petscconf_poison.h> 6df4397b0SStefano Zampini #include <petscfix.h> 73edda6a2SJed Brown #include <stddef.h> 8df4397b0SStefano Zampini 9ac09b921SBarry Smith /* SUBMANSEC = Sys */ 10ac09b921SBarry Smith 11df4397b0SStefano Zampini /*MC 12df4397b0SStefano Zampini PetscErrorCode - datatype used for return error code from almost all PETSc functions 13df4397b0SStefano Zampini 14df4397b0SStefano Zampini Level: beginner 15df4397b0SStefano Zampini 16db781477SPatrick Sanan .seealso: `PetscCall()`, `SETERRQ()` 17df4397b0SStefano Zampini M*/ 18df4397b0SStefano Zampini typedef int PetscErrorCode; 19df4397b0SStefano Zampini 20df4397b0SStefano Zampini /*MC 21df4397b0SStefano Zampini 22df4397b0SStefano Zampini PetscClassId - A unique id used to identify each PETSc class. 23df4397b0SStefano Zampini 24df4397b0SStefano Zampini Notes: 25*87497f52SBarry Smith Use `PetscClassIdRegister()` to obtain a new value for a new class being created. Usually 26df4397b0SStefano Zampini XXXInitializePackage() calls it for each class it defines. 27df4397b0SStefano Zampini 28df4397b0SStefano Zampini Developer Notes: 29*87497f52SBarry Smith Internal integer stored in the `_p_PetscObject` data structure. 30*87497f52SBarry Smith These are all computed by an offset from the lowest one, `PETSC_SMALLEST_CLASSID`. 31df4397b0SStefano Zampini 32df4397b0SStefano Zampini Level: developer 33df4397b0SStefano Zampini 34db781477SPatrick Sanan .seealso: `PetscClassIdRegister()`, `PetscLogEventRegister()`, `PetscHeaderCreate()` 35df4397b0SStefano Zampini M*/ 36df4397b0SStefano Zampini typedef int PetscClassId; 37df4397b0SStefano Zampini 38df4397b0SStefano Zampini /*MC 39df4397b0SStefano Zampini PetscMPIInt - datatype used to represent 'int' parameters to MPI functions. 40df4397b0SStefano Zampini 41df4397b0SStefano Zampini Level: intermediate 42df4397b0SStefano Zampini 43df4397b0SStefano Zampini Notes: 44*87497f52SBarry Smith This is always a 32 bit integer, sometimes it is the same as `PetscInt`, but if PETSc was built with --with-64-bit-indices but 45*87497f52SBarry Smith standard C/Fortran integers are 32 bit then this is NOT the same as `PetscInt`; it remains 32 bit. 46df4397b0SStefano Zampini 47*87497f52SBarry Smith `PetscMPIIntCast`(a,&b) checks if the given `PetscInt` a will fit in a `PetscMPIInt`, if not it 48*87497f52SBarry Smith generates a `PETSC_ERR_ARG_OUTOFRANGE` error. 49df4397b0SStefano Zampini 50db781477SPatrick Sanan .seealso: `PetscBLASInt`, `PetscInt`, `PetscMPIIntCast()` 51df4397b0SStefano Zampini 52df4397b0SStefano Zampini M*/ 53df4397b0SStefano Zampini typedef int PetscMPIInt; 54df4397b0SStefano Zampini 55df4397b0SStefano Zampini /*MC 563edda6a2SJed Brown PetscSizeT - datatype used to represent sizes in memory (like size_t) 573edda6a2SJed Brown 583edda6a2SJed Brown Level: intermediate 593edda6a2SJed Brown 603edda6a2SJed Brown Notes: 613edda6a2SJed Brown This is equivalent to size_t, but defined for consistency with Fortran, which lacks a native equivalent of size_t. 623edda6a2SJed Brown 63db781477SPatrick Sanan .seealso: `PetscInt`, `PetscInt64`, `PetscCount` 643edda6a2SJed Brown 653edda6a2SJed Brown M*/ 663edda6a2SJed Brown typedef size_t PetscSizeT; 673edda6a2SJed Brown 683edda6a2SJed Brown /*MC 6982a78a4eSJed Brown PetscCount - signed datatype used to represent counts 7082a78a4eSJed Brown 7182a78a4eSJed Brown Level: intermediate 7282a78a4eSJed Brown 7382a78a4eSJed Brown Notes: 7482a78a4eSJed Brown This is equivalent to ptrdiff_t, but defined for consistency with Fortran, which lacks a native equivalent of ptrdiff_t. 7582a78a4eSJed Brown 76*87497f52SBarry Smith Use `PetscCount_FMT` to format with `PetscPrintf()`, `printf()`, and related functions. 7782a78a4eSJed Brown 78db781477SPatrick Sanan .seealso: `PetscInt`, `PetscInt64`, `PetscSizeT` 7982a78a4eSJed Brown 8082a78a4eSJed Brown M*/ 8182a78a4eSJed Brown typedef ptrdiff_t PetscCount; 8263a3b9bcSJacob Faibussowitsch #define PetscCount_FMT "td" 8382a78a4eSJed Brown 8482a78a4eSJed Brown /*MC 85df4397b0SStefano Zampini PetscEnum - datatype used to pass enum types within PETSc functions. 86df4397b0SStefano Zampini 87df4397b0SStefano Zampini Level: intermediate 88df4397b0SStefano Zampini 89db781477SPatrick Sanan .seealso: `PetscOptionsGetEnum()`, `PetscOptionsEnum()`, `PetscBagRegisterEnum()` 90df4397b0SStefano Zampini M*/ 91df4397b0SStefano Zampini typedef enum { ENUM_DUMMY } PetscEnum; 92df4397b0SStefano Zampini 93df4397b0SStefano Zampini typedef short PetscShort; 94df4397b0SStefano Zampini typedef char PetscChar; 95df4397b0SStefano Zampini typedef float PetscFloat; 96df4397b0SStefano Zampini 97df4397b0SStefano Zampini /*MC 98df4397b0SStefano Zampini PetscInt - PETSc type that represents an integer, used primarily to 99df4397b0SStefano 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. 100df4397b0SStefano Zampini 101df4397b0SStefano Zampini Notes: 102*87497f52SBarry Smith 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. 103df4397b0SStefano Zampini 104df4397b0SStefano Zampini Level: beginner 105df4397b0SStefano Zampini 106db781477SPatrick Sanan .seealso: `PetscBLASInt`, `PetscMPIInt`, `PetscReal`, `PetscScalar`, `PetscComplex`, `PetscInt`, `MPIU_REAL`, `MPIU_SCALAR`, `MPIU_COMPLEX`, `MPIU_INT` 107df4397b0SStefano Zampini M*/ 108df4397b0SStefano Zampini 109df4397b0SStefano Zampini #if defined(PETSC_HAVE_STDINT_H) 110df4397b0SStefano Zampini # include <stdint.h> 111df4397b0SStefano Zampini #endif 112df4397b0SStefano Zampini #if defined (PETSC_HAVE_INTTYPES_H) 113df4397b0SStefano Zampini # if !defined(__STDC_FORMAT_MACROS) 114df4397b0SStefano Zampini # define __STDC_FORMAT_MACROS /* required for using PRId64 from c++ */ 115df4397b0SStefano Zampini # endif 116df4397b0SStefano Zampini # include <inttypes.h> 117df4397b0SStefano Zampini # if !defined(PRId64) 118df4397b0SStefano Zampini # define PRId64 "ld" 119df4397b0SStefano Zampini # endif 120df4397b0SStefano Zampini #endif 121df4397b0SStefano Zampini 122df4397b0SStefano 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 */ 123df4397b0SStefano Zampini typedef int64_t PetscInt64; 124df4397b0SStefano Zampini #elif (PETSC_SIZEOF_LONG_LONG == 8) 125df4397b0SStefano Zampini typedef long long PetscInt64; 126df4397b0SStefano Zampini #elif defined(PETSC_HAVE___INT64) 127df4397b0SStefano Zampini typedef __int64 PetscInt64; 128df4397b0SStefano Zampini #else 129df4397b0SStefano Zampini # error "cannot determine PetscInt64 type" 130df4397b0SStefano Zampini #endif 131df4397b0SStefano Zampini 132df4397b0SStefano Zampini #if defined(PETSC_USE_64BIT_INDICES) 133df4397b0SStefano Zampini typedef PetscInt64 PetscInt; 134df4397b0SStefano Zampini #else 135df4397b0SStefano Zampini typedef int PetscInt; 136df4397b0SStefano Zampini #endif 137df4397b0SStefano Zampini 138c93fae50SJacob Faibussowitsch #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 */ 139c93fae50SJacob Faibussowitsch # define MPIU_INT64 MPI_INT64_T 140c93fae50SJacob Faibussowitsch # define PetscInt64_FMT PRId64 141c93fae50SJacob Faibussowitsch #elif (PETSC_SIZEOF_LONG_LONG == 8) 142c93fae50SJacob Faibussowitsch # define MPIU_INT64 MPI_LONG_LONG_INT 143c93fae50SJacob Faibussowitsch # define PetscInt64_FMT "lld" 144c93fae50SJacob Faibussowitsch #elif defined(PETSC_HAVE___INT64) 145c93fae50SJacob Faibussowitsch # define MPIU_INT64 MPI_INT64_T 146c93fae50SJacob Faibussowitsch # define PetscInt64_FMT "ld" 147c93fae50SJacob Faibussowitsch #else 148c93fae50SJacob Faibussowitsch # error "cannot determine PetscInt64 type" 149c93fae50SJacob Faibussowitsch #endif 150c93fae50SJacob Faibussowitsch 151df4397b0SStefano Zampini /*MC 152df4397b0SStefano Zampini PetscBLASInt - datatype used to represent 'int' parameters to BLAS/LAPACK functions. 153df4397b0SStefano Zampini 154df4397b0SStefano Zampini Notes: 155*87497f52SBarry Smith Usually this is the same as `PetscIn`t, but if PETSc was built with --with-64-bit-indices but 156*87497f52SBarry Smith standard C/Fortran integers are 32 bit then this may not be the same as `PetscInt`, 157*87497f52SBarry Smith except on some BLAS/LAPACK implementations that support 64 bit integers see the notes below. 158df4397b0SStefano Zampini 159*87497f52SBarry Smith `PetscErrorCode` `PetscBLASIntCast`(a,&b) checks if the given `PetscInt` a will fit in a `PetscBLASInt`, if not it 160*87497f52SBarry Smith generates a `PETSC_ERR_ARG_OUTOFRANGE` error 161df4397b0SStefano Zampini 162df4397b0SStefano Zampini Installation Notes: 163cdc6ee08SBarry Smith ./configure automatically determines the size of the integers used by BLAS/LAPACK except when --with-batch is used 164cdc6ee08SBarry 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 165cdc6ee08SBarry Smith 166cdc6ee08SBarry Smith MATLAB ships with BLAS and LAPACK that use 64 bit integers, for example if you run ./configure with, the option 167df4397b0SStefano Zampini --with-blaslapack-lib=[/Applications/MATLAB_R2010b.app/bin/maci64/libmwblas.dylib,/Applications/MATLAB_R2010b.app/bin/maci64/libmwlapack.dylib] 168df4397b0SStefano Zampini 169cdc6ee08SBarry 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 170cdc6ee08SBarry Smith against the 64 bit version, otherwise it use the 32 bit version 171df4397b0SStefano Zampini 172cdc6ee08SBarry 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 173df4397b0SStefano Zampini 174df4397b0SStefano Zampini External packages such as hypre, ML, SuperLU etc do not provide any support for passing 64 bit integers to BLAS/LAPACK so cannot 175cdc6ee08SBarry 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 176cdc6ee08SBarry Smith these external libraries while using 64 bit integer BLAS/LAPACK. 177df4397b0SStefano Zampini 178df4397b0SStefano Zampini Level: intermediate 179df4397b0SStefano Zampini 180db781477SPatrick Sanan .seealso: `PetscMPIInt`, `PetscInt`, `PetscBLASIntCast()` 181df4397b0SStefano Zampini 182df4397b0SStefano Zampini M*/ 183df4397b0SStefano Zampini #if defined(PETSC_HAVE_64BIT_BLAS_INDICES) 184c93fae50SJacob Faibussowitsch # define PetscBLASInt_FMT PetscInt64_FMT 185df4397b0SStefano Zampini typedef PetscInt64 PetscBLASInt; 186df4397b0SStefano Zampini #else 187c93fae50SJacob Faibussowitsch # define PetscBLASInt_FMT "d" 188df4397b0SStefano Zampini typedef int PetscBLASInt; 189df4397b0SStefano Zampini #endif 190df4397b0SStefano Zampini 191eee0c0a6SToby Isaac /*MC 192eee0c0a6SToby Isaac PetscCuBLASInt - datatype used to represent 'int' parameters to cuBLAS/cuSOLVER functions. 193eee0c0a6SToby Isaac 194eee0c0a6SToby Isaac Notes: 195eee0c0a6SToby Isaac As of this writing PetscCuBLASInt is always the system `int`. 196eee0c0a6SToby Isaac 197*87497f52SBarry Smith `PetscErrorCode` `PetscCuBLASIntCast`(a,&b) checks if the given `PetscInt` a will fit in a `PetscCuBLASInt`, if not it 198*87497f52SBarry Smith generates a `PETSC_ERR_ARG_OUTOFRANGE` error 199eee0c0a6SToby Isaac 200eee0c0a6SToby Isaac Level: intermediate 201eee0c0a6SToby Isaac 202db781477SPatrick Sanan .seealso: `PetscBLASInt`, `PetscMPIInt`, `PetscInt`, `PetscCuBLASIntCast()` 203eee0c0a6SToby Isaac 204eee0c0a6SToby Isaac M*/ 205eee0c0a6SToby Isaac typedef int PetscCuBLASInt; 206eee0c0a6SToby Isaac 207df4397b0SStefano Zampini /*E 208b94d7dedSBarry Smith PetscBool - Logical variable. Actually an enum in C and a logical in Fortran. 209df4397b0SStefano Zampini 210df4397b0SStefano Zampini Level: beginner 211df4397b0SStefano Zampini 212df4397b0SStefano Zampini Developer Note: 213*87497f52SBarry Smith 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 214df4397b0SStefano Zampini boolean values. It is not easy to have a simple macro that that will work properly in all circumstances with all three mechanisms. 215df4397b0SStefano Zampini 216b94d7dedSBarry Smith .seealso: `PETSC_TRUE`, `PETSC_FALSE`, `PetscNot()`, `PetscBool3` 217df4397b0SStefano Zampini E*/ 218df4397b0SStefano Zampini typedef enum { PETSC_FALSE,PETSC_TRUE } PetscBool; 219df4397b0SStefano Zampini 220b94d7dedSBarry Smith /*E 221b94d7dedSBarry Smith PetscBool3 - Ternary logical variable. Actually an enum in C and a 4 byte integer in Fortran. 222b94d7dedSBarry Smith 223b94d7dedSBarry Smith Level: beginner 224b94d7dedSBarry Smith 225*87497f52SBarry Smith Note: 226b94d7dedSBarry Smith Should not be used with the if (flg) or if (!flg) syntax. 227b94d7dedSBarry Smith 228b94d7dedSBarry Smith .seealso: `PETSC_TRUE`, `PETSC_FALSE`, `PetscNot()`, `PETSC_BOOL3_TRUE`, `PETSC_BOOL3_FALSE`, `PETSC_BOOL3_UKNOWN` 229b94d7dedSBarry Smith E*/ 230b94d7dedSBarry Smith typedef enum { PETSC_BOOL3_FALSE,PETSC_BOOL3_TRUE, PETSC_BOOL3_UNKNOWN = -1 } PetscBool3; 231b94d7dedSBarry Smith 232b94d7dedSBarry Smith #define PetscBool3ToBool(a) ((a) == PETSC_BOOL3_TRUE ? PETSC_TRUE : PETSC_FALSE) 233b94d7dedSBarry Smith #define PetscBoolToBool3(a) ((a) == PETSC_TRUE ? PETSC_BOOL3_TRUE : PETSC_BOOL3_FALSE) 234b94d7dedSBarry Smith 235df4397b0SStefano Zampini /*MC 236*87497f52SBarry Smith PetscReal - PETSc type that represents a real number version of `PetscScalar` 237df4397b0SStefano Zampini 238df4397b0SStefano Zampini Notes: 239*87497f52SBarry Smith For MPI calls that require datatypes, use `MPIU_REAL` as the datatype for `PetscReal` and `MPIU_SUM`, `MPIU_MAX`, etc. for operations. 240*87497f52SBarry Smith They will automatically work correctly regardless of the size of `PetscReal`. 241df4397b0SStefano Zampini 242*87497f52SBarry Smith See `PetscScalar` for details on how to ./configure the size of `PetscReal`. 243df4397b0SStefano Zampini 244df4397b0SStefano Zampini Level: beginner 245df4397b0SStefano Zampini 246db781477SPatrick Sanan .seealso: `PetscScalar`, `PetscComplex`, `PetscInt`, `MPIU_REAL`, `MPIU_SCALAR`, `MPIU_COMPLEX`, `MPIU_INT` 247df4397b0SStefano Zampini M*/ 248df4397b0SStefano Zampini 249df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 250df4397b0SStefano Zampini typedef float PetscReal; 251df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL_DOUBLE) 252df4397b0SStefano Zampini typedef double PetscReal; 253df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 254df4397b0SStefano Zampini # if defined(__cplusplus) 255df4397b0SStefano Zampini extern "C" { 256df4397b0SStefano Zampini # endif 257df4397b0SStefano Zampini # include <quadmath.h> 258df4397b0SStefano Zampini # if defined(__cplusplus) 259df4397b0SStefano Zampini } 260df4397b0SStefano Zampini # endif 261df4397b0SStefano Zampini typedef __float128 PetscReal; 262df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FP16) 263df4397b0SStefano Zampini typedef __fp16 PetscReal; 264df4397b0SStefano Zampini #endif /* PETSC_USE_REAL_* */ 265df4397b0SStefano Zampini 266df4397b0SStefano Zampini /*MC 267*87497f52SBarry Smith PetscComplex - PETSc type that represents a complex number with precision matching that of `PetscReal`. 268df4397b0SStefano Zampini 269df4397b0SStefano Zampini Synopsis: 270df4397b0SStefano Zampini #include <petscsys.h> 271df4397b0SStefano Zampini PetscComplex number = 1. + 2.*PETSC_i; 272df4397b0SStefano Zampini 273df4397b0SStefano Zampini Notes: 274*87497f52SBarry Smith For MPI calls that require datatypes, use `MPIU_COMPLEX` as the datatype for `PetscComplex` and `MPIU_SUM` etc for operations. 275*87497f52SBarry Smith They will automatically work correctly regardless of the size of `PetscComplex`. 276df4397b0SStefano Zampini 277*87497f52SBarry Smith See PetscScalar for details on how to ./configure the size of `PetscReal` 278df4397b0SStefano Zampini 279df4397b0SStefano Zampini Complex numbers are automatically available if PETSc was able to find a working complex implementation 280df4397b0SStefano Zampini 281*87497f52SBarry Smith Petsc has a 'fix' for complex numbers to support expressions such as std::complex<PetscReal> + `PetscInt`, which are not supported by the standard 282a966371cSJunchao Zhang C++ library, but are convenient for petsc users. If the C++ compiler is able to compile code in petsccxxcomplexfix.h (This is checked by 283a966371cSJunchao Zhang configure), we include petsccxxcomplexfix.h to provide this convenience. 284a966371cSJunchao Zhang 285*87497f52SBarry Smith If the fix causes conflicts, or one really does not want this fix for a particular C++ file, one can define `PETSC_SKIP_CXX_COMPLEX_FIX` 286a966371cSJunchao Zhang at the beginning of the C++ file to skip the fix. 287a966371cSJunchao Zhang 288df4397b0SStefano Zampini Level: beginner 289df4397b0SStefano Zampini 290db781477SPatrick Sanan .seealso: `PetscReal`, `PetscScalar`, `PetscComplex`, `PetscInt`, `MPIU_REAL`, `MPIU_SCALAR`, `MPIU_COMPLEX`, `MPIU_INT`, `PETSC_i` 291df4397b0SStefano Zampini M*/ 292df4397b0SStefano Zampini #if !defined(PETSC_SKIP_COMPLEX) 2937a19d461SSatish Balay # if defined(PETSC_CLANGUAGE_CXX) 2947a19d461SSatish Balay # if !defined(PETSC_USE_REAL___FP16) && !defined(PETSC_USE_REAL___FLOAT128) 2957a19d461SSatish Balay # if defined(__cplusplus) && defined(PETSC_HAVE_CXX_COMPLEX) /* enable complex for library code */ 296df4397b0SStefano Zampini # define PETSC_HAVE_COMPLEX 1 2977a19d461SSatish Balay # elif !defined(__cplusplus) && defined(PETSC_HAVE_C99_COMPLEX) && defined(PETSC_HAVE_CXX_COMPLEX) /* User code only - conditional on libary code complex support */ 2987a19d461SSatish Balay # define PETSC_HAVE_COMPLEX 1 2997a19d461SSatish Balay # endif 300450fc7c9SSatish Balay # elif defined(PETSC_USE_REAL___FLOAT128) && defined(PETSC_HAVE_C99_COMPLEX) 301450fc7c9SSatish Balay # define PETSC_HAVE_COMPLEX 1 3027a19d461SSatish Balay # endif 3037a19d461SSatish Balay # else /* !PETSC_CLANGUAGE_CXX */ 3047a19d461SSatish Balay # if !defined(PETSC_USE_REAL___FP16) 3057a19d461SSatish Balay # if !defined(__cplusplus) && defined(PETSC_HAVE_C99_COMPLEX) /* enable complex for library code */ 3067a19d461SSatish Balay # define PETSC_HAVE_COMPLEX 1 3077a19d461SSatish Balay # elif defined(__cplusplus) && defined(PETSC_HAVE_C99_COMPLEX) && defined(PETSC_HAVE_CXX_COMPLEX) /* User code only - conditional on libary code complex support */ 3087a19d461SSatish Balay # define PETSC_HAVE_COMPLEX 1 3097a19d461SSatish Balay # endif 3107a19d461SSatish Balay # endif 3117a19d461SSatish Balay # endif /* PETSC_CLANGUAGE_CXX */ 3127a19d461SSatish Balay #endif /* !PETSC_SKIP_COMPLEX */ 3137a19d461SSatish Balay 3147a19d461SSatish Balay #if defined(PETSC_HAVE_COMPLEX) 3157a19d461SSatish Balay #if defined(__cplusplus) /* C++ complex support */ 31611d22bbfSJunchao Zhang /* Locate a C++ complex template library */ 31711d22bbfSJunchao Zhang #if defined(PETSC_DESIRE_KOKKOS_COMPLEX) /* Defined in petscvec_kokkos.hpp for *.kokkos.cxx files */ 31811d22bbfSJunchao Zhang #define petsccomplexlib Kokkos 31911d22bbfSJunchao Zhang #include <Kokkos_Complex.hpp> 320022afdc5SJed Brown #elif defined(__CUDACC__) || defined(__HIPCC__) 321df4397b0SStefano Zampini #define petsccomplexlib thrust 322df4397b0SStefano Zampini #include <thrust/complex.h> 323450fc7c9SSatish Balay #elif defined(PETSC_USE_REAL___FLOAT128) 324450fc7c9SSatish Balay #include <complex.h> 325df4397b0SStefano Zampini #else 326df4397b0SStefano Zampini #define petsccomplexlib std 327df4397b0SStefano Zampini #include <complex> 328df4397b0SStefano Zampini #endif 32911d22bbfSJunchao Zhang 33011d22bbfSJunchao Zhang /* Define PetscComplex based on the precision */ 331df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 332df4397b0SStefano Zampini typedef petsccomplexlib::complex<float> PetscComplex; 333df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL_DOUBLE) 334df4397b0SStefano Zampini typedef petsccomplexlib::complex<double> PetscComplex; 335df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 336450fc7c9SSatish Balay typedef __complex128 PetscComplex; 33711d22bbfSJunchao Zhang #endif 33811d22bbfSJunchao Zhang 339a966371cSJunchao Zhang /* Include a PETSc C++ complex 'fix'. Check PetscComplex manual page for details */ 340a966371cSJunchao Zhang #if defined(PETSC_HAVE_CXX_COMPLEX_FIX) && !defined(PETSC_SKIP_CXX_COMPLEX_FIX) 34139829747SLisandro Dalcin #include <petsccxxcomplexfix.h> 34211d22bbfSJunchao Zhang #endif 3437a19d461SSatish Balay #else /* c99 complex support */ 344df4397b0SStefano Zampini #include <complex.h> 345df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) || defined(PETSC_USE_REAL___FP16) 346df4397b0SStefano Zampini typedef float _Complex PetscComplex; 347df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL_DOUBLE) 348df4397b0SStefano Zampini typedef double _Complex PetscComplex; 349df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 350df4397b0SStefano Zampini typedef __complex128 PetscComplex; 351df4397b0SStefano Zampini #endif /* PETSC_USE_REAL_* */ 3527a19d461SSatish Balay #endif /* !__cplusplus */ 3537a19d461SSatish Balay #endif /* PETSC_HAVE_COMPLEX */ 354df4397b0SStefano Zampini 355df4397b0SStefano Zampini /*MC 356df4397b0SStefano Zampini PetscScalar - PETSc type that represents either a double precision real number, a double precision 357df4397b0SStefano Zampini complex number, a single precision real number, a __float128 real or complex or a __fp16 real - if the code is configured 358df4397b0SStefano Zampini with --with-scalar-type=real,complex --with-precision=single,double,__float128,__fp16 359df4397b0SStefano Zampini 360df4397b0SStefano Zampini Notes: 361*87497f52SBarry Smith For MPI calls that require datatypes, use `MPIU_SCALAR` as the datatype for `PetscScalar` and `MPIU_SUM`, etc for operations. They will automatically work correctly regardless of the size of `PetscScalar`. 362df4397b0SStefano Zampini 363df4397b0SStefano Zampini Level: beginner 364df4397b0SStefano Zampini 365db781477SPatrick Sanan .seealso: `PetscReal`, `PetscComplex`, `PetscInt`, `MPIU_REAL`, `MPIU_SCALAR`, `MPIU_COMPLEX`, `MPIU_INT`, `PetscRealPart()`, `PetscImaginaryPart()` 366df4397b0SStefano Zampini M*/ 367df4397b0SStefano Zampini 3687a19d461SSatish Balay #if defined(PETSC_USE_COMPLEX) && defined(PETSC_HAVE_COMPLEX) 369df4397b0SStefano Zampini typedef PetscComplex PetscScalar; 370df4397b0SStefano Zampini #else /* PETSC_USE_COMPLEX */ 371df4397b0SStefano Zampini typedef PetscReal PetscScalar; 372df4397b0SStefano Zampini #endif /* PETSC_USE_COMPLEX */ 373df4397b0SStefano Zampini 374df4397b0SStefano Zampini /*E 375*87497f52SBarry Smith PetscCopyMode - Determines how an array or `PetscObject` passed to certain functions is copied or retained by the aggregate `PetscObject` 376df4397b0SStefano Zampini 377df4397b0SStefano Zampini Level: beginner 378df4397b0SStefano Zampini 3795d80c0bfSVaclav Hapla For the array input: 380*87497f52SBarry Smith $ `PETSC_COPY_VALUES` - the array values are copied into new space, the user is free to reuse or delete the passed in array 381*87497f52SBarry Smith $ `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 382df4397b0SStefano Zampini $ delete the array. The array MUST have been obtained with PetscMalloc(). Hence this mode cannot be used in Fortran. 383*87497f52SBarry Smith $ `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 3845d80c0bfSVaclav Hapla $ the array but the user must delete the array after the object is destroyed. 3855d80c0bfSVaclav Hapla 3865d80c0bfSVaclav Hapla For the PetscObject input: 387*87497f52SBarry Smith $ `PETSC_COPY_VALUES` - the input `PetscObject` is cloned into the aggregate `PetscObject`; the user is free to reuse/modify the input `PetscObject` without side effects. 388*87497f52SBarry Smith $ `PETSC_OWN_POINTER` - the input `PetscObject` is referenced by pointer (with reference count), thus should not be modified by the user. (Modification may cause errors or unintended side-effects in this or a future version of PETSc.) 389*87497f52SBarry Smith For either case above, the input `PetscObject` should be destroyed by the user when no longer needed (the aggregate object increases its reference count). 390*87497f52SBarry Smith $ `PETSC_USE_POINTER` - invalid for `PetscObject` inputs. 391df4397b0SStefano Zampini 392df4397b0SStefano Zampini E*/ 393df4397b0SStefano Zampini typedef enum {PETSC_COPY_VALUES, PETSC_OWN_POINTER, PETSC_USE_POINTER} PetscCopyMode; 394df4397b0SStefano Zampini 395df4397b0SStefano Zampini /*MC 396*87497f52SBarry Smith PETSC_FALSE - False value of `PetscBool` 397df4397b0SStefano Zampini 398df4397b0SStefano Zampini Level: beginner 399df4397b0SStefano Zampini 400df4397b0SStefano Zampini Note: 401df4397b0SStefano Zampini Zero integer 402df4397b0SStefano Zampini 403*87497f52SBarry Smith .seealso: `PetscBool`, `PetscBool3`, `PETSC_TRUE` 404df4397b0SStefano Zampini M*/ 405df4397b0SStefano Zampini 406df4397b0SStefano Zampini /*MC 407*87497f52SBarry Smith PETSC_TRUE - True value of `PetscBool` 408df4397b0SStefano Zampini 409df4397b0SStefano Zampini Level: beginner 410df4397b0SStefano Zampini 411df4397b0SStefano Zampini Note: 412df4397b0SStefano Zampini Nonzero integer 413df4397b0SStefano Zampini 414*87497f52SBarry Smith .seealso: `PetscBool`, `PetscBool3`, `PETSC_FALSE` 415df4397b0SStefano Zampini M*/ 416df4397b0SStefano Zampini 417df4397b0SStefano Zampini /*MC 418df4397b0SStefano Zampini PetscLogDouble - Used for logging times 419df4397b0SStefano Zampini 420df4397b0SStefano Zampini Notes: 421df4397b0SStefano Zampini Contains double precision numbers that are not used in the numerical computations, but rather in logging, timing etc. 422df4397b0SStefano Zampini 4236b3bf505SPatrick Sanan Level: developer 4246b3bf505SPatrick Sanan 425df4397b0SStefano Zampini M*/ 426df4397b0SStefano Zampini typedef double PetscLogDouble; 427df4397b0SStefano Zampini 428df4397b0SStefano Zampini /*E 429df4397b0SStefano Zampini PetscDataType - Used for handling different basic data types. 430df4397b0SStefano Zampini 431df4397b0SStefano Zampini Level: beginner 432df4397b0SStefano Zampini 433df4397b0SStefano Zampini Notes: 434*87497f52SBarry Smith Use of this should be avoided if one can directly use `MPI_Datatype` instead. 435df4397b0SStefano Zampini 436*87497f52SBarry Smith `PETSC_INT` is the datatype for a `PetscInt`, regardless of whether it is 4 or 8 bytes. 437*87497f52SBarry Smith `PETSC_REAL`, `PETSC_COMPLEX` and `PETSC_SCALAR` are the datatypes for `PetscReal`, `PetscComplex` and `PetscScalar`, regardless of their sizes. 438277f51e8SBarry Smith 439*87497f52SBarry Smith Developer Notes: 440df4397b0SStefano Zampini It would be nice if we could always just use MPI Datatypes, why can we not? 441df4397b0SStefano Zampini 442*87497f52SBarry Smith If you change any values in `PetscDatatype` make sure you update their usage in 443fa131761SBarry Smith share/petsc/matlab/PetscBagRead.m and share/petsc/matlab/@PetscOpenSocket/read/write.m 44491e29162SBarry Smith 445277f51e8SBarry Smith TODO: Add PETSC_INT32 and remove use of improper PETSC_ENUM 446277f51e8SBarry Smith 447db781477SPatrick Sanan .seealso: `PetscBinaryRead()`, `PetscBinaryWrite()`, `PetscDataTypeToMPIDataType()`, 448db781477SPatrick Sanan `PetscDataTypeGetSize()` 449df4397b0SStefano Zampini 450df4397b0SStefano Zampini E*/ 4512ebc710fSVaclav Hapla typedef enum {PETSC_DATATYPE_UNKNOWN = 0, 4522ebc710fSVaclav Hapla PETSC_DOUBLE = 1, PETSC_COMPLEX = 2, PETSC_LONG = 3, PETSC_SHORT = 4, PETSC_FLOAT = 5, 45349de050aSVaclav Hapla PETSC_CHAR = 6, PETSC_BIT_LOGICAL = 7, PETSC_ENUM = 8, PETSC_BOOL = 9, PETSC___FLOAT128 = 10, 454ea127267SVaclav Hapla PETSC_OBJECT = 11, PETSC_FUNCTION = 12, PETSC_STRING = 13, PETSC___FP16 = 14, PETSC_STRUCT = 15, 4559e3e4c22SLisandro Dalcin PETSC_INT = 16, PETSC_INT64 = 17} PetscDataType; 456df4397b0SStefano Zampini 457df4397b0SStefano Zampini #if defined(PETSC_USE_REAL_SINGLE) 458df4397b0SStefano Zampini # define PETSC_REAL PETSC_FLOAT 4595117d392SLisandro Dalcin #elif defined(PETSC_USE_REAL_DOUBLE) 4605117d392SLisandro Dalcin # define PETSC_REAL PETSC_DOUBLE 461df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FLOAT128) 462df4397b0SStefano Zampini # define PETSC_REAL PETSC___FLOAT128 463df4397b0SStefano Zampini #elif defined(PETSC_USE_REAL___FP16) 464df4397b0SStefano Zampini # define PETSC_REAL PETSC___FP16 465df4397b0SStefano Zampini #else 466df4397b0SStefano Zampini # define PETSC_REAL PETSC_DOUBLE 467df4397b0SStefano Zampini #endif 4685117d392SLisandro Dalcin 4695117d392SLisandro Dalcin #if defined(PETSC_USE_COMPLEX) 4705117d392SLisandro Dalcin # define PETSC_SCALAR PETSC_COMPLEX 4715117d392SLisandro Dalcin #else 4725117d392SLisandro Dalcin # define PETSC_SCALAR PETSC_REAL 4735117d392SLisandro Dalcin #endif 4745117d392SLisandro Dalcin 475df4397b0SStefano Zampini #define PETSC_FORTRANADDR PETSC_LONG 476df4397b0SStefano Zampini 477df4397b0SStefano Zampini /*S 478df4397b0SStefano Zampini PetscToken - 'Token' used for managing tokenizing strings 479df4397b0SStefano Zampini 480df4397b0SStefano Zampini Level: intermediate 481df4397b0SStefano Zampini 482db781477SPatrick Sanan .seealso: `PetscTokenCreate()`, `PetscTokenFind()`, `PetscTokenDestroy()` 483df4397b0SStefano Zampini S*/ 484df4397b0SStefano Zampini typedef struct _p_PetscToken* PetscToken; 485df4397b0SStefano Zampini 486df4397b0SStefano Zampini /*S 487*87497f52SBarry Smith PetscObject - any PETSc object, `PetscViewer`, `Mat`, `Vec`, `KSP` etc 488df4397b0SStefano Zampini 489df4397b0SStefano Zampini Level: beginner 490df4397b0SStefano Zampini 491*87497f52SBarry Smith Notes: 492df4397b0SStefano Zampini This is the base class from which all PETSc objects are derived from. 493df4397b0SStefano Zampini 494*87497f52SBarry Smith In certain situations one can cast an object, for example a `Vec`, to a `PetscObject` with (`PetscObject`)vec 495*87497f52SBarry Smith 496db781477SPatrick Sanan .seealso: `PetscObjectDestroy()`, `PetscObjectView()`, `PetscObjectGetName()`, `PetscObjectSetName()`, `PetscObjectReference()`, `PetscObjectDereference()` 497df4397b0SStefano Zampini S*/ 498df4397b0SStefano Zampini typedef struct _p_PetscObject* PetscObject; 499df4397b0SStefano Zampini 500df4397b0SStefano Zampini /*MC 501*87497f52SBarry Smith PetscObjectId - unique integer Id for a `PetscObject` 502df4397b0SStefano Zampini 503df4397b0SStefano Zampini Level: developer 504df4397b0SStefano Zampini 505*87497f52SBarry Smith Note: 506*87497f52SBarry Smith Unlike pointer values, object ids are never reused so one may save a `PetscObjectId` and compare it to one obtained later from a `PetscObject` to determine 507*87497f52SBarry Smith if the objects are the same. Never compare two object pointer values. 508df4397b0SStefano Zampini 509db781477SPatrick Sanan .seealso: `PetscObjectState`, `PetscObjectGetId()` 510df4397b0SStefano Zampini M*/ 511df4397b0SStefano Zampini typedef PetscInt64 PetscObjectId; 512df4397b0SStefano Zampini 513df4397b0SStefano Zampini /*MC 514*87497f52SBarry Smith PetscObjectState - integer state for a `PetscObject` 515df4397b0SStefano Zampini 516df4397b0SStefano Zampini Level: developer 517df4397b0SStefano Zampini 518df4397b0SStefano Zampini Notes: 519df4397b0SStefano Zampini Object state is always-increasing and (for objects that track state) can be used to determine if an object has 520df4397b0SStefano 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. 521df4397b0SStefano Zampini 522db781477SPatrick Sanan .seealso: `PetscObjectId`, `PetscObjectStateGet()`, `PetscObjectStateIncrease()`, `PetscObjectStateSet()` 523df4397b0SStefano Zampini M*/ 524df4397b0SStefano Zampini typedef PetscInt64 PetscObjectState; 525df4397b0SStefano Zampini 526df4397b0SStefano Zampini /*S 527df4397b0SStefano Zampini PetscFunctionList - Linked list of functions, possibly stored in dynamic libraries, accessed 528df4397b0SStefano Zampini by string name 529df4397b0SStefano Zampini 530df4397b0SStefano Zampini Level: advanced 531df4397b0SStefano Zampini 532db781477SPatrick Sanan .seealso: `PetscFunctionListAdd()`, `PetscFunctionListDestroy()` 533df4397b0SStefano Zampini S*/ 534df4397b0SStefano Zampini typedef struct _n_PetscFunctionList *PetscFunctionList; 535df4397b0SStefano Zampini 536df4397b0SStefano Zampini /*E 537df4397b0SStefano Zampini PetscFileMode - Access mode for a file. 538df4397b0SStefano Zampini 539df4397b0SStefano Zampini Level: beginner 540df4397b0SStefano Zampini 541*87497f52SBarry Smith $ `FILE_MODE_UNDEFINED` - initial invalid value 542*87497f52SBarry Smith $ `FILE_MODE_READ` - open a file at its beginning for reading 543*87497f52SBarry Smith $ `FILE_MODE_WRITE` - open a file at its beginning for writing (will create if the file does not exist) 544*87497f52SBarry Smith $ `FILE_MODE_APPEND` - open a file at end for writing 545*87497f52SBarry Smith $ `FILE_MODE_UPDATE` - open a file for updating, meaning for reading and writing 546*87497f52SBarry Smith $ `FILE_MODE_APPEND_UPDATE` - open a file for updating, meaning for reading and writing, at the end 547df4397b0SStefano Zampini 548db781477SPatrick Sanan .seealso: `PetscViewerFileSetMode()` 549df4397b0SStefano Zampini E*/ 5507e4fd573SVaclav Hapla typedef enum {FILE_MODE_UNDEFINED=-1, FILE_MODE_READ=0, FILE_MODE_WRITE, FILE_MODE_APPEND, FILE_MODE_UPDATE, FILE_MODE_APPEND_UPDATE} PetscFileMode; 551df4397b0SStefano Zampini 552df4397b0SStefano Zampini typedef void* PetscDLHandle; 553df4397b0SStefano Zampini typedef enum {PETSC_DL_DECIDE=0,PETSC_DL_NOW=1,PETSC_DL_LOCAL=2} PetscDLMode; 554df4397b0SStefano Zampini 555df4397b0SStefano Zampini /*S 5567243573dSPierre Jolivet PetscObjectList - Linked list of PETSc objects, each accessible by string name 557df4397b0SStefano Zampini 558df4397b0SStefano Zampini Level: developer 559df4397b0SStefano Zampini 560*87497f52SBarry Smith Note: 561*87497f52SBarry Smith Used by `PetscObjectCompose()` and `PetscObjectQuery()` 562df4397b0SStefano Zampini 563db781477SPatrick Sanan .seealso: `PetscObjectListAdd()`, `PetscObjectListDestroy()`, `PetscObjectListFind()`, `PetscObjectCompose()`, `PetscObjectQuery()`, `PetscFunctionList` 564df4397b0SStefano Zampini S*/ 565df4397b0SStefano Zampini typedef struct _n_PetscObjectList *PetscObjectList; 566df4397b0SStefano Zampini 567df4397b0SStefano Zampini /*S 568df4397b0SStefano Zampini PetscDLLibrary - Linked list of dynamics libraries to search for functions 569df4397b0SStefano Zampini 570df4397b0SStefano Zampini Level: advanced 571df4397b0SStefano Zampini 572db781477SPatrick Sanan .seealso: `PetscDLLibraryOpen()` 573df4397b0SStefano Zampini S*/ 574df4397b0SStefano Zampini typedef struct _n_PetscDLLibrary *PetscDLLibrary; 575df4397b0SStefano Zampini 576df4397b0SStefano Zampini /*S 577df4397b0SStefano Zampini PetscContainer - Simple PETSc object that contains a pointer to any required data 578df4397b0SStefano Zampini 579df4397b0SStefano Zampini Level: advanced 580df4397b0SStefano Zampini 581*87497f52SBarry Smith Note: 582*87497f52SBarry Smith This is useful to attach arbitrary data to a `PetscObject` with `PetscObjectCompose()` and `PetscObjectQuery()` 583*87497f52SBarry Smith 584*87497f52SBarry Smith .seealso: `PetscObject`, `PetscContainerCreate()`, `PetscObjectCompose()`, `PetscObjectQuery()` 585df4397b0SStefano Zampini S*/ 586df4397b0SStefano Zampini typedef struct _p_PetscContainer* PetscContainer; 587df4397b0SStefano Zampini 588df4397b0SStefano Zampini /*S 589df4397b0SStefano Zampini PetscRandom - Abstract PETSc object that manages generating random numbers 590df4397b0SStefano Zampini 591df4397b0SStefano Zampini Level: intermediate 592df4397b0SStefano Zampini 593db781477SPatrick Sanan .seealso: `PetscRandomCreate()`, `PetscRandomGetValue()`, `PetscRandomType` 594df4397b0SStefano Zampini S*/ 595df4397b0SStefano Zampini typedef struct _p_PetscRandom* PetscRandom; 596df4397b0SStefano Zampini 597df4397b0SStefano Zampini /* 598df4397b0SStefano Zampini In binary files variables are stored using the following lengths, 599df4397b0SStefano Zampini regardless of how they are stored in memory on any one particular 600df4397b0SStefano Zampini machine. Use these rather then sizeof() in computing sizes for 601df4397b0SStefano Zampini PetscBinarySeek(). 602df4397b0SStefano Zampini */ 603df4397b0SStefano Zampini #define PETSC_BINARY_INT_SIZE (32/8) 604df4397b0SStefano Zampini #define PETSC_BINARY_FLOAT_SIZE (32/8) 605df4397b0SStefano Zampini #define PETSC_BINARY_CHAR_SIZE (8/8) 606df4397b0SStefano Zampini #define PETSC_BINARY_SHORT_SIZE (16/8) 607df4397b0SStefano Zampini #define PETSC_BINARY_DOUBLE_SIZE (64/8) 608df4397b0SStefano Zampini #define PETSC_BINARY_SCALAR_SIZE sizeof(PetscScalar) 609df4397b0SStefano Zampini 610df4397b0SStefano Zampini /*E 611*87497f52SBarry Smith PetscBinarySeekType - argument to `PetscBinarySeek()` 612df4397b0SStefano Zampini 613df4397b0SStefano Zampini Level: advanced 614df4397b0SStefano Zampini 615db781477SPatrick Sanan .seealso: `PetscBinarySeek()`, `PetscBinarySynchronizedSeek()` 616df4397b0SStefano Zampini E*/ 617df4397b0SStefano Zampini typedef enum {PETSC_BINARY_SEEK_SET = 0,PETSC_BINARY_SEEK_CUR = 1,PETSC_BINARY_SEEK_END = 2} PetscBinarySeekType; 618df4397b0SStefano Zampini 619df4397b0SStefano Zampini /*E 620df4397b0SStefano Zampini PetscBuildTwoSidedType - algorithm for setting up two-sided communication 621df4397b0SStefano Zampini 622*87497f52SBarry Smith $ `PETSC_BUILDTWOSIDED_ALLREDUCE` - classical algorithm using an MPI_Allreduce with 623df4397b0SStefano Zampini $ a buffer of length equal to the communicator size. Not memory-scalable due to 624df4397b0SStefano Zampini $ the large reduction size. Requires only MPI-1. 625*87497f52SBarry Smith $ `PETSC_BUILDTWOSIDED_IBARRIER` - nonblocking algorithm based on MPI_Issend and MPI_Ibarrier. 626df4397b0SStefano Zampini $ Proved communication-optimal in Hoefler, Siebert, and Lumsdaine (2010). Requires MPI-3. 627*87497f52SBarry Smith $ `PETSC_BUILDTWOSIDED_REDSCATTER` - similar to above, but use more optimized function 628df4397b0SStefano Zampini $ that only communicates the part of the reduction that is necessary. Requires MPI-2. 629df4397b0SStefano Zampini 630df4397b0SStefano Zampini Level: developer 631df4397b0SStefano Zampini 632db781477SPatrick Sanan .seealso: `PetscCommBuildTwoSided()`, `PetscCommBuildTwoSidedSetType()`, `PetscCommBuildTwoSidedGetType()` 633df4397b0SStefano Zampini E*/ 634df4397b0SStefano Zampini typedef enum { 635df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_NOTSET = -1, 636df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_ALLREDUCE = 0, 637df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_IBARRIER = 1, 638df4397b0SStefano Zampini PETSC_BUILDTWOSIDED_REDSCATTER = 2 639df4397b0SStefano Zampini /* Updates here must be accompanied by updates in finclude/petscsys.h and the string array in mpits.c */ 640df4397b0SStefano Zampini } PetscBuildTwoSidedType; 641df4397b0SStefano Zampini 642aa2d33c1SMatthew G. Knepley /* NOTE: If you change this, you must also change the values in src/vec/f90-mod/petscvec.h */ 643df4397b0SStefano Zampini /*E 644df4397b0SStefano Zampini InsertMode - Whether entries are inserted or added into vectors or matrices 645df4397b0SStefano Zampini 646df4397b0SStefano Zampini Level: beginner 647df4397b0SStefano Zampini 648db781477SPatrick Sanan .seealso: `VecSetValues()`, `MatSetValues()`, `VecSetValue()`, `VecSetValuesBlocked()`, 649db781477SPatrick Sanan `VecSetValuesLocal()`, `VecSetValuesBlockedLocal()`, `MatSetValuesBlocked()`, 650db781477SPatrick Sanan `MatSetValuesBlockedLocal()`, `MatSetValuesLocal()`, `VecScatterBegin()`, `VecScatterEnd()` 651df4397b0SStefano Zampini E*/ 6520c91829fSJunchao 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; 653df4397b0SStefano Zampini 654df4397b0SStefano Zampini /*MC 655df4397b0SStefano Zampini INSERT_VALUES - Put a value into a vector or matrix, overwrites any previous value 656df4397b0SStefano Zampini 657df4397b0SStefano Zampini Level: beginner 658df4397b0SStefano Zampini 659db781477SPatrick Sanan .seealso: `InsertMode`, `VecSetValues()`, `MatSetValues()`, `VecSetValue()`, `VecSetValuesBlocked()`, 660db781477SPatrick Sanan `VecSetValuesLocal()`, `VecSetValuesBlockedLocal()`, `MatSetValuesBlocked()`, `ADD_VALUES`, 661db781477SPatrick Sanan `MatSetValuesBlockedLocal()`, `MatSetValuesLocal()`, `VecScatterBegin()`, `VecScatterEnd()`, `MAX_VALUES` 662df4397b0SStefano Zampini 663df4397b0SStefano Zampini M*/ 664df4397b0SStefano Zampini 665df4397b0SStefano Zampini /*MC 666df4397b0SStefano Zampini ADD_VALUES - Adds a value into a vector or matrix, if there previously was no value, just puts the 667df4397b0SStefano Zampini value into that location 668df4397b0SStefano Zampini 669df4397b0SStefano Zampini Level: beginner 670df4397b0SStefano Zampini 671db781477SPatrick Sanan .seealso: `InsertMode`, `VecSetValues()`, `MatSetValues()`, `VecSetValue()`, `VecSetValuesBlocked()`, 672db781477SPatrick Sanan `VecSetValuesLocal()`, `VecSetValuesBlockedLocal()`, `MatSetValuesBlocked()`, `INSERT_VALUES`, 673db781477SPatrick Sanan `MatSetValuesBlockedLocal()`, `MatSetValuesLocal()`, `VecScatterBegin()`, `VecScatterEnd()`, `MAX_VALUES` 674df4397b0SStefano Zampini 675df4397b0SStefano Zampini M*/ 676df4397b0SStefano Zampini 677df4397b0SStefano Zampini /*MC 678df4397b0SStefano Zampini MAX_VALUES - Puts the maximum of the scattered/gathered value and the current value into each location 679df4397b0SStefano Zampini 680df4397b0SStefano Zampini Level: beginner 681df4397b0SStefano Zampini 682db781477SPatrick Sanan .seealso: `InsertMode`, `VecScatterBegin()`, `VecScatterEnd()`, `ADD_VALUES`, `INSERT_VALUES` 683df4397b0SStefano Zampini 684df4397b0SStefano Zampini M*/ 685df4397b0SStefano Zampini 686421aa1e7SJunchao Zhang /*MC 687421aa1e7SJunchao Zhang MIN_VALUES - Puts the minimal of the scattered/gathered value and the current value into each location 688421aa1e7SJunchao Zhang 689421aa1e7SJunchao Zhang Level: beginner 690421aa1e7SJunchao Zhang 691db781477SPatrick Sanan .seealso: `InsertMode`, `VecScatterBegin()`, `VecScatterEnd()`, `ADD_VALUES`, `INSERT_VALUES` 692421aa1e7SJunchao Zhang 693421aa1e7SJunchao Zhang M*/ 694421aa1e7SJunchao Zhang 695df4397b0SStefano Zampini /*S 696df4397b0SStefano Zampini PetscSubcomm - A decomposition of an MPI communicator into subcommunicators 697df4397b0SStefano Zampini 698df4397b0SStefano Zampini Notes: 699*87497f52SBarry Smith After a call to `PetscSubcommSetType()`, `PetscSubcommSetTypeGeneral()`, or `PetscSubcommSetFromOptions()` one may call 700*87497f52SBarry Smith $ `PetscSubcommChild()` returns the associated subcommunicator on this process 701*87497f52SBarry Smith $ `PetscSubcommContiguousParent()` returns a parent communitor but with all child of the same subcommunicator having contiguous rank 702df4397b0SStefano Zampini 703df4397b0SStefano Zampini Sample Usage: 704*87497f52SBarry Smith .vb 705*87497f52SBarry Smith `PetscSubcommCreate()` 706*87497f52SBarry Smith `PetscSubcommSetNumber()` 707*87497f52SBarry Smith `PetscSubcommSetType`(`PETSC_SUBCOMM_INTERLACED`); 708*87497f52SBarry Smith ccomm = `PetscSubcommChild()` 709*87497f52SBarry Smith `PetscSubcommDestroy()` 710*87497f52SBarry Smith .ve 711df4397b0SStefano Zampini 712df4397b0SStefano Zampini Level: advanced 713df4397b0SStefano Zampini 714df4397b0SStefano Zampini Notes: 715*87497f52SBarry Smith $ `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 716*87497f52SBarry Smith $ `PETSC_SUBCOMM_CONTIGUOUS` - each new communicator contains a set of process with contiguous ranks in the original MPI communicator 717*87497f52SBarry Smith $ `PETSC_SUBCOMM_INTERLACED` - each new communictor contains a set of processes equally far apart in rank from the others in that new communicator 718df4397b0SStefano Zampini 719df4397b0SStefano Zampini Example: Consider a communicator with six processes split into 3 subcommunicators. 720*87497f52SBarry Smith $ `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 721*87497f52SBarry Smith $ `PETSC_SUBCOMM_INTERLACED` - the first communicator contains rank 0,3, the second 1,4 and the third 2,5 722df4397b0SStefano Zampini 723*87497f52SBarry Smith Developer Note: 724*87497f52SBarry Smith This is used in objects such as `PCREDUNDANT` to manage the subcommunicators on which the redundant computations 725df4397b0SStefano Zampini are performed. 726df4397b0SStefano Zampini 727db781477SPatrick Sanan .seealso: `PetscSubcommCreate()`, `PetscSubcommSetNumber()`, `PetscSubcommSetType()`, `PetscSubcommView()`, `PetscSubcommSetFromOptions()` 728df4397b0SStefano Zampini 729df4397b0SStefano Zampini S*/ 730df4397b0SStefano Zampini typedef struct _n_PetscSubcomm* PetscSubcomm; 731df4397b0SStefano Zampini typedef enum {PETSC_SUBCOMM_GENERAL=0,PETSC_SUBCOMM_CONTIGUOUS=1,PETSC_SUBCOMM_INTERLACED=2} PetscSubcommType; 732df4397b0SStefano Zampini 733df4397b0SStefano Zampini /*S 734df4397b0SStefano Zampini PetscHeap - A simple class for managing heaps 735df4397b0SStefano Zampini 736df4397b0SStefano Zampini Level: intermediate 737df4397b0SStefano Zampini 738db781477SPatrick Sanan .seealso: `PetscHeapCreate()`, `PetscHeapAdd()`, `PetscHeapPop()`, `PetscHeapPeek()`, `PetscHeapStash()`, `PetscHeapUnstash()`, `PetscHeapView()`, `PetscHeapDestroy()` 739df4397b0SStefano Zampini S*/ 740df4397b0SStefano Zampini typedef struct _PetscHeap *PetscHeap; 741df4397b0SStefano Zampini 742df4397b0SStefano Zampini typedef struct _n_PetscShmComm* PetscShmComm; 743df4397b0SStefano Zampini typedef struct _n_PetscOmpCtrl* PetscOmpCtrl; 744df4397b0SStefano Zampini 745df4397b0SStefano Zampini /*S 746df4397b0SStefano Zampini PetscSegBuffer - a segmented extendable buffer 747df4397b0SStefano Zampini 748df4397b0SStefano Zampini Level: developer 749df4397b0SStefano Zampini 750db781477SPatrick Sanan .seealso: `PetscSegBufferCreate()`, `PetscSegBufferGet()`, `PetscSegBufferExtract()`, `PetscSegBufferDestroy()` 751df4397b0SStefano Zampini S*/ 752df4397b0SStefano Zampini typedef struct _n_PetscSegBuffer *PetscSegBuffer; 753df4397b0SStefano Zampini 754df4397b0SStefano Zampini typedef struct _n_PetscOptionsHelpPrinted *PetscOptionsHelpPrinted; 755df4397b0SStefano Zampini #endif 756