| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex81a.c | 6 #define Q 5 /* everything is hardwired for a 5x5 MatNest for now */ macro 12 Mat array[Q * Q], A, a; in main() 14 IS rows[Q]; in main() 25 for (i = 0; i < Q * Q; ++i) array[i] = NULL; in main() 26 for (i = 0; i < Q; ++i) { in main() 28 …(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main() 30 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main() 32 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main() 34 PetscCall(MatAssemblyBegin(array[(Q + 1) * i], MAT_FINAL_ASSEMBLY)); in main() 35 PetscCall(MatAssemblyEnd(array[(Q + 1) * i], MAT_FINAL_ASSEMBLY)); in main() [all …]
|
| H A D | ex81.c | 5 #define Q 5 /* everything is hardwired for a 5x5 MatNest for now */ macro 11 Mat array[Q * Q], A, a; in main() 13 IS rows[Q]; in main() 14 PetscInt i, j, *outer, M, N, found = Q; in main() 29 …(j == Q, PETSC_COMM_WORLD, PETSC_ERR_USER, "Sum of outer fieldsplit sizes (%" PetscInt_FMT ") grea… in main() 33 for (i = 0; i < Q * Q; ++i) array[i] = NULL; in main() 34 for (i = 0; i < Q; ++i) { in main() 36 …(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main() 38 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main() 40 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main() [all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex54.c | 73 Mat Q, C, V, A, B; in main() local 85 PetscCall(MatCreate(PETSC_COMM_WORLD, &Q)); in main() 86 PetscCall(MatSetType(Q, MATDENSE)); in main() 87 PetscCall(fill(Q, v)); in main() 89 PetscCall(MatCreateVecs(Q, &a, NULL)); in main() 90 PetscCall(MatCreateNormalHermitian(Q, &C)); in main() 96 PetscCall(KSPSetOperators(QRsolver, Q, Q)); in main() 97 PetscCall(MatViewFromOptions(Q, NULL, "-sys_view")); in main() 106 PetscCall(MatConvert(Q, MATAIJ, MAT_INPLACE_MATRIX, &Q)); in main() 107 PetscCall(MatCreateNormalHermitian(Q, &C)); in main() [all …]
|
| H A D | ex84.c | 6 PetscErrorCode LSCLoadOperators(Mat *A, Mat *Q, Mat *L, Vec *rhs, IS *velocity, IS *pressure) in LSCLoadOperators() argument 14 PetscCall(MatCreate(PETSC_COMM_WORLD, Q)); in LSCLoadOperators() 27 PetscCall(MatLoad(*Q, viewer)); in LSCLoadOperators() 35 Mat A, Q, L = NULL; in port_lsd_bfbt() local 50 if (commute_lsc) PetscCall(LSCLoadOperators(&A, &Q, &L, &b, &isu, &isp)); in port_lsd_bfbt() 51 else PetscCall(LSCLoadOperators(&A, &Q, NULL, &b, &isu, &isp)); in port_lsd_bfbt() 71 PetscCall(PetscObjectCompose((PetscObject)lsc_pc_pmat, "LSC_Qscale", (PetscObject)Q)); in port_lsd_bfbt() 78 PetscCall(MatDestroy(&Q)); in port_lsd_bfbt()
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex42.c | 36 CEED_QFUNCTION(Mass)(PetscCtx ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const … in Mass() argument 41 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) v[i] = qdata[i] * u[i]; in Mass() 70 CEED_QFUNCTION(SetupMassGeoCube)(PetscCtx ctx, const CeedInt Q, const CeedScalar *const *in, CeedSc… in SetupMassGeoCube() argument 75 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in SetupMassGeoCube() 79 {J[i + Q * 0], J[i + Q * 3]}, in SetupMassGeoCube() 80 {J[i + Q * 1], J[i + Q * 4]}, in SetupMassGeoCube() 81 {J[i + Q * 2], J[i + Q * 5]} in SetupMassGeoCube() 98 …qdata[i + Q * 0] = (dxdX[0][0] * dxdX[1][1] - dxdX[1][0] * dxdX[0][1]) * w[i]; /* det J * weight */ in SetupMassGeoCube() 132 CEED_QFUNCTION(SetupMassGeoSphere)(PetscCtx ctx, const CeedInt Q, const CeedScalar *const *in, Ceed… in SetupMassGeoSphere() argument 137 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in SetupMassGeoSphere() [all …]
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | hybrid_hexwedge.msh | 9 …Y�@�@Y�@�@@�@Y�{�G��@��Q�>�@Y���(\�… 11 �@�Q��o�@Y��ǥ@��(\�o�@Y�@�@p�@Y���@p�@… 13 …@�@Y�$��@Y�%@�@��@Y�&��Q��@R��+��@Y�'… 14 …Q�7�@�V�@Y�N���(\9�@�G����@Y�O|�@p�@Y�PX�@p�… 16 �@�Q��o�@��ǥ@��(\�o�@�@�@p�@���@p�@… 18 …@�@���@�@�@��@���Q��@R��+��@�… 19 …m�@���@X�@�|�@@�@��Q�7�@�V�@… 34 …v��7CeR|���8CQe����9Q…
|
| H A D | mesh-3d-box-innersphere_bin.msh | 9 (g�������w�ɿUDo�?�Q��ɐ�5�/7�Ƴ� 10 �TDo�?�Q��ɐ�Z�/7�Ƴ?fZu#��?�J 12 �*��?�Q���k�?U L�\ο8_[f���?�R�}�˿��a��ͤ�9�O�$A��?(Jeܱ��?gz�j-E�?:����ƿ�6|zƿ�� 13 ��;��*��ɿ��~8���?[�Y�h�?<�̒Ր�?��M��ƿf��ږ�?=L�Ɛ��ÿ%|��v��?�-Q�`��?>+� 17 …yO�?+�>'�?0�4?����O�0d��?�? �g�If������[��P����r��o��ww��?Z�9�����Q�… 19 …?�@S>�K�?���_/���?I�C`�v�?s>GPm%ӿ��v���)�?�xI!��?���۵�?��������?�Q 9���?tS� |�?���… 26 …QRVQWRQSX…
|
| H A D | gmsh-3d-binary-32.msh | 32 GHIJKLMNOPQRS������?@������?@������?���… 50 (P !PQQ'… 54 …QOOGGROKOS… 56 …W^UNU_WOV_UPUZV2(Q iR… 60 …Q��[7l~�{>|�>{}�o{}�o{… 62 ;Dk|Wk\|G/R�HQc�7.d~…
|
| H A D | surfacesphere_bin.msh | 9 …7���?7�gm=o�俙�ׁ��?��V����8�����V�w_r�`,�?�Fr46ֿ9�!"J�ݿ�W�%���?�^Q块�:8'���SϿ�����?�… 11 t��?Y��]2ʿB���a��ǿr/��Q�?������C��Ϡ�<ۿ�5G���?��,ڼ�D�l!�����y�!��?��i��=ѿE� 13 �?֏o ��H�>�pGο����Q�?ޡ�|���I��[%H��.k݉��"=�Ċ�J_�h�ܿŧI���I��+�߿KĈtR���F����ؿ�0~… 14 …�<�-�$пӗ�ͥ9�A�7��xٿN#��`�ӿ=3�0�LqW/��Oq�M��W��ɌC�ʿ1Nfœ�տPv"ӥ�:�y"�N�9^��ſQV�h 16 …*%����ʿ�ԤϽ�ſG��bN��Ug p�&�Y�Wٗ���l!�Vi���E�� F#[$ӿ&����W�2�Q�俠^jP`�����ҿX)�t… 18 …�U�h��Wٗ�?���l!�?�*p�&�ix#[$�?�����?q��E��j�8^jP`�?0�����?Q�Q��k%��8���?�lS8k�… 21 …h�\f��?���l!�?d p�&��?Z�Wٗ�?�&����?h���E��? F#[$�?������?�2�Q��?�^jP`�?�%�lS8k… 22 '�b濾�M1ӿ9je����?���~Q�㿿ڍw��ie�{�?�7$��3�пEs}˴�����&m#�?�wВ_�����Aο��� ��?��… 23 t�ǿ�������?��a��ǿ�/��Qǿ���0���?�[LX'Ŀ2�'�XTs,ڼ�? 25 ӿ�J��|���?���pGοb���Q忷4wG�{�?u�W��e�j��}�?�B�>U�2�?��R��v�6f���?�>���:�?H/��mٿ�… [all …]
|
| H A D | square_periodic_bin.msh | 11 …#����:�?$���e1��?%�4L��@Q%�P��@&0�… 14 …�jx2#@ϧ���?B��흝�?ͧ���?C�C1�@���I��@D�Q�;8@R�� r@… 15 …a��?N��Ֆ�:�?d�T���?O����@�V�B�@Pz�8��}#@ ����5�?Q 17 …W6��)�!@Č�w�6@XEMDU��?���O@Y˹���\�?p� Q�!@Zbz!@… 34 …Q$@2TIA+_=B… 35 …%_DXY6Y#QBZ… 38 …Q�.PA�1�…
|
| H A D | gmsh-3d-binary-64.msh | 30 GHIJKLMNOPQR… 44 … !PQQ… 48 …QOOJJR… 50 …X[TPT\X2(Q … 54 …Q��u`z��`>F��…
|
| H A D | square_bin_physnames.msh | 28 …���?�\�����?O`�����?r^�8���?P���k#�?T�c���?Q�_�8���?����k#�?… 30 …v��Uץ�?f>$��?w�i��?��Uץ�?x7A$��?pT�Q��?ykpT�Q��?p� 45 …O IP vHQ wI… 46 …Qo $F^p G_q … 48 …Q;y� P=w� Ox&�…
|
| H A D | gmsh-3d-binary.msh2 | 45 �?P��<���?P�f�x���?�f�x���?Q�MUUU��?�������?R�������?… 55 …!PQQ… 59 …QOOJ… 61 …O(WYVP(VZWQ2 j… 65 …Q��dn`z|�d`>F|…
|
| /petsc/include/ |
| H A D | petscfe.h | 176 …CEED_QFUNCTION(PlexQFunction##fname)(PetscCtx ctx, const CeedInt Q, const CeedScalar *const *in, C… 183 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) \ 187 const CeedScalar x[2] = {qdata[i + Q * 1], qdata[i + Q * 2]}; \ 189 {qdata[i + Q * 3], qdata[i + Q * 5]}, \ 190 {qdata[i + Q * 4], qdata[i + Q * 6]} \ 192 …ar u_x[2] = {invJ[0][0] * du[i + Q * 0] + invJ[1][0] * du[i + Q * 1], invJ[0][1] * du[i + Q * 0] +… 201 dv[i + Q * 0] = qdata[i + Q * 0] * (invJ[0][0] * f1[0] + invJ[0][1] * f1[1]); \ 202 dv[i + Q * 1] = qdata[i + Q * 0] * (invJ[1][0] * f1[0] + invJ[1][1] * f1[1]); \ 203 v[i] = qdata[i + Q * 0] * f0[0]; \
|
| /petsc/src/dm/impls/plex/kokkos/ |
| H A D | plexlocalizationletkf.kokkos.cxx | 69 …scInt n_obs_vertex, const PetscInt n_obs_local, const PetscInt n_dof, Vec Vecxyz[3], Mat H, Mat *Q) in DMPlexGetLETKFLocalizationMatrix() argument 78 PetscAssertPointer(Q, 6); in DMPlexGetLETKFLocalizationMatrix() 113 PetscCall(MatCreate(comm, Q)); in DMPlexGetLETKFLocalizationMatrix() 114 PetscCall(MatSetSizes(*Q, n_vert_local, n_obs_local, PETSC_DETERMINE, n_obs_global)); in DMPlexGetLETKFLocalizationMatrix() 115 PetscCall(MatSetType(*Q, MATAIJ)); in DMPlexGetLETKFLocalizationMatrix() 116 PetscCall(MatSeqAIJSetPreallocation(*Q, n_obs_vertex, NULL)); in DMPlexGetLETKFLocalizationMatrix() 117 PetscCall(MatMPIAIJSetPreallocation(*Q, n_obs_vertex, NULL, n_obs_vertex, NULL)); in DMPlexGetLETKFLocalizationMatrix() 118 PetscCall(MatSetFromOptions(*Q)); in DMPlexGetLETKFLocalizationMatrix() 119 PetscCall(MatSetUp(*Q)); in DMPlexGetLETKFLocalizationMatrix() 121 …PetscCall(PetscInfo((PetscObject)*Q, "Computing LETKF localization matrix: %" PetscInt_FMT " verti… in DMPlexGetLETKFLocalizationMatrix() [all …]
|
| /petsc/src/dm/interface/ |
| H A D | dmceed.c | 91 CEED_QFUNCTION(Geometry2D)(PetscCtx ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *… in Geometry2D() argument 96 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in Geometry2D() 99 {Jac[i + Q * 0], Jac[i + Q * 2]}, in Geometry2D() 100 {Jac[i + Q * 1], Jac[i + Q * 3]} in Geometry2D() 104 qdata[i + Q * 0] = det * w[i]; in Geometry2D() 105 qdata[i + Q * 1] = x[i + Q * 0]; in Geometry2D() 106 qdata[i + Q * 2] = x[i + Q * 1]; in Geometry2D() 107 qdata[i + Q * 3] = J[1][1] / det; in Geometry2D() 108 qdata[i + Q * 4] = -J[1][0] / det; in Geometry2D() 109 qdata[i + Q * 5] = -J[0][1] / det; in Geometry2D() [all …]
|
| /petsc/src/ksp/ksp/impls/lcd/ |
| H A D | lcd.c | 13 PetscCall(VecDuplicateVecs(ksp->work[0], restart + 1, &lcd->Q)); in KSPSetUp_LCD() 78 PetscCall(KSP_PCApply(ksp, Z, lcd->Q[it])); in KSPSolve_LCD() 83 PetscCall(VecDot(lcd->P[it], lcd->Q[it], &den)); in KSPSolve_LCD() 87 PetscCall(VecAXPY(R, -alfa, lcd->Q[it])); in KSPSolve_LCD() 102 PetscCall(KSP_PCApply(ksp, Z, lcd->Q[it + 1])); in KSPSolve_LCD() 105 PetscCall(VecDot(lcd->P[j], lcd->Q[it + 1], &num)); in KSPSolve_LCD() 107 PetscCall(VecDot(lcd->P[j], lcd->Q[j], &den)); in KSPSolve_LCD() 110 PetscCall(VecAXPY(lcd->Q[it + 1], beta, lcd->Q[j])); in KSPSolve_LCD() 130 if (lcd->Q) PetscCall(VecDestroyVecs(lcd->restart + 1, &lcd->Q)); in KSPReset_LCD()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex11.h | 228 CEED_QFUNCTION(PhysicsRiemann_SW_Rusanov_CEED)(PetscCtx ctx, CeedInt Q, const CeedScalar *const in[… in PhysicsRiemann_SW_Rusanov_CEED() argument 239 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in PhysicsRiemann_SW_Rusanov_CEED() 241 const CeedScalar qL[3] = {xL[i + Q * 0], xL[i + Q * 1], xL[i + Q * 2]}; in PhysicsRiemann_SW_Rusanov_CEED() 242 const CeedScalar qR[3] = {xR[i + Q * 0], xR[i + Q * 1], xR[i + Q * 2]}; in PhysicsRiemann_SW_Rusanov_CEED() 243 const CeedScalar n[2] = {geom[i + Q * 0], geom[i + Q * 1]}; in PhysicsRiemann_SW_Rusanov_CEED() 247 PetscPrintf(PETSC_COMM_SELF, "Face %d Normal\n", (int)info[i + Q * 0]); in PhysicsRiemann_SW_Rusanov_CEED() 249 PetscPrintf(PETSC_COMM_SELF, "Cell %d Element Residual: left state\n", (int)info[i + Q * 1]); in PhysicsRiemann_SW_Rusanov_CEED() 251 PetscPrintf(PETSC_COMM_SELF, "Cell %d Element Residual: right state\n", (int)info[i + Q * 2]); in PhysicsRiemann_SW_Rusanov_CEED() 256 cL[i + Q * j] = -flux[j] / geom[i + Q * 2]; in PhysicsRiemann_SW_Rusanov_CEED() 257 cR[i + Q * j] = flux[j] / geom[i + Q * 3]; in PhysicsRiemann_SW_Rusanov_CEED() [all …]
|
| /petsc/src/ksp/ksp/impls/tfqmr/ |
| H A D | tfqmr.c | 16 Vec X, B, V, P, R, RP, T, T1, Q, U, D, AUQ; in KSPSolve_TFQMR() local 25 Q = ksp->work[4]; in KSPSolve_TFQMR() 71 PetscCall(VecWAXPY(Q, -a, V, U)); /* q <- u - a v */ in KSPSolve_TFQMR() 72 PetscCall(VecWAXPY(T, 1.0, U, Q)); /* t <- u + q */ in KSPSolve_TFQMR() 88 PetscCall(VecAYPX(D, cf, Q)); in KSPSolve_TFQMR() 109 PetscCall(VecWAXPY(U, b, Q, R)); /* u <- r + b q */ in KSPSolve_TFQMR() 110 PetscCall(VecAXPY(Q, b, P)); in KSPSolve_TFQMR() 111 PetscCall(VecWAXPY(P, b, Q, U)); /* p <- u + b(q + b p) */ in KSPSolve_TFQMR() 112 PetscCall(KSP_PCApplyBAorAB(ksp, P, V, Q)); /* v <- K p */ in KSPSolve_TFQMR()
|
| /petsc/src/ksp/ksp/impls/cgs/ |
| H A D | cgs.c | 19 Vec X, B, V, P, R, RP, T, Q, U, AUQ; in KSPSolve_CGS() local 35 Q = ksp->work[4]; in KSPSolve_CGS() 91 PetscCall(VecWAXPY(Q, -a, V, U)); /* q <- u - a v */ in KSPSolve_CGS() 92 PetscCall(VecWAXPY(T, 1.0, U, Q)); /* t <- u + q */ in KSPSolve_CGS() 115 PetscCall(VecWAXPY(U, b, Q, R)); /* u <- r + b q */ in KSPSolve_CGS() 116 PetscCall(VecAXPY(Q, b, P)); in KSPSolve_CGS() 117 PetscCall(VecWAXPY(P, b, Q, U)); /* p <- u + b(q + b p) */ in KSPSolve_CGS() 118 PetscCall(KSP_PCApplyBAorAB(ksp, P, V, Q)); /* v <- K p */ in KSPSolve_CGS()
|
| /petsc/src/ksp/ksp/impls/cg/pipelcg/ |
| H A D | pipelcg.c | 20 Vec *Q; /* Q vectors (auxiliary bases) */ member 54 PetscCall(VecDuplicateVecs(plcg->p, 3 * (l - 1) + 1, &plcg->Q)); in KSPSetUp_PIPELCG() 71 PetscCall(VecDestroyVecs(3 * (l - 1) + 1, &plcg->Q)); in KSPReset_PIPELCG() 140 Vec *Z = plcg->Z, *U = plcg->U, *V = plcg->V, *Q = plcg->Q; in KSPSolve_InnerLoop_PIPELCG() local 165 if (it < l - 1) PetscCall(VecCopy(Z[l - it - 1], Q[3 * it])); in KSPSolve_InnerLoop_PIPELCG() 186 for (j = 0; j < l - 1; ++j) PetscCall(VecScale(Q[3 * j], 1.0 / beta)); in KSPSolve_InnerLoop_PIPELCG() 244 PetscCall(VecCopy(Q[0], V[0])); in KSPSolve_InnerLoop_PIPELCG() 258 temp = Q[3 * j + 2]; in KSPSolve_InnerLoop_PIPELCG() 259 for (i = 2; i > 0; i--) Q[3 * j + i] = Q[3 * j + i - 1]; in KSPSolve_InnerLoop_PIPELCG() 260 Q[3 * j] = temp; in KSPSolve_InnerLoop_PIPELCG() [all …]
|
| /petsc/src/snes/impls/al/ |
| H A D | al.c | 210 static PetscErrorCode SNESNewtonALComputeFunction_NEWTONAL(SNES snes, Vec X, Vec Q) in SNESNewtonALComputeFunction_NEWTONAL() argument 220 PetscCall(VecZeroEntries(Q)); in SNESNewtonALComputeFunction_NEWTONAL() 224 …PetscCallBack("SNES callback NewtonAL tangent load function", (*computealfunction)(snes, X, Q, ctx… in SNESNewtonALComputeFunction_NEWTONAL() 235 PetscCall(VecAXPY(Q, 1, al->vec_rhs_orig)); in SNESNewtonALComputeFunction_NEWTONAL() 260 PetscErrorCode SNESNewtonALComputeFunction(SNES snes, Vec X, Vec Q) in SNESNewtonALComputeFunction() argument 265 PetscValidHeaderSpecific(Q, VEC_CLASSID, 3); in SNESNewtonALComputeFunction() 267 PetscCheckSameComm(snes, 1, Q, 3); in SNESNewtonALComputeFunction() 269 PetscCall(PetscLogEventBegin(SNES_NewtonALEval, snes, X, Q, 0)); in SNESNewtonALComputeFunction() 270 PetscTryMethod(snes, "SNESNewtonALComputeFunction_C", (SNES, Vec, Vec), (snes, X, Q)); in SNESNewtonALComputeFunction() 271 PetscCall(PetscLogEventEnd(SNES_NewtonALEval, snes, X, Q, 0)); in SNESNewtonALComputeFunction() [all …]
|
| /petsc/src/ksp/ksp/impls/bcgs/pipebcgs/ |
| H A D | pipebcgs.c | 19 Vec X, B, S, R, RP, Y, Q, P2, Q2, R2, S2, W, Z, W2, Z2, T, V; in KSPSolve_PIPEBCGS() local 31 Q = ksp->work[4]; in KSPSolve_PIPEBCGS() 104 PetscCall(VecWAXPY(Q, -alpha, S, R)); /* q <- r - alpha s */ in KSPSolve_PIPEBCGS() 108 PetscCall(VecDotBegin(Q, Y, &d1)); /* d1 <- (q,y) */ in KSPSolve_PIPEBCGS() 111 PetscCall(PetscCommSplitReductionBegin(PetscObjectComm((PetscObject)Q))); in KSPSolve_PIPEBCGS() 115 PetscCall(VecDotEnd(Q, Y, &d1)); in KSPSolve_PIPEBCGS() 120 PetscCall(VecDot(Q, Q, &d1)); in KSPSolve_PIPEBCGS() 137 PetscCall(VecWAXPY(R, -omega, Y, Q)); /* r <- q - omega y */ in KSPSolve_PIPEBCGS()
|
| /petsc/src/ksp/ksp/impls/cr/pipecr/ |
| H A D | pipecr.c | 25 Vec X, B, Z, P, W, Q, U, M, N; in KSPSolve_PIPECR() local 40 Q = ksp->work[5]; in KSPSolve_PIPECR() 102 PetscCall(VecCopy(M, Q)); /* q <- m */ in KSPSolve_PIPECR() 108 PetscCall(VecAYPX(Q, beta, M)); /* q <- m + beta * q */ in KSPSolve_PIPECR() 112 PetscCall(VecAXPY(U, -alpha, Q)); /* u <- u - alpha * q */ in KSPSolve_PIPECR()
|
| /petsc/src/dm/impls/swarm/ |
| H A D | swarmpic_plex.c | 25 PetscDualSpace Q; in private_PetscFECreateDefault_scalar_pk1() local 43 PetscCall(PetscDualSpaceCreate(PetscObjectComm((PetscObject)dm), &Q)); in private_PetscFECreateDefault_scalar_pk1() 44 PetscCall(PetscDualSpaceSetType(Q, PETSCDUALSPACELAGRANGE)); in private_PetscFECreateDefault_scalar_pk1() 47 PetscCall(PetscDualSpaceSetDM(Q, K)); in private_PetscFECreateDefault_scalar_pk1() 49 PetscCall(PetscDualSpaceSetNumComponents(Q, Nc)); in private_PetscFECreateDefault_scalar_pk1() 50 PetscCall(PetscDualSpaceSetOrder(Q, order)); in private_PetscFECreateDefault_scalar_pk1() 51 PetscCall(PetscDualSpaceLagrangeSetTensor(Q, tensor)); in private_PetscFECreateDefault_scalar_pk1() 53 PetscCall(PetscDualSpaceSetType(Q, PETSCDUALSPACELAGRANGE)); in private_PetscFECreateDefault_scalar_pk1() 54 PetscCall(PetscDualSpaceSetUp(Q)); in private_PetscFECreateDefault_scalar_pk1() 61 PetscCall(PetscFESetDualSpace(*fem, Q)); in private_PetscFECreateDefault_scalar_pk1() [all …]
|