Lines Matching refs:dim
847 PetscInt dim; in DMPlexP1FieldCreate_Private() local
852 PetscCall(DMGetDimension(dm, &dim)); in DMPlexP1FieldCreate_Private()
855 PetscCall(PetscFECreateLagrange(comm, dim, size, PETSC_TRUE, 1, PETSC_DETERMINE, &fe)); in DMPlexP1FieldCreate_Private()
955 static void identity(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscI… in identity() argument
1034 static PetscErrorCode LAPACKsyevFail(PetscInt dim, PetscScalar Mpos[]) in LAPACKsyevFail() argument
1040 for (i = 0; i < dim; ++i) { in LAPACKsyevFail()
1043 for (j = 0; j < dim; ++j) { in LAPACKsyevFail()
1044 …if (j < dim - 1) PetscCall(PetscPrintf(PETSC_COMM_SELF, "%15.8e, ", (double)PetscAbsScalar(Mpos[i … in LAPACKsyevFail()
1045 … else PetscCall(PetscPrintf(PETSC_COMM_SELF, "%15.8e", (double)PetscAbsScalar(Mpos[i * dim + j]))); in LAPACKsyevFail()
1047 if (i < dim - 1) PetscCall(PetscPrintf(PETSC_COMM_SELF, "]\n")); in LAPACKsyevFail()
1053 static PetscErrorCode DMPlexMetricModify_Private(PetscInt dim, PetscReal h_min, PetscReal h_max, Pe… in DMPlexMetricModify_Private() argument
1060 PetscCall(PetscMalloc2(dim * dim, &Mpos, dim, &eigs)); in DMPlexMetricModify_Private()
1063 for (i = 0; i < dim; ++i) { in DMPlexMetricModify_Private()
1064 Mpos[i * dim + i] = Mp[i * dim + i]; in DMPlexMetricModify_Private()
1065 for (j = i + 1; j < dim; ++j) { in DMPlexMetricModify_Private()
1066 Mpos[i * dim + j] = 0.5 * (Mp[i * dim + j] + Mp[j * dim + i]); in DMPlexMetricModify_Private()
1067 Mpos[j * dim + i] = Mpos[i * dim + j]; in DMPlexMetricModify_Private()
1072 if (dim == 1) { in DMPlexMetricModify_Private()
1081 PetscCall(PetscBLASIntCast(5 * dim, &lwork)); in DMPlexMetricModify_Private()
1082 PetscCall(PetscMalloc1(5 * dim, &work)); in DMPlexMetricModify_Private()
1087 PetscCall(PetscBLASIntCast(dim, &nb)); in DMPlexMetricModify_Private()
1092 PetscCall(PetscMalloc1(3 * dim, &rwork)); in DMPlexMetricModify_Private()
1100 for (i = 0; i < dim; ++i) { in DMPlexMetricModify_Private()
1101 Mpos[i * dim + i] = Mp[i * dim + i]; in DMPlexMetricModify_Private()
1102 for (j = i + 1; j < dim; ++j) { in DMPlexMetricModify_Private()
1103 Mpos[i * dim + j] = 0.5 * (Mp[i * dim + j] + Mp[j * dim + i]); in DMPlexMetricModify_Private()
1104 Mpos[j * dim + i] = Mpos[i * dim + j]; in DMPlexMetricModify_Private()
1107 PetscCall(LAPACKsyevFail(dim, Mpos)); in DMPlexMetricModify_Private()
1117 for (i = 0; i < dim; ++i) { in DMPlexMetricModify_Private()
1124 for (i = 0; i < dim; ++i) { in DMPlexMetricModify_Private()
1130 for (i = 0; i < dim; ++i) { in DMPlexMetricModify_Private()
1131 for (j = 0; j < dim; ++j) { in DMPlexMetricModify_Private()
1132 Mp[i * dim + j] = 0.0; in DMPlexMetricModify_Private()
1133 for (k = 0; k < dim; ++k) Mp[i * dim + j] += Mpos[k * dim + i] * eigs[k] * Mpos[k * dim + j]; in DMPlexMetricModify_Private()
1168 PetscInt dim, vStart, vEnd, v; in DMPlexMetricEnforceSPD() local
1176 PetscCall(DMGetDimension(dm, &dim)); in DMPlexMetricEnforceSPD()
1208 else nrow = dim; in DMPlexMetricEnforceSPD()
1226 static void detMFunc(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscI… in detMFunc() argument
1230 f0[0] = PetscPowScalar(u[0], p / (2.0 * p + dim)); in detMFunc()
1266 PetscInt dim, Nd, vStart, vEnd, v, i; in DMPlexMetricNormalize() local
1275 PetscCall(DMGetDimension(dm, &dim)); in DMPlexMetricNormalize()
1279 else Nd = dim * dim; in DMPlexMetricNormalize()
1314 factGlob = PetscPowReal(target / realIntegral, 2.0 / dim); in DMPlexMetricNormalize()
1332 met[0] *= factGlob * PetscPowReal(PetscAbsScalar(det[0]), -1.0 / (2 * p + dim)); in DMPlexMetricNormalize()
1339 else nrow = dim; in DMPlexMetricNormalize()
1347 fact = factGlob * PetscPowReal(PetscAbsScalar(detM[0]), -1.0 / (2 * p + dim)); in DMPlexMetricNormalize()
1466 static PetscErrorCode DMPlexMetricIntersection_Private(PetscInt dim, PetscScalar M1[], PetscScalar … in DMPlexMetricIntersection_Private() argument
1474 if (dim == 1) { in DMPlexMetricIntersection_Private()
1480 PetscCall(PetscMalloc4(dim * dim, &evecs, dim * dim, &sqrtM1, dim * dim, &isqrtM1, dim, &evals)); in DMPlexMetricIntersection_Private()
1481 for (i = 0; i < dim; ++i) { in DMPlexMetricIntersection_Private()
1482 for (j = 0; j < dim; ++j) evecs[i * dim + j] = M1[i * dim + j]; in DMPlexMetricIntersection_Private()
1488 PetscCall(PetscBLASIntCast(5 * dim, &lwork)); in DMPlexMetricIntersection_Private()
1489 PetscCall(PetscMalloc1(5 * dim, &work)); in DMPlexMetricIntersection_Private()
1495 PetscCall(PetscBLASIntCast(dim, &nb)); in DMPlexMetricIntersection_Private()
1500 PetscCall(PetscMalloc1(3 * dim, &rwork)); in DMPlexMetricIntersection_Private()
1508 PetscCall(LAPACKsyevFail(dim, M1)); in DMPlexMetricIntersection_Private()
1514 for (i = 0; i < dim; ++i) { in DMPlexMetricIntersection_Private()
1515 for (k = 0; k < dim; ++k) { in DMPlexMetricIntersection_Private()
1516 sqrtM1[i * dim + k] = 0.0; in DMPlexMetricIntersection_Private()
1517 isqrtM1[i * dim + k] = 0.0; in DMPlexMetricIntersection_Private()
1518 for (j = 0; j < dim; ++j) { in DMPlexMetricIntersection_Private()
1520 sqrtM1[i * dim + k] += evecs[j * dim + i] * sqrtj * evecs[j * dim + k]; in DMPlexMetricIntersection_Private()
1521 isqrtM1[i * dim + k] += evecs[j * dim + i] * (1.0 / sqrtj) * evecs[j * dim + k]; in DMPlexMetricIntersection_Private()
1527 for (i = 0; i < dim; ++i) { in DMPlexMetricIntersection_Private()
1528 for (l = 0; l < dim; ++l) { in DMPlexMetricIntersection_Private()
1529 evecs[i * dim + l] = 0.0; in DMPlexMetricIntersection_Private()
1530 for (j = 0; j < dim; ++j) { in DMPlexMetricIntersection_Private()
1531 …for (k = 0; k < dim; ++k) evecs[i * dim + l] += isqrtM1[j * dim + i] * M2[j * dim + k] * isqrtM1[k… in DMPlexMetricIntersection_Private()
1541 PetscCall(PetscMalloc1(3 * dim, &rwork)); in DMPlexMetricIntersection_Private()
1549 for (i = 0; i < dim; ++i) { in DMPlexMetricIntersection_Private()
1550 for (l = 0; l < dim; ++l) { in DMPlexMetricIntersection_Private()
1551 evecs[i * dim + l] = 0.0; in DMPlexMetricIntersection_Private()
1552 for (j = 0; j < dim; ++j) { in DMPlexMetricIntersection_Private()
1553 …for (k = 0; k < dim; ++k) evecs[i * dim + l] += isqrtM1[j * dim + i] * M2[j * dim + k] * isqrtM1[k… in DMPlexMetricIntersection_Private()
1557 PetscCall(LAPACKsyevFail(dim, evecs)); in DMPlexMetricIntersection_Private()
1563 for (i = 0; i < dim; ++i) evals[i] = PetscMax(evals[i], 1.0); in DMPlexMetricIntersection_Private()
1566 for (i = 0; i < dim; ++i) { in DMPlexMetricIntersection_Private()
1567 for (m = 0; m < dim; ++m) { in DMPlexMetricIntersection_Private()
1568 M2[i * dim + m] = 0.0; in DMPlexMetricIntersection_Private()
1569 for (j = 0; j < dim; ++j) { in DMPlexMetricIntersection_Private()
1570 for (k = 0; k < dim; ++k) { in DMPlexMetricIntersection_Private()
1571 …for (l = 0; l < dim; ++l) M2[i * dim + m] += sqrtM1[j * dim + i] * evecs[j * dim + k] * evals[k] *… in DMPlexMetricIntersection_Private()
1636 PetscInt dim, vStart, vEnd, nrow; in DMPlexMetricIntersection() local
1639 PetscCall(DMGetDimension(dm, &dim)); in DMPlexMetricIntersection()
1641 else nrow = dim; in DMPlexMetricIntersection()