| /petsc/src/tao/quadratic/impls/bqpip/ |
| H A D | bqpip.c | 10 static PetscErrorCode QPIPComputeResidual(TAO_BQPIP *qp, Tao tao) in QPIPComputeResidual() argument
12 PetscReal dtmp = 1.0 - qp->psteplength; in QPIPComputeResidual()
17 PetscCall(VecScale(qp->R3, dtmp)); in QPIPComputeResidual()
18 PetscCall(VecScale(qp->R5, dtmp)); in QPIPComputeResidual()
19 qp->pinfeas = dtmp * qp->pinfeas; in QPIPComputeResidual()
21 PetscCall(VecCopy(qp->S, tao->gradient)); in QPIPComputeResidual()
22 PetscCall(VecAXPY(tao->gradient, -1.0, qp->Z)); in QPIPComputeResidual()
24 PetscCall(MatMult(tao->hessian, tao->solution, qp->RHS)); in QPIPComputeResidual()
25 PetscCall(VecScale(qp->RHS, -1.0)); in QPIPComputeResidual()
26 PetscCall(VecAXPY(qp->RHS, -1.0, qp->C)); in QPIPComputeResidual()
[all …]
|
| /petsc/src/ts/utils/dmplexlandau/tutorials/output/ |
| H A D | ex1_re.out | 49 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=5.508970e-02, cell 0, qp 1/9
50 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=4.379594e-02, cell 0, qp 2/9
51 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=3.250658e-02, cell 0, qp 3/9
52 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=2.594060e-02, cell 4, qp 1/9
53 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=1.468071e-02, cell 4, qp 2/9
54 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=3.676391e-03, cell 4, qp 3/9
55 [0] <vec:seq> adaptToleranceFEM(): 1) Found another inner r=3.676391e-03, cell 8, qp 1/9, d=0.00…
59 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=5.508970e-02, cell 0, qp 1/9
60 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=4.379594e-02, cell 0, qp 2/9
61 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=3.250658e-02, cell 0, qp 3/9
[all …]
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex51.c | 407 PetscReal z1, z, xm, xl, q, qp, Ln, scale; in leggaulob() local
426 qAndLEvaluation(n, z, &q, &qp, &Ln); in leggaulob()
428 z = z1 - q / qp; /* Newton's method. */ in leggaulob()
430 qAndLEvaluation(n, z, &q, &qp, &Ln); in leggaulob()
438 qAndLEvaluation(n, 0.0, &q, &qp, &Ln); in leggaulob()
448 static void qAndLEvaluation(PetscInt n, PetscReal x, PetscReal *q, PetscReal *qp, PetscReal *Ln) in qAndLEvaluation() argument
480 *qp = Lnp1p - Lnm1p; in qAndLEvaluation()
|
| /petsc/src/dm/dt/space/impls/subspace/ |
| H A D | spacesubspace.c | 316 const PetscReal *qp; in PetscSpaceSetUp_Subspace() local
320 PetscCall(PetscQuadratureGetData(q, NULL, NULL, &qNp, &qp, &qw)); in PetscSpaceSetUp_Subspace()
329 … for (j = 0; j < subDim; j++) allPoints[origDim * offset + i] += Jx[i * subDim + j] * qp[j]; in PetscSpaceSetUp_Subspace()
332 for (i = 0; i < PetscMin(subDim, origDim); i++) allPoints[origDim * offset + i] += qp[i]; in PetscSpaceSetUp_Subspace()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex11.c | 161 static void PhysicsRiemann_Advect(PetscInt dim, PetscInt Nf, const PetscReal *qp, const PetscReal *… in PhysicsRiemann_Advect() argument
173 wind[0] = -qp[1]; in PhysicsRiemann_Advect()
174 wind[1] = qp[0]; in PhysicsRiemann_Advect()
182 comp2[i] = qp[i] * qp[i]; in PhysicsRiemann_Advect()
186 wind[0] = -qp[1]; in PhysicsRiemann_Advect()
187 wind[1] = qp[0]; in PhysicsRiemann_Advect()
|
| H A D | ex11.h | 170 static void PhysicsRiemann_SW_Rusanov(PetscInt dim, PetscInt Nf, const PetscReal *qp, const PetscRe… in PhysicsRiemann_SW_Rusanov() argument
270 static void PhysicsRiemann_SW_HLL(PetscInt dim, PetscInt Nf, const PetscReal *qp, const PetscReal *… in PhysicsRiemann_SW_HLL() argument
784 static void PhysicsRiemann_Euler_Godunov(PetscInt dim, PetscInt Nf, const PetscReal *qp, const Pets… in PhysicsRiemann_Euler_Godunov() argument
|
| H A D | ex18.c | 286 static void riemann_advection(PetscInt dim, PetscInt Nf, const PetscReal *qp, const PetscReal *n, c… in riemann_advection() argument
294 static void riemann_coupled_advection(PetscInt dim, PetscInt Nf, const PetscReal *qp, const PetscRe… in riemann_coupled_advection() argument
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 291 PetscScalar qp; in SymBroydenCompactDenseUpdateArrays() local
305 PetscCallBLAS("BLASdot", qp = BLASdot_(&m, q0, &ione, p0, &ione)); in SymBroydenCompactDenseUpdateArrays()
306 PetscCallBLAS("BLASdot", qp += BLASdot_(&m, q1, &ione, p1, &ione)); in SymBroydenCompactDenseUpdateArrays()
308 sBis = StB0S[i + i * ldasbs] + qp; in SymBroydenCompactDenseUpdateArrays()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexsubmesh.c | 1336 PetscInt coneSize, supportSize, qf, qn, qp, e; in DMPlexConstructCohesiveCells_Internal() local
1375 for (e = 0, qn = 0, qp = 0, qf = 0; e < supportSize; ++e) { in DMPlexConstructCohesiveCells_Internal()
1381 if ((val == 1) || (val == -(shift + 1))) ++qp; in DMPlexConstructCohesiveCells_Internal()
1386 PetscCall(DMPlexSetSupportSize(sdm, splitp, qp + 1)); in DMPlexConstructCohesiveCells_Internal()
1395 for (e = 0, qn = 0, qp = 0, qf = 0; e < supportSize; ++e) { in DMPlexConstructCohesiveCells_Internal()
1401 if ((val == dim - 1) || (val == -(shift + dim - 1))) ++qp; in DMPlexConstructCohesiveCells_Internal()
1406 PetscCall(DMPlexSetSupportSize(sdm, splitp, qp + 1)); in DMPlexConstructCohesiveCells_Internal()
1480 PetscInt coneSize, supportSize, q, qf, qn, qp, v, e, s; in DMPlexConstructCohesiveCells_Internal() local
1618 for (e = 0, qp = 0; e < supportSize; ++e) { in DMPlexConstructCohesiveCells_Internal()
1625 supportNew[qp++] = edge + pMaxNew[dep + 1]; in DMPlexConstructCohesiveCells_Internal()
[all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 3379 PetscInt dim, qdim, d, Na, o, Nq, q, qp; in PetscQuadratureComputePermutations() local
3401 for (qp = 0; qp < Nq; ++qp) { in PetscQuadratureComputePermutations()
3404 for (d = 0; d < dim; ++d) diff += PetscAbsReal(txq[d] - xq[qp * dim + d]); in PetscQuadratureComputePermutations()
3407 …PetscCheck(qp < Nq, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Transformed quad point %" PetscInt_FMT " doe… in PetscQuadratureComputePermutations()
3408 idx[q] = qp; in PetscQuadratureComputePermutations()
|
| /petsc/doc/ |
| H A D | petsc.bib | 31563 @Article{ kanzow.qi:qp-free,
|