Home
last modified time | relevance | path

Searched refs:W (Results 1 – 25 of 137) sorted by relevance

123456

/petsc/src/mat/impls/baij/seq/
H A Dbaij.h276 …PetscErrorCode PetscKernel_A_gets_A_times_B_2(PetscScalar *A, const PetscScalar *B, PetscScalar *W) in PetscKernel_A_gets_A_times_B_2() argument
278 PetscCall(PetscArraycpy(W, A, 4)); in PetscKernel_A_gets_A_times_B_2()
279 A[0] = W[0] * B[0] + W[2] * B[1]; in PetscKernel_A_gets_A_times_B_2()
280 A[1] = W[1] * B[0] + W[3] * B[1]; in PetscKernel_A_gets_A_times_B_2()
281 A[2] = W[0] * B[2] + W[2] * B[3]; in PetscKernel_A_gets_A_times_B_2()
282 A[3] = W[1] * B[2] + W[3] * B[3]; in PetscKernel_A_gets_A_times_B_2()
316 …PetscErrorCode PetscKernel_A_gets_A_times_B_3(PetscScalar *A, const PetscScalar *B, PetscScalar *W) in PetscKernel_A_gets_A_times_B_3() argument
318 PetscCall(PetscArraycpy(W, A, 9)); in PetscKernel_A_gets_A_times_B_3()
319 A[0] = W[0] * B[0] + W[3] * B[1] + W[6] * B[2]; in PetscKernel_A_gets_A_times_B_3()
320 A[1] = W[1] * B[0] + W[4] * B[1] + W[7] * B[2]; in PetscKernel_A_gets_A_times_B_3()
[all …]
/petsc/src/sys/classes/draw/impls/x/
H A Dximpl.h46 static inline void PetscDrawXiSetPixVal(PetscDraw_X *W, PetscDrawXiPixVal pix) in PetscDrawXiSetPixVal() argument
48 if (W->gc.cur_pix != pix) { in PetscDrawXiSetPixVal()
49 XSetForeground(W->disp, W->gc.set, pix); in PetscDrawXiSetPixVal()
50 W->gc.cur_pix = pix; in PetscDrawXiSetPixVal()
55 …#define PetscDrawXiValidColor(W, color) PetscCheck((color) >= 0 && (color) < PETSC_DRAW_MAXCOLOR, … argument
57 #define PetscDrawXiValidColor(W, color) \ argument
62 #define PetscDrawXiSetColor(W, color) \ argument
64 PetscDrawXiValidColor(W, color); \
65 PetscDrawXiSetPixVal(W, (W)->cmapping[(color)]); \
/petsc/src/snes/linesearch/impls/basic/
H A Dlinesearchbasic.c7 Vec X, F, Y, W; in SNESLineSearchApply_Basic() local
12 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, NULL)); in SNESLineSearchApply_Basic()
22 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Basic()
23 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Basic()
26 PetscCall(SNESLineSearchPostCheck(linesearch, X, Y, W, &changed_y, &changed_w)); in SNESLineSearchApply_Basic()
28 if (!changed_w) PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Basic()
29 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Basic()
32 PetscCall((*linesearch->ops->snesfunc)(snes, W, F)); in SNESLineSearchApply_Basic()
38 PetscCall(VecNormBegin(W, NORM_2, &linesearch->xnorm)); in SNESLineSearchApply_Basic()
40 PetscCall(VecNormEnd(W, NORM_2, &linesearch->xnorm)); in SNESLineSearchApply_Basic()
[all …]
/petsc/src/snes/linesearch/impls/cp/
H A Dlinesearchcp.c7 Vec X, Y, F, W; in SNESLineSearchApply_CP() local
16 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, NULL)); in SNESLineSearchApply_CP()
47 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_CP()
48 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_CP()
49 PetscCall((*linesearch->ops->snesfunc)(snes, W, F)); in SNESLineSearchApply_CP()
51 PetscCall((*linesearch->ops->vidirderiv)(snes, F, W, Y, &fty)); in SNESLineSearchApply_CP()
102 PetscCall(VecWAXPY(W, -0.5 * (lambda + lambda_old), Y, X)); in SNESLineSearchApply_CP()
103 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_CP()
104 PetscCall((*linesearch->ops->snesfunc)(snes, W, F)); in SNESLineSearchApply_CP()
106 PetscCall((*linesearch->ops->vidirderiv)(snes, F, W, Y, &fty_mid1)); in SNESLineSearchApply_CP()
[all …]
/petsc/src/vec/vec/tests/
H A Dex34.c8 Vec V, W; in main() local
89 PetscCall(VecDuplicate(V, &W)); in main()
90 PetscCall(VecCopy(V, W)); in main()
98 PetscCall(VecNorm(W, NORM_1, &nrm1)); in main()
99 PetscCall(VecNorm(W, NORM_2, &nrm2)); in main()
106 PetscCall(VecCopy(V, W)); in main()
114 PetscCall(VecNorm(W, NORM_1, &nrm1)); in main()
115 PetscCall(VecNorm(W, NORM_2, &nrm2)); in main()
138 PetscCall(VecNorm(W, NORM_1, &nrm2)); in main()
142 PetscCall(VecNorm(W, NORM_INFINITY, &nrm4)); in main()
[all …]
H A Dex38.c8 Vec X, Y, Z, W; in main() local
55 PetscCall(VecGetSubVector(X, is2, &W)); in main()
58 PetscCall(VecView(W, viewer)); in main()
65 PetscCall(VecGetArray(W, &w)); in main()
67 PetscCall(VecRestoreArray(W, &w)); in main()
70 PetscCall(VecRestoreSubVector(X, is2, &W)); in main()
/petsc/src/snes/linesearch/impls/secant/
H A Dlinesearchsecant.c11 Vec W; in SNESLineSearchApply_Secant() local
25 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, NULL)); in SNESLineSearchApply_Secant()
49 PetscCall(VecWAXPY(W, -lambda_mid, Y, X)); in SNESLineSearchApply_Secant()
50 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Secant()
53 PetscCall((*linesearch->ops->snesfunc)(snes, W, F)); in SNESLineSearchApply_Secant()
56 PetscCall((*linesearch->ops->vinorm)(snes, F, W, &fnrm_mid)); in SNESLineSearchApply_Secant()
62 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Secant()
63 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Secant()
64 PetscCall((*linesearch->ops->snesfunc)(snes, W, F)); in SNESLineSearchApply_Secant()
67 PetscCall((*linesearch->ops->vinorm)(snes, F, W, &fnrm)); in SNESLineSearchApply_Secant()
[all …]
/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
39 W = ksp->work[4]; in KSPSolve_PIPECR()
48 PetscCall(KSP_MatMult(ksp, Amat, X, W)); /* w <- b - Ax */ in KSPSolve_PIPECR()
49 PetscCall(VecAYPX(W, -1.0, B)); in KSPSolve_PIPECR()
51 PetscCall(VecCopy(B, W)); /* w <- b (x is 0) */ in KSPSolve_PIPECR()
53 PetscCall(KSP_PCApply(ksp, W, U)); /* u <- Bw */ in KSPSolve_PIPECR()
59 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECR()
63 PetscCall(KSP_MatMult(ksp, Amat, U, W)); in KSPSolve_PIPECR()
77 PetscCall(KSP_PCApply(ksp, W, M)); /* m <- Bw */ in KSPSolve_PIPECR()
80 PetscCall(VecDotBegin(W, U, &gamma)); in KSPSolve_PIPECR()
[all …]
/petsc/share/petsc/datafiles/meshes/
H A Dgmsh-3d-binary-32.msh36 ���T���U���V���W���X���Y���Z���[���\���]���^���_���`���^������?�������?^������?������?�������?���…
47 …>���?����>�?�>���?W���>�? �>���?W���>�?�������?�������?�������?���>�?���>�?���>�?X���>…
56 …W���6���W���\������7������Z������8������Z���X���9���X���]������:������Z������;���V���Z����…
60 …j���c�����������d���l���7���~�������y���]���n�����������E���x���q���}�������W���\���D���|�������S�…
H A Dgmsh-3d-binary-64.msh33 �������T�������U�������V�������W�������X�������Y�������Z�������[�������\�������]�������^�������_���…
41 …X���>�?�>���?�������?�������?�������?X���>�? �>���?W���>�?���>�?���>�?���>�? �>���?W�…
50 …W�������9�������W�������_��������������:��������������\��������������;�������X�������\��������…
54 …W�������}���������������U�������^�������D�������{���������������E�������x�������o�������|���������…
H A Dmesh-3d-box-innersphere_bin.msh17 …��������0����RͿ0����R�?U����������0����R�?0����RͿV����������>��<���?�������<W������������������<��…
26 …W������������������������W���T������������������������R���W�������������������…
54 …���U������ ���p������������������V������ ������p���������������W������ ������o����…
89 ����������B����������������������)���=���'���W�������������������������C���
96 …W����������������������V���r���;��������������������������2���E���?������������������…
H A Dsurfacesphere_bin.msh9 …�������07���?7����gm=o�俙�ׁ��?��V����8��������V�w_r�`,�?�Fr46ֿ9����!"J�ݿ�W�%���?�^Q块�:���8'���S…
14 �ݿ��p?�|ӿM����<�-�$пӗ�ͥ9�A�7��xٿN���#��`�ӿ=3�0�L׿qW/��O���q�M��W��ɌC�ʿ1Nfœ�տP���v"ӥ�:׿�y"�N�9^�…
16 …��}n���ǿ�h�\f��zĥt�VǿT���*%����ʿ�ԤϽ�ſG��bN��U���g p�&�׿Y��Wٗ���l!�V���i���E�� F#[$ӿ&����W����…
18 …W�e����n�\f��?4��t�V�?�n���ǿf�����Ͻ��?u��bN��?�����ʿg���“�%�c�?��a���?3����U�h�����Wٗ�?���l!�…
21 …?q�M��W�?�ɌC��?����8^����?v"ӥ�:�?�y"�N�?����H��bN��?,%�����?�ԤϽ��?���� �a���?ٕ���U�? G�%�c�?��…
25 ӿ����J��|���?���pGοb���Q忷���4wG�{�?u�W��e�j��}�?����B�>U�2�?��R��v�6f���?����>��� :�?H/��mٿ�…
35 …I���V���L������������F���������V���I���K������������G������������W���������������H�…
H A Dgmsh-3d-binary.msh245 …���?�������?A������?U���^������?�������?^������?V���������?�������?������?W���`������?�������?`�…
51 …��>�?�>���?�>���?�����>���?����>�?�>���?���� �>���?W���>�?X���>�?����X���>�? �>���?W���…
61 …W���Z���������������<������(������`���W������������=������(���%���_���$������������>�…
65 …*���X���D�������������������d���T���]���y����������������������d���W���>���`���|�������…
/petsc/src/snes/linesearch/impls/bt/
H A Dlinesearchbt.c61 Vec X, F, Y, W, G; in SNESLineSearchApply_BT() local
75 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, &G)); in SNESLineSearchApply_BT()
102 PetscCall(VecCopy(X, W)); in SNESLineSearchApply_BT()
123 PetscCall(MatMult(jac, Y, W)); in SNESLineSearchApply_BT()
124 PetscCall(VecDotRealPart(F, W, &initslope)); in SNESLineSearchApply_BT()
131 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_BT()
132 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_BT()
141 PetscCall(SNESComputeObjective(snes, W, &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()
[all …]
/petsc/src/snes/linesearch/impls/bisection/
H A Dlinesearchbisection.c7 Vec X, F, Y, W, G; in SNESLineSearchApply_Bisection() local
17 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, &G)); in SNESLineSearchApply_Bisection()
43 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Bisection()
44 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Bisection()
45 PetscCall((*linesearch->ops->snesfunc)(snes, W, G)); in SNESLineSearchApply_Bisection()
53 PetscCall((*linesearch->ops->vidirderiv)(snes, G, W, Y, &fty)); in SNESLineSearchApply_Bisection()
144 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Bisection()
145 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Bisection()
146 PetscCall((*linesearch->ops->snesfunc)(snes, W, G)); in SNESLineSearchApply_Bisection()
154 PetscCall((*linesearch->ops->vidirderiv)(snes, G, W, Y, &fty)); in SNESLineSearchApply_Bisection()
[all …]
/petsc/src/ts/tests/
H A Dex13.c21 Vec W, W2, Wdot; in main() local
31 PetscCall(VecCreate(PETSC_COMM_WORLD, &W)); in main()
32 PetscCall(VecSetSizes(W, 1, PETSC_DECIDE)); in main()
33 PetscCall(VecSetUp(W)); in main()
34 PetscCall(VecDuplicate(W, &Wdot)); in main()
35 PetscCall(VecDuplicate(W, &W2)); in main()
59 PetscCall(VecSet(W, func(p, TT[i]))); in main()
61 PetscCall(TSTrajectorySet(tj, ts, II[i], TT[i], W)); in main()
67 PetscCall(TSTrajectoryGetVecs(tj, ts, PETSC_DECIDE, &testtime, W, Wdot)); in main()
68 PetscCall(VecGetArrayRead(W, &aW)); in main()
[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
36 W = ksp->work[9]; in KSPSolve_PIPEBCGS()
78 PetscCall(KSP_MatMult(ksp, pc->mat, R2, W)); /* w <- A r2 */ in KSPSolve_PIPEBCGS()
81 PetscCall(VecDotBegin(W, RP, &d2)); /* d2 <- (w,rp) */ in KSPSolve_PIPEBCGS()
82 PetscCall(PetscCommSplitReductionBegin(PetscObjectComm((PetscObject)W))); in KSPSolve_PIPEBCGS()
83 PetscCall(KSP_PCApply(ksp, W, W2)); /* w2 <- K w */ in KSPSolve_PIPEBCGS()
85 PetscCall(VecDotEnd(W, RP, &d2)); in KSPSolve_PIPEBCGS()
95 PetscCall(VecCopy(W, S)); /* s <- w */ in KSPSolve_PIPEBCGS()
100 …PetscCall(VecAXPBYPCZ(S, 1.0, -beta * omega, beta, W, Z)); /* s <- beta * s + w - beta * ome… in KSPSolve_PIPEBCGS()
106 PetscCall(VecWAXPY(Y, -alpha, Z, W)); /* y <- w - alpha z */ in KSPSolve_PIPEBCGS()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecg/
H A Dpipecg.c25 Vec X, B, Z, P, W, Q, U, M, N, R, S; in KSPSolve_PIPECG() local
39 W = ksp->work[4]; in KSPSolve_PIPECG()
61 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECG()
67 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECG()
73 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECG()
79 PetscCall(KSP_MatMult(ksp, Amat, U, W)); in KSPSolve_PIPECG()
99 PetscCall(VecDotBegin(W, U, &delta)); in KSPSolve_PIPECG()
102 PetscCall(KSP_PCApply(ksp, W, M)); /* m <- Bw */ in KSPSolve_PIPECG()
111 PetscCall(VecDotEnd(W, U, &delta)); in KSPSolve_PIPECG()
129 PetscCall(VecCopy(W, S)); /* s <- w */ in KSPSolve_PIPECG()
[all …]
/petsc/src/ksp/pc/impls/deflation/
H A Ddeflation.c198 static PetscErrorCode PCDeflationSetSpace_Deflation(PC pc, Mat W, PetscBool transpose) in PCDeflationSetSpace_Deflation() argument
204 PetscCall(PetscObjectReference((PetscObject)W)); in PCDeflationSetSpace_Deflation()
207 def->Wt = W; in PCDeflationSetSpace_Deflation()
209 PetscCall(MatDestroy(&def->W)); in PCDeflationSetSpace_Deflation()
210 def->W = W; in PCDeflationSetSpace_Deflation()
238 PetscErrorCode PCDeflationSetSpace(PC pc, Mat W, PetscBool transpose) in PCDeflationSetSpace() argument
242 PetscValidHeaderSpecific(W, MAT_CLASSID, 2); in PCDeflationSetSpace()
244 PetscTryMethod(pc, "PCDeflationSetSpace_C", (PC, Mat, PetscBool), (pc, W, transpose)); in PCDeflationSetSpace()
408 PetscCall(MatMultHermitianTranspose(def->W, r, w1)); /* w1 <- W'*r */ in PCPreSolve_Deflation()
411 PetscCall(MatMult(def->W, w2, r)); /* r <- W*w2 */ in PCPreSolve_Deflation()
[all …]
/petsc/src/snes/impls/gs/
H A Dgssecant.c11 Vec W, G, F; in SNESComputeNGSDefaultSecant() local
28 W = snes->work[0]; in SNESComputeNGSDefaultSecant()
64 PetscCall(VecCopy(X, W)); in SNESComputeNGSDefaultSecant()
65 PetscCall(VecGetArray(W, &wa)); in SNESComputeNGSDefaultSecant()
67 PetscCall(VecRestoreArray(W, &wa)); in SNESComputeNGSDefaultSecant()
69 PetscCall((*func)(snes, W, G, fctx)); in SNESComputeNGSDefaultSecant()
76 PetscCall(VecGetArray(W, &wa)); in SNESComputeNGSDefaultSecant()
99 PetscCall(VecRestoreArray(W, &wa)); in SNESComputeNGSDefaultSecant()
117 PetscCall(VecSwap(X, W)); in SNESComputeNGSDefaultSecant()
/petsc/src/ksp/ksp/impls/cg/pipecgrr/
H A Dpipecgrr.c26 Vec X, B, Z, P, W, Q, U, M, N, R, S; in KSPSolve_PIPECGRR() local
40 W = ksp->work[4]; in KSPSolve_PIPECGRR()
63 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECGRR()
71 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECGRR()
79 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECGRR()
86 PetscCall(KSP_MatMult(ksp, Amat, U, W)); in KSPSolve_PIPECGRR()
125 PetscCall(VecDotBegin(W, U, &delta)); in KSPSolve_PIPECGRR()
136 PetscCall(VecNormBegin(W, NORM_2, &dw)); in KSPSolve_PIPECGRR()
139 PetscCall(KSP_PCApply(ksp, W, M)); /* m <- Bw */ in KSPSolve_PIPECGRR()
148 PetscCall(VecDotEnd(W, U, &delta)); in KSPSolve_PIPECGRR()
[all …]
/petsc/src/snes/linesearch/impls/nleqerr/
H A Dlinesearchnleqerr.c36 Vec X, F, Y, W, G; in SNESLineSearchApply_NLEQERR() local
49 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, &G)); in SNESLineSearchApply_NLEQERR()
78 PetscCall(VecCopy(X, W)); 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()
164 PetscCall(VecWAXPY(W, -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()
174 PetscCall(VecNorm(W, NORM_2, &wnorm)); in SNESLineSearchApply_NLEQERR()
185 PetscCall(VecWAXPY(G, -(1.0 - lambda), Y, W)); in SNESLineSearchApply_NLEQERR()
209 PetscCall(VecAXPY(G, -1.0, W)); in SNESLineSearchApply_NLEQERR()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipeprcg/
H A Dpipeprcg.c44 Vec X, B, R, RT, W, WT, P, S, ST, U, UT, PRTST[3]; in KSPSolve_PIPEPRCG() local
61 W = ksp->work[2]; in KSPSolve_PIPEPRCG()
81 PetscCall(KSP_MatMult(ksp, Amat, RT, W)); /* w <- A rt */ in KSPSolve_PIPEPRCG()
82 PetscCall(KSP_PCApply(ksp, W, WT)); /* wt <- B w */ in KSPSolve_PIPEPRCG()
85 PetscCall(VecCopy(W, S)); /* p <- rt */ in KSPSolve_PIPEPRCG()
140 PetscCall(VecAXPY(W, -alpha, U)); /* w <- w - alpha * u */ in KSPSolve_PIPEPRCG()
143 PetscCall(VecAYPX(S, beta, W)); /* s <- w + beta * s */ in KSPSolve_PIPEPRCG()
162 PetscCall(KSP_MatMult(ksp, Amat, RT, W)); /* w <- A rt */ in KSPSolve_PIPEPRCG()
163 PetscCall(KSP_PCApply(ksp, W, WT)); /* wt <- B w */ in KSPSolve_PIPEPRCG()
/petsc/src/snes/impls/ntrdc/
H A Dntrdc.c274 static PetscErrorCode SNESNewtonTRDCPostCheck(SNES snes, Vec X, Vec Y, Vec W, PetscBool *changed_Y,… in SNESNewtonTRDCPostCheck() argument
282 PetscCall((*tr->postcheck)(snes, X, Y, W, changed_Y, changed_W, tr->postcheckctx)); in SNESNewtonTRDCPostCheck()
298 Vec X, F, Y, G, W, GradF, YNtmp; in SNESSolve_NEWTONTRDC() local
320 W = snes->work[2]; /* temporary vector */ in SNESSolve_NEWTONTRDC()
389 PetscCall(VecStrideSet(W, j, inorms[j])); in SNESSolve_NEWTONTRDC()
400 PetscCall(MatDiagonalScale(jac, NULL, W)); in SNESSolve_NEWTONTRDC()
413 PetscCall(MatMult(jac, GradF, W)); in SNESSolve_NEWTONTRDC()
414 PetscCall(VecDotRealPart(W, W, &gTBg)); /* completes GradF^T J^T J GradF */ in SNESSolve_NEWTONTRDC()
440 PetscCall(VecWAXPY(W, tau, YNtmp, YCtmp)); in SNESSolve_NEWTONTRDC()
441 PetscCall(VecAXPY(W, -tau, YCtmp)); in SNESSolve_NEWTONTRDC()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/rescale/
H A Dsymbrdnrescale.c65 static PetscErrorCode DiagonalUpdate(SymBroydenRescale ldb, Vec D, Vec s, Vec y, Vec V, Vec W, Vec … in DiagonalUpdate() argument
74 PetscCall(VecPointwiseMult(W, D, s)); in DiagonalUpdate()
75 PetscCall(VecDotRealPart(W, s, &stDs)); in DiagonalUpdate()
86 PetscCall(VecPointwiseMult(BFGS, W, W)); in DiagonalUpdate()
98 PetscCall(VecPointwiseMult(DFP, DFP, W)); in DiagonalUpdate()
100 PetscCall(VecCopy(DFP, W)); in DiagonalUpdate()
101 PetscCall(VecConjugate(W)); in DiagonalUpdate()
102 PetscCall(VecAXPY(DFP, 1.0, W)); in DiagonalUpdate()
145 …PetscCall(DiagonalUpdate(ldb, ldb->invDnew, s_last, y_last, ldb->V, ldb->W, ldb->BFGS, ldb->DFP, l… in SymBroydenRescaleUpdateDiagonal()
151 …PetscCall(DiagonalUpdate(ldb, ldb->invDnew, y_last, s_last, ldb->V, ldb->W, ldb->DFP, ldb->BFGS, 1… in SymBroydenRescaleUpdateDiagonal()
[all …]

123456