Home
last modified time | relevance | path

Searched refs:dimC (Results 1 – 4 of 4) sorted by relevance

/petsc/src/dm/field/impls/shell/
H A Ddmfieldshell.c37 PetscInt dimC, dim, numPoints, Nq, p, Nc; in DMFieldShellEvaluateFEDefault() local
44 PetscCall(DMGetCoordinateDim(dm, &dimC)); in DMFieldShellEvaluateFEDefault()
47 PetscCall(PetscMalloc1(dimC * Nq * numPoints, &pfArray)); in DMFieldShellEvaluateFEDefault()
48 for (p = 0; p < numPoints * dimC * Nq; p++) pfArray[p] = geom->v[p]; in DMFieldShellEvaluateFEDefault()
49 …PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)pointIS), dimC, dimC * Nq * numPoints… in DMFieldShellEvaluateFEDefault()
63 for (i = 0; i < dimC; i++) d[i] = sD[(p * Nc + q) * dimC + i]; in DMFieldShellEvaluateFEDefault()
64 for (i = 0; i < dimC; i++) sD[(p * Nc + q) * dimC + i] = 0.; in DMFieldShellEvaluateFEDefault()
65 for (i = 0; i < dimC; i++) { in DMFieldShellEvaluateFEDefault()
66 …for (j = 0; j < dimC; j++) sD[(p * Nc + q) * dimC + i] += geom->J[(p * dimC + j) * dimC + i] * d[j… in DMFieldShellEvaluateFEDefault()
80 for (i = 0; i < dimC; i++) d[i] = rD[(p * Nc + q) * dimC + i]; in DMFieldShellEvaluateFEDefault()
[all …]
/petsc/src/dm/impls/plex/tests/
H A Dex22.c8 PetscInt i, j, dimC, dimR; in testIdentity() local
13 PetscCall(DMGetCoordinateDim(dm, &dimC)); in testIdentity()
16 PetscCall(DMGetWorkArray(dm, dimC * numPoints, MPIU_REAL, &mapped)); in testIdentity()
45 for (j = 0; j < dimC; j++) { in testIdentity()
46 PetscCall(PetscPrintf(PETSC_COMM_SELF, "%+f", (double)mapped[i * dimC + j])); in testIdentity()
47 if (j < dimC - 1) PetscCall(PetscPrintf(PETSC_COMM_SELF, ",")); in testIdentity()
71 for (j = 0; j < dimC; j++) { in testIdentity()
72 …PetscSNPrintfCount(strBuf + offset, BUFSIZ - offset, "%+f", &count, (double)mapped[i * dimC + j])); in testIdentity()
74 if (j < dimC - 1) { in testIdentity()
95 PetscCall(DMRestoreWorkArray(dm, dimC * numPoints, MPIU_REAL, &mapped)); in testIdentity()
[all …]
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c458 PetscInt Nq, Nc, dimC, qNc, N; in DMFieldEvaluateFV_DS() local
468 PetscCall(DMGetCoordinateDim(field->dm, &dimC)); in DMFieldEvaluateFV_DS()
492 if (D) PetscCall(DMGetWorkArray(field->dm, N * dimC, mpitype, &qD)); in DMFieldEvaluateFV_DS()
493 if (H) PetscCall(DMGetWorkArray(field->dm, N * dimC * dimC, mpitype, &qH)); in DMFieldEvaluateFV_DS()
538 for (j = 0; j < Nc * dimC; j++) sD[i * Nc * dimC + j] = 0.; in DMFieldEvaluateFV_DS()
544 for (l = 0; l < dimC; l++) { in DMFieldEvaluateFV_DS()
545 …for (m = 0; m < dim; m++) pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * sqD[((i * Nq… in DMFieldEvaluateFV_DS()
547 …for (l = 0; l < dimC; l++) sD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD… in DMFieldEvaluateFV_DS()
550 for (k = 0; k < Nc * dimC; k++) sD[i * Nc * dimC + k] /= vol; in DMFieldEvaluateFV_DS()
559 for (j = 0; j < Nc * dimC; j++) rD[i * Nc * dimC + j] = 0.; in DMFieldEvaluateFV_DS()
[all …]
/petsc/src/dm/impls/plex/
H A Dplexgeometry.c527 PetscInt degree, dimR, dimC; in DMPlexLocatePoint_Quad_2D_Internal() local
538 PetscCall(DMGetDimension(cdm, &dimC)); in DMPlexLocatePoint_Quad_2D_Internal()
554 for (PetscInt d = 0; d < dimC; d++) pointR[d] = PetscRealPart(point[d]); in DMPlexLocatePoint_Quad_2D_Internal()
555 …PetscCall(DMPlexCoordinatesToReference_FE(cdm, fe, c, 1, pointR, ref, coords, dimC, dimR, 10, &err… in DMPlexLocatePoint_Quad_2D_Internal()
559 PetscReal real[3], inverseError = 0, normPoint = DMPlex_NormD_Internal(dimC, pointR); in DMPlexLocatePoint_Quad_2D_Internal()
562 PetscCall(DMPlexReferenceToCoordinates_FE(cdm, fe, c, 1, ref, real, coords, dimC, dimR)); in DMPlexLocatePoint_Quad_2D_Internal()
563 inverseError = DMPlex_DistRealD_Internal(dimC, real, pointR); in DMPlexLocatePoint_Quad_2D_Internal()
643 PetscInt degree, dimR, dimC; in DMPlexLocatePoint_Hex_3D_Internal() local
654 PetscCall(DMGetDimension(cdm, &dimC)); in DMPlexLocatePoint_Hex_3D_Internal()
670 for (PetscInt d = 0; d < dimC; d++) pointR[d] = PetscRealPart(point[d]); in DMPlexLocatePoint_Hex_3D_Internal()
[all …]