Home
last modified time | relevance | path

Searched refs:fs (Results 1 – 22 of 22) sorted by relevance

/petsc/src/mat/impls/aij/seq/seqcusparse/
H A Daijcusparse.cu204 Mat_SeqAIJCUSPARSETriFactors *fs = static_cast<Mat_SeqAIJCUSPARSETriFactors *>(A->spptr); in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU() local
213 …if (!fs->csrRowPtr) { // Is't the first time to do the setup? Use csrRowPtr sin… in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
229 PetscCallCUDA(cudaMalloc(&fs->csrRowPtr, sizeof(*fs->csrRowPtr) * (m + 1))); in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
230 PetscCallCUDA(cudaMalloc(&fs->csrColIdx, sizeof(*fs->csrColIdx) * Mnz)); in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
231 PetscCallCUDA(cudaMalloc(&fs->csrVal, sizeof(*fs->csrVal) * Mnz)); in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
232 …PetscCallCUDA(cudaMemcpy(fs->csrRowPtr, Mi, sizeof(*fs->csrRowPtr) * (m + 1), cudaMemcpyHostToDevi… in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
233 …PetscCallCUDA(cudaMemcpy(fs->csrColIdx, Mj, sizeof(*fs->csrColIdx) * Mnz, cudaMemcpyHostToDevice)); in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
244 …PetscCallCUSPARSE(cusparseCreateCsr(&fs->spMatDescr_L, m, m, Mnz, fs->csrRowPtr, fs->csrColIdx, fs… in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
245 …PetscCallCUSPARSE(cusparseSpMatSetAttribute(fs->spMatDescr_L, CUSPARSE_SPMAT_FILL_MODE, &fillMode,… in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
246 …PetscCallCUSPARSE(cusparseSpMatSetAttribute(fs->spMatDescr_L, CUSPARSE_SPMAT_DIAG_TYPE, &diagType,… in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
[all …]
/petsc/src/mat/impls/aij/seq/seqhipsparse/
H A Daijhipsparse.hip.cxx1249 Mat_SeqAIJHIPSPARSETriFactors *fs = (Mat_SeqAIJHIPSPARSETriFactors *)fact->spptr; in MatSolve_SeqAIJHIPSPARSE_ILU0() local
1260 PetscCallHIPSPARSE(hipsparseDnVecSetValues(fs->dnVecDescr_X, (void *)barray)); in MatSolve_SeqAIJHIPSPARSE_ILU0()
1261 PetscCallHIPSPARSE(hipsparseDnVecSetValues(fs->dnVecDescr_Y, fs->Y)); in MatSolve_SeqAIJHIPSPARSE_ILU0()
1263 …PetscCallHIPSPARSE(hipsparseSpSV_solve(fs->handle, HIPSPARSE_OPERATION_NON_TRANSPOSE, &PETSC_HIPSP… in MatSolve_SeqAIJHIPSPARSE_ILU0()
1264 …fs->dnVecDescr_X, fs->dnVecDescr_Y, hipsparse_scalartype, HIPSPARSE_SPSV_ALG_DEFAULT, fs->spsvDesc… in MatSolve_SeqAIJHIPSPARSE_ILU0()
1266 …PetscCallHIPSPARSE(hipsparseSpSV_solve(fs->handle, HIPSPARSE_OPERATION_NON_TRANSPOSE, &PETSC_HIPSP… in MatSolve_SeqAIJHIPSPARSE_ILU0()
1267 …fs->dnVecDescr_X, fs->dnVecDescr_Y, hipsparse_scalartype, HIPSPARSE_SPSV_ALG_DEFAULT, fs->spsvDesc… in MatSolve_SeqAIJHIPSPARSE_ILU0()
1270 PetscCallHIPSPARSE(hipsparseDnVecSetValues(fs->dnVecDescr_X, xarray)); in MatSolve_SeqAIJHIPSPARSE_ILU0()
1272 …PetscCallHIPSPARSE(hipsparseSpSV_solve(fs->handle, HIPSPARSE_OPERATION_NON_TRANSPOSE, &PETSC_HIPSP… in MatSolve_SeqAIJHIPSPARSE_ILU0()
1273 …fs->dnVecDescr_Y, fs->dnVecDescr_X, hipsparse_scalartype, HIPSPARSE_SPSV_ALG_DEFAULT, fs->spsvDesc… in MatSolve_SeqAIJHIPSPARSE_ILU0()
[all …]
/petsc/src/ts/tutorials/multirate/
H A Dex6.c144 hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf); in FVSample_2WaySplit()
155 } else if (i < ctx->fs) { in FVSample_2WaySplit()
165 xi = ctx->xmin + ctx->sf * hs + (ctx->fs - ctx->sf) * hf + 0.5 * hs + (i - ctx->fs) * hs; in FVSample_2WaySplit()
192 hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf); in SolutionErrorNorms_2WaySplit()
196 if (i < ctx->sf || i > ctx->fs - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_2WaySplit()
208 PetscInt i, j, k, Mx, dof, xs, xm, sf = ctx->sf, fs = ctx->fs; in FVRHSFunction_2WaySplit() local
219 hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf); in FVRHSFunction_2WaySplit()
261 (*ctx->limit2)(&info, cjmpL, cjmpR, ctx->sf, ctx->fs, i, ctx->cslope); in FVRHSFunction_2WaySplit()
298 } else if (i > sf && i < fs) { /* fast region */ in FVRHSFunction_2WaySplit()
310 } else if (i == fs) { /* interface between the fast region and the slow region */ in FVRHSFunction_2WaySplit()
[all …]
H A Dex4.c488 hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf); in FVSample_2WaySplit()
499 } else if (i < ctx->fs) { in FVSample_2WaySplit()
509 xi = ctx->xmin + ctx->sf * hs + (ctx->fs - ctx->sf) * hf + 0.5 * hs + (i - ctx->fs) * hs; in FVSample_2WaySplit()
536 hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf); in SolutionErrorNorms_2WaySplit()
540 if (i < ctx->sf || i > ctx->fs - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_2WaySplit()
552 PetscInt i, j, k, Mx, dof, xs, xm, sf = ctx->sf, fs = ctx->fs; in FVRHSFunction_2WaySplit() local
563 hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf); in FVRHSFunction_2WaySplit()
630 (*ctx->limit2)(&info, cjmpL, cjmpR, ctx->sf, ctx->fs, i, ctx->cslope); in FVRHSFunction_2WaySplit()
667 } else if (i > sf && i < fs) { /* fast region */ in FVRHSFunction_2WaySplit()
679 } else if (i == fs) { /* interface between the fast region and the slow region */ in FVRHSFunction_2WaySplit()
[all …]
H A Dfinitevolume1d.c143 … PetscScalar *jL, const PetscScalar *jR, const PetscInt sf, const PetscInt fs, PetscInt n, PetscSc… in Limit2_Upwind() argument
148 … PetscScalar *jL, const PetscScalar *jR, const PetscInt sf, const PetscInt fs, PetscInt n, PetscSc… in Limit2_LaxWendroff() argument
157 } else if (n > sf && n < fs - 1) { /* fast components */ in Limit2_LaxWendroff()
159 } else if (n == fs - 1) { /* fast component next to slow components */ in Limit2_LaxWendroff()
161 } else if (n == fs) { /* slow component next to fast components */ in Limit2_LaxWendroff()
167 … PetscScalar *jL, const PetscScalar *jR, const PetscInt sf, const PetscInt fs, PetscInt n, PetscSc… in Limit2_BeamWarming() argument
176 } else if (n > sf && n < fs - 1) { in Limit2_BeamWarming()
178 } else if (n == fs - 1) { in Limit2_BeamWarming()
180 } else if (n == fs) { in Limit2_BeamWarming()
186 … PetscScalar *jL, const PetscScalar *jR, const PetscInt sf, const PetscInt fs, PetscInt n, PetscSc… in Limit2_Fromm() argument
[all …]
H A Dex7.c71 PetscInt sf, fs; /* slow-fast and fast-slow interfaces */ member
172 PetscInt i, j, Mx, dof, xs, xm, sf = ctx->sf, fs = ctx->fs; in FVRHSFunction() local
204 if (i < sf || i > fs + 1) { in FVRHSFunction()
253 } else if (i > sf + 1 && i < fs) { in FVRHSFunction()
269 } else if (i == fs) { in FVRHSFunction()
285 } else if (i == fs + 1) { in FVRHSFunction()
321 PetscInt i, j, Mx, dof, xs, xm, islow = 0, sf = ctx->sf, fs = ctx->fs; in FVRHSFunctionslow() local
383 } else if (i == fs) { in FVRHSFunctionslow()
397 } else if (i == fs + 1) { in FVRHSFunctionslow()
414 } else if (i > fs + 1) { in FVRHSFunctionslow()
[all …]
H A Dfinitevolume1d.h130 PetscInt sf, fs; /* slow-fast and fast-slow interfaces */ member
/petsc/src/dm/impls/plex/exodusii/
H A Dplexexodusii2.c1159 PetscInt i, j, fs, fsSize; in DMView_PlexExodusII() local
1171 for (fs = 0; fs < num_fs; ++fs) { in DMView_PlexExodusII()
1172 PetscCall(DMLabelGetStratumIS(fsLabel, fsIdx[fs], &stratumIS)); in DMView_PlexExodusII()
1180 for (fs = 0; fs < num_fs; ++fs) { in DMView_PlexExodusII()
1181 PetscCall(DMLabelGetStratumIS(fsLabel, fsIdx[fs], &stratumIS)); in DMView_PlexExodusII()
1185 PetscCallExternal(ex_put_set_param, exo->exoid, EX_SIDE_SET, fsIdx[fs], fsSize, 0); in DMView_PlexExodusII()
1187 if (fs < num_fs - 1) elem_ind[fs + 1] = elem_ind[fs] + fsSize; in DMView_PlexExodusII()
1193 elem_list[elem_ind[fs] + i] = points[2] + 1; in DMView_PlexExodusII()
1198 …PetscCall(DMPlexGetTransitiveClosure(dm, elem_list[elem_ind[fs] + i] - 1, PETSC_TRUE, &numPoints, … in DMView_PlexExodusII()
1223 side_list[elem_ind[fs] + i] = j; in DMView_PlexExodusII()
[all …]
/petsc/src/ksp/pc/impls/gamg/
H A Dclassical.c190 PetscInt fn, cn, fs, fe, cs, ce, i, j, ncols, col, row_f, row_c, cmax = 0, idx, noff; in PCGAMGProlongator_Classical_Direct() local
205 PetscCall(MatGetOwnershipRange(A, &fs, &fe)); in PCGAMGProlongator_Classical_Direct()
206 fn = fe - fs; in PCGAMGProlongator_Classical_Direct()
256 for (i = fs; i < fe; i++) { in PCGAMGProlongator_Classical_Direct()
257 Amax_pos[i - fs] = 0.; in PCGAMGProlongator_Classical_Direct()
258 Amax_neg[i - fs] = 0.; in PCGAMGProlongator_Classical_Direct()
261 …if ((PetscRealPart(-rval[j]) > Amax_neg[i - fs]) && i != rcol[j]) Amax_neg[i - fs] = PetscAbsScala… in PCGAMGProlongator_Classical_Direct()
262 …if ((PetscRealPart(rval[j]) > Amax_pos[i - fs]) && i != rcol[j]) Amax_pos[i - fs] = PetscAbsScalar… in PCGAMGProlongator_Classical_Direct()
308 row_f = i + fs; in PCGAMGProlongator_Classical_Direct()
385 row_f = i + fs; in PCGAMGProlongator_Classical_Direct()
[all …]
/petsc/src/snes/tutorials/
H A Dbuild.zig20 … var path = std.fs.path.join(std.heap.c_allocator, &[_][] const u8 { PETSC_DIR,PETSC_ARCH,"lib"});
24 path = std.fs.path.join(std.heap.c_allocator, &[_][] const u8 { PETSC_DIR,"include"});
26 … path = std.fs.path.join(std.heap.c_allocator, &[_][] const u8 { PETSC_DIR,PETSC_ARCH,"include"});
H A Dex18.c138 …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
184 fs = ds * (t0 - ts); in FormFunction()
207 fs = ds * (t0 - ts); in FormFunction()
208 } else fs = zero; in FormFunction()
233 fs = ds * (t0 - ts); in FormFunction()
234 } else fs = zero; in FormFunction()
255 fs = zero; in FormFunction()
277 fs = ds * (t0 - ts); in FormFunction()
282 f[j][i] = -hydhx * (fe - fw) - hxdhy * (fn - fs); in FormFunction()
/petsc/include/petsc/private/
H A Ddmimpl.h520 const PetscSection fs = dm->localSection->field[field]; in DMGetGlobalFieldOffset_Private() local
524 const PetscInt lfoff = fs->atlasOff[point - s->pStart]; in DMGetGlobalFieldOffset_Private()
525 const PetscInt fdof = fs->atlasDof[point - s->pStart]; in DMGetGlobalFieldOffset_Private()
526 const PetscInt fcdof = fs->bc ? fs->bc->atlasDof[point - fs->bc->pStart] : 0; in DMGetGlobalFieldOffset_Private()
/petsc/systems/Apple/OSX/bin/
H A Dmakedmg29 hdiutil create -srcfolder tempsource -volname PETSc -fs HFS+ -fsargs "-c c=64,a=16,e=16" -format UD…
/petsc/src/ts/utils/
H A Dtsconvest.c206 PetscSection s, fs; in PetscConvEstGetConvRateTS_Spatial_Private() local
211 PetscCall(PetscSectionGetField(s, f, &fs)); in PetscConvEstGetConvRateTS_Spatial_Private()
212 PetscCall(PetscSectionGetConstrainedStorageSize(fs, &lsize)); in PetscConvEstGetConvRateTS_Spatial_Private()
/petsc/systems/Apple/iOS/bin/
H A Dmakedmg34 hdiutil create -srcfolder tempsource -volname PETSc -fs HFS+ -fsargs "-c c=64,a=16,e=16" -format UD…
/petsc/src/snes/tests/
H A Dex20.c138 …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
187 fs = ds * (t0 - ts); in FormFunction()
220 fs = ds * (t0 - ts); in FormFunction()
222 fs = zero; in FormFunction()
268 fs = ds * (t0 - ts); in FormFunction()
270 fs = zero; in FormFunction()
312 fs = zero; in FormFunction()
352 fs = ds * (t0 - ts); in FormFunction()
389 fs = ds * (t0 - ts); in FormFunction()
418 fs = ds * (t0 - ts); in FormFunction()
[all …]
/petsc/src/snes/utils/
H A Dconvest.c405 PetscSection s, fs; in PetscConvEstGetConvRateSNES_Private() local
410 PetscCall(PetscSectionGetField(s, f, &fs)); in PetscConvEstGetConvRateSNES_Private()
411 PetscCall(PetscSectionGetConstrainedStorageSize(fs, &lsize)); in PetscConvEstGetConvRateSNES_Private()
/petsc/src/dm/impls/plex/cgns/
H A Dplexcgns2.c2239 …for (PetscInt fs = 0; fs < num_fs_local; ++fs) PetscCall(DMLabelGetStratumIS(fsLabel, fsValuesLoca… in DMView_PlexCGNS() local
2248 for (PetscInt fs = 0; fs < num_fs_local; ++fs) PetscCall(ISDestroy(&fsPoints[fs])); in DMView_PlexCGNS() local
2307 for (PetscInt fs = 0; fs < num_fs_global; ++fs) { in DMView_PlexCGNS() local
2309 const PetscInt fsID = fsValuesGlobal[fs]; in DMView_PlexCGNS()
/petsc/src/ts/tutorials/
H A Dex18.c897 PetscSection fs; in MonitorFunctionals() local
907 PetscCall(PetscSectionClone(s, &fs)); in MonitorFunctionals()
908 PetscCall(PetscSectionSetFieldName(fs, 0, NULL)); in MonitorFunctionals()
909 PetscCall(PetscSectionSetFieldName(fs, 1, name)); in MonitorFunctionals()
910 PetscCall(DMSetLocalSection(fdm[f], fs)); in MonitorFunctionals()
911 PetscCall(PetscSectionDestroy(&fs)); in MonitorFunctionals()
/petsc/src/dm/impls/plex/
H A Dplexsubmesh.c920 PetscInt d, depth = 0, nleaves, loc, Ng, numFS, fs, fStart, fEnd, ghostCell, cEnd, c; in DMPlexConstructGhostCells_Internal() local
933 for (fs = 0; fs < numFS; ++fs) { in DMPlexConstructGhostCells_Internal()
938 PetscCall(DMLabelGetStratumIS(label, values[fs], &faceIS)); in DMPlexConstructGhostCells_Internal()
969 for (fs = 0; fs < numFS; ++fs) { in DMPlexConstructGhostCells_Internal()
974 PetscCall(DMLabelGetStratumIS(label, values[fs], &faceIS)); in DMPlexConstructGhostCells_Internal()
996 for (fs = 0; fs < numFS; ++fs) { in DMPlexConstructGhostCells_Internal()
1001 PetscCall(DMLabelGetStratumIS(label, values[fs], &faceIS)); in DMPlexConstructGhostCells_Internal()
/petsc/src/dm/dt/dualspace/impls/lagrange/
H A Ddspacelagrange.c1754 PetscInt fs[2]; in DMPlexPointIsTensor_Internal_Given() local
1759 fs[0] = f; in DMPlexPointIsTensor_Internal_Given()
1760 fs[1] = f2; in DMPlexPointIsTensor_Internal_Given()
1761 PetscCall(DMPlexGetMeet(dm, 2, fs, &meetSize, &meet)); in DMPlexPointIsTensor_Internal_Given()
1763 PetscCall(DMPlexRestoreMeet(dm, 2, fs, &meetSize, &meet)); in DMPlexPointIsTensor_Internal_Given()
/petsc/share/petsc/datafiles/meshes/
H A Dtestcase3D.cas4433 (vof/compressive/regime-discr/fs-limiter 1.5)
4557 (mp/fs-drag/viscous/shear? #t)
4558 (mp/fs-drag/tau-max 1000000000000000.)
4559 (mp/fs-drag/tau-min 1e-15)
4560 (mp/fs-drag/cd-max 100000.)
4561 (mp/fs-drag/cd-min 0.01)
4565 (mp/flow-regime/aiad/fs-grad-frac 0.1)
4566 (mp/flow-regime/aiad/ncells-fs 5)
9674 (beach-fs-level . 0)