Lines Matching refs:Nc
186 …" PetscInt_FMT "\n", PetscAbsInt(sp->k), sp->k < 0 ? "(stored in dual form) " : "", sp->Nc, pdim)); in PetscDualSpaceView_ASCII()
188 …Printf(v, "Dual space with %" PetscInt_FMT " components, size %" PetscInt_FMT "\n", sp->Nc, pdim)); in PetscDualSpaceView_ASCII()
308 …ponents", "The number of components", "PetscDualSpaceSetNumComponents", sp->Nc, &sp->Nc, NULL, 1)); in PetscDualSpaceSetFromOptions()
454 s->Nc = 1; in PetscDualSpaceCreate()
498 (*spNew)->Nc = sp->Nc; in PetscDualSpaceDuplicate()
616 PetscErrorCode PetscDualSpaceGetNumComponents(PetscDualSpace sp, PetscInt *Nc) in PetscDualSpaceGetNumComponents() argument
620 PetscAssertPointer(Nc, 2); in PetscDualSpaceGetNumComponents()
621 *Nc = sp->Nc; in PetscDualSpaceGetNumComponents()
636 PetscErrorCode PetscDualSpaceSetNumComponents(PetscDualSpace sp, PetscInt Nc) in PetscDualSpaceSetNumComponents() argument
641 sp->Nc = Nc; in PetscDualSpaceSetNumComponents()
1137 …etscDualSpace sp, PetscInt f, PetscReal time, PetscFEGeom *cgeom, PetscInt Nc, PetscErrorCode (*fu… in PetscDualSpaceApplyDefault() argument
1154 …= Nc, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_SIZ, "The quadrature components %" PetscInt_… in PetscDualSpaceApplyDefault()
1155 PetscCall(DMGetWorkArray(dm, Nc, MPIU_SCALAR, &val)); in PetscDualSpaceApplyDefault()
1162 PetscCall((*func)(dE, time, x, Nc, val, ctx)); in PetscDualSpaceApplyDefault()
1164 PetscCall((*func)(dE, time, &cgeom->v[dE * q], Nc, val, ctx)); in PetscDualSpaceApplyDefault()
1166 for (c = 0; c < Nc; ++c) *value += val[c] * weights[q * Nc + c]; in PetscDualSpaceApplyDefault()
1168 PetscCall(DMRestoreWorkArray(dm, Nc, MPIU_SCALAR, &val)); in PetscDualSpaceApplyDefault()
1299 PetscInt Nc, nrows, ncols; in PetscDualSpaceCreateAllDataDefault() local
1305 PetscCall(PetscDualSpaceGetNumComponents(sp, &Nc)); in PetscDualSpaceCreateAllDataDefault()
1323 ncols = numPoints * Nc; in PetscDualSpaceCreateAllDataDefault()
1325 PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF, nrows, ncols, maxNumPoints * Nc, NULL, &A)); in PetscDualSpaceCreateAllDataDefault()
1333 PetscCheck(fnc == Nc, PETSC_COMM_SELF, PETSC_ERR_PLIB, "functional component mismatch"); in PetscDualSpaceCreateAllDataDefault()
1335 for (i = 0; i < Np * Nc; i++) PetscCall(MatSetValue(A, f, offset * Nc, w[i], INSERT_VALUES)); in PetscDualSpaceCreateAllDataDefault()
1413 PetscInt Nc; in PetscDualSpaceCreateInteriorDataDefault() local
1432 PetscCall(PetscDualSpaceGetNumComponents(sp, &Nc)); in PetscDualSpaceCreateInteriorDataDefault()
1449 nnz[f] = Np * Nc; in PetscDualSpaceCreateInteriorDataDefault()
1453 PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF, spdim0, numPoints * Nc, 0, nnz, &imat)); in PetscDualSpaceCreateInteriorDataDefault()
1471 … for (i = 0; i < Np * Nc; i++) PetscCall(MatSetValue(imat, f, matoffset + i, w[i], INSERT_VALUES)); in PetscDualSpaceCreateInteriorDataDefault()
1473 matoffset += Np * Nc; in PetscDualSpaceCreateInteriorDataDefault()
1564 …DualSpace sp, PetscInt f, PetscReal time, PetscFVCellGeom *cgeom, PetscInt Nc, PetscErrorCode (*fu… in PetscDualSpaceApplyFVM() argument
1579 …= Nc, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_SIZ, "The quadrature components %" PetscInt_… in PetscDualSpaceApplyFVM()
1580 PetscCall(DMGetWorkArray(dm, Nc, MPIU_SCALAR, &val)); in PetscDualSpaceApplyFVM()
1583 PetscCall((*func)(dimEmbed, time, cgeom->centroid, Nc, val, ctx)); in PetscDualSpaceApplyFVM()
1584 for (c = 0; c < Nc; ++c) *value += val[c] * weights[q * Nc + c]; in PetscDualSpaceApplyFVM()
1586 PetscCall(DMRestoreWorkArray(dm, Nc, MPIU_SCALAR, &val)); in PetscDualSpaceApplyFVM()
1884 …Type trans, PetscBool isInverse, PetscFEGeom *fegeom, PetscInt Nv, PetscInt Nc, PetscScalar vals[]) in PetscDualSpaceTransform() argument
1903 for (c = 0; c < Nc; c++) vals[v * Nc + c] *= Jstar[0]; in PetscDualSpaceTransform()
1906 …for (c = 0; c < Nc; c += 2) DMPlex_Mult2DReal_Internal(Jstar, 1, &vals[v * Nc + c], &vals[v * Nc +… in PetscDualSpaceTransform()
1909 …for (c = 0; c < Nc; c += 3) DMPlex_Mult3DReal_Internal(Jstar, 1, &vals[v * Nc + c], &vals[v * Nc +… in PetscDualSpaceTransform()
1940 …Type trans, PetscBool isInverse, PetscFEGeom *fegeom, PetscInt Nv, PetscInt Nc, PetscScalar vals[]) in PetscDualSpaceTransformGradient() argument
1953 for (c = 0; c < Nc; ++c) { in PetscDualSpaceTransformGradient()
1956 vals[(v * Nc + c) * dim] *= fegeom->invJ[0]; in PetscDualSpaceTransformGradient()
1959 …Plex_MultTranspose2DReal_Internal(fegeom->invJ, 1, &vals[(v * Nc + c) * dim], &vals[(v * Nc + c) *… in PetscDualSpaceTransformGradient()
1962 …Plex_MultTranspose3DReal_Internal(fegeom->invJ, 1, &vals[(v * Nc + c) * dim], &vals[(v * Nc + c) *… in PetscDualSpaceTransformGradient()
1971 …for (c = 0; c < Nc; ++c) DMPlex_MultTransposeReal_Internal(fegeom->invJ, dim, dE, 1, &vals[(v * Nc… in PetscDualSpaceTransformGradient()
1975 if (Nc == 1 || Nc != dim) PetscFunctionReturn(PETSC_SUCCESS); in PetscDualSpaceTransformGradient()
1985 …DMPlex_MultTranspose2DReal_Internal(fegeom->J, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + … in PetscDualSpaceTransformGradient()
1988 …DMPlex_MultTranspose3DReal_Internal(fegeom->J, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + … in PetscDualSpaceTransformGradient()
2000 …Plex_MultTranspose2DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim +… in PetscDualSpaceTransformGradient()
2003 …Plex_MultTranspose3DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim +… in PetscDualSpaceTransformGradient()
2018 … DMPlex_Mult2DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + d]); in PetscDualSpaceTransformGradient()
2021 … DMPlex_Mult3DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + d]); in PetscDualSpaceTransformGradient()
2026 for (c = 0; c < Nc; ++c) vals[(v * Nc + c) * dim + d] *= fegeom->detJ[0]; in PetscDualSpaceTransformGradient()
2034 … DMPlex_Mult2DReal_Internal(fegeom->J, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + d]); in PetscDualSpaceTransformGradient()
2037 … DMPlex_Mult3DReal_Internal(fegeom->J, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + d]); in PetscDualSpaceTransformGradient()
2042 for (c = 0; c < Nc; ++c) vals[(v * Nc + c) * dim + d] /= fegeom->detJ[0]; in PetscDualSpaceTransformGradient()
2073 …Type trans, PetscBool isInverse, PetscFEGeom *fegeom, PetscInt Nv, PetscInt Nc, PetscScalar vals[]) in PetscDualSpaceTransformHessian() argument
2086 for (c = 0; c < Nc; ++c) { in PetscDualSpaceTransformHessian()
2089 vals[(v * Nc + c) * dim * dim] *= PetscSqr(fegeom->invJ[0]); in PetscDualSpaceTransformHessian()
2092 …DMPlex_PTAP2DReal_Internal(fegeom->invJ, &vals[(v * Nc + c) * dim * dim], &vals[(v * Nc + c) * dim… in PetscDualSpaceTransformHessian()
2095 …DMPlex_PTAP3DReal_Internal(fegeom->invJ, &vals[(v * Nc + c) * dim * dim], &vals[(v * Nc + c) * dim… in PetscDualSpaceTransformHessian()
2104 …for (c = 0; c < Nc; ++c) DMPlex_PTAPReal_Internal(fegeom->invJ, dim, dE, &vals[(v * Nc + c) * dE *… in PetscDualSpaceTransformHessian()
2108 if (Nc == 1 || Nc != dim) PetscFunctionReturn(PETSC_SUCCESS); in PetscDualSpaceTransformHessian()
2142 …acePullback(PetscDualSpace dsp, PetscFEGeom *fegeom, PetscInt Nq, PetscInt Nc, PetscScalar pointEv… in PetscDualSpacePullback() argument
2168 PetscCall(PetscDualSpaceTransform(dsp, trans, PETSC_TRUE, fegeom, Nq, Nc, pointEval)); in PetscDualSpacePullback()
2194 …Pushforward(PetscDualSpace dsp, PetscFEGeom *fegeom, PetscInt Nq, PetscInt Nc, PetscScalar pointEv… in PetscDualSpacePushforward() argument
2220 PetscCall(PetscDualSpaceTransform(dsp, trans, PETSC_FALSE, fegeom, Nq, Nc, pointEval)); in PetscDualSpacePushforward()
2246 …ardGradient(PetscDualSpace dsp, PetscFEGeom *fegeom, PetscInt Nq, PetscInt Nc, PetscScalar pointEv… in PetscDualSpacePushforwardGradient() argument
2272 PetscCall(PetscDualSpaceTransformGradient(dsp, trans, PETSC_FALSE, fegeom, Nq, Nc, pointEval)); in PetscDualSpacePushforwardGradient()
2298 …wardHessian(PetscDualSpace dsp, PetscFEGeom *fegeom, PetscInt Nq, PetscInt Nc, PetscScalar pointEv… in PetscDualSpacePushforwardHessian() argument
2324 PetscCall(PetscDualSpaceTransformHessian(dsp, trans, PETSC_FALSE, fegeom, Nq, Nc, pointEval)); in PetscDualSpacePushforwardHessian()