Home
last modified time | relevance | path

Searched refs:Q (Results 1 – 25 of 79) sorted by relevance

1234

/petsc/src/ksp/ksp/tutorials/
H A Dex81a.c6 #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 Dex81.c5 #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 Dex54.c73 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 Dex84.c6 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 Dex42.c36 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 Dhybrid_hexwedge.msh9 …����Y�����������������@�@������Y���������@�@�����@�@������Y����{�G��@��Q�>�@������Y������(\�…
11 �@�Q��o�@������Y���������ǥ@��(\�o�@������Y���������@�@�����p�@������Y�����������@�����p�@������…
13 …��@�@������Y�$������������������@������Y�%��������@�@�������@������Y�&�����Q��@R��+��@������Y�'�…
14Q�7�@�����V�@������Y�N������(\9�@�G����@������Y�O��������|�@�����p�@������Y�P��������X�@�����p�…
16 �@�Q��o�@�����������������ǥ@��(\�o�@�����������������@�@�����p�@�������������������@�����p�@������…
18 …��@�@���������������������������@�����������������@�@�������@��������������Q��@R��+��@����������…
19 …����m�@�������������������@�����X�@�����������������|�@�����@�@�������������Q�7�@�����V�@�������…
34 …��v�����������7��� ���C���e���R���|���������������8���C������Q���e�������������������9���Q������…
H A Dmesh-3d-box-innersphere_bin.msh9 (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 …���Q���R���V���������������������Q���W���R���������������������Q���S���X����������…
H A Dgmsh-3d-binary-32.msh32 ���G���H���I���J���K���L���M���N���O���P���Q���R���S���������?��������@������?@������?�����������…
50 ��� ����� �����������(���������P��� ������!���P��������� ���Q������������Q���'�����…
54Q��� ���������O������������O���G������G���R������������O���������K���O������������S�…
56 …W���^���U���N���U���_���W���O���V���_���U���P���U���Z���V������2������(���Q������ ���i���R����…
60Q�����������[���7���l���~����������{���>���|�������>���{������}�������o������{���}�������o���{�…
62 ��;���D���k���|�����W���k���\���|�����G���/���R���������H���Q���c���������7���.���d���~���…
H A Dsurfacesphere_bin.msh9 …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�L׿qW/��O���q�M��W��ɌC�ʿ1Nfœ�տP���v"ӥ�:׿�y"�N�9^���ſQ���V�h
16 …�*%����ʿ�ԤϽ�ſG��bN��U���g p�&�׿Y��Wٗ���l!�V���i���E�� F#[$ӿ&����W����2�Q�俠^jP`�����ҿX���)�t…
18 …�U�h�����Wٗ�?���l!�?�*p�&�׿i���x#[$�?�����?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�'�X￱���Ts,ڼ�?
25 ӿ����J��|���?���pGοb���Q忷���4wG�{�?u�W��e�j��}�?����B�>U�2�?��R��v�6f���?����>��� :�?H/��mٿ�…
[all …]
H A Dsquare_periodic_bin.msh11 …��������#���������������:�?��������$��������������e1��?��������%����4L��@Q%�P��@��������&���0�…
14 …�jx2#@ϧ͹���?��������B�����흝�?ͧ͹���?��������C����C1�@���I��@��������D����Q�;8@R�� r@�������…
15 …a��?��������N�����Ֆ�:�?d�T���?��������O��������@�V�B�@��������P���z�8��}#@ ����5�?��������Q���
17 …������W���6��)�!@Č�w�6@��������X���EMDU��?���O@��������Y���˹���\�?p� Q�!@��������Z�����bz!@…
34Q���$������������@����������2���T���I������������A����������+���_���=������������B�����…
35 …%���_���D������������X�������������Y���6������������Y����������#���Q���B������������Z�…
38Q�����������������������.���P���A�����������������������������1���������������������…
H A Dgmsh-3d-binary-64.msh30 �������G�������H�������I�������J�������K�������L�������M�������N�������O�������P�������Q�������R���…
44 …���� ��������������!�������P��������������������� �������Q����������������������������Q�����…
48Q������� ���������������������O����������������������������O�������J��������������J�������R�…
50 …X�������[�������T�������P�������T�������\�������X����������2������(�������Q�������������� �����…
54Q�����������������������u�������`�������z�����������������������`�������>�������F�����������������…
H A Dsquare_bin_physnames.msh28 …���?�\�����?��������O���`�����?r^�8���?��������P������k#�?T�c���?��������Q����_�8���?����k#�?��…
30 …�����v�����Uץ�?f>޲$��?��������w���� i��?��Uץ�?��������x���7A޲$��?pT�Q��?��������y���kpT�Q��?p�
45 …���O��� ���������I���������������P��� ���������v���H������������Q��� ���������w���I�…
46Q������������o��� ������$���F���^������������p��� ������ ���G���_������������q��� ����…
48Q���;���y���������������� ������P���=���w���������������� ������O���x���&��������������…
H A Dgmsh-3d-binary.msh245 �����?��������P��<���?P����f�x���?���������f�x���?Q����MUUU��?���������������?R����������?��������…
55 …��������������!���P������������������������ ���Q���������������������������Q���…
59Q��� ������������������������O���������������������������O���J����������������…
61 …���O������(���W���Y���V������������P������(���V���Z���W������������Q������2������ ���j�…
65Q�����������������������d���n���`���z���|�������������������d���`���>���F���|����������…
/petsc/include/
H A Dpetscfe.h176 …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 Dplexlocalizationletkf.kokkos.cxx69 …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 Ddmceed.c91 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 Dlcd.c13 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 Dex11.h228 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 Dtfqmr.c16 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 Dcgs.c19 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 Dpipelcg.c20 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 Dal.c210 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 Dpipebcgs.c19 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 Dpipecr.c25 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 Dswarmpic_plex.c25 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 …]

1234