Home
last modified time | relevance | path

Searched refs:rho1 (Results 1 – 7 of 7) sorted by relevance

/petsc/src/ksp/ksp/impls/symmlq/
H A Dsymmlq.c18 PetscScalar c = 1.0, cold = 1.0, s = 0.0, sold = 0.0, coold, soold, rho0, rho1, rho2, rho3; in KSPSolve_SYMMLQ() local
150 rho1 = PetscSqrtScalar(rho0 * rho0 + beta * beta); /* gamma */ in KSPSolve_SYMMLQ()
155 c = rho0 / rho1; in KSPSolve_SYMMLQ()
156 s = beta / rho1; in KSPSolve_SYMMLQ()
158 if (ksp->its == 1) ceta = beta1 / rho1; in KSPSolve_SYMMLQ()
159 else ceta = -(rho2 * ceta_old + rho3 * ceta_oold) / rho1; in KSPSolve_SYMMLQ()
/petsc/src/ksp/ksp/impls/bcgs/qmrcgs/
H A Dqmrcgs.c18 PetscScalar eta, rho1, rho2, alpha, eta2, omega, beta, cf, cf1, uu; in KSPSolve_QMRCGS() local
85 PetscCall(VecDot(RP, RP, &rho1)); in KSPSolve_QMRCGS()
100 if (rho1 == 0) { in KSPSolve_QMRCGS()
106 alpha = rho1 / rho2; in KSPSolve_QMRCGS()
191 beta = (alpha * rho2) / (omega * rho1); in KSPSolve_QMRCGS()
193 rho1 = rho2; in KSPSolve_QMRCGS()
/petsc/src/ksp/ksp/impls/bcgsl/
H A Dbcgsl.c15 PetscScalar rho0, rho1; in KSPSolve_BCGSL() local
101 PetscCall(VecDot(VVR[j], VRT, &rho1)); in KSPSolve_BCGSL()
102 KSPCheckDot(ksp, rho1); in KSPSolve_BCGSL()
103 if (rho1 == 0.0) { in KSPSolve_BCGSL()
107 beta = alpha * (rho1 / rho0); in KSPSolve_BCGSL()
108 rho0 = rho1; in KSPSolve_BCGSL()
122 alpha = rho1 / sigma; in KSPSolve_BCGSL()
/petsc/src/ksp/ksp/impls/minres/
H A Dminres.c481 PetscScalar rho0, rho1, rho2, rho3, dp = 0.0; in KSPSolve_MINRES_OLD() local
596 rho1 = PetscSqrtScalar(rho0 * rho0 + beta * beta); in KSPSolve_MINRES_OLD()
609 c = rho0 / rho1; in KSPSolve_MINRES_OLD()
610 s = beta / rho1; in KSPSolve_MINRES_OLD()
618 PetscCall(VecAXPBY(W, 1.0 / rho1, 0.0, U)); in KSPSolve_MINRES_OLD()
621 PetscScalar alphas[] = {-rho2 / rho1, -rho3 / rho1}; in KSPSolve_MINRES_OLD()
/petsc/src/ts/adapt/impls/dsp/
H A Dadaptdsp.c131 PetscReal rho1 = dsp->hratio[0]; in TSAdaptChoose_DSP() local
140 rho0 *= PetscPowReal(rho1, -a2); in TSAdaptChoose_DSP()
/petsc/src/ts/tutorials/
H A Dex11.h582 …lar *uxm, PetscScalar *px, PetscScalar *utx, PetscScalar *ubx, PetscScalar *gam, PetscScalar *rho1) in riemannsolver() argument
602 *rho1 = *rho1r; in riemannsolver()
606 *rho1 = *rho1l; in riemannsolver()
666 *rho1 = *rho1r; in riemannsolver()
670 *rho1 = *rho1l; in riemannsolver()
/petsc/src/dm/impls/stag/tutorials/
H A Dex4.c30 PetscScalar eta1, eta2, rho1, rho2, gy, eta_characteristic; member
329 return PetscRealPart(x) < (ctx->xmax - ctx->xmin) / 2.0 ? ctx->rho1 : ctx->rho2; in GetRho_layers()
345 return (xx * xx > 0.15 * 0.15 || yy * yy > 0.15 * 0.15) ? ctx->rho1 : ctx->rho2; in GetRho_sinker_box2()
364 …f_width || PetscAbsReal(yy) > half_width || PetscAbsReal(zz) > half_width) ? ctx->rho1 : ctx->rho2; in GetRho_sinker_box3()
384 return (xx * xx + yy * yy + zz * zz > half_width * half_width) ? ctx->rho1 : ctx->rho2; in GetRho_sinker_sphere3()
412 return ctx->rho1 + ctx->rho2 * PetscExpScalar(-20.0 * (xx * xx + yy * yy + zz * zz)); in GetRho_blob3()
469 ctx->rho1 = 0.0; in CtxCreateAndSetFromOptions()
478 ctx->rho1 = 3200; in CtxCreateAndSetFromOptions()