| /petsc/src/mat/impls/aij/seq/ |
| H A D | symtranspose.c | 15 PetscInt an = A->cmap->N, am = A->rmap->N; in MatTransposeSymbolic_SeqAIJ() local 20 PetscCall(PetscCalloc1(an + 1, &ati)); in MatTransposeSymbolic_SeqAIJ() 27 for (i = 0; i < an; i++) ati[i + 1] += ati[i]; in MatTransposeSymbolic_SeqAIJ() 30 PetscCall(PetscMalloc1(an, &atfill)); in MatTransposeSymbolic_SeqAIJ() 31 PetscCall(PetscArraycpy(atfill, ati, an)); in MatTransposeSymbolic_SeqAIJ() 43 …PetscCall(MatCreateSeqAIJWithArrays(PetscObjectComm((PetscObject)A), an, am, ati, atj, NULL, &At)); in MatTransposeSymbolic_SeqAIJ() 50 at->maxnz = ati[an]; in MatTransposeSymbolic_SeqAIJ() 60 PetscInt an = A->cmap->N, am = A->rmap->N; in MatTranspose_SeqAIJ() local 80 PetscCall(PetscCalloc1(an + 1, &ati)); in MatTranspose_SeqAIJ() 86 for (i = 0; i < an; i++) ati[i + 1] += ati[i]; in MatTranspose_SeqAIJ() [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex18.c | 138 …PetscScalar t0, tn, ts, te, tw, an, as, ae, aw, dn, ds, de, dw, fn = 0.0, fs = 0.0, fe = 0.0, fw… in FormFunction() local 187 an = 0.5 * (t0 + tn); in FormFunction() 188 dn = PetscPowScalar(an, beta); in FormFunction() 212 an = 0.5 * (t0 + tn); in FormFunction() 213 dn = PetscPowScalar(an, beta); in FormFunction() 238 an = 0.5 * (t0 + tn); in FormFunction() 239 dn = PetscPowScalar(an, beta); in FormFunction() 258 an = 0.5 * (t0 + tn); in FormFunction() 259 dn = PetscPowScalar(an, beta); in FormFunction() 297 PetscScalar dn, ds, de, dw, an, as, ae, aw, bn, bs, be, bw, gn, gs, ge, gw; in FormJacobian() local [all …]
|
| H A D | ex70.c | 386 PetscScalar an = s->hx / s->hy, anb = s->hx / (s->hy / 2); in StokesStencilLaplacian() local 392 vals[0] = -(ae + awb + asb + an); in StokesStencilLaplacian() 396 vals[2] = an; in StokesStencilLaplacian() 410 vals[1] = -(aeb + aw + asb + an); in StokesStencilLaplacian() 412 vals[2] = an; in StokesStencilLaplacian() 426 vals[1] = -(ae + awb + as + an); in StokesStencilLaplacian() 430 vals[3] = an; in StokesStencilLaplacian() 438 vals[2] = -(aeb + aw + as + an); in StokesStencilLaplacian() 440 vals[3] = an; in StokesStencilLaplacian() 446 vals[1] = -(ae + aw + asb + an); in StokesStencilLaplacian() [all …]
|
| /petsc/src/mat/utils/ |
| H A D | multequal.c | 14 PetscInt am, an, bm, bn, k; in MatMultEqual_Private() local 29 PetscCall(MatGetLocalSize(A, &am, &an)); in MatMultEqual_Private() 31 … an == bn, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Mat A,Mat B: local dim %" PetscInt_FMT " %" PetscI… in MatMultEqual_Private() 121 PetscInt am, an, bm, bn, cm, cn, k; in MatMatMultEqual_Private() local 138 PetscCall(MatGetLocalSize(A, &am, &an)); in MatMatMultEqual_Private() 142 PetscInt tt = an; in MatMatMultEqual_Private() 143 an = am; in MatMatMultEqual_Private() 151 …an == bm && am == cm && bn == cn, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Mat A, B, C local dim %" Pe… in MatMatMultEqual_Private() 442 PetscInt i, am, an, bm, bn, cm, cn; in MatProjMultEqual_Private() local 446 PetscCall(MatGetLocalSize(A, &am, &an)); in MatProjMultEqual_Private() [all …]
|
| /petsc/doc/developers/ |
| H A D | mrmanagement.md | 5 At any given time, at least one of the {any}`sec_core_developers` is an 19 …- Add an appropriate label, or label ``workflow::Waiting-on-Submitter`` and ask the submitter/assi… 20 * - Without an assignee 21 - Assign the submitter (if the MR is from a fork, also list an appropriate developer)
|
| H A D | style.md | 8 concrete realization of an abstract base class. For example, 21 letters. When they consist of several complete words, there is an 27 - have an appended `_Private` (for example, `StashValues_Private`) 29 - have an appended `_Subtype` (for example, `MatMultSeq_AIJ`). 41 function name, an underscore, and the name of the implementation, for 52 08. Options database keys are lower case, have an underscore between 57 an underscore in them unless they refer to another package that uses 58 an underscore, for example, `MATSOLVERSUPERLU_DIST`. 77 …to `PetscInt` from `PetscCount`. Since, in most configurations an array of `PetscCount` requires t… 78 of an array of `PetscInt` most index arrays (such as `ISGetIndices()` use `PetscInt`, [all …]
|
| /petsc/doc/changes/ |
| H A D | 322.md | 20 - Change `MPIU_Allreduce()` to always returns an MPI error code that should be checked with `PetscC… 94 preconditioners based on an alternating direction iteration, e.g., setting $M$ for 112 - Add support for `PETSC_DETERMINE` as an argument to `KSPSetTolerances()` to set the parameter bac… 118 - Add support for `PETSC_DETERMINE` as an argument to `SNESSetTolerances()` to set the parameter ba… 137 - Add support for `PETSC_DETERMINE` as an argument to `TSSetTolerances()` to set the parameter back… 139 - Add support for `PETSC_DETERMINE` as an argument to `TSSetMaxSteps()` and `TSSetMaxTime()` 151 - Add support for `PETSC_DETERMINE` as an argument to `TaoSetTolerances()` and `TaoSetConstraintTol… 160 - Add an additional return vector to `DMCreateMassMatrixLumped()` to retrieve the local mass lumping 171 - Add an argument to `DMPlexLabelCohesiveComplete()` in order to change behavior at surface boundary 213 - Add `PETSC_NULL_ENUM` to be used instead of `PETSC_NULL_INTEGER` when a pointer to an `enum` is e… [all …]
|
| H A D | 2028.md | 23 - `TSSetMonitor()` now takes an additional final optional argument 30 - Changed calling sequence of SNES monitoring routines to add an 33 - Added additional final argument to `SNESSetMonitor()` an 52 - Added additional final argument to `KSPSetMonitor()` an optional 125 PetscTruth * instead of an int\* 128 - Changed these routines to return an error code. extern int 142 - PetscFree() now always returns an error code that may be checked 144 user has no access to an error condition in free(), PetscFree() 167 Also you can put anything in \{anything} as long as it is an
|
| H A D | 2024.md | 60 - S`NESSetConvergenceHistory()` takes an additional argument; 78 - Added support for an alternative way to compute the 84 accepts an optional user-defined context variable. Likewise, 110 - Added an additional argument to `KSPSetResidualHistory()`; see 133 - Added an additional argument to `MatCopy()`; 135 - Added an additional argument to `MatDuplicate()`, 152 the manual change your current code to work, you only need add an 169 - Changed the `VecType` variable type from an enum type to a 188 - `DrawTensorContour()` has changed; its final argument is now an 200 - Changed `ViewerType` from an enum to a char \*. [all …]
|
| H A D | 2029.md | 23 - SNESGetJacobian() now has an additional optional final argument 25 - SNESGetFunction() now has an additional optional final argument 42 - MGSetLevels() now takes an additional optional final argument that 43 is an array of communicators to use on the SLES on each level 137 src/contrib/oberman/laplacian_q1 for an example of its use
|
| H A D | 231.md | 48 - -ksp_view, -ksp_monitor now take an optional filename argument 81 - -snes-view, -snes_monitor now take an optional filename argument 86 - -ts_view, -ts_monitor now take an optional filename argument 95 the linear system matrix can optionally fill also an alternate
|
| H A D | 2017.md | 66 - Added an interface to the SPAI preconditioner implementation of 75 an error if a new nonzero is generated in a sparse matrix. 86 routines linked into an application code. This can decrease the 157 pierr is an additional integer error code passed to the error 167 petsc/src/snes/examples/tutorials/ex3.c for an example. 176 - Now the option -log_summary takes [filename] as an optional 198 be used in place of the flag CONF. For an example of usage, see
|
| /petsc/src/sys/utils/ |
| H A D | sorti.c | 1038 PetscErrorCode PetscMergeIntArray(PetscInt an, const PetscInt aI[], PetscInt bn, const PetscInt bI[… in PetscMergeIntArray() argument 1044 PetscCall(PetscMalloc1(an + bn, L)); in PetscMergeIntArray() 1048 while (ak < an && bk < bn) { in PetscMergeIntArray() 1064 if (ak < an) { in PetscMergeIntArray() 1065 PetscCall(PetscArraycpy(L_ + k, aI + ak, an - ak)); in PetscMergeIntArray() 1066 k += (an - ak); in PetscMergeIntArray() 1103 PetscErrorCode PetscMergeIntArrayPair(PetscInt an, const PetscInt aI[], const PetscInt aJ[], PetscI… in PetscMergeIntArrayPair() argument 1110 n_ = an + bn; in PetscMergeIntArrayPair() 1117 while (ak < an && bk < bn) { in PetscMergeIntArrayPair() 1130 if (ak < an) { in PetscMergeIntArrayPair() [all …]
|
| /petsc/src/tao/leastsquares/tutorials/matlab/ |
| H A D | README | 4 24 - fails on destroy with an object in an invalid state; succeeds with nf_const = 50
|
| /petsc/doc/manualpages/doctext/ |
| H A D | doctextcommon.txt | 21 _Given an operation `A-B` (e.g. `A` = `write`, `B` = `read`) on an object or memory 31 memory regions may execute in an arbitrary order and offer no guarantee of sequential 38 considered to be in an undefined state. A well-defined state can only be restored by
|
| /petsc/src/snes/tests/ |
| H A D | ex20.c | 138 …PetscScalar t0, tn, ts, te, tw, an, as, ae, aw, dn, ds, de, dw, fn = 0.0, fs = 0.0, fe = 0.0, f… in FormFunction() local 190 an = 0.5 * (t0 + tn); in FormFunction() 191 dn = PetscPowScalar(an, beta); in FormFunction() 227 an = 0.5 * (t0 + tn); in FormFunction() 228 dn = PetscPowScalar(an, beta); in FormFunction() 275 an = 0.5 * (t0 + tn); in FormFunction() 276 dn = PetscPowScalar(an, beta); in FormFunction() 315 an = 0.5 * (t0 + tn); in FormFunction() 316 dn = PetscPowScalar(an, beta); in FormFunction() 392 an = 0.5 * (t0 + tn); in FormFunction() [all …]
|
| /petsc/src/ts/tests/output/ |
| H A D | ex26_arkimex_explicit_stage.out | 3 …]PETSC ERROR: This scheme requires an identity mass matrix, however the TSIFunctionFn you provided…
|
| /petsc/doc/manual/ |
| H A D | dmstag.md | 23 distributed in parallel over the ranks of an `MPI_Comm`. It is, however, 42 giving a rectangular local subdomain on each rank. This is extended by an element-wise stencil width 43 of *ghost* elements to create an atlas of overlapping patches. 70 strata have different numbers of points (e.g. in 1D there is an "extra" vertex on the right). This … 77 …representation of the unknown. The 4 local representations are shown, with an (elementwise) stenci… 84 The following snippet give an example of this usage. 124 Internally, each point is canonically associated with an element (top-level point (cell)). For purp… 125 regular-blocked storage, an element is grouped with lower-dimensional points left of, below ("down"… 132 with an element is 147 …) unknown in the global representation. See {numref}`figure_dmstag_numbering_local` for an example. [all …]
|
| H A D | getting_started.md | 26 Each class consists of an abstract interface (simply a set of calling 27 sequences corresponding to an abstract base class in C++) and an implementation for each algorithm … 157 `[-display name]` - start debugger only on encountering an error 174 argument `file` optionally indicates an alternative name for the PETSc 178 for runtime customization. The final argument, `help`, is an optional 206 {cite}`using-mpi` provides an excellent introduction to the 222 Most PETSc functions return a `PetscErrorCode`, an integer 223 indicating whether an error occurred during the call. The error code 224 is set to be nonzero if an error has been detected; otherwise, it is 226 return value, while for the Fortran version, each PETSc routine has an integer error variable as [all …]
|
| H A D | fortran.md | 133 - One cannot pass a scalar variable to a function expecting an array, {any}`sec_passarray`. 135 - One must pass pointers to arrays for arguments that output an array, for example `PetscScalar, po… 151 In the Fortran version, each PETSc routine has as its final argument an 153 nonzero if an error has been detected; otherwise, it is zero. For 176 a single value one can use the Fortran [] notation to pass the scalar as an array, for example 207 Where the code expects an array, then use `PETSC_NULL_XXX_ARRAY`. For example: 225 Finally when a subroutine returns a `PetscObject` through an argument, to check if it is `NULL` you… 253 For PETSc routine arguments that return an array of `PetscInt`, `PetscScalar`, `PetscReal` or of PE… 254 one passes in a pointer to an array and the PETSc routine returns an array containing the values. F… 378 Certain PETSc functions return a context in an argument, for example, `SNESGetApplicationContext()`… [all …]
|
| /petsc/share/petsc/saws/ |
| H A D | documentation | 49 …o loaded. Some are retrieved from the Argonne servers or other servers, so an internet connection … 55 …an object called "matInfo" since using a tree structure to store just a few options would be overk… 57 In Javascript, let's say we have an object 'Car' with property 'color'. 64 For arrays, deleting an element works the same. So, it is important to remember that the array leng…
|
| /petsc/src/mat/tests/ |
| H A D | ex104.c | 12 PetscInt i, M = 10, N = 5, j, nrows, ncols, am, an, rstart, rend; in main() local 134 PetscCall(MatGetLocalSize(A, &am, &an)); in main() 163 PetscInt am, an, aM, aN; in main() local 165 PetscCall(MatGetLocalSize(A, &am, &an)); in main() 167 …PetscCall(MatCreateDense(PetscObjectComm((PetscObject)A), PETSC_DECIDE, an, aM + 10, aN, NULL, &B)… in main()
|
| H A D | ex122.c | 12 PetscInt nza, am, an; in main() local 33 PetscCall(MatGetLocalSize(A, &am, &an)); in main()
|
| /petsc/src/sys/tests/linter/output/ |
| H A D | testValidFunctionDocStrings.patch | 283 -+ alpha, an alpha 284 ++ alpha - an alpha 289 + alpha, an alpha 412 . asdasd - an arg 415 . asdasd - an arg 487 + foo - the input description for an in-out param 493 + foo - the input description for an in-out param 504 - foo - the output description for an in-out param 509 -- foo - the output description for an in-out param 510 +. foo - the output description for an in-out param
|
| /petsc/ |
| H A D | CONTRIBUTING | 13 act of submitting a pull request or patch (with or without an explicit 14 Signed-off-by tag) will be understood as an affirmation of the 26 of my knowledge, is covered under an appropriate open-source
|