| /petsc/src/mat/impls/aij/mpi/strumpack/ |
| H A D | strumpack.c | 14 STRUMPACK_SparseSolver *S = (STRUMPACK_SparseSolver *)A->data; in MatDestroy_STRUMPACK() local 18 PetscStackCallExternalVoid("STRUMPACK_destroy", STRUMPACK_destroy(S)); in MatDestroy_STRUMPACK() 53 STRUMPACK_SparseSolver *S = (STRUMPACK_SparseSolver *)F->data; in MatSTRUMPACKSetReordering_STRUMPACK() local 56 …PetscStackCallExternalVoid("STRUMPACK_set_reordering_method", STRUMPACK_set_reordering_method(*S, … in MatSTRUMPACKSetReordering_STRUMPACK() 61 STRUMPACK_SparseSolver *S = (STRUMPACK_SparseSolver *)F->data; in MatSTRUMPACKGetReordering_STRUMPACK() local 64 …RUMPACK_reordering_method", *reordering = (MatSTRUMPACKReordering)STRUMPACK_reordering_method(*S)); in MatSTRUMPACKGetReordering_STRUMPACK() 118 STRUMPACK_SparseSolver *S = (STRUMPACK_SparseSolver *)F->data; in MatSTRUMPACKSetColPerm_STRUMPACK() local 121 …PetscStackCallExternalVoid("STRUMPACK_set_matching", STRUMPACK_set_matching(*S, cperm ? STRUMPACK_… in MatSTRUMPACKSetColPerm_STRUMPACK() 126 STRUMPACK_SparseSolver *S = (STRUMPACK_SparseSolver *)F->data; in MatSTRUMPACKGetColPerm_STRUMPACK() local 129 …PetscStackCallExternalVoid("STRUMPACK_matching", *cperm = (PetscBool)(STRUMPACK_matching(*S) != ST… in MatSTRUMPACKGetColPerm_STRUMPACK() [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex203.c | 25 Mat A, S; in main() local 45 PetscCall(MatCreateShell(PETSC_COMM_WORLD, 2, 2, 2, 2, user, &S)); in main() 46 PetscCall(MatShellSetOperation(S, MATOP_GET_DIAGONAL, (PetscErrorCodeFn *)MatGetDiagonal_User)); in main() 47 PetscCall(MatSetUp(S)); in main() 49 PetscCall(MatShift(S, 42)); in main() 50 PetscCall(MatGetDiagonal(S, Y)); in main() 52 PetscCall(MatDiagonalSet(S, X, ADD_VALUES)); in main() 53 PetscCall(MatGetDiagonal(S, Y)); in main() 55 PetscCall(MatScale(S, 42)); in main() 56 PetscCall(MatGetDiagonal(S, Y)); in main() [all …]
|
| H A D | ex218.c | 43 Mat S; in main() local 56 PetscCall(MatCreateShell(PETSC_COMM_WORLD, 2, 2, 2, 2, user, &S)); in main() 57 PetscCall(MatSetUp(S)); in main() 58 PetscCall(MatShellSetOperation(S, MATOP_MULT, (PetscErrorCodeFn *)MatMult_User)); in main() 59 …PetscCall(MatShellSetOperation(S, MATOP_MULT_TRANSPOSE, (PetscErrorCodeFn *)MatMultTranspose_User)… in main() 61 PetscCall(MatShellTestMult(S, MyFunction, base, user, NULL)); in main() 62 PetscCall(MatShellTestMultTranspose(S, MyFunction, base, user, NULL)); in main() 66 PetscCall(MatDestroy(&S)); in main()
|
| H A D | ex69.c | 6 static PetscErrorCode MatMult_S(Mat S, Vec x, Vec y) in MatMult_S() argument 11 PetscCall(MatShellGetContext(S, &A)); in MatMult_S() 18 static PetscErrorCode MatMultTranspose_S(Mat S, Vec x, Vec y) in MatMultTranspose_S() argument 23 PetscCall(MatShellGetContext(S, &A)); in MatMultTranspose_S() 34 Mat A, B, C, S; in main() local 120 PetscCall(MatCreateShell(PetscObjectComm((PetscObject)v), nloc, nloc, n, n, A, &S)); in main() 121 PetscCall(MatShellSetOperation(S, MATOP_MULT, (PetscErrorCodeFn *)MatMult_S)); in main() 122 … PetscCall(MatShellSetOperation(S, MATOP_MULT_TRANSPOSE, (PetscErrorCodeFn *)MatMultTranspose_S)); in main() 123 PetscCall(MatShellSetVecType(S, vtype)); in main() 126 S = A; in main() [all …]
|
| H A D | ex221.c | 65 Mat A, S; in main() local 162 PetscCall(MatCreateShell(PETSC_COMM_WORLD, m, n, M, N, user, &S)); in main() 163 PetscCall(MatShellSetOperation(S, MATOP_MULT, (PetscErrorCodeFn *)MatMult_User)); in main() 164 …PetscCall(MatShellSetOperation(S, MATOP_MULT_TRANSPOSE, (PetscErrorCodeFn *)MatMultTranspose_User)… in main() 165 …if (cong) PetscCall(MatShellSetOperation(S, MATOP_GET_DIAGONAL, (PetscErrorCodeFn *)MatGetDiagonal… in main() 166 PetscCall(MatShellSetOperation(S, MATOP_COPY, (PetscErrorCodeFn *)MatCopy_User)); in main() 167 PetscCall(MatShellSetOperation(S, MATOP_DESTROY, (PetscErrorCodeFn *)MatDestroy_User)); in main() 176 PetscCall(MatMultAddEqual(A, S, 10, &flg)); in main() 178 PetscCall(MatMultTransposeAddEqual(A, S, 10, &flg)); in main() 180 PetscCall(MatMultHermitianTransposeAddEqual(A, S, 10, &flg)); in main() [all …]
|
| H A D | ex177.c | 9 PetscScalar *S, *T; in main() local 32 PetscCall(PetscMalloc2(p * q, &S, p * q, &T)); in main() 33 for (i = 0; i < p * q; i++) S[i] = 0; in main() 38 S[i + p * j] = ((PetscReal)((i + 1) * (j + 1))) / ((PetscReal)(p + q)); in main() 46 PetscCall(MatCreateKAIJ(A, p, q, S, T, &TA)); in main() 90 PetscCall(MatCreateKAIJ(A, p, q, S, NULL, &TA)); in main() 120 PetscCall(MatCreateKAIJ(A, p, q, S, T, &TA)); in main() 147 if (i == j) S[i + j * p] = T[i + j * p] = 2.0; in main() 148 else S[i + j * p] = T[i + j * p] = 0.0; in main() 151 PetscCall(MatCreateKAIJ(A, p, q, S, T, &TA)); in main() [all …]
|
| /petsc/src/ksp/ksp/utils/schurm/ |
| H A D | schurm.c | 175 PetscErrorCode MatCreateSchurComplement(Mat A00, Mat Ap00, Mat A01, Mat A10, Mat A11, Mat *S) in MatCreateSchurComplement() argument 179 PetscCall(MatCreate(PetscObjectComm((PetscObject)A00), S)); in MatCreateSchurComplement() 180 PetscCall(MatSetType(*S, MATSCHURCOMPLEMENT)); in MatCreateSchurComplement() 181 PetscCall(MatSchurComplementSetSubMatrices(*S, A00, Ap00, A01, A10, A11)); in MatCreateSchurComplement() 217 PetscErrorCode MatSchurComplementSetSubMatrices(Mat S, Mat A00, Mat Ap00, Mat A01, Mat A10, Mat A11) in MatSchurComplementSetSubMatrices() argument 219 Mat_SchurComplement *Na = (Mat_SchurComplement *)S->data; in MatSchurComplementSetSubMatrices() 223 PetscCall(PetscObjectTypeCompare((PetscObject)S, MATSCHURCOMPLEMENT, &isschur)); in MatSchurComplementSetSubMatrices() 225 …PetscCheck(!S->assembled, PetscObjectComm((PetscObject)S), PETSC_ERR_ARG_WRONGSTATE, "Use MatSchur… in MatSchurComplementSetSubMatrices() 244 PetscCall(MatSetSizes(S, A10->rmap->n, A01->cmap->n, A10->rmap->N, A01->cmap->N)); in MatSchurComplementSetSubMatrices() 255 PetscCall(MatSetUp(S)); in MatSchurComplementSetSubMatrices() [all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex21.c | 107 Mat A, S = NULL, Sexplicit = NULL, Sp, B, C; in main() local 124 …PetscCall(MatGetSchurComplement(A, is0, is0, is1, is1, MAT_INITIAL_MATRIX, &S, ainv_type, MAT_IGNO… in main() 125 PetscCall(MatSetFromOptions(S)); in main() 126 PetscCall(MatComputeOperator(S, MATAIJ, &Sexplicit)); in main() 130 PetscCall(MatSchurComplementSetAinvType(S, MAT_SCHUR_COMPLEMENT_AINV_FULL)); in main() 131 PetscCall(MatSchurComplementGetPmat(S, MAT_INITIAL_MATRIX, &Sp)); in main() 134 PetscCall(MatSchurComplementSetAinvType(S, MAT_SCHUR_COMPLEMENT_AINV_DIAG)); in main() 142 PetscCall(MatMatMult(S, B, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &C)); in main() 148 PetscCall(MatDestroy(&S)); in main() 153 …PetscCall(MatGetSchurComplement(A, is1, is1, is0, is0, MAT_INITIAL_MATRIX, &S, ainv_type, MAT_IGNO… in main() [all …]
|
| H A D | ex82.c | 11 Vec b, x, *S; in main() local 29 PetscCall(VecDuplicateVecs(b, N, &S)); in main() 42 …PetscCall(KSPComputeRitz(ksp, harmonic ? PETSC_FALSE : PETSC_TRUE, PETSC_TRUE, &Na, S, tetar, teta… in main() 56 PetscCall(VecView(S[i], PETSC_VIEWER_STDOUT_WORLD)); in main() 60 PetscCall(VecView(S[i + 1], PETSC_VIEWER_STDOUT_WORLD)); in main() 69 PetscCall(VecDestroyVecs(N, &S)); in main()
|
| /petsc/lib/petsc/bin/maint/abi-compliance-checker/modules/Internals/ |
| H A D | Basic.pm | 517 my $S = $_[0]; 519 if(index($S, '@')==-1 520 and index($S, '$')==-1) { 521 return ($S, "", ""); 524 if($S=~/\A([^\@\$\?]+)([\@\$]+)([^\@\$]+)\Z/) { 528 return ($S, "", ""); 669 my $S = substr($Params, $Pos, 1); 670 if(defined $B{$S}) { 671 $B{$S} += 1; 673 if($S eq "," and [all …]
|
| H A D | Path.pm | 84 my $S = "/"; 86 $S = "\\"; 88 return join($S, @_);
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | mesh-3d-box-innersphere_bin.msh | 11 …��#�?^|Uȩ�?���~Ӈ�?2� ��0��d�;ſ�������?3��V�����jrϐ�¿Z�B���ɿ4ӮS$(�?.�70+��0����?5… 17 …S�οK��_M'ڶ�q� ��[�?��t;��?L<Mv�d�?�%�����?��.��?M-�7�E?��i��7�?aa[`�?N�Y�yO�?… 19 �g�?�-�)��?�(${9��?|���}��?�@S>�K�?���_/���?I�C`�v�?s>GPm%ӿ��v���)�?�xI!��?���۵�?��… 24 …TSS… 25 XSTXV… 26 …SYYUW… 52 fjcS
|
| H A D | square_periodic_bin.msh | 8 $@t���:�?$@�?#����?$@�<7*�S�? 11 …hz���@ �K���R�?!�<7*�S�?"… 17 …S�#j�9�@�pz��C@T,h�="@�ȶu�V@U\T+�m@,Z[7*��?V�,]t��?M�… 34 …Q3OXR6Y8S 35 …VaXb&Sbc… 36 …Sj/=fk'aRl… 38 …SE�c�3Lg�…
|
| H A D | square_bin_physnames.msh | 26 ��?12)��1��?˄�ʼnS�?2E�킊�?��4T���?3��4T���?/�KH���?4�7-�… 28 …k#�?T�c���?Q�_�8���?����k#�?R�k��[��?�b���T�?Sc���T�?�(��HM�?… 33 …?�L��O�?�o���D�?����S�.�?s��v�?����v�?�����E�?�C… 44 …Sg7< Rf6= Ui<… 45 …Q wIR %xFS 46 …S?o} T@p~ UAq… 48 …S�?� T�@� U�A�…
|
| H A D | surfacesphere_bin.msh | 9 …=o�俙�ׁ��?��V����8�����V�w_r�`,�?�Fr46ֿ9�!"J�ݿ�W�%���?�^Q块�:8'���SϿ�����?��V4�;f… 15 \�ٿ���+�?Ŀ"P"�T��Rڕ���U� 16 G�%�c���a��ĿS�}n���ǿ�h�\f��zĥt�VǿT*%����ʿ�ԤϽ�ſG��bN��Ug p�&�Y�Wٗ���l!�Vi���E�� F… 19 …|��w��ѿ��U9���?��¿v�?}D;���㿏e�G���?:�S����?~P}���Iݿ����.��?������?���=�鿴I��r��?�"w… 21 …䨟K��]��ƿ"�!�@�?�#�c�BX�� ��X"Ŀ���o��?�J��3Ŀ(��U���Ze�?��E���ǿ�p�06�ƿS@�����?�� 22 …e��?R]��,^п�s%�4�ٿ���V4��?�����SϿ�����꿫���RE�?1`V����'�����ǿ�Y�{P�N�?���:�q���ſ�… 24 &Р�<ۿ@%G��濲imi��=�?`�!���濴�x�!�俳V�p���?{S�������@V�ѿ�1ۦr��?y�vNݿ%���Hqҿ����o ��?�#… 26 …?�81¶�&`0�1i�?�9�S����?B;���㿎e�G���?���¿v�?��w��ѿ��U9���?�������?P}���Iݿ����.��?��"… 35 …R.<RO=S…
|
| /petsc/src/mat/impls/aij/seq/aijsell/ |
| H A D | aijsell.c | 12 Mat S; /* The SELL formatted "shadow" matrix. */ member 43 PetscCall(MatDestroy(&aijsell->S)); in MatConvert_SeqAIJSELL_SeqAIJ() 62 PetscCall(MatDestroy(&aijsell->S)); in MatDestroy_SeqAIJSELL() 86 if (aijsell->S && aijsell->state == state) { in MatSeqAIJSELL_build_shadow() 92 if (aijsell->S) { in MatSeqAIJSELL_build_shadow() 93 PetscCall(MatConvert_SeqAIJ_SeqSELL(A, MATSEQSELL, MAT_REUSE_MATRIX, &aijsell->S)); in MatSeqAIJSELL_build_shadow() 95 PetscCall(MatConvert_SeqAIJ_SeqSELL(A, MATSEQSELL, MAT_INITIAL_MATRIX, &aijsell->S)); in MatSeqAIJSELL_build_shadow() 115 aijsell_dest->S = NULL; in MatDuplicate_SeqAIJSELL() 153 PetscCall(MatMult_SeqSELL(aijsell->S, xx, yy)); in MatMult_SeqAIJSELL() 163 PetscCall(MatMultTranspose_SeqSELL(aijsell->S, xx, yy)); in MatMultTranspose_SeqAIJSELL() [all …]
|
| /petsc/src/ksp/ksp/impls/tsirm/ |
| H A D | tsirm.c | 6 Mat A, S; member 21 PetscCall(MatCreate(PetscObjectComm((PetscObject)ksp), &tsirm->S)); in KSPSetUp_TSIRM() 22 …PetscCall(MatSetSizes(tsirm->S, tsirm->Iend - tsirm->Istart, PETSC_DECIDE, tsirm->size, tsirm->siz… in KSPSetUp_TSIRM() 23 PetscCall(MatSetType(tsirm->S, MATDENSE)); in KSPSetUp_TSIRM() 24 PetscCall(MatSetUp(tsirm->S)); in KSPSetUp_TSIRM() 27 PetscCall(MatCreateVecs(tsirm->S, &tsirm->Alpha, &tsirm->r)); in KSPSetUp_TSIRM() 77 …PetscCall(MatSetValues(tsirm->S, tsirm->Iend - tsirm->Istart, ind_row, 1, &col, array, INSERT_VALU… in KSPSolve_TSIRM() 89 PetscCall(MatAssemblyBegin(tsirm->S, MAT_FINAL_ASSEMBLY)); in KSPSolve_TSIRM() 90 PetscCall(MatAssemblyEnd(tsirm->S, MAT_FINAL_ASSEMBLY)); in KSPSolve_TSIRM() 92 PetscCall(MatMatMult(tsirm->A, tsirm->S, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &AS)); in KSPSolve_TSIRM() [all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | brdn.c | 53 LMBasis S = NULL, Y = NULL; in BroydenKernel_Recursive() local 61 PetscCall(MatLMVMGetUpdatedBasis(B, S_t, &S, NULL, NULL)); in BroydenKernel_Recursive() 65 PetscCall(LMBasisGetWorkVec(S, &G)); in BroydenKernel_Recursive() 73 PetscCall(LMBasisGetVecRead(S, i, &s_i)); in BroydenKernel_Recursive() 83 PetscCall(LMBasisRestoreVecRead(S, i, &s_i)); in BroydenKernel_Recursive() 92 PetscCall(LMBasisRestoreWorkVec(S, &G)); in BroydenKernel_Recursive() 117 LMBasis S, Y; in BroydenKernelHermitianTranspose_Recursive() local 121 PetscCall(MatLMVMGetUpdatedBasis(B, S_t, &S, NULL, NULL)); in BroydenKernelHermitianTranspose_Recursive() 128 PetscCall(LMBasisGetVecRead(S, i, &s_i)); in BroydenKernelHermitianTranspose_Recursive() 138 PetscCall(LMBasisRestoreVecRead(S, i, &s_i)); in BroydenKernelHermitianTranspose_Recursive() [all …]
|
| /petsc/src/binding/petsc4py/docs/source/ |
| H A D | citing.rst | 13 * S. Balay, S. Abhyankar, M. Adams, S. Benson, J. Brown, 16 D. Karpeyev, D. Kaushik, M. Knepley, F. Kong, S. Kruger, 19 S. Zampini, H. Zhang, and H. Zhang, J. Zhang,
|
| /petsc/src/ksp/ksp/impls/cg/pipeprcg/ |
| H A D | pipeprcg.c | 44 Vec X, B, R, RT, W, WT, P, S, ST, U, UT, PRTST[3]; in KSPSolve_PIPEPRCG() local 64 S = ksp->work[5]; in KSPSolve_PIPEPRCG() 85 PetscCall(VecCopy(W, S)); /* p <- rt */ in KSPSolve_PIPEPRCG() 92 PetscCall(VecDotBegin(P, S, mu_p)); in KSPSolve_PIPEPRCG() 93 PetscCall(VecDotBegin(ST, S, gamma_p)); in KSPSolve_PIPEPRCG() 96 PetscCall(VecDotEnd(P, S, mu_p)); /* mu <- (p,s) */ in KSPSolve_PIPEPRCG() 97 PetscCall(VecDotEnd(ST, S, gamma_p)); /* gamma <- (st,s) */ in KSPSolve_PIPEPRCG() 138 PetscCall(VecAXPY(R, -alpha, S)); /* r <- r - alpha * s */ in KSPSolve_PIPEPRCG() 143 PetscCall(VecAYPX(S, beta, W)); /* s <- w + beta * s */ in KSPSolve_PIPEPRCG() 152 PetscCall(VecMDotBegin(S, 3, PRTST, mudelgam)); in KSPSolve_PIPEPRCG() [all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/dfp/ |
| H A D | dfp.c | 17 LMBasis S, Y; in DFPKernel_Recursive() local 25 PetscCall(MatLMVMGetUpdatedBasis(B, S_t, &S, NULL, NULL)); in DFPKernel_Recursive() 29 PetscCall(LMBasisGetWorkVec(S, &G)); in DFPKernel_Recursive() 42 PetscCall(LMBasisGetVecRead(S, i, &s_i)); in DFPKernel_Recursive() 44 PetscCall(LMBasisRestoreVecRead(S, i, &s_i)); in DFPKernel_Recursive() 51 PetscCall(LMBasisGetVecRead(S, i, &s_i)); in DFPKernel_Recursive() 53 PetscCall(LMBasisRestoreVecRead(S, i, &s_i)); in DFPKernel_Recursive() 63 PetscCall(LMBasisRestoreWorkVec(S, &G)); in DFPKernel_Recursive() 81 LMBasis S, Y; in DFPKernel_CompactDense() local 84 PetscCall(MatLMVMGetUpdatedBasis(B, S_t, &S, NULL, NULL)); in DFPKernel_CompactDense() [all …]
|
| /petsc/src/mat/impls/kaij/ |
| H A D | kaij.c | 86 PetscErrorCode MatKAIJGetS(Mat A, PetscInt *m, PetscInt *n, PetscScalar *S[]) in MatKAIJGetS() argument 93 if (S) *S = b->S; in MatKAIJGetS() 117 PetscErrorCode MatKAIJGetSRead(Mat A, PetscInt *m, PetscInt *n, const PetscScalar *S[]) in MatKAIJGetSRead() argument 124 if (S) *S = b->S; in MatKAIJGetSRead() 145 PetscErrorCode MatKAIJRestoreS(Mat A, PetscScalar *S[]) in MatKAIJRestoreS() argument 148 if (S) *S = NULL; in MatKAIJRestoreS() 170 PetscErrorCode MatKAIJRestoreSRead(Mat A, const PetscScalar *S[]) in MatKAIJRestoreSRead() argument 173 if (S) *S = NULL; in MatKAIJRestoreSRead() 346 PetscErrorCode MatKAIJSetS(Mat A, PetscInt p, PetscInt q, const PetscScalar S[]) in MatKAIJSetS() argument 351 PetscCall(PetscFree(a->S)); in MatKAIJSetS() [all …]
|
| /petsc/src/binding/petsc4py/ |
| H A D | DESCRIPTION.rst | 36 * S. Balay, S. Abhyankar, M. Adams, S. Benson, J. Brown, 39 D. Karpeyev, D. Kaushik, M. Knepley, F. Kong, S. Kruger, 42 S. Zampini, H. Zhang, and H. Zhang, J. Zhang,
|
| /petsc/src/ksp/ksp/impls/bcgs/fbcgs/ |
| H A D | fbcgs.c | 19 Vec X, B, V, P, R, RP, T, S, P2, S2; in KSPSolve_FBCGS() local 32 S = ksp->work[4]; in KSPSolve_FBCGS() 92 PetscCall(VecWAXPY(S, -alpha, V, R)); /* s <- r - alpha v */ in KSPSolve_FBCGS() 94 PetscCall(KSP_PCApply(ksp, S, S2)); /* s2 <- K s */ in KSPSolve_FBCGS() 97 PetscCall(VecDotNorm2(S, T, &d1, &d2)); in KSPSolve_FBCGS() 100 PetscCall(VecDot(S, S, &d1)); in KSPSolve_FBCGS() 120 PetscCall(VecWAXPY(R, -omega, T, S)); /* r <- s - omega t */ in KSPSolve_FBCGS()
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_mat_schur.py | 21 S = PETSc.Mat().createSchurComplement(A00, A00, A01, A10, A11) 22 M, N = S.getSize() 25 m, n = S.getLocalSize() 28 A00_dup, A00p_dup, A01_dup, A10_dup, A11_dup = S.getSchurComplementSubMatrices()
|