Home
last modified time | relevance | path

Searched refs:R (Results 1 – 25 of 172) sorted by relevance

1234567

/petsc/src/snes/tutorials/output/
H A Dex20_2d_p1_gmg_fcycle_rate.out1 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.698164 R^2 0.5911…
2 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.69301 R^2 0.541848
3 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.73243 R^2 0.805662
4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.696497 R^2 0.7674…
5 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.700915 R^2 0.6791…
6 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.760929 R^2 0.5625…
7 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.728701 R^2 0.8390…
8 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.742781 R^2 0.6120…
9 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.694817 R^2 0.7814…
10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.700569 R^2 0.7003…
[all …]
H A Dex20_2d_p1_gmg_vcycle_rate.out1 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.774356 R^2 0.9846…
2 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.613993 R^2 0.9107…
3 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.689662 R^2 0.5206…
4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.688226 R^2 0.5129…
5 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.729655 R^2 0.8322…
6 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.712464 R^2 0.8676…
7 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.773593 R^2 0.9848…
8 … Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.6147 R^2 0.91186
9 … Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.6897 R^2 0.520871
10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.688241 R^2 0.5129…
[all …]
H A Dex20_2d_p1_gmg_vcycle_adapt_rate.out1 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.774356 R^2 0.9846…
2 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.719169 R^2 0.6944…
3 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.765723 R^2 0.6269…
4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.767943 R^2 0.5237…
5 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.729265 R^2 0.84772
6 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.72737 R^2 0.943032
7 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.770778 R^2 0.9834…
8 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.718876 R^2 0.7020…
9 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.76585 R^2 0.625507
10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.768003 R^2 0.5232…
[all …]
H A Dex20_2d_p1_gmg_vcycle_cr.out1 … mg_levels_1_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.617436 R^2 0.947533
2 … mg_levels_2_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.604551 R^2 0.962706
3 …Linear mg_levels_3_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.625511 R^2 0…
4 … mg_levels_1_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.617436 R^2 0.947533
5 … mg_levels_2_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.604551 R^2 0.962706
6 …Linear mg_levels_3_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.625511 R^2 0…
7 … mg_levels_1_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.617436 R^2 0.947533
8 … mg_levels_2_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.604551 R^2 0.962706
9 …Linear mg_levels_3_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.625511 R^2 0…
10 … mg_levels_1_cr_ solve converged due to CONVERGED_ITS iterations 5 error rate 0.617436 R^2 0.947533
[all …]
H A Dex20_2d_p1_scalable_rate.out1 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.785382 R^2 0.9647…
2 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.648211 R^2 0.9771…
3 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.736109 R^2 0.9663…
4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.656634 R^2 0.9768…
5 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.764618 R^2 0.9750…
6 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.644732 R^2 0.9831…
7 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.776005 R^2 0.9734…
8 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.636778 R^2 0.9840…
9 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.75345 R^2 0.943944
10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.647287 R^2 0.9823…
[all …]
/petsc/src/vec/vec/tests/output/
H A Dex37_5.out43 Testing Wnorms 2: E? 0 A? 0 R? 0
44 Testing Wnorms 2: E? 0 A? 0 R? 1
45 Testing Wnorms 2: E? 0 A? 1 R? 0
46 Testing Wnorms 2: E? 0 A? 1 R? 1
47 Testing Wnorms 2: E? 1 A? 0 R? 0
48 Testing Wnorms 2: E? 1 A? 0 R? 1
49 Testing Wnorms 2: E? 1 A? 1 R? 0
50 Testing Wnorms 2: E? 1 A? 1 R? 1
51 Testing Wnorms inf: E? 0 A? 0 R? 0
52 Testing Wnorms inf: E? 0 A? 0 R? 1
[all …]
H A Dex37_4.out43 Testing Wnorms 2: E? 0 A? 0 R? 0
44 Testing Wnorms 2: E? 0 A? 0 R? 1
45 Testing Wnorms 2: E? 0 A? 1 R? 0
46 Testing Wnorms 2: E? 0 A? 1 R? 1
47 Testing Wnorms 2: E? 1 A? 0 R? 0
48 Testing Wnorms 2: E? 1 A? 0 R? 1
49 Testing Wnorms 2: E? 1 A? 1 R? 0
50 Testing Wnorms 2: E? 1 A? 1 R? 1
51 Testing Wnorms inf: E? 0 A? 0 R? 0
52 Testing Wnorms inf: E? 0 A? 0 R? 1
[all …]
H A Dex37_1.out43 Testing Wnorms 2: E? 0 A? 0 R? 0
44 Testing Wnorms 2: E? 0 A? 0 R? 1
45 Testing Wnorms 2: E? 0 A? 1 R? 0
46 Testing Wnorms 2: E? 0 A? 1 R? 1
47 Testing Wnorms 2: E? 1 A? 0 R? 0
48 Testing Wnorms 2: E? 1 A? 0 R? 1
49 Testing Wnorms 2: E? 1 A? 1 R? 0
50 Testing Wnorms 2: E? 1 A? 1 R? 1
51 Testing Wnorms inf: E? 0 A? 0 R? 0
52 Testing Wnorms inf: E? 0 A? 0 R? 1
[all …]
H A Dex37_3.out51 Testing Wnorms 2: E? 0 A? 0 R? 0
52 Testing Wnorms 2: E? 0 A? 0 R? 1
53 Testing Wnorms 2: E? 0 A? 1 R? 0
54 Testing Wnorms 2: E? 0 A? 1 R? 1
55 Testing Wnorms 2: E? 1 A? 0 R? 0
56 Testing Wnorms 2: E? 1 A? 0 R? 1
57 Testing Wnorms 2: E? 1 A? 1 R? 0
58 Testing Wnorms 2: E? 1 A? 1 R? 1
59 Testing Wnorms inf: E? 0 A? 0 R? 0
60 Testing Wnorms inf: E? 0 A? 0 R? 1
[all …]
/petsc/src/mat/tests/
H A Dex264.c22 Mat block[4], H, R, C, M; in main() local
37 PetscCall(MatCreate(PETSC_COMM_WORLD, &R)); in main()
38 PetscCall(MatSetSizes(R, PETSC_DECIDE, PETSC_DECIDE, n, n)); in main()
39 PetscCall(MatSetFromOptions(R)); in main()
45 PetscCall(MatGetOwnershipRange(R, &Istart, &Iend)); in main()
47 if (i > 1) PetscCall(MatSetValue(R, i, i - 2, a, INSERT_VALUES)); in main()
48 if (i > 0) PetscCall(MatSetValue(R, i, i - 1, b, INSERT_VALUES)); in main()
49 PetscCall(MatSetValue(R, i, i, c, INSERT_VALUES)); in main()
50 if (i < n - 1) PetscCall(MatSetValue(R, i, i + 1, PetscConj(b), INSERT_VALUES)); in main()
51 if (i < n - 2) PetscCall(MatSetValue(R, i, i + 2, PetscConj(a), INSERT_VALUES)); in main()
[all …]
H A Dex161.c8 Mat A, R, C, C_dense, C_sparse, Rt_dense, P, PtAP; in main() local
51 PetscCall(MatCreate(PETSC_COMM_SELF, &R)); in main()
52 PetscCall(MatSetSizes(R, 2, 4, 2, 4)); in main()
53 PetscCall(MatSetType(R, MATSEQAIJ)); in main()
54 PetscCall(MatSetFromOptions(R)); in main()
55 PetscCall(MatSetUp(R)); in main()
58 PetscCall(MatSetValues(R, 1, &row, 1, &col, &one, ADD_VALUES)); in main()
61 PetscCall(MatSetValues(R, 1, &row, 1, &col, &one, ADD_VALUES)); in main()
65 PetscCall(MatSetValues(R, 1, &row, 1, &col, &one, ADD_VALUES)); in main()
68 PetscCall(MatSetValues(R, 1, &row, 1, &col, &one, ADD_VALUES)); in main()
[all …]
/petsc/src/ksp/ksp/tests/
H A Dex87.c24 Mat H, R, C, block[4]; in main() local
56 PetscCall(MatCreate(PETSC_COMM_WORLD, &R)); in main()
57 PetscCall(MatSetSizes(R, PETSC_DECIDE, PETSC_DECIDE, n, n)); in main()
58 PetscCall(MatSetFromOptions(R)); in main()
64 PetscCall(MatGetOwnershipRange(R, &Istart, &Iend)); in main()
66 if (i > 1) PetscCall(MatSetValue(R, i, i - 2, a, INSERT_VALUES)); in main()
67 if (i > 0) PetscCall(MatSetValue(R, i, i - 1, b, INSERT_VALUES)); in main()
68 PetscCall(MatSetValue(R, i, i, c, INSERT_VALUES)); in main()
69 if (i < n - 1) PetscCall(MatSetValue(R, i, i + 1, PetscConj(b), INSERT_VALUES)); in main()
70 if (i < n - 2) PetscCall(MatSetValue(R, i, i + 2, PetscConj(a), INSERT_VALUES)); in main()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dmatrart.c25 PetscErrorCode MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart(Mat A, Mat R, PetscReal fill, Mat C) in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart() argument
38 PetscCall(MatTransposeSymbolic(R, &P)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
42 PetscCall(MatSetBlockSizes(C, R->rmap->bs, R->rmap->bs)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
94 …t_FMT ")/(m*ncolors)=%g\n", RARt_dense->rmap->n, RARt_dense->cmap->n, R->cmap->n, R->rmap->n, c->n… in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
103 static PetscErrorCode MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense(Mat R, Mat A, Mat B, Mat RAB, Pet… in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense() argument
105 Mat_SeqAIJ *a = (Mat_SeqAIJ *)A->data, *r = (Mat_SeqAIJ *)R->data; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
112 PetscInt *rj, rm = R->rmap->n, dm = RAB->rmap->n, dn = RAB->cmap->n; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
118 …scCheck(am == R->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Number columns in R %" PetscInt_FMT… in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
119 …R->rmap->n == RAB->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Number rows in RAB %" PetscInt_FM… in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
133 PetscCall(MatMatMultNumeric_SeqAIJ_SeqDense(R, AB_den, RAB)); in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
[all …]
/petsc/share/petsc/datafiles/meshes/
H A Dmesh-3d-box-innersphere_bin.msh12 �*��?�Q���k�?U L�\ο8���_[f���?�R�}��˿��a��ͤ�9����O�$A��?(Jeܱ��?gz�j-E�?:�������׿ƿ �6|zƿ��
17 …R�����������}��}��?�}��}��?S�����������b���̿�b���̿T����������0����RͿ0����R�?U����������0����R�?0…
H A Dsquare_periodic_bin.msh9 …��$@�K���R�?�����������������$@hz����@�����������������$@4t���@�����������������$@����B� @��…
11 …����R�I*�@����������������������p@�����������������������B� @�������������������4t���@����…
14 …͹���?��������B�����흝�?ͧ͹���?��������C����C1�@���I��@��������D����Q�;8@R�� r@��������E���:6…
16 �e.�<�?"����5�?��������R���� �q
18 @&=R��@��������a��������@�oqr@��������b���/f��_@�#����@��������c����R�41@����
34 …���P����������*���I���;������������Q����������3���O���X������������R����������6���Y���8�…
36 …U���S������������j����������/���=���f������������k����������'���a���R������������l�����…
37 …)���L���:��������������������������8���Y�����������������������?���R���a��������������…
/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmresdeflation.c22 PetscInt i, L, R, j; in KSPAGMRESQuickSort() local
35 R = fin[i] - 1; in KSPAGMRESQuickSort()
36 if (L < R) { in KSPAGMRESQuickSort()
42 while (L < R) { in KSPAGMRESQuickSort()
43 abs_val = PetscSqrtReal(val_r[R] * val_r[R] + val_i[R] * val_i[R]); in KSPAGMRESQuickSort()
44 while (abs_val >= abs_pivot && L < R) { in KSPAGMRESQuickSort()
45 R--; in KSPAGMRESQuickSort()
46 abs_val = PetscSqrtReal(val_r[R] * val_r[R] + val_i[R] * val_i[R]); in KSPAGMRESQuickSort()
48 if (L < R) { in KSPAGMRESQuickSort()
49 val_r[L] = val_r[R]; in KSPAGMRESQuickSort()
[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
35 R = ksp->work[0]; in KSPSolve_PIPECG()
49 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_PIPECG()
50 PetscCall(VecAYPX(R, -1.0, B)); in KSPSolve_PIPECG()
52 PetscCall(VecCopy(B, R)); /* r <- b (x is 0) */ in KSPSolve_PIPECG()
55 PetscCall(KSP_PCApply(ksp, R, U)); /* u <- Br */ in KSPSolve_PIPECG()
65 PetscCall(VecNormBegin(R, NORM_2, &dp)); /* dp <- r'*r = e'*A'*A*e */ in KSPSolve_PIPECG()
66 PetscCall(PetscCommSplitReductionBegin(PetscObjectComm((PetscObject)R))); in KSPSolve_PIPECG()
68 PetscCall(VecNormEnd(R, NORM_2, &dp)); in KSPSolve_PIPECG()
71 PetscCall(VecDotBegin(R, U, &gamma)); /* gamma <- u'*r */ in KSPSolve_PIPECG()
[all …]
/petsc/src/sys/classes/draw/utils/
H A Dcmap.c55 static PetscErrorCode PetscDrawCmap_Hue(int mapsize, unsigned char R[], unsigned char G[], unsigned… in PetscDrawCmap_Hue() argument
65 R[i] = (unsigned char)(255 * PetscMin(r, 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
118 R[i] = (unsigned char)(255 * PetscMin(r, 1.0)); in PetscDrawCmap_Jet()
125 static PetscErrorCode PetscDrawCmap_Hot(int mapsize, unsigned char R[], unsigned char G[], unsigned… in PetscDrawCmap_Hot() argument
154 R[i] = (unsigned char)(255 * PetscMin(r, 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/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
27 R = ksp->work[0]; in KSPSolve_PIPEBCGS()
54 PetscCall(VecCopy(B, R)); in KSPSolve_PIPEBCGS()
55 PetscCall(VecAXPY(R, -1.0, Q2)); in KSPSolve_PIPEBCGS()
57 PetscCall(VecCopy(B, R)); in KSPSolve_PIPEBCGS()
61 if (ksp->normtype != KSP_NORM_NONE) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_PIPEBCGS()
73 PetscCall(VecCopy(R, RP)); /* rp <- r */ in KSPSolve_PIPEBCGS()
75 PetscCall(VecDotBegin(R, RP, &rho)); /* rho <- (r,rp) */ in KSPSolve_PIPEBCGS()
76 PetscCall(PetscCommSplitReductionBegin(PetscObjectComm((PetscObject)R))); in KSPSolve_PIPEBCGS()
77 PetscCall(KSP_PCApply(ksp, R, R2)); /* r2 <- K r */ in KSPSolve_PIPEBCGS()
[all …]
/petsc/src/dm/dt/tests/
H A Dex2.c9 PetscReal src_points[1000], target_points[1000], *R; in main() local
39 PetscCall(PetscMalloc1((nsrc_points - 1) * (ntarget_points - 1), &R)); in main()
41 …cDTReconstructPoly(degrees[i], nsrc_points - 1, src_points, ntarget_points - 1, target_points, R)); in main()
43 if (PetscAbs(R[j]) < 10 * PETSC_MACHINE_EPSILON) R[j] = 0; in main()
47 … PetscCall(PetscRealView(nsrc_points - 1, R + j * (nsrc_points - 1), PETSC_VIEWER_STDOUT_WORLD)); in main()
50 PetscCall(PetscFree(R)); in main()
/petsc/src/ts/tutorials/autodiff/adolc-utils/
H A Dsparse.cxx78 …tscErrorCode GenerateSeedMatrixPlusRecovery(ISColoring iscoloring, PetscScalar **S, PetscScalar *R) in GenerateSeedMatrixPlusRecovery() argument
91 R[indices[j]] = colour; in GenerateSeedMatrixPlusRecovery()
113 …etRecoveryMatrix(PetscScalar **S, unsigned int **sparsity, PetscInt m, PetscInt p, PetscScalar **R) in GetRecoveryMatrix() argument
120 R[i][colour] = -1.; in GetRecoveryMatrix()
124 R[i][colour] = j; in GetRecoveryMatrix()
147 PetscErrorCode RecoverJacobian(Mat A, InsertMode mode, PetscInt m, PetscInt p, PetscScalar **R, Pet… in RecoverJacobian() argument
152 PetscInt j = (PetscInt)R[i][colour]; in RecoverJacobian()
177 PetscErrorCode RecoverJacobianLocal(Mat A, InsertMode mode, PetscInt m, PetscInt p, PetscScalar **R… in RecoverJacobianLocal() argument
182 PetscInt j = (PetscInt)R[i][colour]; in RecoverJacobianLocal()
205 PetscErrorCode RecoverDiagonal(Vec diag, InsertMode mode, PetscInt m, PetscScalar *R, PetscScalar *… in RecoverDiagonal() argument
[all …]
/petsc/src/ksp/ksp/impls/cg/cgne/
H A Dcgne.c46 Vec X, B, Z, R, P, T; in KSPSolve_CGNE() local
61 R = ksp->work[0]; in KSPSolve_CGNE()
79 PetscCall(KSP_MatMultTranspose(ksp, Amat, P, R)); in KSPSolve_CGNE()
80 PetscCall(VecAYPX(R, -1.0, T)); in KSPSolve_CGNE()
82 PetscCall(VecCopy(T, R)); /* r <- b (x is 0) */ in KSPSolve_CGNE()
85 PetscCall(KSP_PCApplyTranspose(ksp, R, T)); in KSPSolve_CGNE()
87 PetscCall(KSP_PCApply(ksp, R, T)); in KSPSolve_CGNE()
94 PetscCall(VecNorm(R, NORM_2, &dp)); /* dp <- r'*r */ in KSPSolve_CGNE()
96 PetscCall(VecXDot(Z, R, &beta)); in KSPSolve_CGNE()
109 PetscCall(VecXDot(Z, R, &beta)); /* beta <- r'z */ in KSPSolve_CGNE()
[all …]
/petsc/src/ksp/ksp/impls/cg/
H A Dcg.c125 Vec X, B, Z, R, P, W; in KSPSolve_CG() local
139 R = ksp->work[0]; in KSPSolve_CG()
154 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_CG()
156 PetscCall(VecAYPX(R, -1.0, B)); in KSPSolve_CG()
162 PetscCall(VecCopy(B, R)); /* r <- b (x is 0) */ in KSPSolve_CG()
166 PetscCall(VecFlag(R, ksp->reason == KSP_DIVERGED_PC_FAILED)); in KSPSolve_CG()
170 PetscCall(KSP_PCApply(ksp, R, Z)); /* z <- Br */ in KSPSolve_CG()
175 PetscCall(VecNorm(R, NORM_2, &dp)); /* dp <- r'*r = e'*A'*A*e */ in KSPSolve_CG()
179 PetscCall(KSP_PCApply(ksp, R, Z)); /* z <- Br */ in KSPSolve_CG()
180 PetscCall(VecXDot(Z, R, &beta)); /* beta <- z'*r */ in KSPSolve_CG()
[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
59 R = ksp->work[0]; in KSPSolve_PIPEPRCG()
74 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_PIPEPRCG()
75 PetscCall(VecAYPX(R, -1.0, B)); in KSPSolve_PIPEPRCG()
77 PetscCall(VecCopy(B, R)); /* r <- b */ in KSPSolve_PIPEPRCG()
80 PetscCall(KSP_PCApply(ksp, R, RT)); /* rt <- Br */ in KSPSolve_PIPEPRCG()
91 PetscCall(VecDotBegin(RT, R, &nu)); in KSPSolve_PIPEPRCG()
95 PetscCall(VecDotEnd(RT, R, &nu)); /* nu <- (rt,r) */ in KSPSolve_PIPEPRCG()
110 PetscCall(VecNormBegin(R, NORM_2, &dp)); in KSPSolve_PIPEPRCG()
111 PetscCall(PetscCommSplitReductionBegin(PetscObjectComm((PetscObject)R))); in KSPSolve_PIPEPRCG()
[all …]
/petsc/src/ksp/ksp/impls/symmlq/
H A Dsymmlq.c21 Vec X, B, R, Z, U, V, W, UOLD, VOLD, Wbar; in KSPSolve_SYMMLQ() local
32 R = ksp->work[0]; in KSPSolve_SYMMLQ()
50 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - A*x */ in KSPSolve_SYMMLQ()
51 PetscCall(VecAYPX(R, -1.0, B)); in KSPSolve_SYMMLQ()
53 PetscCall(VecCopy(B, R)); /* r <- b (x is 0) */ in KSPSolve_SYMMLQ()
56 PetscCall(KSP_PCApply(ksp, R, Z)); /* z <- B*r */ in KSPSolve_SYMMLQ()
57 PetscCall(VecDot(R, Z, &dp)); /* dp = r'*z; */ in KSPSolve_SYMMLQ()
77 PetscCall(VecCopy(R, V)); /* v <- r; */ in KSPSolve_SYMMLQ()
103 PetscCall(VecCopy(R, V)); in KSPSolve_SYMMLQ()
120 PetscCall(KSP_MatMult(ksp, Amat, U, R)); /* r <- Amat*u; */ in KSPSolve_SYMMLQ()
[all …]

1234567