| /petsc/src/vec/is/sf/tests/ |
| H A D | ex4.c | 11 Vec a, b, ba; in main() local 77 PetscCall(VecCreateSeq(PETSC_COMM_WORLD, maxleafB + 1, &ba)); in main() 95 PetscCall(VecGetArray(ba, &arrayW)); in main() 100 PetscCall(VecRestoreArray(ba, &arrayW)); in main() 103 PetscCall(VecView(ba, NULL)); in main() 109 PetscCall(VecGetArray(ba, &arrayW)); in main() 113 PetscCall(VecRestoreArray(ba, &arrayW)); in main() 116 PetscCall(VecView(ba, NULL)); in main() 118 PetscCall(VecDestroy(&ba)); in main()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 566 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local 650 diag = ba + i * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N() 651 u = ba + ili * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N() 660 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N() 669 u = ba + j * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N() 686 diag = ba + k * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N() 697 u = ba + j * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N() 739 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local 782 diag = ba + i * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() 783 u = ba + ili * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() [all …]
|
| H A D | sbaijfact4.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local 55 diag = ba + i * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() 56 u = ba + ili * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() 86 PetscCall(PetscArraycpy(ba + ili * 9, uik, 9)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() 95 u = ba + j * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() 122 diag = ba + k * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() 132 u = ba + j * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
|
| H A D | sbaijfact3.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local 93 diag = ba + i * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3() 94 u = ba + ili * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3() 124 PetscCall(PetscArraycpy(ba + ili * 9, uik, 9)); in MatCholeskyFactorNumeric_SeqSBAIJ_3() 133 u = ba + j * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3() 160 diag = ba + k * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3() 170 u = ba + j * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
|
| H A D | sbaijfact5.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() local 56 diag = ba + i * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() 57 u = ba + ili * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() 103 PetscCall(PetscArraycpy(ba + ili * 16, uik, 16)); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() 112 u = ba + j * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() 147 diag = ba + k * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() 161 u = ba + j * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
|
| H A D | sbaijfact6.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_4() local 94 diag = ba + i * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4() 95 u = ba + ili * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4() 141 PetscCall(PetscArraycpy(ba + ili * 16, uik, 16)); in MatCholeskyFactorNumeric_SeqSBAIJ_4() 150 u = ba + j * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4() 185 diag = ba + k * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4() 200 u = ba + j * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
|
| H A D | sbaijfact8.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local 55 d = ba + i * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() 56 u = ba + ili * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() 122 PetscCall(PetscArraycpy(ba + ili * 25, uik, 25)); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() 131 u = ba + j * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() 176 d = ba + k * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() 186 u = ba + j * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
|
| H A D | sbaijfact7.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local 93 d = ba + i * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5() 94 u = ba + ili * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5() 160 PetscCall(PetscArraycpy(ba + ili * 25, uik, 25)); in MatCholeskyFactorNumeric_SeqSBAIJ_5() 169 u = ba + j * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5() 214 d = ba + k * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5() 224 u = ba + j * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
|
| H A D | sbaijfact10.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local 63 d = ba + i * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 64 u = ba + ili * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 228 PetscCall(PetscArraycpy(ba + ili * 36, uik, 36)); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 237 u = ba + j * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 332 d = ba + k * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 342 u = ba + j * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
|
| H A D | sbaijfact9.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local 95 d = ba + i * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6() 96 u = ba + ili * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6() 223 PetscCall(PetscArraycpy(ba + ili * 36, uik, 36)); in MatCholeskyFactorNumeric_SeqSBAIJ_6() 232 u = ba + j * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6() 327 d = ba + k * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6() 337 u = ba + j * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
|
| H A D | sbaijfact12.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local 59 d = ba + i * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() 60 u = ba + ili * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() 228 PetscCall(PetscArraycpy(ba + ili * 49, uik, 49)); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() 237 u = ba + j * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() 359 d = ba + k * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() 369 u = ba + j * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
|
| H A D | sbaijfact11.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local 96 d = ba + i * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7() 97 u = ba + ili * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7() 265 PetscCall(PetscArraycpy(ba + ili * 49, uik, 49)); in MatCholeskyFactorNumeric_SeqSBAIJ_7() 274 u = ba + j * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7() 396 d = ba + k * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7() 406 u = ba + j * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
|
| /petsc/src/mat/tests/ |
| H A D | ex303.c | 12 const PetscScalar *ba; in main() local 56 PetscCall(MatSeqAIJGetArrayRead(B, &ba)); in main() 62 …i[i + 1] - bi[i]; j++) PetscCall(MatSetValue(B2, i, garray[bj[bi[i] + j]], ba[bi[i] + j], INSERT_V… in main() 68 PetscCall(MatSeqAIJRestoreArrayRead(B, &ba)); in main()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 235 MatScalar *ba = b->a, *aa = a->a; in MatILUFactorNumeric_SeqBAIJ_2_inplace() local 277 pv = ba + 4 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_2_inplace() 306 pv = ba + 4 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_inplace() 318 w = ba + 4 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_2_inplace() 346 MatScalar *ba = b->a, *aa = a->a; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() local 384 pv = ba + 4 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() 413 pv = ba + 4 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() 431 w = ba + 4 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() 623 MatScalar *ba = b->a, *aa = a->a; in MatILUFactorNumeric_SeqBAIJ_1_inplace() local 650 pv = ba + diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_1_inplace() [all …]
|
| H A D | baijfact4.c | 16 MatScalar *ba = b->a, *aa = a->a, *pv, *v, *rtmp, *multiplier, *v_work, *pc, *w; in MatILUFactorNumeric_SeqBAIJ_N_inplace() local 52 pv = ba + bs2 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_N_inplace() 63 pv = ba + bs2 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
|
| H A D | baijfact13.c | 22 MatScalar *ba = b->a, *aa = a->a; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local 74 pv = ba + 9 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_3_inplace() 127 pv = ba + 9 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace() 144 w = ba + 9 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace() 282 MatScalar *ba = b->a, *aa = a->a; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace() local 331 pv = ba + 9 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace() 385 pv = ba + 9 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace() 402 w = ba + 9 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace()
|
| H A D | baijfact11.c | 24 MatScalar *ba = b->a, *aa = a->a; in MatILUFactorNumeric_SeqBAIJ_4_inplace() local 92 pv = ba + 16 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_4_inplace() 177 pv = ba + 16 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace() 201 w = ba + 16 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace() 355 MatScalar *ba = b->a, *aa = a->a; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace() local 420 pv = ba + 16 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace() 504 pv = ba + 16 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace() 528 w = ba + 16 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace()
|
| H A D | baijfact9.c | 26 MatScalar *ba = b->a, work[25]; in MatILUFactorNumeric_SeqBAIJ_5_inplace() local 123 pv = ba + 25 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_5_inplace() 246 pv = ba + 25 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace() 283 w = ba + 25 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace() 428 MatScalar *ba = b->a, *aa = a->a, work[25]; in MatILUFactorNumeric_SeqBAIJ_5_NaturalOrdering_inplace() local 506 pv = ba + 25 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_5_NaturalOrdering_inplace() 628 pv = ba + 25 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_NaturalOrdering_inplace() 661 w = ba + 25 * diag_offset[i]; in MatILUFactorNumeric_SeqBAIJ_5_NaturalOrdering_inplace()
|
| /petsc/src/mat/impls/aij/seq/bas/ |
| H A D | basfactor.c | 101 MatScalar *ba = b->a; in MatCholeskyFactorNumeric_SeqAIJ_Bas() local 110 PetscCall(PetscFree(ba)); in MatCholeskyFactorNumeric_SeqAIJ_Bas() 144 PetscCall(spbas_matrix_to_crs(matrix_L, &ba, &bi, &bj)); in MatCholeskyFactorNumeric_SeqAIJ_Bas() 147 b->a = ba; in MatCholeskyFactorNumeric_SeqAIJ_Bas()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 229 MatScalar *ba; in MatLUFactorNumeric_SeqAIJ() local 238 PetscCall(MatSeqAIJGetArrayWrite(B, &ba)); in MatLUFactorNumeric_SeqAIJ() 294 pv = ba + bdiag[row]; in MatLUFactorNumeric_SeqAIJ() 299 pv = ba + bdiag[row + 1] + 1; in MatLUFactorNumeric_SeqAIJ() 311 pv = ba + bi[i]; in MatLUFactorNumeric_SeqAIJ() 320 pv = ba + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ() 335 pv = ba + bdiag[i]; in MatLUFactorNumeric_SeqAIJ() 355 PetscCall(MatSeqAIJRestoreArrayWrite(B, &ba)); in MatLUFactorNumeric_SeqAIJ() 408 MatScalar *ba; in MatLUFactorNumeric_SeqAIJ_inplace() local 417 PetscCall(MatSeqAIJGetArrayWrite(B, &ba)); in MatLUFactorNumeric_SeqAIJ_inplace() [all …]
|
| H A D | matmatmult.c | 258 const PetscScalar *aa, *ba, *baj; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted() local 262 PetscCall(MatSeqAIJGetArrayRead(B, &ba)); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted() 290 baj = PetscSafePointerPlusOffset(ba, bi[brow]); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted() 315 PetscCall(MatSeqAIJRestoreArrayRead(B, &ba)); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted() 329 const PetscScalar *aa, *ba, *baj; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable() local 334 PetscCall(MatSeqAIJGetArrayRead(B, &ba)); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable() 353 baj = ba + bi[brow]; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable() 376 PetscCall(MatSeqAIJRestoreArrayRead(B, &ba)); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable() 1304 MatScalar *aa = a->a, *aval, *ba = b->a, *bval, *ca, *cval; in MatMatTransposeMultNumeric_SeqAIJ_SeqAIJ() local 1349 bval = ba + bi[brow]; in MatMatTransposeMultNumeric_SeqAIJ_SeqAIJ() [all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mmaij.c | 153 const PetscScalar *ba; in MatDisAssemble_MPIAIJ() local 182 PetscCall(MatSeqAIJGetArrayRead(B, &ba)); in MatDisAssemble_MPIAIJ() 186 v = ba[ct++]; in MatDisAssemble_MPIAIJ() 190 PetscCall(MatSeqAIJRestoreArrayRead(B, &ba)); in MatDisAssemble_MPIAIJ()
|
| H A D | mpiaij.c | 478 …MatSeqXAIJReallocateAIJ(B, bm, 1, nrow2, row, col, rmax2, ba, bi, bj, rp2, ap2, bimax, nonew, MatS… 496 PetscScalar *aa, *ba; in MatSetValuesRow_MPIAIJ() local 508 PetscCall(MatSeqAIJGetArray(mat->B, &ba)); in MatSetValuesRow_MPIAIJ() 509 PetscCall(PetscArraycpy(ba + b->i[row], v, l)); in MatSetValuesRow_MPIAIJ() 510 PetscCall(MatSeqAIJRestoreArray(mat->B, &ba)); in MatSetValuesRow_MPIAIJ() 522 PetscCall(MatSeqAIJGetArray(mat->B, &ba)); in MatSetValuesRow_MPIAIJ() 523 …PetscCall(PetscArraycpy(ba + b->i[row] + l, v + l + a->i[row + 1] - a->i[row], b->i[row + 1] - b->… in MatSetValuesRow_MPIAIJ() 524 PetscCall(MatSeqAIJRestoreArray(mat->B, &ba)); in MatSetValuesRow_MPIAIJ() 545 MatScalar *aa, *ba; in MatSetValues_MPIAIJ() local 552 PetscCall(MatSeqAIJGetArray(B, &ba)); in MatSetValues_MPIAIJ() [all …]
|
| /petsc/src/mat/impls/aij/seq/umfpack/ |
| H A D | umfpack.c | 114 const PetscScalar *ba; in MatSolve_UMFPACK_Private() local 131 PetscCall(VecGetArrayRead(b, &ba)); in MatSolve_UMFPACK_Private() 134 …(uflag, ai, aj, (PetscReal *)av, NULL, (PetscReal *)xa, NULL, (PetscReal *)ba, NULL, lu->Numeric, … in MatSolve_UMFPACK_Private() 136 …status = umfpack_UMF_wsolve(uflag, ai, aj, av, xa, ba, lu->Numeric, lu->Control, lu->Info, lu->Wi,… in MatSolve_UMFPACK_Private() 144 PetscCall(VecRestoreArrayRead(b, &ba)); in MatSolve_UMFPACK_Private()
|
| /petsc/src/mat/impls/aij/mpi/crl/ |
| H A D | mcrl.c | 44 PetscScalar *aa = Aij->a, *ba = Bij->a, *acols, *array; in MatMPIAIJCRL_create_aijcrl() local 63 acols[j * m + i] = *ba++; in MatMPIAIJCRL_create_aijcrl()
|