| /petsc/src/ts/tutorials/power_grid/output/ |
| H A D | ex9opt_2.out | 4 Testing hand-coded Gradient, if (for double precision runs) ||G - Gfd||/||G|| is 6 ||Gfd|| 0.0021849, ||G|| = 0.00218604, angle cosine = (Gfd'G)/||Gfd||||G|| = 1. 7 2-norm ||G - Gfd||/||G|| = 0.000520807, ||G - Gfd|| = 1.1385e-06 8 max-norm ||G - Gfd||/||G|| = 0.000520807, ||G - Gfd|| = 1.1385e-06 11 ||Gfd|| 1., ||G|| = 1., angle cosine = (Gfd'G)/||Gfd||||G|| = 1. 12 2-norm ||G - Gfd||/||G|| = 2.22045e-16, ||G - Gfd|| = 2.22045e-16 13 max-norm ||G - Gfd||/||G|| = 2.22045e-16, ||G - Gfd|| = 2.22045e-16 15 ||Gfd|| 0.0736232, ||G|| = 0.0736226, angle cosine = (Gfd'G)/||Gfd||||G|| = 1. 16 2-norm ||G - Gfd||/||G|| = 8.2006e-06, ||G - Gfd|| = 6.03754e-07 17 max-norm ||G - Gfd||/||G|| = 8.2006e-06, ||G - Gfd|| = 6.03754e-07 [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex170.c | 44 Mat G; in CreateGraph() local 47 PetscCall(MatCreate(comm, &G)); in CreateGraph() 53 PetscCall(MatSetSizes(G, PETSC_DETERMINE, PETSC_DETERMINE, 5, 5)); in CreateGraph() 54 PetscCall(MatSetUp(G)); in CreateGraph() 55 PetscCall(MatCreateVecs(G, &D, NULL)); in CreateGraph() 57 PetscCall(MatDiagonalSet(G, D, INSERT_VALUES)); in CreateGraph() 65 PetscCall(MatSetSizes(G, PETSC_DETERMINE, PETSC_DETERMINE, 5, 5)); in CreateGraph() 66 PetscCall(MatSetFromOptions(G)); in CreateGraph() 67 PetscCall(MatSeqAIJSetPreallocation(G, 2, NULL)); in CreateGraph() 68 PetscCall(MatSetUp(G)); in CreateGraph() [all …]
|
| H A D | ex145.c | 7 Mat A, F, B, X, C, Aher, G; in main() local 157 PetscCall(MatCreate(PETSC_COMM_WORLD, &G)); in main() 158 PetscCall(MatSetSizes(G, m, n, PETSC_DECIDE, PETSC_DECIDE)); in main() 159 PetscCall(MatSetType(G, MATELEMENTAL)); in main() 160 PetscCall(MatSetFromOptions(G)); in main() 161 PetscCall(MatSetUp(G)); in main() 162 PetscCall(MatAssemblyBegin(G, MAT_FINAL_ASSEMBLY)); in main() 163 PetscCall(MatAssemblyEnd(G, MAT_FINAL_ASSEMBLY)); in main() 164 PetscCall(MatCopy(Aher, G, SAME_NONZERO_PATTERN)); in main() 167 PetscCall(MatCholeskyFactor(G, 0, 0)); in main() [all …]
|
| H A D | ex245.c | 7 Mat A, F, B, X, C, Aher, G; in main() local 138 PetscCall(MatDuplicate(Aher, MAT_COPY_VALUES, &G)); in main() 141 PetscCall(MatCholeskyFactor(G, 0, 0)); in main() 144 PetscCall(MatView(G, PETSC_VIEWER_STDOUT_WORLD)); in main() 148 PetscCall(MatSolve(G, b, x)); in main() 149 PetscCall(MatMatSolve(G, B, X)); in main() 150 PetscCall(MatDestroy(&G)); in main() 153 PetscCall(MatGetFactor(Aher, MATSOLVERSCALAPACK, MAT_FACTOR_CHOLESKY, &G)); in main() 154 PetscCall(MatCholeskyFactorSymbolic(G, Aher, 0, NULL)); in main() 155 PetscCall(MatCholeskyFactorNumeric(G, Aher, NULL)); in main() [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | therm_meth_ethanol.dat | 3 H 120186H 1 G 0300.00 5000.00 1000.00 1 7 H2 121286H 2 G 0300.00 5000.00 1000.00 1 11 O 120186O 1 G 0300.00 5000.00 1000.00 1 15 O2 121386O 2 G 0300.00 5000.00 1000.00 1 19 OH 7/13/ 0 RUCICH 1O 1 0 0G 300.000 5000.000 1710.000 01 23 H2O 20387H 2O 1 G 0300.00 5000.00 1000.00 1 27 N2 121286N 2 G 0300.00 5000.00 1000.00 1 31 HO2 L 5/89H 1O 2 G 200.0 3500.0 1000.0 1 35 H2O2 120186H 2O 2 G 0300.00 5000.00 1000.00 1 39 AR 120186AR 1 G 0300.00 5000.00 1000.00 1 [all …]
|
| H A D | ex53.c | 149 const PetscReal G = PetscRealPart(constants[0]); in f0_quadratic_linear_u() local 154 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in f0_quadratic_linear_u() 157 for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * t; in f0_quadratic_linear_u() 158 f0[dim - 1] -= 2.0 * lambda + 4.0 * G - alpha * t; in f0_quadratic_linear_u() 231 const PetscReal G = PetscRealPart(constants[0]); in f0_quadratic_trig_u() local 236 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in f0_quadratic_trig_u() 239 for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * PetscCosReal(t); in f0_quadratic_trig_u() 240 f0[dim - 1] -= 2.0 * lambda + 4.0 * G - alpha * PetscCosReal(t); in f0_quadratic_trig_u() 344 const PetscReal G = PetscRealPart(constants[0]); in f0_trig_linear_u() local 349 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in f0_trig_linear_u() [all …]
|
| /petsc/src/ksp/ksp/impls/cg/pipelcg/ |
| H A D | pipelcg.c | 6 #define G(i, j) (plcg->G[((j) * (2 * l + 1)) + (shift((i), (j)))]) macro 7 #define G_noshift(i, j) (plcg->G[((j) * (2 * l + 1)) + (i)]) 22 PetscScalar *G; /* such that Z = VG (band matrix)*/ member 181 beta = PetscSqrtReal(PetscRealPart(G(0, 0))); in KSPSolve_InnerLoop_PIPELCG() 182 G(0, 0) = 1.0; in KSPSolve_InnerLoop_PIPELCG() 192 …for (j = PetscMax(0, it - 3 * l + 1); j <= it - 2 * l; j++) G(j, it - l + 1) = G(it - 2 * l + 1, j… in KSPSolve_InnerLoop_PIPELCG() 198 for (j = PetscMax(it - 3 * l + 1, 0); j <= it - l + 1; ++j) G(j, it - l + 1) *= invbeta2; in KSPSolve_InnerLoop_PIPELCG() 203 …for (k = PetscMax(it - 3 * l + 1, 0); k <= j - 1; ++k) sum_dummy = sum_dummy + G(k, j) * G(k, it -… in KSPSolve_InnerLoop_PIPELCG() 204 G(j, it - l + 1) = (G(j, it - l + 1) - sum_dummy) / G(j, j); in KSPSolve_InnerLoop_PIPELCG() 208 …Max(it - 3 * l + 1, 0); k <= it - l; ++k) sum_dummy = sum_dummy + G(k, it - l + 1) * G(k, it - l +… in KSPSolve_InnerLoop_PIPELCG() [all …]
|
| /petsc/src/sys/classes/draw/utils/ |
| H A D | cmap.c | 55 static PetscErrorCode PetscDrawCmap_Hue(int mapsize, unsigned char R[], unsigned char G[], unsigned… in PetscDrawCmap_Hue() argument 66 G[i] = (unsigned char)(255 * PetscMin(g, 1.0)); in PetscDrawCmap_Hue() 72 static PetscErrorCode PetscDrawCmap_Gray(int mapsize, unsigned char R[], unsigned char G[], unsigne… in PetscDrawCmap_Gray() argument 75 …for (int i = 0; i < mapsize; i++) R[i] = G[i] = B[i] = (unsigned char)((255.0 * i) / (mapsize - 1)… in PetscDrawCmap_Gray() 79 static PetscErrorCode PetscDrawCmap_Jet(int mapsize, unsigned char R[], unsigned char G[], unsigned… in PetscDrawCmap_Jet() argument 119 G[i] = (unsigned char)(255 * PetscMin(g, 1.0)); in PetscDrawCmap_Jet() 125 static PetscErrorCode PetscDrawCmap_Hot(int mapsize, unsigned char R[], unsigned char G[], unsigned… in PetscDrawCmap_Hot() argument 155 G[i] = (unsigned char)(255 * PetscMin(g, 1.0)); in PetscDrawCmap_Hot() 161 static PetscErrorCode PetscDrawCmap_Bone(int mapsize, unsigned char R[], unsigned char G[], unsigne… in PetscDrawCmap_Bone() argument 164 (void)PetscDrawCmap_Hot(mapsize, R, G, B); in PetscDrawCmap_Bone() [all …]
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex85.c | 7 Mat G, A, B, C, D; in main() local 18 PetscCall(MatCreate(PETSC_COMM_WORLD, &G)); in main() 19 PetscCall(MatSetSizes(G, m, n, M, N)); in main() 20 PetscCall(MatSetFromOptions(G)); in main() 21 PetscCall(PetscObjectTypeCompare((PetscObject)G, MATNEST, &isnest)); in main() 22 PetscCall(MatGetOwnershipRange(G, &rStart, NULL)); in main() 23 PetscCall(MatGetOwnershipRangeColumn(G, &cStart, NULL)); in main() 70 PetscCall(MatNestSetSubMats(G, 2, NULL, 2, NULL, submat)); in main() 84 PetscCall(MatSetLocalToGlobalMapping(G, rowMap, colMap)); in main() 96 PetscCall(MatGetLocalSubMatrix(G, row, col, &A)); in main() [all …]
|
| /petsc/src/snes/linesearch/impls/nleqerr/ |
| H A D | linesearchnleqerr.c | 36 Vec X, F, Y, W, G; in SNESLineSearchApply_NLEQERR() local 49 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, &G)); in SNESLineSearchApply_NLEQERR() 79 PetscCall(VecCopy(F, G)); in SNESLineSearchApply_NLEQERR() 99 PetscCall(VecWAXPY(G, +1.0, Y, W)); /* bar_delta_x - delta_x; +1 because Y is -delta_x */ in SNESLineSearchApply_NLEQERR() 100 PetscCall(VecNormBegin(G, NORM_2, &gnorm)); in SNESLineSearchApply_NLEQERR() 101 PetscCall(VecNormEnd(G, NORM_2, &gnorm)); in SNESLineSearchApply_NLEQERR() 147 PetscCall(VecWAXPY(G, -lambda, Y, X)); in SNESLineSearchApply_NLEQERR() 165 PetscCall(SNESComputeFunction(snes, W, G)); in SNESLineSearchApply_NLEQERR() 168 PetscCall(KSPSolve(snes->ksp, G, W)); in SNESLineSearchApply_NLEQERR() 185 PetscCall(VecWAXPY(G, -(1.0 - lambda), Y, W)); in SNESLineSearchApply_NLEQERR() [all …]
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | gmsh-3d-binary-64.msh | 30 GHIJKLMNOPQR… 48 …GLH"HMG#ILG$G… 50 …*E*^UF+[,G,[… 54 …{b��-}~�G|…
|
| H A D | square_periodic_bin.msh | 11 …?�@(���w�� @`��d�Z�?)ބb[���?9]C,>�?*ߙv�G@�<�l.@+… 14 …8@R�� r@E:6���v@k��,I@F�r�Zd�"@�uj����?G����~�?�)�����?… 19 �!@dx��#�&@b��on@e��п@G�JG���?f 34 …G*;WHX<I(b@… 35 …Xb&Sbc)G0d… 36 …Gg%DEh68`i… 38 …G�,DR�,SE�…
|
| H A D | gmsh-3d-binary-32.msh | 32 GHIJKLMNOPQRS������?@������?@������?���… 54 …GGROKOSK"L!!LI… 56 …&]XC*_+DW_*E*\WF+_,G,_VHV… 62 …G/R�HQc�7.d~<Eq}Xn]�…
|
| H A D | square_bin_physnames.msh | 26 …z���JQ�?�X��G>�?C��03���?R���JQ�?D2O�����?��a�r��?E�Bdj]�?��03… 30 …r�C��^��?f������?se������?��-IM�?t�(ffff�?�G��^��?u �-… 44 …E /XDF E1YG W9b… 45 GN HO IP … 46 …GT .B�U /D�V … 47 �� ^F�� _G��… 48 …U�A� .F�� /G��…
|
| H A D | hybrid_hexwedge.msh | 9 …@�@Y�@�@Y�@�@@�@Y�{�G��@��Q�>�@Y�… 14 …G��@@@Y�HX�@@�@Y�I@�@��@Y�JR���9�@… 19 …�|�@@�@��Q�7�@�V�@����(\9�@�G����@�… 31 …����&-,%����G_Tv��� G 34 …Gd�}��6GVd�v��7CeR|���8…
|
| H A D | surfacesphere_bin.msh | 11 t��?Y��]2ʿB���a��ǿr/��Q�?������C��Ϡ�<ۿ�5G���?��,ڼ�D�l!�����y�!��?��i��=ѿE� 12 �����V,�@V��?:%�p���F���vNݿ2���Hq�?��ۦr��G>�#���鿻�Y� 16 G�%�c���a��ĿS�}n���ǿ�h�\f��zĥt�VǿT*%����ʿ�ԤϽ�ſG��bN��Ug p�&�Y�Wٗ���l!�Vi���E�� F… 19 …�Ze�?"(��U�?{�81¶�&`0�1i�?�~/A�?|��w��ѿ��U9���?��¿v�?}D;���㿏e�G���?:�S����?~P}���… 21 …��?�8^���?v"ӥ�:�?�y"�N�?�H��bN��?,%�����?�ԤϽ��?��a���?ٕ���U�?G�%�c�?�ĥt�V�?�}n�… 24 &Р�<ۿ@%G��濲imi��=�?`�!���濴�x�!�俳V�p���?{S�������@V�ѿ�1ۦr��?y�vNݿ%���Hqҿ����o ��?�#… 26 …e�?�x��w���?�rn��ǿV{m���?��~/A�?�81¶�&`0�1i�?�9�S����?B;���㿎e�G���?���¿v�?��w��ѿ… 34 …,FE7GEG7…
|
| /petsc/src/binding/petsc4py/demo/legacy/taosolve/ |
| H A D | rosenbrock.py | 43 def formGradient(self, tao, x, G): argument 47 G.zeroEntries() 51 G[2*i] = -4*alpha*t1*x[2*i] - 2*t2; 52 G[2*i+1] = 2*alpha*t1; 54 def formObjGrad(self, tao, x, G): argument 59 G.zeroEntries() 64 G[2*i] = -4*alpha*t1*x[2*i] - 2*t2; 65 G[2*i+1] = 2*alpha*t1;
|
| /petsc/src/tao/interface/ |
| H A D | taosolver_fg.c | 127 PetscErrorCode TaoComputeGradient(Tao tao, Vec X, Vec G) in TaoComputeGradient() argument 134 PetscValidHeaderSpecific(G, VEC_CLASSID, 3); in TaoComputeGradient() 136 PetscCheckSameComm(tao, 1, G, 3); in TaoComputeGradient() 139 PetscCall(PetscLogEventBegin(TAO_GradientEval, tao, X, G, NULL)); in TaoComputeGradient() 140 … PetscCallBack("Tao callback gradient", (*tao->ops->computegradient)(tao, X, G, tao->user_gradP)); in TaoComputeGradient() 141 PetscCall(PetscLogEventEnd(TAO_GradientEval, tao, X, G, NULL)); in TaoComputeGradient() 144 PetscCall(PetscLogEventBegin(TAO_ObjGradEval, tao, X, G, NULL)); in TaoComputeGradient() 145 …jective/gradient", (*tao->ops->computeobjectiveandgradient)(tao, X, &dummy, G, tao->user_objgradP)… in TaoComputeGradient() 146 PetscCall(PetscLogEventEnd(TAO_ObjGradEval, tao, X, G, NULL)); in TaoComputeGradient() 151 PetscCall(TaoTestGradient(tao, X, G)); in TaoComputeGradient() [all …]
|
| /petsc/src/mat/graphops/color/utils/ |
| H A D | weights.c | 7 Mat G = mc->mat; in MatColoringCreateLexicalWeights() local 10 PetscCall(MatGetOwnershipRange(G, &s, &e)); in MatColoringCreateLexicalWeights() 20 Mat G = mc->mat; in MatColoringCreateRandomWeights() local 26 PetscCall(MatGetOwnershipRange(G, &s, &e)); in MatColoringCreateRandomWeights() 35 PetscErrorCode MatColoringGetDegrees(Mat G, PetscInt distance, PetscInt *degrees) in MatColoringGetDegrees() argument 51 PetscCall(MatGetOwnershipRange(G, &s, &e)); in MatColoringGetDegrees() 53 PetscCall(ISCreateStride(PetscObjectComm((PetscObject)G), n, s, 1, &ris)); in MatColoringGetDegrees() 54 PetscCall(MatIncreaseOverlap(G, 1, &ris, distance)); in MatColoringGetDegrees() 56 PetscCall(MatCreateSubMatrices(G, 1, &ris, &ris, MAT_INITIAL_MATRIX, &lGs)); in MatColoringGetDegrees() 59 …PetscCheck(isSEQAIJ, PetscObjectComm((PetscObject)G), PETSC_ERR_SUP, "Requires an MPI/SEQAIJ Matri… in MatColoringGetDegrees() [all …]
|
| /petsc/src/mat/graphops/color/impls/jp/ |
| H A D | jp.c | 32 Mat G = mc->mat, dG, oG; in MCJPGreatestWeight_Private() local 46 PetscCall(PetscObjectTypeCompare((PetscObject)G, MATSEQAIJ, &isSeq)); in MCJPGreatestWeight_Private() 47 PetscCall(PetscObjectTypeCompare((PetscObject)G, MATMPIAIJ, &isMPI)); in MCJPGreatestWeight_Private() 48 …PetscCheck(isSeq || isMPI, PetscObjectComm((PetscObject)G), PETSC_ERR_ARG_WRONGSTATE, "MatColoring… in MCJPGreatestWeight_Private() 55 aij = (Mat_MPIAIJ *)G->data; in MCJPGreatestWeight_Private() 67 PetscCall(MatGetLayouts(G, &layout, NULL)); in MCJPGreatestWeight_Private() 72 dG = G; in MCJPGreatestWeight_Private() 132 Mat G = mc->mat, dG, oG; in MCJPInitialLocalColor_Private() local 148 PetscCall(MatGetOwnershipRange(G, &s, &e)); in MCJPInitialLocalColor_Private() 150 PetscCall(PetscObjectBaseTypeCompare((PetscObject)G, MATSEQAIJ, &isSeq)); in MCJPInitialLocalColor_Private() [all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/dfp/ |
| H A D | dfp.c | 19 Vec G = X; in DFPKernel_Recursive() local 29 PetscCall(LMBasisGetWorkVec(S, &G)); in DFPKernel_Recursive() 30 PetscCall(VecCopy(X, G)); in DFPKernel_Recursive() 38 PetscCall(VecDot(G, y_i, &yitx)); in DFPKernel_Recursive() 43 PetscCall(VecAXPY(G, -a, s_i)); in DFPKernel_Recursive() 46 PetscCall(MatLMVMApplyJ0Mode(mode)(B, G, BX)); in DFPKernel_Recursive() 63 PetscCall(LMBasisRestoreWorkVec(S, &G)); in DFPKernel_Recursive() 121 Vec G = X; in DFPKernel_Dense() local 137 PetscCall(LMBasisGetWorkVec(Y, &G)); in DFPKernel_Dense() 138 PetscCall(VecCopy(X, G)); in DFPKernel_Dense() [all …]
|
| /petsc/src/snes/linesearch/impls/bt/ |
| H A D | linesearchbt.c | 61 Vec X, F, Y, W, G; in SNESLineSearchApply_BT() local 75 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, &G)); in SNESLineSearchApply_BT() 103 PetscCall(VecCopy(F, G)); in SNESLineSearchApply_BT() 143 PetscCall((*linesearch->ops->snesfunc)(snes, W, G)); in SNESLineSearchApply_BT() 146 PetscCall((*linesearch->ops->vinorm)(snes, G, W, &gnorm)); in SNESLineSearchApply_BT() 148 PetscCall(VecNorm(G, NORM_2, &gnorm)); in SNESLineSearchApply_BT() 207 PetscCall((*linesearch->ops->snesfunc)(snes, W, G)); in SNESLineSearchApply_BT() 210 PetscCall((*linesearch->ops->vinorm)(snes, G, W, &gnorm)); in SNESLineSearchApply_BT() 212 PetscCall(VecNorm(G, NORM_2, &gnorm)); in SNESLineSearchApply_BT() 286 PetscCall((*linesearch->ops->snesfunc)(snes, W, G)); in SNESLineSearchApply_BT() [all …]
|
| /petsc/lib/petsc/bin/maint/abi-compliance-checker/modules/Internals/ |
| H A D | Basic.pm | 534 my $G = undef; 537 $G = "macos"; 540 $G = "bsd"; 543 $G = "beos"; 546 $G = "symbian"; 549 $G = "windows"; 552 $G = "solaris"; 556 $G = "linux"; 559 return $G;
|
| /petsc/src/sys/classes/draw/impls/x/ |
| H A D | xcolor.c | 41 unsigned char G[PETSC_DRAW_MAXCOLOR - PETSC_DRAW_BASIC_COLORS]; in PetscDrawSetUpColormap_Shared() local 61 PetscCall(PetscDrawUtilitySetCmap(NULL, ncolors, R, G, B)); in PetscDrawSetUpColormap_Shared() 64 colordef.green = (unsigned short)(G[i] << 8); in PetscDrawSetUpColormap_Shared() 71 gCpalette[k][1] = G[i]; in PetscDrawSetUpColormap_Shared() 91 unsigned char G[PETSC_DRAW_MAXCOLOR - PETSC_DRAW_BASIC_COLORS]; in PetscDrawSetUpColormap_Private() local 128 PetscCall(PetscDrawUtilitySetCmap(NULL, ncolors, R, G, B)); in PetscDrawSetUpColormap_Private() 131 colordef.green = (unsigned short)(G[i] << 8); in PetscDrawSetUpColormap_Private() 151 gCpalette[k][1] = G[i]; in PetscDrawSetUpColormap_Private()
|
| /petsc/src/snes/impls/gs/ |
| H A D | gssecant.c | 11 Vec W, G, F; in SNESComputeNGSDefaultSecant() local 29 G = snes->work[1]; in SNESComputeNGSDefaultSecant() 69 PetscCall((*func)(snes, W, G, fctx)); in SNESComputeNGSDefaultSecant() 71 if (B) PetscCall(VecAXPY(G, -1.0, B)); in SNESComputeNGSDefaultSecant() 79 PetscCall(VecGetArrayRead(G, &ga)); in SNESComputeNGSDefaultSecant() 98 PetscCall(VecRestoreArrayRead(G, &ga)); in SNESComputeNGSDefaultSecant() 118 PetscCall(VecSwap(F, G)); in SNESComputeNGSDefaultSecant()
|