| /petsc/src/ts/tutorials/output/ |
| H A D | ex24_2.out | 1 0 t= 0.0e+00 dt= 1.0e-05 f= 2.0e+06 df= 2.6e+05 it=( 0, 0) 2 1 t= 1.0e-05 dt= 1.1e-05 f= 1.5e+06 df= 2.1e+05 it=( 3, 3) 3 2 t= 2.1e-05 dt= 1.5e-05 f= 1.2e+06 df= 1.7e+05 it=( 3, 3) 4 3 t= 3.6e-05 dt= 2.0e-05 f= 8.5e+05 df= 1.4e+05 it=( 3, 3) 5 4 t= 5.6e-05 dt= 2.9e-05 f= 6.0e+05 df= 1.0e+05 it=( 3, 3) 6 5 t= 8.6e-05 dt= 4.4e-05 f= 4.0e+05 df= 7.7e+04 it=( 3, 3) 7 6 t= 1.3e-04 dt= 6.9e-05 f= 2.4e+05 df= 5.4e+04 it=( 3, 3) 8 7 t= 2.0e-04 dt= 1.2e-04 f= 1.4e+05 df= 3.5e+04 it=( 3, 3) 9 8 t= 3.2e-04 dt= 2.1e-04 f= 6.9e+04 df= 2.1e+04 it=( 3, 3) 10 9 t= 5.3e-04 dt= 4.3e-04 f= 3.0e+04 df= 1.1e+04 it=( 4, 4) [all …]
|
| H A D | ex24_1.out | 1 0 t=0.0000e+00 dt=1.0000e-03 f=1.3005e+05 df=7.5724e+04 it=( 0, 0) 2 1 t=1.0000e-03 dt=1.1000e-03 f=2.5359e+04 df=2.2799e+04 it=( 1, 1) 3 2 t=2.1000e-03 dt=3.9248e-03 f=5.5935e+03 df=7.0287e+03 it=( 1, 1) 4 3 t=6.0248e-03 dt=1.4459e-02 f=1.5992e+03 df=2.0987e+03 it=( 1, 1) 5 4 t=2.0484e-02 dt=5.0883e-02 f=6.4634e+02 df=6.5600e+02 it=( 1, 1) 6 5 t=7.1367e-02 dt=6.6206e-02 f=1.9371e+02 df=5.5459e+02 it=( 1, 1) 7 6 t=1.3757e-01 dt=2.0251e-01 f=5.3877e+01 df=1.9944e+02 it=( 1, 1) 8 7 t=3.4008e-01 dt=6.1534e-01 f=1.0661e+01 df=7.2201e+01 it=( 1, 1) 9 8 t=9.5542e-01 dt=3.6083e+00 f=2.1794e+00 df=1.3544e+01 it=( 1, 1) 10 9 t=4.5637e+00 dt=1.8537e+01 f=1.5598e+00 df=2.9000e+00 it=( 1, 1) [all …]
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | wedge_cylinder.cas | 28 4 2 3 2f 2e 3 a 29 4 3 4 30 2f 4 b 38 4 3 2f 39 9 a b 40 4 35 7 37 36 8 f 51 4 35 7 d 3f f 3a 53 4 7 37 41 d f 10 57 4 8 9 f e 11 3b 58 4 9 39 43 f 11 12 59 4 9 a 10 f 12 3c 60 4 3f d 41 40 f 16 [all …]
|
| H A D | square.cas | 115 3e d 1f 0 124 11 f 6 0 125 f 20 2e 0 139 2e 2f 14 0 140 2f 51 4a 0 148 4 3f 1e 1d 149 3e 3f 1d 22 150 3e c 22 1f 151 d c 1f 5 166 11 10 2f 6 [all …]
|
| H A D | testcase3D.cas | 11 (sg-swirl? #f) 13 (rf-energy? #f) 14 (rp-inviscid? #f) 15 (rp-lam? #f) 16 (rp-atm? #f) 17 (rp-hvac? #f) 18 (rp-les? #f) 19 (rp-les-sgske? #f) 20 (rp-sa? #f) 21 (rp-des? #f) [all …]
|
| /petsc/include/petsc/private/ |
| H A D | petsclegacycupmblas.h | 9 …lasXaxpy(a, b, c, d, e, f, g) cublasCaxpy((a), (b), (cuComplex *)(c), (cuComp… argument 11 …ne cublasXdotu(a, b, c, d, e, f, g) cublasCdotu((a), (b), (cuComplex *)(c), (… argument 12 …ine cublasXdot(a, b, c, d, e, f, g) cublasCdotc((a), (b), (cuComplex *)(c), … argument 13 …ne cublasXswap(a, b, c, d, e, f) cublasCswap((a), (b), (cuComplex *)(c), (… argument 17 …ne cublasXgemv(a, b, c, d, e, f, g, h, i, j, k, l) cublasCgemv((a), (b), (c), (d), (cuComple… argument 18 …#define cublasXgemm(a, b, c, d, e, f, g, h, i, j, k, l, m, n) cublasCgemm((a), (b), (c), (d), (e),… argument 19 …#define cublasXgeam(a, b, c, d, e, f, g, h, i, j, k, l, m) cublasCgeam((a), (b), (c), (d), (e),… argument 20 #define cublasXgemvStridedBatched(a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p) \ argument 21 …cublasCgemvStridedBatched((a), (b), (c), (d), (const cuComplex *)(e), (const cuComplex *)(f), (g),… 23 …a, b, c, d, e, f, g) cublasZaxpy((a), (b), (cuDoubleComplex *)(c), (cuDoubleC… argument [all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dtds.c | 113 PetscInt Nf, numConstants, f; in PetscDSView_Ascii() local 122 for (f = 0; f < Nf; ++f) { in PetscDSView_Ascii() 130 PetscCall(PetscDSGetDiscretization(ds, f, &obj)); in PetscDSView_Ascii() 143 …((PetscObject)ds), PETSC_ERR_ARG_WRONG, "Unknown discretization type for field %" PetscInt_FMT, f); in PetscDSView_Ascii() 146 if (ds->implicit[f]) PetscCall(PetscViewerASCIIPrintf(viewer, " (implicit)")); in PetscDSView_Ascii() 152 PetscCall(PetscViewerASCIIPrintf(viewer, " %" PetscInt_FMT "-jet", ds->jetDegree[f])); in PetscDSView_Ascii() 164 if (b->field != f) continue; in PetscDSView_Ascii() 213 …for (f = 0; f < numConstants; ++f) PetscCall(PetscViewerASCIIPrintf(viewer, "%g\n", (double)PetscR… in PetscDSView_Ascii() 366 PetscInt dim, dimEmbed, NbMax = 0, NcMax = 0, NqMax = 0, NsMax = 1, f; in PetscDSSetUp() local 384 for (f = 0; f < Nf; ++f) { in PetscDSSetUp() [all …]
|
| /petsc/src/mat/tests/ |
| H A D | mmio.c | 18 FILE *f; in mm_read_unsymmetric_sparse() local 25 if ((f = fopen(fname, "r")) == NULL) return -1; in mm_read_unsymmetric_sparse() 27 if (mm_read_banner(f, &matcode) != 0) { in mm_read_unsymmetric_sparse() 41 if (mm_read_mtx_crd_size(f, &M, &N, &nz) != 0) { in mm_read_unsymmetric_sparse() 65 if (fscanf(f, "%d %d %lg\n", &ia[i], &ja[i], &val[i]) != 3) { in mm_read_unsymmetric_sparse() 72 fclose(f); in mm_read_unsymmetric_sparse() 86 int mm_read_banner(FILE *f, MM_typecode *matcode) in mm_read_banner() argument 98 if (fgets(line, MM_MAX_LINE_LENGTH, f) == NULL) return MM_PREMATURE_EOF; in mm_read_banner() 140 int mm_write_mtx_crd_size(FILE *f, int M, int N, int nz) in mm_write_mtx_crd_size() argument 142 if (fprintf(f, "%d %d %d\n", M, N, nz) < 0) return MM_COULD_NOT_WRITE_FILE; in mm_write_mtx_crd_size() [all …]
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex12.c | 41 PetscPointFn *f[4] = {f0, f1, f2, f3}; in TestSetIndex() local 61 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, i, f[i], 0, NULL)); in TestSetIndex() 62 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, j, f[j], 0, NULL)); in TestSetIndex() 63 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, k, f[k], 0, NULL)); in TestSetIndex() 64 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, l, f[l], 0, NULL)); in TestSetIndex() 65 PetscCall(CheckResidual(wf, key, 4, f, 0, NULL)); in TestSetIndex() 79 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, 0, NULL, i, f[i])); in TestSetIndex() 80 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, 0, NULL, j, f[j])); in TestSetIndex() 81 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, 0, NULL, k, f[k])); in TestSetIndex() 82 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, 0, NULL, l, f[l])); in TestSetIndex() [all …]
|
| H A D | ex3.c | 116 PetscInt f; in main() local 125 for (f = 0; f < 14; ++f) { in main() 129 if ((f == 6) || (f == 7) || (f == 9) || (f == 11)) continue; in main() 131 if (f == 8) continue; in main() 133 …PetscCall(PetscDTTanhSinhIntegrate(funcs[f], bounds[f * 2 + 0], bounds[f * 2 + 1], digits, NULL, &… in main() 134 …Real(integral - analytic[f]) > PetscMax(epsilon, PetscPowRealInt(10.0, -digits)) || PetscIsInfOrNa… in main() 135 …nc%2" PetscInt_FMT " is wrong: %g (%g)\n", f + 1, (double)integral, (double)PetscAbsReal(integral … in main() 139 for (f = 0; f < 14; ++f) { in main() 142 …PetscCall(PetscDTTanhSinhIntegrateMPFR(funcs[f], bounds[f * 2 + 0], bounds[f * 2 + 1], digits, NUL… in main() 143 …f]) > PetscPowRealInt(10.0, -digits)) PetscCall(PetscPrintf(PETSC_COMM_SELF, "The integral of func… in main()
|
| /petsc/src/dm/impls/plex/tests/output/ |
| H A D | ex36_0.out | 13 rank 0 f 25 cells 0,1 normal 1.,0.,0. centroid 0.125,0.5,0.5 14 rank 0 f 26 cells 1,2 normal 1.,0.,0. centroid 0.25,0.5,0.5 15 rank 0 f 39 cells 2,3 normal 1.,0.,0. centroid 0.375,0.5,0.5 16 rank 1 f 29 cells 0,3 normal -1.,-0.,-0. centroid 0.375,0.5,0.5 17 rank 1 f 30 cells 0,1 normal 1.,0.,0. centroid 0.5,0.5,0.5 18 rank 1 f 31 cells 1,2 normal 1.,0.,0. centroid 0.625,0.5,0.5 19 rank 1 f 49 cells 2,4 normal 1.,0.,0. centroid 0.75,0.5,0.5 20 rank 2 f 19 cells 0,2 normal -1.,-0.,-0. centroid 0.75,0.5,0.5 21 rank 2 f 20 cells 0,1 normal 1.,0.,0. centroid 0.875,0.5,0.5 35 rank 0 f 25 cells 0,1 normal 1.,0.,0. centroid 0.125,0.5,0.5 [all …]
|
| /petsc/src/mat/impls/aij/seq/seqcusparse/ |
| H A D | cusparsematimpl.h | 23 …parseXcsrilu02_bufferSize(a, b, c, d, e, f, g, h, i) cusparseCcsrilu02_bufferSize(a, b, c, d, (cu… argument 24 …usparseXcsrilu02_analysis(a, b, c, d, e, f, g, h, i, j) cusparseCcsrilu02_analysis(a, b, c, d, (cu… argument 25 …#define cusparseXcsrilu02(a, b, c, d, e, f, g, h, i, j) cusparseCcsrilu02(a, b, c, d, (cu… argument 26 …sparseXcsric02_bufferSize(a, b, c, d, e, f, g, h, i) cusparseCcsric02_bufferSize(a, b, c, d, (cu… argument 27 …cusparseXcsric02_analysis(a, b, c, d, e, f, g, h, i, j) cusparseCcsric02_analysis(a, b, c, d, (cu… argument 28 …#define cusparseXcsric02(a, b, c, d, e, f, g, h, i, j) cusparseCcsric02(a, b, c, d, (cuC… argument 32 …seXcsrilu02_bufferSize(a, b, c, d, e, f, g, h, i) cusparseZcsrilu02_bufferSize(a, b, c, d, (cuDou… argument 33 …arseXcsrilu02_analysis(a, b, c, d, e, f, g, h, i, j) cusparseZcsrilu02_analysis(a, b, c, d, (cuDou… argument 34 …fine cusparseXcsrilu02(a, b, c, d, e, f, g, h, i, j) cusparseZcsrilu02(a, b, c, d, (cuDou… argument 35 …rseXcsric02_bufferSize(a, b, c, d, e, f, g, h, i) cusparseZcsric02_bufferSize(a, b, c, d, (cuDou… argument [all …]
|
| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | ex6.c | 39 PetscInt i, j, f, c; in LoadData2D() local 49 for (f = 0; f < user->Nf; ++f) { in LoadData2D() 50 PetscInt ioff = i * user->k[f], joff = j * user->k[f]; in LoadData2D() 52 for (kj = 0; kj <= user->k[f]; ++kj) { in LoadData2D() 53 for (ki = 0; ki <= user->k[f]; ++ki) { in LoadData2D() 54 …for (c = 0; c < user->Nc[f]; ++c) closure[o++] = ((kj + joff) * (Ni * user->k[f] + 1) + ki + ioff)… in LoadData2D() 67 PetscInt i, j, k, f, c; in LoadData3D() local 78 for (f = 0; f < user->Nf; ++f) { in LoadData3D() 79 PetscInt ioff = i * user->k[f], joff = j * user->k[f], koff = k * user->k[f]; in LoadData3D() 81 for (kk = 0; kk <= user->k[f]; ++kk) { in LoadData3D() [all …]
|
| /petsc/src/dm/impls/swarm/ |
| H A D | data_bucket.c | 91 PetscInt f; in DMSwarmDataBucketDestroy() local 95 for (f = 0; f < db->nfields; ++f) PetscCall(DMSwarmDataFieldDestroy(&db->field[f])); in DMSwarmDataBucketDestroy() 105 PetscInt f; in DMSwarmDataBucketQueryForActiveFields() local 109 for (f = 0; f < db->nfields; ++f) { in DMSwarmDataBucketQueryForActiveFields() 110 if (db->field[f]->active) { in DMSwarmDataBucketQueryForActiveFields() 239 …PetscInt current_allocated, current_used, new_used, new_unused, new_buffer, new_allocated, f, end; in DMSwarmDataBucketSetSizes() local 259 for (f = 0; f < db->nfields; f++) { in DMSwarmDataBucketSetSizes() 260 PetscCall(DMSwarmDataFieldSetSize(db->field[f], new_allocated)); in DMSwarmDataBucketSetSizes() 261 PetscCall(DMSwarmDataFieldZeroBlock(db->field[f], current_allocated, new_allocated)); in DMSwarmDataBucketSetSizes() 269 … for (f = 0; f < db->nfields; ++f) PetscCall(DMSwarmDataFieldSetSize(db->field[f], new_allocated)); in DMSwarmDataBucketSetSizes() [all …]
|
| /petsc/src/mat/impls/aij/seq/seqhipsparse/ |
| H A D | hipsparsematimpl.h | 33 …arseXcsrilu02_bufferSize(a, b, c, d, e, f, g, h, i) hipsparseCcsrilu02_bufferSize(a, b, c, d, (hi… argument 34 …sparseXcsrilu02_analysis(a, b, c, d, e, f, g, h, i, j) hipsparseCcsrilu02_analysis(a, b, c, d, (hi… argument 35 …efine hipsparseXcsrilu02(a, b, c, d, e, f, g, h, i, j) hipsparseCcsrilu02(a, b, c, d, (hi… argument 36 …parseXcsric02_bufferSize(a, b, c, d, e, f, g, h, i) hipsparseCcsric02_bufferSize(a, b, c, d, (hi… argument 37 …psparseXcsric02_analysis(a, b, c, d, e, f, g, h, i, j) hipsparseCcsric02_analysis(a, b, c, d, (hi… argument 38 …define hipsparseXcsric02(a, b, c, d, e, f, g, h, i, j) hipsparseCcsric02(a, b, c, d, (hi… argument 42 …eXcsrilu02_bufferSize(a, b, c, d, e, f, g, h, i) hipsparseZcsrilu02_bufferSize(a, b, c, d, (hipDo… argument 43 …rseXcsrilu02_analysis(a, b, c, d, e, f, g, h, i, j) hipsparseZcsrilu02_analysis(a, b, c, d, (hipDo… argument 44 …ne hipsparseXcsrilu02(a, b, c, d, e, f, g, h, i, j) hipsparseZcsrilu02(a, b, c, d, (hipDo… argument 45 …seXcsric02_bufferSize(a, b, c, d, e, f, g, h, i) hipsparseZcsric02_bufferSize(a, b, c, d, (hipDo… argument [all …]
|
| /petsc/src/snes/utils/ |
| H A D | convest.c | 148 PetscInt Nf, f, Nds, s; in PetscConvEstSetUp() local 155 for (f = 0; f < Nf; ++f) ce->initGuess[f] = zero_private; in PetscConvEstSetUp() 167 for (f = 0; f < dsNf; ++f) { in PetscConvEstSetUp() 168 const PetscInt field = fields[f]; in PetscConvEstSetUp() 173 …f = 0; f < Nf; ++f) PetscCheck(ce->exactSol[f], PetscObjectComm((PetscObject)ce), PETSC_ERR_ARG_WR… in PetscConvEstSetUp() 217 PetscInt f; in PetscConvEstMonitorDefault() local 227 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstMonitorDefault() 228 if (f > 0) PetscCall(PetscPrintf(comm, ", ")); in PetscConvEstMonitorDefault() 229 PetscCall(PetscPrintf(comm, "%7" PetscInt_FMT, dofs[f])); in PetscConvEstMonitorDefault() 235 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstMonitorDefault() [all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexsection.c | 7 PetscInt depth, Nf, f, pStart, pEnd; in DMPlexCreateSectionFields() local 15 for (f = 0; f < Nf; ++f) { in DMPlexCreateSectionFields() 19 PetscCall(DMGetField(dm, f, NULL, &obj)); in DMPlexCreateSectionFields() 22 isFE[f] = PETSC_TRUE; in DMPlexCreateSectionFields() 24 isFE[f] = PETSC_FALSE; in DMPlexCreateSectionFields() 32 for (f = 0; f < Nf; ++f) { in DMPlexCreateSectionFields() 33 PetscCall(PetscSectionSetFieldComponents(*section, f, numComp[f])); in DMPlexCreateSectionFields() 34 if (isFE[f]) { in DMPlexCreateSectionFields() 41 PetscCall(DMGetField(dm, f, NULL, (PetscObject *)&fe)); in DMPlexCreateSectionFields() 75 PetscCall(PetscSectionSetFieldSym(*section, f, sym)); in DMPlexCreateSectionFields() [all …]
|
| H A D | plexproject.c | 73 PetscInt coordDim, Nf, *Nc, f, spDim, d, v, tp; in DMProjectPoint_Func_Private() local 84 for (f = 0, v = 0, tp = 0; f < Nf; ++f) { in DMProjectPoint_Func_Private() 85 void *const ctx = ctxs ? ctxs[f] : NULL; in DMProjectPoint_Func_Private() 88 if (!sp[f]) continue; in DMProjectPoint_Func_Private() 89 PetscCall(PetscDSGetCohesive(ds, f, &cohesive)); in DMProjectPoint_Func_Private() 90 PetscCall(PetscDualSpaceGetDimension(sp[f], &spDim)); in DMProjectPoint_Func_Private() 91 if (funcs[f]) { in DMProjectPoint_Func_Private() 92 if (isFE[f]) { in DMProjectPoint_Func_Private() 100 PetscCall(PetscDualSpaceGetDM(sp[f], &rdm)); in DMProjectPoint_Func_Private() 101 PetscCall(PetscDualSpaceGetAllData(sp[f], &allPoints, NULL)); in DMProjectPoint_Func_Private() [all …]
|
| /petsc/src/sys/tests/linter/output/ |
| H A D | testValidPointers.patch | 4 …ters(void *a, char *b, PetscInt *c, PetscMPIInt *d, PetscInt *e, PetscBool *f, PetscScalar *g, Pet… 5 …ters(void *a, char *b, PetscInt *c, PetscMPIInt *d, PetscInt *e, PetscBool *f, PetscScalar *g, Pet… 35 PetscAssertPointer(f, 7); 41 PetscAssertPointer(f, 7); 46 - PetscAssertPointer(f, 7); 47 + PetscAssertPointer(f, 6); 52 PetscAssertPointer(f, 7); 58 PetscAssertPointer(f, 7); 67 -void testValidPointers2(void *a, char *b, PetscInt *c, PetscMPIInt *d, PetscInt *e, PetscBool *f, … 68 …ers2(void *a, char *b, PetscInt *c, PetscMPIInt *d, PetscInt *e, PetscBool *f, PetscScalar *g, Pet… [all …]
|
| H A D | testTypes.patch | 4 -PetscErrorCode testTypes(PetscRandom r, PetscViewer v, PetscObject o, PetscFunctionList f) 5 +static PetscErrorCode testTypes(PetscRandom r, PetscViewer v, PetscObject o, PetscFunctionList f) 14 PetscCheckSameComm(o, -2, f, -2); 20 PetscCheckSameComm(o, -2, f, -2); 21 PetscCheckSameTypeAndComm(r, -3, f, -3); 25 - PetscCheckSameComm(o, -2, f, -2); 26 + PetscCheckSameComm(o, 3, f, 4); 27 PetscCheckSameTypeAndComm(r, -3, f, -3); 31 PetscCheckSameComm(o, -2, f, -2); 32 - PetscCheckSameTypeAndComm(r, -3, f, -3); [all …]
|
| /petsc/lib/petsc/bin/ |
| H A D | petsc_tas_analysis.py | 39 …print(f'\nconfigureTAS.py defaultGraphs was updated with path {cmdLineArgs.setDefaultGraphDir[0]}\… 52 … print(f'\nconfigureTAS.py defaultData was updated with path {cmdLineArgs.setDefaultFileDir[0]}\n') 68 …print(f'\nconfigureTAS.py was updated with\n\talias: {listToAdd[counter]}\n\tpath {listToAdd[count… 70 …print(f'\nconfigureTAS.py was NOT updated with \n\talias: {listToAdd[counter]}\n\tpath {listToAdd[… 72 …print(f'\nWhen using the command line option to add or edit an alias path pair, both must be inclu… 73 f'Your input was:') 83 print(f'\n{alias} was successfully removed from configureTAS.py\n') 85 print(f'\n{alias} was not successful removed from configureTAS.py\n') 87 print(f'\n{alias} was not found in configureTAS.py\nList of valid aliases is:\n') 132 print(f'\nalias: {alias}\nalready has path: {config.filePath[alias]}\n\n' [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex31.c | 73 PetscScalar *f; in RHSFunction_Hull1972A1() local 78 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972A1() 79 f[0] = -y[0]; in RHSFunction_Hull1972A1() 81 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972A1() 103 PetscScalar *f; in IFunction_Hull1972A1() local 107 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972A1() 108 f[0] = -y[0]; in IFunction_Hull1972A1() 110 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972A1() 136 PetscScalar *f; in RHSFunction_Hull1972A2() local 140 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972A2() [all …]
|
| /petsc/src/dm/tutorials/ |
| H A D | ex13f90aux.F90 | 21 function dfdt_vdp(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n, f) 28 PetscReal, dimension(n, ib1:ibn, jb1:jbn, kb1:kbn), intent(inout) :: f 32 dfdt_vdp(1, :, :, :) = f(2, 1, 1, 1) 33 dfdt_vdp(2, :, :, :) = mu*(one - f(1, 1, 1, 1)**2)*f(2, 1, 1, 1) - f(1, 1, 1, 1) 46 function dfdt(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n, f) 50 PetscReal, dimension(n, ib1:ibn, jb1:jbn, kb1:kbn), intent(inout) :: f 68 subroutine petsc_to_local(da, vec, array, f, dof, stw) argument 74 PetscReal, intent(inout), dimension(:, 1 - stw:, 1 - stw:, 1 - stw:) :: f 78 call transform_petsc_us(array, f, stw) 80 subroutine transform_petsc_us(array, f, stw) argument [all …]
|
| /petsc/src/dm/interface/ |
| H A D | dmi.c | 91 for (PetscInt f = 0, off = 0; f < numFields; ++f) { in PetscSectionSelectFields_Private() local 94 if (numComps[f] < 0) continue; in PetscSectionSelectFields_Private() 95 PetscCall(PetscSectionGetFieldComponents(s, f, &Nc)); in PetscSectionSelectFields_Private() 96 …f] <= Nc, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Field %" PetscInt_FMT ": Number of selected … in PetscSectionSelectFields_Private() 97 …f]; ++c, ++off) PetscCheck(comps[off] < Nc, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Field %" P… in PetscSectionSelectFields_Private() 98 bs += numComps[f]; in PetscSectionSelectFields_Private() 101 for (PetscInt f = 0; f < numFields; ++f) { in PetscSectionSelectFields_Private() local 104 PetscCall(PetscSectionGetFieldComponents(s, fields[f], &Nc)); in PetscSectionSelectFields_Private() 117 for (PetscInt f = 0; f < numFields; ++f) { in PetscSectionSelectFields_Private() local 120 PetscCall(PetscSectionGetFieldComponents(s, f, &Nc)); in PetscSectionSelectFields_Private() [all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex8.c | 16 PetscErrorCode f(PetscReal t, Vec UV, Vec F) in f() function 19 PetscScalar *f; in f() local 27 PetscCall(VecGetArrayWrite(F, &f)); in f() 28 for (i = 0; i < n; i++) f[i] = u[i] + v[i]; in f() 30 PetscCall(VecRestoreArrayWrite(F, &f)); in f() 41 PetscScalar *f; in F() local 49 PetscCall(VecGetArrayWrite(F, &f)); in F() 50 f = f + n; in F() 51 for (i = 0; i < n; i++) f[i] = u[i] - v[i]; in F() 52 f = f - n; in F() [all …]
|