| /petsc/src/sys/classes/draw/utils/ |
| H A D | zoom.c | 20 PetscReal dpause, xc, yc, scale = 1.0, w, h, xr, xl, yr, yl, xmin, xmax, ymin, ymax; in PetscDrawZoom() local 53 scale = 0.5; in PetscDrawZoom() 56 scale = 2.0; in PetscDrawZoom() 59 scale = 8 / 10.; in PetscDrawZoom() 62 scale = 10 / 8.; in PetscDrawZoom() 65 scale = 1.0; in PetscDrawZoom() 67 xl = scale * (xl + w - xc) + xc - w * scale; in PetscDrawZoom() 68 xr = scale * (xr - w - xc) + xc + w * scale; in PetscDrawZoom() 69 yl = scale * (yl + h - yc) + yc - h * scale; in PetscDrawZoom() 70 yr = scale * (yr - h - yc) + yc + h * scale; in PetscDrawZoom() [all …]
|
| /petsc/src/ksp/ksp/impls/rich/ |
| H A D | richscale.c | 19 PetscErrorCode KSPRichardsonSetScale(KSP ksp, PetscReal scale) in KSPRichardsonSetScale() argument 23 PetscValidLogicalCollectiveReal(ksp, scale, 2); in KSPRichardsonSetScale() 24 PetscTryMethod(ksp, "KSPRichardsonSetScale_C", (KSP, PetscReal), (ksp, scale)); in KSPRichardsonSetScale() 51 PetscErrorCode KSPRichardsonSetSelfScale(KSP ksp, PetscBool scale) in KSPRichardsonSetSelfScale() argument 55 PetscValidLogicalCollectiveBool(ksp, scale, 2); in KSPRichardsonSetSelfScale() 56 PetscTryMethod(ksp, "KSPRichardsonSetSelfScale_C", (KSP, PetscBool), (ksp, scale)); in KSPRichardsonSetSelfScale()
|
| H A D | rich.c | 23 PetscScalar scale, rdot; in KSPSolve_Richardson() local 60 …if (exists && maxit > 0 && richardsonP->scale == 1.0 && (ksp->converged == KSPConvergedDefault || … in KSPSolve_Richardson() 92 scale = rdot / abr; in KSPSolve_Richardson() 93 PetscCall(PetscInfo(ksp, "Self-scale factor %g\n", (double)PetscRealPart(scale))); in KSPSolve_Richardson() 94 PetscCall(VecAXPY(x, scale, z)); /* x <- x + scale z */ in KSPSolve_Richardson() 95 PetscCall(VecAXPY(r, -scale, w)); /* r <- r - scale*Az */ in KSPSolve_Richardson() 96 PetscCall(VecAXPY(z, -scale, y)); /* z <- z - scale*y */ in KSPSolve_Richardson() 114 PetscCall(VecAXPY(x, richardsonP->scale, z)); /* x <- x + scale z */ in KSPSolve_Richardson() 158 … PetscCall(PetscViewerASCIIPrintf(viewer, " damping factor=%g\n", (double)richardsonP->scale)); in KSPView_Richardson() 172 …Real("-ksp_richardson_scale", "damping factor", "KSPRichardsonSetScale", rich->scale, &tmp, &flg)); in KSPSetFromOptions_Richardson() [all …]
|
| /petsc/src/mat/tests/output/ |
| H A D | ex140_1.out | 18 Pass: scale 46 Pass: scale 75 Pass: scale 103 Pass: scale 132 Pass: scale 160 Pass: scale 189 Pass: scale 217 Pass: scale 246 Pass: scale 274 Pass: scale [all …]
|
| /petsc/src/sys/classes/random/impls/random123/ |
| H A D | random123.c | 43 PetscReal scale = 1. / ((PetscReal)UINT64_MAX + 1.); in PetscRandom123Step() local 44 PetscReal shift = .5 * scale; in PetscRandom123Step() 48 ret = r123->result.v[mod] * scale + shift; in PetscRandom123Step() 104 PetscReal scale = 1. / ((PetscReal)UINT64_MAX + 1.); in PetscRandomGetValuesReal_Random123() local 105 PetscReal shift = .5 * scale; in PetscRandomGetValuesReal_Random123() 118 scale *= PetscRealPart(r->width); in PetscRandomGetValuesReal_Random123() 124 vals[0] = r123_copy.result.v[0] * scale + shift; in PetscRandomGetValuesReal_Random123() 125 vals[1] = r123_copy.result.v[1] * scale + shift; in PetscRandomGetValuesReal_Random123() 126 vals[2] = r123_copy.result.v[2] * scale + shift; in PetscRandomGetValuesReal_Random123() 127 vals[3] = r123_copy.result.v[3] * scale + shift; in PetscRandomGetValuesReal_Random123()
|
| /petsc/src/dm/tutorials/ |
| H A D | ex26.c | 14 PetscReal *zeros, *weights, scale, h, sigma = 1 / sqrt(2), g = 0, mu = 0; in main() local 33 PetscCall(PetscDTFactorial(n, &scale)); in main() 34 scale = exp2(n - 1) * scale * PetscSqrtReal(PETSC_PI) / (n * n); in main() 38 weights[s] = scale / (h * h); in main()
|
| /petsc/src/ksp/pc/impls/lsc/ |
| H A D | lsc.c | 6 Vec Avec0, Avec1, Svec0, scale; member 27 if (lsc->scalediag) PetscCall(VecDuplicate(lsc->Avec0, &lsc->scale)); in PCLSCAllocate_Private() 67 if (lsc->scale) { in PCSetUp_LSC() 69 PetscCall(MatGetDiagonal(Qscale, lsc->scale)); in PCSetUp_LSC() 70 PetscCall(VecReciprocal(lsc->scale)); in PCSetUp_LSC() 75 if (lsc->scale) { in PCSetUp_LSC() 78 PetscCall(MatDiagonalScale(CAdiaginv, NULL, lsc->scale)); in PCSetUp_LSC() 127 if (lsc->scale) PetscCall(VecPointwiseMult(lsc->Avec0, lsc->Avec0, lsc->scale)); in PCApply_LSC() 129 if (lsc->scale) PetscCall(VecPointwiseMult(lsc->Avec1, lsc->Avec1, lsc->scale)); in PCApply_LSC() 148 if (lsc->scale) PetscCall(VecDestroy(&lsc->scale)); in PCReset_LSC()
|
| /petsc/src/snes/interface/noise/ |
| H A D | snesdnest.c | 26 static double scale; in SNESNoise_dnest_() local 158 scale = 0.; in SNESNoise_dnest_() 166 d__2 = scale; in SNESNoise_dnest_() 168 scale = PetscMax(d__2, d__3); in SNESNoise_dnest_() 172 if (scale == 0.) stdv = 0.; in SNESNoise_dnest_() 178 d__1 = fval[i__] / scale; in SNESNoise_dnest_() 181 stdv = scale * PetscSqrtScalar(stdv / (*nf - j)); in SNESNoise_dnest_()
|
| /petsc/src/ksp/pc/impls/jacobi/ |
| H A D | jacobi.c | 63 PetscReal scale; /* for scaling rowl1 off-diagonals */ member 113 j->scale = flg; in PCJacobiSetRowl1Scale_Jacobi() 122 *flg = j->scale; in PCJacobiGetRowl1Scale_Jacobi() 238 if (jac->scale != 1.0) { in PCSetUp_Jacobi() 243 PetscCall(VecScale(diag, jac->scale)); // scale off-diag in PCSetUp_Jacobi() 442 …of off-diagonal elements for rowl1", "PCJacobiSetRowl1Scale", jac->scale, &jac->scale, NULL, 0.0, … in PCSetFromOptions_Jacobi() 458 PetscReal scale; in PCView_Jacobi() local 463 PetscCall(PCJacobiGetRowl1Scale(pc, &scale)); in PCView_Jacobi() 465 …te value of entries" : "", !fixdiag ? ", not checking null diagonal entries" : "", (double)scale)); in PCView_Jacobi() 530 jac->scale = 1.0; in PCCreate_Jacobi() [all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/tests/ |
| H A D | ex2.c | 66 PetscReal err, scale; in main() local 86 PetscCall(VecNorm(v_sym, NORM_2, &scale)); in main() 88 …cInfo(diag, "Diagonal Broyden error %g, relative error %g\n", (double)err, (double)(err / scale))); in main() 89 …err <= PETSC_SMALL * scale, comm, PETSC_ERR_PLIB, "Diagonal Broyden error %g, relative error %g", … in main()
|
| H A D | ex1.c | 29 PetscReal scale; in HermitianTransposeTest() local 45 scale = PetscMax(x_norm * Bhf_norm, Bx_norm * f_norm); in HermitianTransposeTest() 46 …ject)B, "Hermitian transpose error %g, relative error %g \n", (double)err, (double)(err / scale))); in HermitianTransposeTest() 47 …PetscCheck(err <= PETSC_SMALL * scale, PetscObjectComm((PetscObject)B), PETSC_ERR_PLIB, "Hermitian… in HermitianTransposeTest() 59 PetscReal scale; in InverseTest() local 71 scale = PetscMax(x_norm, Bx_norm); in InverseTest() 72 …Info((PetscObject)B, "Inverse error %g, relative error %g\n", (double)err, (double)(err / scale))); in InverseTest() 73 …PetscCheck(err <= 100.0 * PETSC_SMALL * scale, PetscObjectComm((PetscObject)B), PETSC_ERR_PLIB, "I… in InverseTest() 86 PetscReal scale; in IsHermitianTest() local 103 scale = PetscMax(x_norm * By_norm, Bx_norm * y_norm); in IsHermitianTest() [all …]
|
| /petsc/src/dm/impls/plex/tutorials/output/ |
| H A D | ex11_tetrahedron.out | 8 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 29 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 50 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 71 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 92 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 113 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 134 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 155 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 176 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 197 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] [all …]
|
| H A D | ex11_triangle.out | 8 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 28 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 48 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 68 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 88 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 108 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont]
|
| H A D | ex11_quadrilateral.out | 8 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 30 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 52 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 74 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 96 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 118 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 140 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 162 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont]
|
| H A D | ex11_pyramid.out | 8 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 32 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 56 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 80 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 104 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 128 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 152 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 176 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont]
|
| H A D | ex11_triangular_prism.out | 8 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 34 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 60 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 86 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 112 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 138 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 164 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 190 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 216 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 242 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] [all …]
|
| H A D | ex11_tensor_triangular_prism.out | 8 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 34 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 60 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 86 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 112 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 138 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 164 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 190 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 216 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 242 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] [all …]
|
| H A D | ex11_alfeld_tetrahedron.out | 8 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 37 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 66 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 95 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 124 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 153 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 182 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 211 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 240 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] 269 \begin{tikzpicture}[scale = 1.,font=\fontsize{8}{8}\selectfont] [all …]
|
| H A D | ex11_segment.out | 8 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont] 22 \begin{tikzpicture}[scale = 0.5,font=\fontsize{8}{8}\selectfont]
|
| /petsc/src/ksp/pc/impls/mat/tests/ |
| H A D | ex1.c | 8 PetscReal err, scale; in TestVecEquality() local 15 PetscCall(VecNorm(x, NORM_INFINITY, &scale)); in TestVecEquality() 16 …PetscCheck(err <= PETSC_SMALL * scale, PetscObjectComm((PetscObject)x), PETSC_ERR_PLIB, "PC operat… in TestVecEquality() 24 PetscReal err, scale; in TestMatEquality() local 32 PetscCall(MatNorm(x, NORM_FROBENIUS, &scale)); in TestMatEquality() 33 …PetscCheck(err < PETSC_SMALL * m * n * scale, PetscObjectComm((PetscObject)x), PETSC_ERR_PLIB, "PC… in TestMatEquality()
|
| /petsc/src/mat/impls/normal/ |
| H A D | normm.c | 31 PetscScalar shift, scale; in MatCreateSubMatrices_Normal() local 36 …PetscCall(MatShellGetScalingShifts(mat, &shift, &scale, (Vec *)MAT_SHELL_NOT_ALLOWED, (Vec *)MAT_S… in MatCreateSubMatrices_Normal() 49 PetscCall(MatScale((*submat)[M], scale)); in MatCreateSubMatrices_Normal() 62 PetscScalar shift, scale; in MatPermute_Normal() local 66 …PetscCall(MatShellGetScalingShifts(A, &shift, &scale, (Vec *)MAT_SHELL_NOT_ALLOWED, (Vec *)MAT_SHE… in MatPermute_Normal() 76 PetscCall(MatScale(*B, scale)); in MatPermute_Normal() 225 PetscScalar scale, shift; in MatConvert_Normal_AIJ() local 230 …PetscCall(MatShellGetScalingShifts(A, &shift, &scale, &dshift, &left, &right, (Mat *)MAT_SHELL_NOT… in MatConvert_Normal_AIJ() 248 PetscCall(MatScale(*newmat, scale)); in MatConvert_Normal_AIJ() 276 PetscScalar *array, scale; in MatProductNumeric_Normal_Dense() local [all …]
|
| H A D | normmh.c | 14 PetscScalar shift, scale; in MatCreateSubMatrices_NormalHermitian() local 19 …PetscCall(MatShellGetScalingShifts(mat, &shift, &scale, (Vec *)MAT_SHELL_NOT_ALLOWED, (Vec *)MAT_S… in MatCreateSubMatrices_NormalHermitian() 32 PetscCall(MatScale((*submat)[M], scale)); in MatCreateSubMatrices_NormalHermitian() 45 PetscScalar shift, scale; in MatPermute_NormalHermitian() local 49 …PetscCall(MatShellGetScalingShifts(A, &shift, &scale, (Vec *)MAT_SHELL_NOT_ALLOWED, (Vec *)MAT_SHE… in MatPermute_NormalHermitian() 59 PetscCall(MatScale(*B, scale)); in MatPermute_NormalHermitian() 209 PetscScalar scale, shift; in MatConvert_NormalHermitian_AIJ() local 214 …PetscCall(MatShellGetScalingShifts(A, &shift, &scale, &dshift, &left, &right, (Mat *)MAT_SHELL_NOT… in MatConvert_NormalHermitian_AIJ() 238 PetscCall(MatScale(*newmat, scale)); in MatConvert_NormalHermitian_AIJ()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex30.c | 257 Vec scale; in main() local 261 PetscCall(VecDuplicate(x, &scale)); in main() 262 PetscCall(VecGetOwnershipRange(scale, &start, &end)); in main() 263 …for (j = start; j < end; j++) PetscCall(VecSetValue(scale, j, ((PetscReal)(j + 1)) / ((PetscReal)n… in main() 264 PetscCall(VecAssemblyBegin(scale)); in main() 265 PetscCall(VecAssemblyEnd(scale)); in main() 266 PetscCall(PCSetDiagonalScale(pc, scale)); in main() 267 PetscCall(VecDestroy(&scale)); in main()
|
| /petsc/src/ts/adapt/interface/ |
| H A D | tsadapt.c | 609 PetscErrorCode TSAdaptSetScaleSolveFailed(TSAdapt adapt, PetscReal scale) in TSAdaptSetScaleSolveFailed() argument 613 PetscValidLogicalCollectiveReal(adapt, scale, 2); in TSAdaptSetScaleSolveFailed() 614 …PetscCheck(scale > 0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Scale factor %g must be positive… in TSAdaptSetScaleSolveFailed() 615 …PetscCheck(scale <= 1, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Scale factor %g must be less th… in TSAdaptSetScaleSolveFailed() 616 adapt->scale_solve_failed = scale; in TSAdaptSetScaleSolveFailed() 635 PetscErrorCode TSAdaptGetScaleSolveFailed(TSAdapt adapt, PetscReal *scale) in TSAdaptGetScaleSolveFailed() argument 639 if (scale) PetscAssertPointer(scale, 2); in TSAdaptGetScaleSolveFailed() 640 if (scale) *scale = adapt->scale_solve_failed; in TSAdaptGetScaleSolveFailed() 743 PetscReal safety, reject_safety, clip[2], scale, hmin, hmax; in TSAdaptSetFromOptions() local 780 …is factor if solve fails", "TSAdaptSetScaleSolveFailed", adapt->scale_solve_failed, &scale, &set)); in TSAdaptSetFromOptions() [all …]
|
| /petsc/src/mat/impls/aij/seq/spqr/ |
| H A D | aijspqr.c | 19 PetscScalar *ca, scale; in MatWrapCholmod_SPQR_seqaij() local 36 …PetscCall(MatShellGetScalingShifts(B, (PetscScalar *)MAT_SHELL_NOT_ALLOWED, &scale, (Vec *)MAT_SHE… in MatWrapCholmod_SPQR_seqaij() 150 PetscCall(VecScale(X, chol->scale)); in MatSolve_SPQR() 174 PetscCall(MatScale(X, chol->scale)); in MatMatSolve_SPQR() 255 PetscScalar scale; in MatQRFactorNumeric_SPQR() local 257 …PetscCall(MatShellGetScalingShifts(A, (PetscScalar *)MAT_SHELL_NOT_ALLOWED, &scale, (Vec *)MAT_SHE… in MatQRFactorNumeric_SPQR() 260 chol->scale = 1.0 / scale; in MatQRFactorNumeric_SPQR() 264 chol->scale = 1.0; in MatQRFactorNumeric_SPQR()
|