Home
last modified time | relevance | path

Searched refs:dp (Results 1 – 25 of 43) sorted by relevance

12

/petsc/src/ksp/ksp/impls/cr/
H A Dcr.c13 PetscReal dp; in KSPSolve_CR() local
46 PetscCall(VecNormBegin(RT, NORM_2, &dp)); /* dp <- RT'*RT */ in KSPSolve_CR()
48 PetscCall(VecNormEnd(RT, NORM_2, &dp)); /* dp <- RT'*RT */ in KSPSolve_CR()
49 KSPCheckNorm(ksp, dp); in KSPSolve_CR()
51 …dp = 0.0; /* meaningless value that is passed to monitor and convergen… in KSPSolve_CR()
54 PetscCall(VecNormBegin(R, NORM_2, &dp)); /* dp <- R'*R */ in KSPSolve_CR()
56 PetscCall(VecNormEnd(R, NORM_2, &dp)); /* dp <- RT'*RT */ in KSPSolve_CR()
57 KSPCheckNorm(ksp, dp); in KSPSolve_CR()
60 dp = PetscSqrtReal(PetscAbsScalar(btop)); /* dp = sqrt(R,AR) */ in KSPSolve_CR()
69 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_CR()
[all …]
/petsc/src/ksp/ksp/impls/cg/groppcg/
H A Dgroppcg.c27 PetscReal dp = 0.0; in KSPSolve_GROPPCG() local
65 PetscCall(VecNorm(z, NORM_2, &dp)); /* dp <- z'*z = e'*A'*B'*B*A'*e' */ in KSPSolve_GROPPCG()
69 PetscCall(VecNorm(r, NORM_2, &dp)); /* dp <- r'*r = e'*A'*A*e */ in KSPSolve_GROPPCG()
73 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- r'*z = r'*B*r = e'*A'*B*A*e */ in KSPSolve_GROPPCG()
76 dp = 0.0; in KSPSolve_GROPPCG()
81 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_GROPPCG()
82 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_GROPPCG()
83 ksp->rnorm = dp; in KSPSolve_GROPPCG()
84 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); /* test for convergence */ in KSPSolve_GROPPCG()
105 PetscCall(VecNormBegin(r, NORM_2, &dp)); in KSPSolve_GROPPCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecg/
H A Dpipecg.c24 PetscReal dp = 0.0; in KSPSolve_PIPECG() local
59 PetscCall(VecNormBegin(U, NORM_2, &dp)); /* dp <- u'*u = e'*A'*B'*B*A'*e' */ in KSPSolve_PIPECG()
62 PetscCall(VecNormEnd(U, NORM_2, &dp)); in KSPSolve_PIPECG()
65 PetscCall(VecNormBegin(R, NORM_2, &dp)); /* dp <- r'*r = e'*A'*A*e */ in KSPSolve_PIPECG()
68 PetscCall(VecNormEnd(R, NORM_2, &dp)); in KSPSolve_PIPECG()
76 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- r'*u = r'*B*r = e'*A'*B*A*e */ in KSPSolve_PIPECG()
80 dp = 0.0; in KSPSolve_PIPECG()
85 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_PIPECG()
86 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_PIPECG()
87 ksp->rnorm = dp; in KSPSolve_PIPECG()
[all …]
/petsc/src/ksp/ksp/impls/cgs/
H A Dcgs.c20 PetscReal dp = 0.0; in KSPSolve_CGS() local
45 PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_CGS()
46 KSPCheckNorm(ksp, dp); in KSPSolve_CGS()
47 if (ksp->normtype == KSP_NORM_NATURAL) dp *= dp; in KSPSolve_CGS()
48 } else dp = 0.0; in KSPSolve_CGS()
52 ksp->rnorm = dp; in KSPSolve_CGS()
54 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_CGS()
55 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_CGS()
56 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_CGS()
99 dp = PetscAbsScalar(rho); in KSPSolve_CGS()
[all …]
/petsc/src/ksp/ksp/impls/bicg/
H A Dbicg.c15 PetscReal dp; in KSPSolve_BiCG() local
44 PetscCall(VecNorm(Zr, NORM_2, &dp)); /* dp <- z'*z */ in KSPSolve_BiCG()
46 PetscCall(VecNorm(Rr, NORM_2, &dp)); /* dp <- r'*r */ in KSPSolve_BiCG()
47 } else dp = 0.0; in KSPSolve_BiCG()
49 KSPCheckNorm(ksp, dp); in KSPSolve_BiCG()
50 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_BiCG()
53 ksp->rnorm = dp; in KSPSolve_BiCG()
55 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_BiCG()
56 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_BiCG()
90 PetscCall(VecNorm(Zr, NORM_2, &dp)); /* dp <- z'*z */ in KSPSolve_BiCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/
H A Dcg.c123 PetscReal dp = 0.0; in KSPSolve_CG() local
171 PetscCall(VecNorm(Z, NORM_2, &dp)); /* dp <- z'*z = e'*A'*B'*B*A*e */ in KSPSolve_CG()
172 KSPCheckNorm(ksp, dp); in KSPSolve_CG()
175 PetscCall(VecNorm(R, NORM_2, &dp)); /* dp <- r'*r = e'*A'*A*e */ in KSPSolve_CG()
176 KSPCheckNorm(ksp, dp); in KSPSolve_CG()
182 dp = PetscSqrtReal(PetscAbsScalar(beta)); /* dp <- r'*z = r'*B*r = e'*A'*B*A*e */ in KSPSolve_CG()
185 dp = 0.0; in KSPSolve_CG()
200 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_CG()
201 PetscCall(KSPMonitor(ksp, ksp->its, dp)); in KSPSolve_CG()
202 ksp->rnorm = dp; in KSPSolve_CG()
[all …]
/petsc/src/ksp/ksp/impls/cr/pipecr/
H A Dpipecr.c24 PetscReal dp = 0.0; in KSPSolve_PIPECR() local
57 PetscCall(VecNormBegin(U, NORM_2, &dp)); /* dp <- u'*u = e'*A'*B'*B*A'*e' */ in KSPSolve_PIPECR()
60 PetscCall(VecNormEnd(U, NORM_2, &dp)); in KSPSolve_PIPECR()
64 dp = 0.0; in KSPSolve_PIPECR()
69 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_PIPECR()
70 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_PIPECR()
71 ksp->rnorm = dp; in KSPSolve_PIPECR()
72 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); /* test for convergence */ in KSPSolve_PIPECR()
79 if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) PetscCall(VecNormBegin(U, NORM_2, &dp)); in KSPSolve_PIPECR()
86 if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) PetscCall(VecNormEnd(U, NORM_2, &dp)); in KSPSolve_PIPECR()
[all …]
/petsc/src/ksp/ksp/impls/cg/cgne/
H A Dcgne.c45 PetscReal dp = 0.0; in KSPSolve_CGNE() local
92 PetscCall(VecNorm(Z, NORM_2, &dp)); /* dp <- z'*z */ in KSPSolve_CGNE()
94 PetscCall(VecNorm(R, NORM_2, &dp)); /* dp <- r'*r */ in KSPSolve_CGNE()
98 dp = PetscSqrtReal(PetscAbsScalar(beta)); in KSPSolve_CGNE()
99 } else dp = 0.0; in KSPSolve_CGNE()
100 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_CGNE()
101 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_CGNE()
102 ksp->rnorm = dp; in KSPSolve_CGNE()
103 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); /* test for convergence */ in KSPSolve_CGNE()
149 PetscCall(VecNorm(Z, NORM_2, &dp)); /* dp <- z'*z */ in KSPSolve_CGNE()
[all …]
/petsc/src/ksp/ksp/impls/symmlq/
H A Dsymmlq.c19 PetscScalar dp = 0.0; in KSPSolve_SYMMLQ() local
57 PetscCall(VecDot(R, Z, &dp)); /* dp = r'*z; */ in KSPSolve_SYMMLQ()
58 KSPCheckDot(ksp, dp); in KSPSolve_SYMMLQ()
59 if (PetscAbsScalar(dp) < symmlq->haptol) { in KSPSolve_SYMMLQ()
60 …(ksp, "Detected happy breakdown %g tolerance %g\n", (double)PetscAbsScalar(dp), (double)symmlq->ha… in KSPSolve_SYMMLQ()
67 if (dp < 0.0) { in KSPSolve_SYMMLQ()
72 dp = PetscSqrtScalar(dp); in KSPSolve_SYMMLQ()
73 beta = dp; /* beta <- sqrt(r'*z) */ in KSPSolve_SYMMLQ()
129 PetscCall(VecDot(R, Z, &dp)); /* dp <- r'*z; */ in KSPSolve_SYMMLQ()
130 KSPCheckDot(ksp, dp); in KSPSolve_SYMMLQ()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecgrr/
H A Dpipecgrr.c24 …PetscReal dp = 0.0, nsi = 0.0, sqn = 0.0, Anorm = 0.0, rnp = 0.0, pnp = 0.0, snp = 0.0, unp = 0.… in KSPSolve_PIPECGRR() local
60 PetscCall(VecNormBegin(U, NORM_2, &dp)); /* dp <- u'*u = e'*A'*B'*B*A'*e' */ in KSPSolve_PIPECGRR()
64 PetscCall(VecNormEnd(U, NORM_2, &dp)); in KSPSolve_PIPECGRR()
68 PetscCall(VecNormBegin(R, NORM_2, &dp)); /* dp <- r'*r = e'*A'*A*e */ in KSPSolve_PIPECGRR()
72 PetscCall(VecNormEnd(R, NORM_2, &dp)); in KSPSolve_PIPECGRR()
83 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- r'*u = r'*B*r = e'*A'*B*A*e */ in KSPSolve_PIPECGRR()
87 dp = 0.0; in KSPSolve_PIPECGRR()
92 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_PIPECGRR()
93 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_PIPECGRR()
94 ksp->rnorm = dp; in KSPSolve_PIPECGRR()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/qmrcgs/
H A Dqmrcgs.c20 PetscReal dp = 0.0, final, tau, tau2, theta, theta2, c, F, NV, vv; in KSPSolve_QMRCGS() local
64 if (ksp->normtype != KSP_NORM_NONE) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_QMRCGS()
67 ksp->rnorm = dp; in KSPSolve_QMRCGS()
69 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_QMRCGS()
70 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_QMRCGS()
71 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_QMRCGS()
79 if (dp == 0.0) { in KSPSolve_QMRCGS()
82 tau = dp; in KSPSolve_QMRCGS()
139 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_QMRCGS()
162 if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i + 2) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_QMRCGS()
[all …]
/petsc/src/tao/linesearch/impls/morethuente/
H A Dmorethuente.c10 …l *dx, PetscReal *sty, PetscReal *fy, PetscReal *dy, PetscReal *stp, PetscReal *fp, PetscReal *dp);
363 …al *dx, PetscReal *sty, PetscReal *fy, PetscReal *dy, PetscReal *stp, PetscReal *fp, PetscReal *dp) in Tao_mcstep() argument
377 sgnd = *dp * (*dx / PetscAbsReal(*dx)); in Tao_mcstep()
387 theta = 3 * (*fx - *fp) / (*stp - *stx) + *dx + *dp; in Tao_mcstep()
389 s = PetscMax(s, PetscAbsReal(*dp)); in Tao_mcstep()
390 gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s)); in Tao_mcstep()
394 q = ((gamma1 - *dx) + gamma1) + *dp; in Tao_mcstep()
410 theta = 3 * (*fx - *fp) / (*stp - *stx) + *dx + *dp; in Tao_mcstep()
412 s = PetscMax(s, PetscAbsReal(*dp)); in Tao_mcstep()
413 gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s)); in Tao_mcstep()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/
H A Dbcgs.c23 PetscReal dp = 0.0, d2; in KSPSolve_BCGS() local
48 PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_BCGS()
49 KSPCheckNorm(ksp, dp); in KSPSolve_BCGS()
53 ksp->rnorm = dp; in KSPSolve_BCGS()
55 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_BCGS()
56 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_BCGS()
57 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_BCGS()
106 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_BCGS()
114 PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_BCGS()
115 KSPCheckNorm(ksp, dp); in KSPSolve_BCGS()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/fbcgs/
H A Dfbcgs.c20 PetscReal dp = 0.0, d2; in KSPSolve_FBCGS() local
56 if (ksp->normtype != KSP_NORM_NONE) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_FBCGS()
59 ksp->rnorm = dp; in KSPSolve_FBCGS()
61 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_FBCGS()
62 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_FBCGS()
63 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_FBCGS()
113 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_FBCGS()
121 if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i + 2) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_FBCGS()
128 ksp->rnorm = dp; in KSPSolve_FBCGS()
130 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_FBCGS()
[all …]
/petsc/src/ksp/ksp/impls/tfqmr/
H A Dtfqmr.c15 PetscReal dp, dpold, w, dpest, tau, psi, cm; in KSPSolve_TFQMR() local
36 PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_TFQMR()
37 KSPCheckNorm(ksp, dp); in KSPSolve_TFQMR()
39 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = dp; in KSPSolve_TFQMR()
54 tau = dp; in KSPSolve_TFQMR()
55 dpold = dp; in KSPSolve_TFQMR()
75 PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_TFQMR()
76 KSPCheckNorm(ksp, dp); in KSPSolve_TFQMR()
78 if (!m) w = PetscSqrtReal(dp * dpold); in KSPSolve_TFQMR()
79 else w = dp; in KSPSolve_TFQMR()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/pipebcgs/
H A Dpipebcgs.c20 PetscReal dp = 0.0; in KSPSolve_PIPEBCGS() local
61 if (ksp->normtype != KSP_NORM_NONE) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_PIPEBCGS()
62 else dp = 0.0; in KSPSolve_PIPEBCGS()
65 ksp->rnorm = dp; in KSPSolve_PIPEBCGS()
67 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_PIPEBCGS()
68 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_PIPEBCGS()
69 PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_PIPEBCGS()
131 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_PIPEBCGS()
144 … KSP_NORM_NONE && ksp->chknorm < i + 2) PetscCall(VecNormBegin(R, NORM_2, &dp)); /* dp <- norm(r) … in KSPSolve_PIPEBCGS()
154 … if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i + 2) PetscCall(VecNormEnd(R, NORM_2, &dp)); in KSPSolve_PIPEBCGS()
[all …]
/petsc/src/ksp/ksp/impls/fcg/
H A Dfcg.c77 PetscReal dp = 0.0; in KSPSolve_FCG() local
110 PetscCall(VecNorm(Z, NORM_2, &dp)); /* dp <- dqrt(z'*z) = sqrt(e'*A'*B'*B*A*e) */ in KSPSolve_FCG()
111 KSPCheckNorm(ksp, dp); in KSPSolve_FCG()
114 PetscCall(VecNorm(R, NORM_2, &dp)); /* dp <- sqrt(r'*r) = sqrt(e'*A'*A*e) */ in KSPSolve_FCG()
115 KSPCheckNorm(ksp, dp); in KSPSolve_FCG()
121 dp = PetscSqrtReal(PetscAbsScalar(s)); /* dp <- sqrt(r'*z) = sqrt(e'*A'*B*A*e) */ in KSPSolve_FCG()
124 dp = 0.0; in KSPSolve_FCG()
131 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_FCG()
132 PetscCall(KSPMonitor(ksp, 0, dp)); in KSPSolve_FCG()
133 ksp->rnorm = dp; in KSPSolve_FCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipeprcg/
H A Dpipeprcg.c43 PetscReal dp = 0.0; in KSPSolve_PIPEPRCG() local
105 PetscCall(VecNormBegin(RT, NORM_2, &dp)); in KSPSolve_PIPEPRCG()
107 PetscCall(VecNormEnd(RT, NORM_2, &dp)); in KSPSolve_PIPEPRCG()
110 PetscCall(VecNormBegin(R, NORM_2, &dp)); in KSPSolve_PIPEPRCG()
112 PetscCall(VecNormEnd(R, NORM_2, &dp)); in KSPSolve_PIPEPRCG()
115 dp = PetscSqrtReal(PetscAbsScalar(nu)); in KSPSolve_PIPEPRCG()
118 dp = 0.0; in KSPSolve_PIPEPRCG()
124 ksp->rnorm = dp; in KSPSolve_PIPEPRCG()
125 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_PIPEPRCG()
126 PetscCall(KSPMonitor(ksp, i, dp)); in KSPSolve_PIPEPRCG()
[all …]
/petsc/src/ksp/ksp/impls/fcg/pipefcg/
H A Dpipefcg.c85 PetscReal dp = 0.0, delta, *eta, *etas; in KSPSolve_PIPEFCG_cycle() local
158 PetscCall(VecNorm(Z, NORM_2, &dp)); /* dp <- sqrt(z'*z) = sqrt(e'*A'*B'*B*A*e) */ in KSPSolve_PIPEFCG_cycle()
161 PetscCall(VecNorm(R, NORM_2, &dp)); /* dp <- sqrt(r'*r) = sqrt(e'*A'*A*e) */ in KSPSolve_PIPEFCG_cycle()
164 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- sqrt(r'*z) = sqrt(e'*A'*B*A*e) */ in KSPSolve_PIPEFCG_cycle()
167 dp = 0.0; in KSPSolve_PIPEFCG_cycle()
174 ksp->rnorm = dp; in KSPSolve_PIPEFCG_cycle()
175 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_PIPEFCG_cycle()
176 PetscCall(KSPMonitor(ksp, ksp->its, dp)); in KSPSolve_PIPEFCG_cycle()
177 PetscCall((*ksp->converged)(ksp, ksp->its, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_PIPEFCG_cycle()
264 PetscReal dp = 0.0; in KSPSolve_PIPEFCG() local
[all …]
/petsc/src/ksp/ksp/impls/minres/
H A Dminres.c481 PetscScalar rho0, rho1, rho2, rho3, dp = 0.0; in KSPSolve_MINRES_OLD() local
526 PetscCall(VecDot(R, Z, &dp)); in KSPSolve_MINRES_OLD()
527 KSPCheckDot(ksp, dp); in KSPSolve_MINRES_OLD()
529 if (PetscRealPart(dp) < minres->haptol && np > minres->haptol) { in KSPSolve_MINRES_OLD()
530 …LED, "Detected indefinite operator %g tolerance %g", (double)PetscRealPart(dp), (double)minres->ha… in KSPSolve_MINRES_OLD()
531 …p, "Detected indefinite operator %g tolerance %g\n", (double)PetscRealPart(dp), (double)minres->ha… in KSPSolve_MINRES_OLD()
543 dp = PetscAbsScalar(dp); in KSPSolve_MINRES_OLD()
544 dp = PetscSqrtScalar(dp); in KSPSolve_MINRES_OLD()
545 beta = dp; /* beta <- sqrt(r'*z) */ in KSPSolve_MINRES_OLD()
583 PetscCall(VecDot(R, Z, &dp)); in KSPSolve_MINRES_OLD()
[all …]
/petsc/src/snes/tutorials/
H A Dex48.c801 const PetscReal pp = phi[l], *dp = dphi[l]; in THIFunctionLocal() local
802 …fn[l]->u += dp[0] * jw * eta * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp… in THIFunctionLocal()
803 …fn[l]->v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp… in THIFunctionLocal()
1004 const PetscReal pp = phi[l], *dp = dphi[l]; in THIJacobianLocal_2D() local
1011 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + pp * jw … in THIJacobianLocal_2D()
1012 Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * eta * 2. * dpl[1] + dp[1] * jw * eta * dpl[0]; in THIJacobianLocal_2D()
1013 Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * eta * 2. * dpl[0] + dp[0] * jw * eta * dpl[1]; in THIJacobianLocal_2D()
1014 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + pp * jw … in THIJacobianLocal_2D()
1016 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d… in THIJacobianLocal_2D()
1017 …Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d… in THIJacobianLocal_2D()
[all …]
/petsc/src/ksp/pc/impls/bjacobi/bjkokkos/
H A Dbjkokkos.kokkos.cxx93 PetscReal dp, dpold, w, dpest, tau, psi, cm, r0; in BJSolve_TFQMR() local
166 r0 = dp = PetscSqrtReal(PetscRealPart(dpi)); in BJSolve_TFQMR()
169 …single(Kokkos::PerTeam(team), [=]() { printf("%3d KSP Residual norm %14.12e\n", 0, (double)dp); }); in BJSolve_TFQMR()
171 if (dp < atol) { in BJSolve_TFQMR()
188 tau = dp; in BJSolve_TFQMR()
189 dpold = dp; in BJSolve_TFQMR()
229 dp = PetscSqrtReal(PetscRealPart(dpi)); in BJSolve_TFQMR()
231 if (!m) w = PetscSqrtReal(dp * dpold); in BJSolve_TFQMR()
232 else w = dp; in BJSolve_TFQMR()
305 dpold = dp; in BJSolve_TFQMR()
[all …]
/petsc/src/vec/vec/impls/mpi/
H A Dpvecimpl.h137 static inline PetscErrorCode VecDotNorm2_MPI_Default(Vec s, Vec t, PetscScalar *dp, PetscScalar *nm… in VecDotNorm2_MPI_Default() argument
140 PetscCall(VecDotNorm2_SeqFn(s, t, dp, nm)); in VecDotNorm2_MPI_Default()
142 PetscScalar sum[] = {*dp, *nm}; in VecDotNorm2_MPI_Default()
145 *dp = sum[0]; in VecDotNorm2_MPI_Default()
/petsc/src/ksp/ksp/impls/cg/pipelcg/
H A Dpipelcg.c143 PetscReal dp = 0.0, tmp = 0.0, beta = 0.0, invbeta2 = 0.0; in KSPSolve_InnerLoop_PIPELCG() local
314 dp = beta; in KSPSolve_InnerLoop_PIPELCG()
324 dp = PetscAbsScalar(zeta); in KSPSolve_InnerLoop_PIPELCG()
326 ksp->rnorm = dp; in KSPSolve_InnerLoop_PIPELCG()
327 PetscCall(KSPLogResidualHistory(ksp, dp)); in KSPSolve_InnerLoop_PIPELCG()
328 PetscCall(KSPMonitor(ksp, ksp->its, dp)); in KSPSolve_InnerLoop_PIPELCG()
329 PetscCall((*ksp->converged)(ksp, ksp->its, dp, &ksp->reason, ksp->cnvP)); in KSPSolve_InnerLoop_PIPELCG()
/petsc/src/ts/tutorials/
H A Dex14.c285 …d4(const PetscReal dphi[][4][2], PetscReal hx, PetscReal hy, const PrmNode pn[4], PrmNode dp[4][2]) in QuadComputeGrad4()
289 PetscScalar (*restrict dpg)[2][PRMNODE_SIZE] = (PetscScalar (*)[2][PRMNODE_SIZE])dp; in QuadComputeGrad4()
803 const PetscReal pp = phi[l], *dp = dphi[l]; in THIFunctionLocal_3D() local
804 …fn[l]->u += dp[0] * jw * eta * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp… in THIFunctionLocal_3D()
805 …fn[l]->v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp… in THIFunctionLocal_3D()
1136 const PetscReal pp = phi[l], *restrict dp = dphi[l]; in THIJacobianLocal_Momentum() local
1143 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + dp[2] * … in THIJacobianLocal_Momentum()
1144 … Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * eta * 2. * dpl[1] + dp[1] * jw * eta * dpl[0]; in THIJacobianLocal_Momentum()
1145 … Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * eta * 2. * dpl[0] + dp[0] * jw * eta * dpl[1]; in THIJacobianLocal_Momentum()
1146 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + dp[2] * … in THIJacobianLocal_Momentum()
[all …]

12