Lines Matching refs:dimR
527 PetscInt degree, dimR, dimC; in DMPlexLocatePoint_Quad_2D_Internal() local
536 PetscCall(DMGetDimension(dm, &dimR)); 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()
562 PetscCall(DMPlexReferenceToCoordinates_FE(cdm, fe, c, 1, ref, real, coords, dimC, dimR)); in DMPlexLocatePoint_Quad_2D_Internal()
643 PetscInt degree, dimR, dimC; in DMPlexLocatePoint_Hex_3D_Internal() local
652 PetscCall(DMGetDimension(dm, &dimR)); in DMPlexLocatePoint_Hex_3D_Internal()
671 …PetscCall(DMPlexCoordinatesToReference_FE(cdm, fe, c, 1, pointR, ref, coords, dimC, dimR, 10, &err… in DMPlexLocatePoint_Hex_3D_Internal()
678 PetscCall(DMPlexReferenceToCoordinates_FE(cdm, fe, c, 1, ref, real, coords, dimC, dimR)); in DMPlexLocatePoint_Hex_3D_Internal()
1985 const PetscInt dimR = 2; in DMPlexComputeRectangleGeometry_Internal() local
2018 PetscReal xi = points[dimR * i], eta = points[dimR * i + 1]; in DMPlexComputeRectangleGeometry_Internal()
2046 for (k = 0; k < dimR; k++) { in DMPlexComputeRectangleGeometry_Internal()
2049 for (l = 0; l < Nv; l++) val += zCoeff[dim * l + j] * extJ[dimR * l + k]; in DMPlexComputeRectangleGeometry_Internal()
2138 const PetscInt dimR = 3; in DMPlexComputeHexahedronGeometry_Internal() local
2159 PetscReal xi = points[dimR * i], eta = points[dimR * i + 1], theta = points[dimR * i + 2]; in DMPlexComputeHexahedronGeometry_Internal()
2208 for (k = 0; k < dimR; k++) { in DMPlexComputeHexahedronGeometry_Internal()
2211 for (l = 0; l < Nv; l++) val += zCoeff[dim * l + j] * extJ[dimR * l + k]; in DMPlexComputeHexahedronGeometry_Internal()
2253 const PetscInt dimR = 3; in DMPlexComputeTriangularPrismGeometry_Internal() local
2307 … const PetscReal xi = points[dimR * i], eta = points[dimR * i + 1], zeta = points[dimR * i + 2]; in DMPlexComputeTriangularPrismGeometry_Internal()
2349 for (k = 0; k < dimR; k++) { in DMPlexComputeTriangularPrismGeometry_Internal()
2352 for (l = 0; l < Nv; l++) val += coeff[dim * l + j] * extJ[dimR * l + k]; in DMPlexComputeTriangularPrismGeometry_Internal()
3492 static PetscErrorCode DMPlexCoordinatesToReference_NewtonUpdate(PetscInt dimC, PetscInt dimR, Petsc… in DMPlexCoordinatesToReference_NewtonUpdate() argument
3497 if (dimC == dimR && dimR <= 3) { in DMPlexCoordinatesToReference_NewtonUpdate()
3501 switch (dimR) { in DMPlexCoordinatesToReference_NewtonUpdate()
3530 for (l = 0; l < dimR; l++) { in DMPlexCoordinatesToReference_NewtonUpdate()
3537 PetscCall(PetscBLASIntCast(dimR, &m)); in DMPlexCoordinatesToReference_NewtonUpdate()
3545 for (l = 0; l < dimR; l++) guess[l] += PetscRealPart(invJ[l]); in DMPlexCoordinatesToReference_NewtonUpdate()
3550 …nts, const PetscReal realCoords[], PetscReal refCoords[], Vec coords, PetscInt dimC, PetscInt dimR) in DMPlexCoordinatesToReference_Tensor() argument
3552 PetscInt coordSize, i, j, k, l, m, maxIts = 7, numV = (1 << dimR); in DMPlexCoordinatesToReference_Tensor()
3560 …ecting at least %" PetscInt_FMT " coordinates, got %" PetscInt_FMT, dimC * (1 << dimR), coordSize); in DMPlexCoordinatesToReference_Tensor()
3561 PetscCall(DMGetWorkArray(dm, 2 * coordSize + dimR + dimC, MPIU_REAL, &cellData)); in DMPlexCoordinatesToReference_Tensor()
3562 PetscCall(DMGetWorkArray(dm, 3 * dimR * dimC, MPIU_SCALAR, &J)); in DMPlexCoordinatesToReference_Tensor()
3566 resNeg = &cellData[2 * coordSize + dimR]; in DMPlexCoordinatesToReference_Tensor()
3567 invJ = &J[dimR * dimC]; in DMPlexCoordinatesToReference_Tensor()
3568 work = &J[2 * dimR * dimC]; in DMPlexCoordinatesToReference_Tensor()
3569 if (dimR == 2) { in DMPlexCoordinatesToReference_Tensor()
3577 } else if (dimR == 3) { in DMPlexCoordinatesToReference_Tensor()
3589 for (i = 0; i < dimR; i++) { in DMPlexCoordinatesToReference_Tensor()
3599 if (i < dimR - 1) { in DMPlexCoordinatesToReference_Tensor()
3605 PetscCall(PetscArrayzero(refCoords, numPoints * dimR)); in DMPlexCoordinatesToReference_Tensor()
3608 PetscReal *guess = &refCoords[dimR * j]; in DMPlexCoordinatesToReference_Tensor()
3612 for (k = 0; k < dimC * dimR; k++) J[k] = 0.; in DMPlexCoordinatesToReference_Tensor()
3615 for (l = 0; l < dimR; l++) { in DMPlexCoordinatesToReference_Tensor()
3622 for (m = 0; m < dimR; m++) { in DMPlexCoordinatesToReference_Tensor()
3634 for (m = 0; m < dimR; m++) J[dimR * l + m] += coeff * extJ[m]; in DMPlexCoordinatesToReference_Tensor()
3644 … PetscCall(DMPlexCoordinatesToReference_NewtonUpdate(dimC, dimR, J, invJ, work, resNeg, guess)); in DMPlexCoordinatesToReference_Tensor()
3647 PetscCall(DMRestoreWorkArray(dm, 3 * dimR * dimC, MPIU_SCALAR, &J)); in DMPlexCoordinatesToReference_Tensor()
3648 PetscCall(DMRestoreWorkArray(dm, 2 * coordSize + dimR + dimC, MPIU_REAL, &cellData)); in DMPlexCoordinatesToReference_Tensor()
3653 …nts, const PetscReal refCoords[], PetscReal realCoords[], Vec coords, PetscInt dimC, PetscInt dimR) in DMPlexReferenceToCoordinates_Tensor() argument
3655 PetscInt coordSize, i, j, k, l, numV = (1 << dimR); in DMPlexReferenceToCoordinates_Tensor()
3662 …ecting at least %" PetscInt_FMT " coordinates, got %" PetscInt_FMT, dimC * (1 << dimR), coordSize); in DMPlexReferenceToCoordinates_Tensor()
3666 if (dimR == 2) { in DMPlexReferenceToCoordinates_Tensor()
3674 } else if (dimR == 3) { in DMPlexReferenceToCoordinates_Tensor()
3686 for (i = 0; i < dimR; i++) { in DMPlexReferenceToCoordinates_Tensor()
3696 if (i < dimR - 1) { in DMPlexReferenceToCoordinates_Tensor()
3704 const PetscReal *guess = &refCoords[dimR * j]; in DMPlexReferenceToCoordinates_Tensor()
3709 for (l = 0; l < dimR; l++) { in DMPlexReferenceToCoordinates_Tensor()
3727 …al realCoords[], PetscReal refCoords[], Vec coords, PetscInt Nc, PetscInt dimR, PetscInt maxIter, … in DMPlexCoordinatesToReference_FE() argument
3749 PetscCall(DMGetWorkArray(dm, pdim * Nc + pdim * Nc * dimR + Nc, MPIU_REAL, &B)); in DMPlexCoordinatesToReference_FE()
3751 resNeg = &D[pdim * Nc * dimR]; in DMPlexCoordinatesToReference_FE()
3752 PetscCall(DMGetWorkArray(dm, 3 * Nc * dimR, MPIU_SCALAR, &J)); in DMPlexCoordinatesToReference_FE()
3753 invJ = &J[Nc * dimR]; in DMPlexCoordinatesToReference_FE()
3754 work = &invJ[Nc * dimR]; in DMPlexCoordinatesToReference_FE()
3755 for (i = 0; i < numPoints * dimR; i++) refCoords[i] = 0.; in DMPlexCoordinatesToReference_FE()
3760 PetscReal *guess = &refCoords[j * dimR], error = 0; in DMPlexCoordinatesToReference_FE()
3763 for (k = 0; k < Nc * dimR; k++) J[k] = 0.; in DMPlexCoordinatesToReference_FE()
3767 for (m = 0; m < dimR; m++) J[l * dimR + m] += modes[k] * D[(k * Nc + l) * dimR + m]; in DMPlexCoordinatesToReference_FE()
3781 PetscCall(DMPlexCoordinatesToReference_NewtonUpdate(Nc, dimR, J, invJ, work, resNeg, guess)); in DMPlexCoordinatesToReference_FE()
3784 PetscCall(DMRestoreWorkArray(dm, 3 * Nc * dimR, MPIU_SCALAR, &J)); in DMPlexCoordinatesToReference_FE()
3785 PetscCall(DMRestoreWorkArray(dm, pdim * Nc + pdim * Nc * dimR + Nc, MPIU_REAL, &B)); in DMPlexCoordinatesToReference_FE()
3792 …oints, const PetscReal refCoords[], PetscReal realCoords[], Vec coords, PetscInt Nc, PetscInt dimR) in DMPlexReferenceToCoordinates_FE() argument
3855 PetscInt dimC, dimR, depth, i, cellHeight, height; in DMPlexCoordinatesToReference() local
3863 PetscCall(DMGetDimension(dm, &dimR)); in DMPlexCoordinatesToReference()
3865 if (dimR <= 0 || dimC <= 0 || numPoints <= 0) PetscFunctionReturn(PETSC_SUCCESS); in DMPlexCoordinatesToReference()
3892 isSimplex = (coneSize == (dimR + 1)) ? PETSC_TRUE : PETSC_FALSE; in DMPlexCoordinatesToReference()
3893 isTensor = (coneSize == ((depth == 1) ? (1 << dimR) : (2 * dimR))) ? PETSC_TRUE : PETSC_FALSE; in DMPlexCoordinatesToReference()
3904 CoordinatesRealToRef(dimC, dimR, x0, v0, invJ, &realCoords[dimC * i], &refCoords[dimR * i]); in DMPlexCoordinatesToReference()
3908 …oordinatesToReference_Tensor(coordDM, cell, numPoints, realCoords, refCoords, coords, dimC, dimR)); in DMPlexCoordinatesToReference()
3911 …sToReference_FE(coordDM, fe, cell, numPoints, realCoords, refCoords, coords, dimC, dimR, 7, NULL)); in DMPlexCoordinatesToReference()
3938 PetscInt dimC, dimR, depth, i, cellHeight, height; in DMPlexReferenceToCoordinates() local
3946 PetscCall(DMGetDimension(dm, &dimR)); in DMPlexReferenceToCoordinates()
3948 if (dimR <= 0 || dimC <= 0 || numPoints <= 0) PetscFunctionReturn(PETSC_SUCCESS); in DMPlexReferenceToCoordinates()
3975 isSimplex = (coneSize == (dimR + 1)) ? PETSC_TRUE : PETSC_FALSE; in DMPlexReferenceToCoordinates()
3976 isTensor = (coneSize == ((depth == 1) ? (1 << dimR) : (2 * dimR))) ? PETSC_TRUE : PETSC_FALSE; in DMPlexReferenceToCoordinates()
3986 CoordinatesRefToReal(dimC, dimR, xi0, v0, J, &refCoords[dimR * i], &realCoords[dimC * i]); in DMPlexReferenceToCoordinates()
3990 …eferenceToCoordinates_Tensor(coordDM, cell, numPoints, refCoords, realCoords, coords, dimC, dimR)); in DMPlexReferenceToCoordinates()
3993 …eferenceToCoordinates_FE(coordDM, fe, cell, numPoints, refCoords, realCoords, coords, dimC, dimR)); in DMPlexReferenceToCoordinates()