| /petsc/share/petsc/datafiles/meshes/ |
| H A D | hybrid_hexwedge.msh | 7 …Y���������@�@��������������Y�������(\��@�������?������Y���������ʥ@��������������Y���������@�@�…
8 ���������@�(\��?�@������Y� ��������ǥ@�z�GA�@������Y� ��������@�@�����@�@������Y�
|
| /petsc/src/mat/utils/ |
| H A D | axpy.c | 3 static PetscErrorCode MatTransposeAXPY_Private(Mat Y, PetscScalar a, Mat X, MatStructure str, Mat T) in MatTransposeAXPY_Private() argument
21 A = Y; in MatTransposeAXPY_Private()
31 A = Y; in MatTransposeAXPY_Private()
43 static PetscErrorCode MatAXPY_BasicWithTypeCompare(Mat Y, PetscScalar a, Mat X, MatStructure str) in MatAXPY_BasicWithTypeCompare() argument
48 PetscCall(MatIsShell(Y, &flg)); in MatAXPY_BasicWithTypeCompare()
52 PetscCall(MatGetOperation(Y, MATOP_AXPY, (PetscErrorCodeFn **)&f)); in MatAXPY_BasicWithTypeCompare()
54 PetscCall((*f)(Y, a, X, str)); in MatAXPY_BasicWithTypeCompare()
59 PetscCall(PetscObjectBaseTypeCompareAny((PetscObject)Y, &flg, MATSEQDENSE, MATMPIDENSE, "")); in MatAXPY_BasicWithTypeCompare()
63 PetscCall(MatAXPY_Dense_Nest(Y, a, X)); in MatAXPY_BasicWithTypeCompare()
71 PetscCall(MatConvert(X, ((PetscObject)Y)->type_name, MAT_INITIAL_MATRIX, &C)); in MatAXPY_BasicWithTypeCompare()
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex31.c | 71 PetscErrorCode RHSFunction_Hull1972A1(TS ts, PetscReal t, Vec Y, Vec F, void *s) in RHSFunction_Hull1972A1() argument
77 PetscCall(VecGetArrayRead(Y, &y)); in RHSFunction_Hull1972A1()
80 PetscCall(VecRestoreArrayRead(Y, &y)); in RHSFunction_Hull1972A1()
85 PetscErrorCode RHSJacobian_Hull1972A1(TS ts, PetscReal t, Vec Y, Mat A, Mat B, void *s) in RHSJacobian_Hull1972A1() argument
92 PetscCall(VecGetArrayRead(Y, &y)); in RHSJacobian_Hull1972A1()
96 PetscCall(VecRestoreArrayRead(Y, &y)); in RHSJacobian_Hull1972A1()
100 PetscErrorCode IFunction_Hull1972A1(TS ts, PetscReal t, Vec Y, Vec Ydot, Vec F, void *s) in IFunction_Hull1972A1() argument
106 PetscCall(VecGetArrayRead(Y, &y)); in IFunction_Hull1972A1()
109 PetscCall(VecRestoreArrayRead(Y, &y)); in IFunction_Hull1972A1()
116 PetscErrorCode IJacobian_Hull1972A1(TS ts, PetscReal t, Vec Y, Vec Ydot, PetscReal a, Mat A, Mat B,… in IJacobian_Hull1972A1() argument
[all …]
|
| H A D | ex36.c | 47 static PetscErrorCode IFunctionImplicit(TS ts, PetscReal t, Vec Y, Vec Ydot, Vec F, PetscCtx ctx) in IFunctionImplicit() argument
54 PetscCall(VecGetArrayRead(Y, &y)); in IFunctionImplicit()
64 PetscCall(VecRestoreArrayRead(Y, &y)); in IFunctionImplicit()
73 static PetscErrorCode IJacobianImplicit(TS ts, PetscReal t, Vec Y, Vec Ydot, PetscReal a, Mat A, Ma… in IJacobianImplicit() argument
80 PetscCall(VecGetArrayRead(Y, &y)); in IJacobianImplicit()
101 PetscCall(VecRestoreArrayRead(Y, &y)); in IJacobianImplicit()
116 Vec Y; /* solution will be stored here */ in main() local
138 PetscCall(MatCreateVecs(A, &Y, NULL)); in main()
140 PetscCall(VecGetArray(Y, &y)); in main()
146 PetscCall(VecRestoreArray(Y, &y)); in main()
[all …]
|
| /petsc/src/vec/vec/tests/ |
| H A D | ex44.c | 13 Vec X, Y; in main() local
61 PetscCall(VecCreate(PETSC_COMM_WORLD, &Y)); in main()
62 PetscCall(VecSetSizes(Y, m, PETSC_DECIDE)); in main()
66 PetscCall(VecSetFromOptions(Y)); in main()
69 PetscCall(VecSet(Y, 1.0)); in main()
99 PetscCall(VecScatterCreate(X, fromISStrided, Y, toISStrided, &vscatSStoSS)); in main()
100 PetscCall(VecScatterBegin(vscatSStoSS, X, Y, addv, mode)); in main()
101 PetscCall(VecScatterEnd(vscatSStoSS, X, Y, addv, mode)); in main()
105 PetscCall(VecScatterCreate(Y, fromISGeneral, X, toISStrided, &vscatSGtoSS)); in main()
106 PetscCall(VecScatterBegin(vscatSGtoSS, Y, X, addv, mode)); in main()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex97.c | 27 Vec Y; in Compare2() local
31 PetscCall(VecDuplicate(X[0], &Y)); in Compare2()
32 PetscCall(VecCopy(X[0], Y)); in Compare2()
33 PetscCall(VecAYPX(Y, -1.0, X[1])); in Compare2()
34 PetscCall(VecNorm(Y, NORM_INFINITY, &norm)); in Compare2()
45 PetscCall(VecView(Y, PETSC_VIEWER_STDOUT_WORLD)); in Compare2()
47 PetscCall(VecDestroy(&Y)); in Compare2()
51 static PetscErrorCode CheckMatrices(Mat A, Mat B, Vec left, Vec right, Vec X, Vec Y, Vec X1, Vec Y1) in CheckMatrices() argument
67 PetscCall(MatMultTranspose(A, Y, rtmp[0])); in CheckMatrices()
68 PetscCall(MatMultTranspose(B, Y, rtmp[1])); in CheckMatrices()
[all …]
|
| H A D | ex205.c | 10 static PetscErrorCode MatMult_User(Mat A, Vec X, Vec Y) in MatMult_User() argument
16 PetscCall(MatMult(user->A, X, Y)); in MatMult_User()
53 Vec X, Y; in main() local
69 PetscCall(VecDuplicate(X, &Y)); in main()
70 PetscCall(VecSetValues(Y, 2, inds, yvals, INSERT_VALUES)); in main()
71 PetscCall(VecAssemblyBegin(Y)); in main()
72 PetscCall(VecAssemblyEnd(Y)); in main()
87 PetscCall(MatDiagonalScale(S1, X, Y)); in main()
89 PetscCall(MatMult(S1, X, Y)); in main()
90 PetscCall(VecView(Y, PETSC_VIEWER_STDOUT_WORLD)); in main()
[all …]
|
| H A D | ex99.c | 27 Vec Y; in Compare2() local
31 PetscCall(VecDuplicate(X[0], &Y)); in Compare2()
32 PetscCall(VecCopy(X[0], Y)); in Compare2()
33 PetscCall(VecAYPX(Y, -1.0, X[1])); in Compare2()
34 PetscCall(VecNorm(Y, NORM_INFINITY, &norm)); in Compare2()
45 PetscCall(VecView(Y, PETSC_VIEWER_STDOUT_WORLD)); in Compare2()
47 PetscCall(VecDestroy(&Y)); in Compare2()
51 static PetscErrorCode CheckMatrices(Mat A, Mat B, Vec left, Vec right, Vec X, Vec Y, Vec X1, Vec Y1) in CheckMatrices() argument
69 PetscCall(MatMultTranspose(A, Y, rtmp[0])); in CheckMatrices()
70 PetscCall(MatMultTranspose(B, Y, rtmp[1])); in CheckMatrices()
[all …]
|
| H A D | ex203.c | 26 Vec X, Y; in main() local
40 PetscCall(VecDuplicate(X, &Y)); in main()
50 PetscCall(MatGetDiagonal(S, Y)); in main()
51 PetscCall(VecView(Y, PETSC_VIEWER_STDOUT_WORLD)); in main()
53 PetscCall(MatGetDiagonal(S, Y)); in main()
54 PetscCall(VecView(Y, PETSC_VIEWER_STDOUT_WORLD)); in main()
56 PetscCall(MatGetDiagonal(S, Y)); in main()
57 PetscCall(VecView(Y, PETSC_VIEWER_STDOUT_WORLD)); in main()
62 PetscCall(VecDestroy(&Y)); in main()
|
| H A D | ex185.c | 8 Vec X, Y; in main() local
17 PetscCall(MatCreateVecs(A, &X, &Y)); in main()
22 PetscCall(MatMult(A, X, Y)); in main()
23 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in main()
46 PetscCall(MatGetDiagonal(A, Y)); in main()
50 PetscCall(MatSolve(Af, X, Y)); in main()
51 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in main()
55 PetscCall(MatSolve(A, X, Y)); in main()
56 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in main()
60 PetscCall(MatSOR(A, X, 1.5, SOR_FORWARD_SWEEP, 0.0, 1, 1, Y)); in main()
[all …]
|
| H A D | ex154.c | 7 Mat F, A, B, X, Y, S; in main() local
90 PetscCall(MatMatMult(A, X, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &Y)); in main()
91 PetscCall(MatViewFromOptions(Y, NULL, "-Y_view")); in main()
92 PetscCall(MatAXPY(Y, -1.0, B, SAME_NONZERO_PATTERN)); in main()
93 PetscCall(MatViewFromOptions(Y, NULL, "-err_view")); in main()
94 PetscCall(MatNorm(Y, NORM_FROBENIUS, &norm)); in main()
97 PetscCall(MatConvert(Y, MATAIJ, MAT_INPLACE_MATRIX, &Y)); in main()
98 PetscCall(MatFilter(Y, PETSC_SMALL, PETSC_TRUE, PETSC_FALSE)); in main()
99 PetscCall(MatViewFromOptions(Y, NULL, "-aij_err_view")); in main()
105 PetscCall(MatDestroy(&Y)); in main()
|
| /petsc/src/tao/constrained/impls/almm/ |
| H A D | almmutils.c | 166 PetscErrorCode TaoALMMGetMultipliers(Tao tao, Vec *Y) in TaoALMMGetMultipliers() argument
170 PetscAssertPointer(Y, 2); in TaoALMMGetMultipliers()
171 PetscUseMethod(tao, "TaoALMMGetMultipliers_C", (Tao, Vec *), (tao, Y)); in TaoALMMGetMultipliers()
175 PetscErrorCode TaoALMMGetMultipliers_Private(Tao tao, Vec *Y) in TaoALMMGetMultipliers_Private() argument
181 *Y = auglag->Y; in TaoALMMGetMultipliers_Private()
208 PetscErrorCode TaoALMMSetMultipliers(Tao tao, Vec Y) in TaoALMMSetMultipliers() argument
212 PetscValidHeaderSpecific(Y, VEC_CLASSID, 2); in TaoALMMSetMultipliers()
213 PetscTryMethod(tao, "TaoALMMSetMultipliers_C", (Tao, Vec), (tao, Y)); in TaoALMMSetMultipliers()
217 PetscErrorCode TaoALMMSetMultipliers_Private(Tao tao, Vec Y) in TaoALMMSetMultipliers_Private() argument
226 if (Y == auglag->Y) PetscFunctionReturn(PETSC_SUCCESS); in TaoALMMSetMultipliers_Private()
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbadbrdn.c | 3 static PetscErrorCode MatMult_LMVMSymBadBrdn_Recursive(Mat B, Vec X, Vec Y) in MatMult_LMVMSymBadBrdn_Recursive() argument
11 PetscCall(DFPKernel_Recursive(B, MATLMVM_MODE_PRIMAL, X, Y)); in MatMult_LMVMSymBadBrdn_Recursive()
13 PetscCall(BFGSKernel_Recursive(B, MATLMVM_MODE_PRIMAL, X, Y)); in MatMult_LMVMSymBadBrdn_Recursive()
15 PetscCall(SymBroydenKernel_Recursive(B, MATLMVM_MODE_PRIMAL, X, Y, PETSC_TRUE)); in MatMult_LMVMSymBadBrdn_Recursive()
20 static PetscErrorCode MatSolve_LMVMSymBadBrdn_Recursive(Mat B, Vec X, Vec Y) in MatSolve_LMVMSymBadBrdn_Recursive() argument
28 PetscCall(BFGSKernel_Recursive(B, MATLMVM_MODE_DUAL, X, Y)); in MatSolve_LMVMSymBadBrdn_Recursive()
30 PetscCall(DFPKernel_Recursive(B, MATLMVM_MODE_DUAL, X, Y)); in MatSolve_LMVMSymBadBrdn_Recursive()
32 PetscCall(SymBroydenKernel_Recursive(B, MATLMVM_MODE_DUAL, X, Y, PETSC_FALSE)); in MatSolve_LMVMSymBadBrdn_Recursive()
37 static PetscErrorCode MatMult_LMVMSymBadBrdn_CompactDense(Mat B, Vec X, Vec Y) in MatMult_LMVMSymBadBrdn_CompactDense() argument
45 PetscCall(DFPKernel_CompactDense(B, MATLMVM_MODE_PRIMAL, X, Y)); in MatMult_LMVMSymBadBrdn_CompactDense()
[all …]
|
| /petsc/src/dm/impls/patch/tests/output/ |
| H A D | ex1_3.out | 4 X range of indices: 0 2, Y range of indices: 0 2
6 X range of indices: 2 4, Y range of indices: 0 2
8 X range of indices: 0 2, Y range of indices: 2 4
10 X range of indices: 2 4, Y range of indices: 2 4
13 X range of indices: 0 2, Y range of indices: 0 2
15 X range of indices: 2 3, Y range of indices: 0 2
17 X range of indices: 0 2, Y range of indices: 2 3
19 X range of indices: 2 3, Y range of indices: 2 3
47 X range of indices: 0 2, Y range of indices: 0 2
49 X range of indices: 2 3, Y range of indices: 0 2
[all …]
|
| /petsc/src/ts/adapt/impls/glee/ |
| H A D | adaptglee.c | 5 Vec Y; member
11 Vec X, Y, E; in TSAdaptChoose_GLEE() local
26 if (!glee->Y && adapt->glee_use_local) { in TSAdaptChoose_GLEE()
27 PetscCall(VecDuplicate(X, &glee->Y)); /*create vector to store previous step global error*/ in TSAdaptChoose_GLEE()
28 …PetscCall(VecZeroEntries(glee->Y)); /*set error to zero on the first step - may not work if erro… in TSAdaptChoose_GLEE()
34 …if (adapt->glee_use_local) PetscCall(VecAXPY(E, -1.0, glee->Y)); /* local error = current error - … in TSAdaptChoose_GLEE()
42 if (!glee->Y) PetscCall(VecDuplicate(X, &glee->Y)); in TSAdaptChoose_GLEE()
43 Y = glee->Y; in TSAdaptChoose_GLEE()
44 PetscCall(TSEvaluateStep(ts, order - 1, Y, NULL)); in TSAdaptChoose_GLEE()
45 PetscCall(TSErrorWeightedNorm(ts, X, Y, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_GLEE()
[all …]
|
| /petsc/src/dm/tests/output/ |
| H A D | ex34_1.out | 4 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 4
6 X range of indices: 4 8, Y range of indices: 0 3, Z range of indices: 0 4
8 X range of indices: 8 11, Y range of indices: 0 3, Z range of indices: 0 4
10 X range of indices: 0 4, Y range of indices: 3 5, Z range of indices: 0 4
12 X range of indices: 4 8, Y range of indices: 3 5, Z range of indices: 0 4
14 X range of indices: 8 11, Y range of indices: 3 5, Z range of indices: 0 4
16 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 4 7
18 X range of indices: 4 8, Y range of indices: 0 3, Z range of indices: 4 7
20 X range of indices: 8 11, Y range of indices: 0 3, Z range of indices: 4 7
22 X range of indices: 0 4, Y range of indices: 3 5, Z range of indices: 4 7
[all …]
|
| /petsc/src/snes/linesearch/impls/cp/ |
| H A D | linesearchcp.c | 7 Vec X, Y, F, W; in SNESLineSearchApply_CP() local
16 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, NULL)); in SNESLineSearchApply_CP()
25 PetscCall(SNESLineSearchPreCheck(linesearch, X, Y, &changed_y)); in SNESLineSearchApply_CP()
30 PetscCall((*linesearch->ops->vidirderiv)(snes, F, X, Y, &fty_old)); in SNESLineSearchApply_CP()
32 PetscCall(VecDot(F, Y, &fty_old)); in SNESLineSearchApply_CP()
47 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_CP()
51 PetscCall((*linesearch->ops->vidirderiv)(snes, F, W, Y, &fty)); in SNESLineSearchApply_CP()
53 PetscCall(VecDot(F, Y, &fty)); in SNESLineSearchApply_CP()
102 PetscCall(VecWAXPY(W, -0.5 * (lambda + lambda_old), Y, X)); in SNESLineSearchApply_CP()
106 PetscCall((*linesearch->ops->vidirderiv)(snes, F, W, Y, &fty_mid1)); in SNESLineSearchApply_CP()
[all …]
|
| /petsc/src/dm/impls/stag/tests/output/ |
| H A D | ex50_dim-3_coords-false.out | 4 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 3
118 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 3
160 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 3
202 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
280 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
310 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
340 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
370 X range of indices: 0 4, Y range of indices: 0 2, Z range of indices: 0 3
448 X range of indices: 0 4, Y range of indices: 0 2, Z range of indices: 0 3
478 X range of indices: 0 4, Y range of indices: 0 2, Z range of indices: 0 3
[all …]
|
| H A D | ex50_dim-2_coords-false.out | 4 X range of indices: 0 4, Y range of indices: 0 3
46 X range of indices: 0 4, Y range of indices: 0 3
64 X range of indices: 0 4, Y range of indices: 0 3
82 X range of indices: 0 4, Y range of indices: 0 2
112 X range of indices: 0 4, Y range of indices: 0 2
126 X range of indices: 0 3, Y range of indices: 0 3
159 X range of indices: 0 3, Y range of indices: 0 3
174 X range of indices: 0 3, Y range of indices: 0 2
198 X range of indices: 0 3, Y range of indices: 0 2
210 X range of indices: 0 3, Y range of indices: 0 2
[all …]
|
| H A D | ex50_dim-2_coords-true.out | 4 X range of indices: 0 4, Y range of indices: 0 3
46 X range of indices: 0 4, Y range of indices: 0 3
64 X range of indices: 0 4, Y range of indices: 0 3
82 X range of indices: 0 4, Y range of indices: 0 2
112 X range of indices: 0 4, Y range of indices: 0 2
126 X range of indices: 0 3, Y range of indices: 0 3
159 X range of indices: 0 3, Y range of indices: 0 3
174 X range of indices: 0 3, Y range of indices: 0 2
198 X range of indices: 0 3, Y range of indices: 0 2
210 X range of indices: 0 3, Y range of indices: 0 2
[all …]
|
| H A D | ex50_dim-3_coords-true.out | 4 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 3
119 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 3
162 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 3
205 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
284 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
315 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
346 X range of indices: 0 4, Y range of indices: 0 3, Z range of indices: 0 2
377 X range of indices: 0 4, Y range of indices: 0 2, Z range of indices: 0 3
456 X range of indices: 0 4, Y range of indices: 0 2, Z range of indices: 0 3
487 X range of indices: 0 4, Y range of indices: 0 2, Z range of indices: 0 3
[all …]
|
| /petsc/src/snes/linesearch/impls/basic/ |
| H A D | linesearchbasic.c | 7 Vec X, F, Y, W; in SNESLineSearchApply_Basic() local
12 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, NULL)); in SNESLineSearchApply_Basic()
19 PetscCall(SNESLineSearchPreCheck(linesearch, X, Y, &changed_y)); in SNESLineSearchApply_Basic()
22 PetscCall(VecWAXPY(W, -lambda, Y, X)); 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()
37 PetscCall(VecNormBegin(Y, NORM_2, &linesearch->ynorm)); in SNESLineSearchApply_Basic()
39 PetscCall(VecNormEnd(Y, NORM_2, &linesearch->ynorm)); in SNESLineSearchApply_Basic()
|
| /petsc/src/vec/vec/tutorials/output/ |
| H A D | ex44_serial.out | 22 Expected Y vector
38 Testing VecConcatenate() for Y = [X[1], X[2], ...]
55 Testing VecConcatenate() for Y = [X[1], X[2], ...] w/o IS
97 Testing VecScatter for Y -> X[0]
103 Testing VecScatter for Y -> X[1]
111 Testing VecScatter for Y -> X[2]
122 Testing VecScatter for X[:] -> Y
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | badbrdn.c | 66 LMBasis Y; in BadBroydenKernel_Recursive_Inner() local
70 PetscCall(MatLMVMGetUpdatedBasis(B, Y_t, &Y, NULL, NULL)); in BadBroydenKernel_Recursive_Inner()
77 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in BadBroydenKernel_Recursive_Inner()
79 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in BadBroydenKernel_Recursive_Inner()
104 LMBasis Y; in BadBroydenRecursiveBasisUpdate() local
112 PetscCall(MatLMVMGetUpdatedBasis(B, Y_t, &Y, NULL, NULL)); in BadBroydenRecursiveBasisUpdate()
137 PetscCall(LMBasisGetVecRead(Y, j, &y_j)); in BadBroydenRecursiveBasisUpdate()
141 PetscCall(LMBasisRestoreVecRead(Y, j, &y_j)); in BadBroydenRecursiveBasisUpdate()
148 PETSC_INTERN PetscErrorCode BadBroydenKernel_Recursive(Mat B, MatLMVMMode mode, Vec X, Vec Y) in BadBroydenKernel_Recursive() argument
153 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, Y)); in BadBroydenKernel_Recursive()
[all …]
|
| /petsc/src/snes/linesearch/impls/bisection/ |
| H A D | linesearchbisection.c | 7 Vec X, F, Y, W, G; in SNESLineSearchApply_Bisection() local
17 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, &G)); in SNESLineSearchApply_Bisection()
25 PetscCall(SNESLineSearchPreCheck(linesearch, X, Y, &changed_y)); in SNESLineSearchApply_Bisection()
28 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in SNESLineSearchApply_Bisection()
36 PetscCall((*linesearch->ops->vidirderiv)(snes, F, X, Y, &fty_left)); in SNESLineSearchApply_Bisection()
38 PetscCall(VecDot(F, Y, &fty_left)); in SNESLineSearchApply_Bisection()
43 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Bisection()
53 PetscCall((*linesearch->ops->vidirderiv)(snes, G, W, Y, &fty)); in SNESLineSearchApply_Bisection()
55 PetscCall(VecDot(G, Y, &fty)); in SNESLineSearchApply_Bisection()
144 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Bisection()
[all …]
|