| /petsc/src/mat/impls/sbaij/seq/ ! |
| H A D | aijsbaij.c | 12 MatScalar *av, *bv; in MatConvert_SeqSBAIJ_SeqAIJ() local 67 av = a->a; in MatConvert_SeqSBAIJ_SeqAIJ() 77 *(bv + rowstart[itmp]) = *(av + k * bs + j); in MatConvert_SeqSBAIJ_SeqAIJ() 82 av += bs2; in MatConvert_SeqSBAIJ_SeqAIJ() 91 *(bv + rowstart[itmp]) = aconj ? PetscConj(*(av + j * bs + k)) : *(av + j * bs + k); in MatConvert_SeqSBAIJ_SeqAIJ() 100 *(bv + rowstart[itmp]) = *(av + k * bs + j); in MatConvert_SeqSBAIJ_SeqAIJ() 105 av += bs2; in MatConvert_SeqSBAIJ_SeqAIJ() 169 MatScalar *av, *bv; in MatConvert_SeqAIJ_SeqSBAIJ() local 208 av = a->a + adiag[i]; in MatConvert_SeqAIJ_SeqSBAIJ() 213 *bv = *av; in MatConvert_SeqAIJ_SeqSBAIJ() [all …]
|
| /petsc/src/binding/petsc4py/test/ ! |
| H A D | test_mat_aij.py | 82 ai, aj, av = self._set_values() 87 ai, aj, av = self._set_values_ijv() 99 ai, aj, av = self._set_values() 104 ai, aj, av = self._set_values_ijv() 128 _, ai, aj, av = self._get_aijv() 129 csr = [ai, aj, av] 146 ai, aj, av = self._set_values() 151 ai, aj, av = self._set_values() 159 ai, aj, av = self._set_values_ijv() 164 ai, aj, av = self._set_values_ijv() [all …]
|
| H A D | test_optdb.py | 107 av = (1, 0, 1) 113 self.opts.setValue('av', tov(av)) 130 self.assertEqual(len(arrayv), len(av)) 131 for v1, v2 in zip(arrayv, av): 137 self.assertEqual(len(arrayv), len(av)) 138 for v1, v2 in zip(arrayv, av):
|
| /petsc/src/mat/impls/h2opus/cuda/ ! |
| H A D | math2opusutils.cu | 43 const PetscScalar *av; in VecSign() local 62 PetscCall(VecCUDAGetArrayRead(v, &av)); in VecSign() 64 SignVector_Functor sign_vector(av, as); in VecSign() 67 PetscCall(VecCUDARestoreArrayRead(v, &av)); in VecSign() 71 PetscCall(VecGetArrayRead(v, &av)); in VecSign() 73 for (i = 0; i < n; i++) as[i] = PetscAbsScalar(av[i]) < 0 ? -1. : 1.; in VecSign() 75 PetscCall(VecRestoreArrayRead(v, &av)); in VecSign()
|
| /petsc/src/mat/impls/aij/seq/umfpack/ ! |
| H A D | umfpack.c | 113 PetscScalar *av = a->a, *xa; in MatSolve_UMFPACK_Private() local 134 …status = umfpack_UMF_wsolve(uflag, ai, aj, (PetscReal *)av, NULL, (PetscReal *)xa, NULL, (PetscRea… 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() 171 PetscScalar *av = a->a; in MatLUFactorNumeric_UMFPACK() local 179 …status = umfpack_UMF_numeric(ai, aj, (double *)av, NULL, lu->Symbolic, &lu->Numeric, lu->Control, … in MatLUFactorNumeric_UMFPACK() 181 status = umfpack_UMF_numeric(ai, aj, av, lu->Symbolic, &lu->Numeric, lu->Control, lu->Info); in MatLUFactorNumeric_UMFPACK() 208 PetscScalar *av = a->a; in MatLUFactorSymbolic_UMFPACK() local 287 …status = umfpack_UMF_qsymbolic(n, m, ai, aj, av, lu->perm_c, &lu->Symbolic, lu->Control, lu->Info); in MatLUFactorSymbolic_UMFPACK() 293 status = umfpack_UMF_symbolic(n, m, ai, aj, av, &lu->Symbolic, lu->Control, lu->Info); in MatLUFactorSymbolic_UMFPACK()
|
| /petsc/src/mat/impls/aij/seq/ ! |
| H A D | symtranspose.c | 63 const MatScalar *aa, *av; in MatTranspose_SeqAIJ() local 76 PetscCall(MatSeqAIJGetArrayRead(A, &av)); in MatTranspose_SeqAIJ() 77 aa = av; in MatTranspose_SeqAIJ() 120 PetscCall(MatSeqAIJRestoreArrayRead(A, &av)); in MatTranspose_SeqAIJ()
|
| /petsc/src/snes/tutorials/output/ ! |
| H A D | ex9_5.out | 2 errors on 5 x 5 grid: av |u-uexact| = 4.531e-02, |u-uexact|_inf = 1.635e-01
|
| H A D | ex9_5_alt.out | 2 errors on 5 x 5 grid: av |u-uexact| = 4.531e-02, |u-uexact|_inf = 1.635e-01
|
| H A D | ex9_4.out | 3 errors on 5 x 5 grid: av |u-uexact| = 4.531e-02, |u-uexact|_inf = 1.635e-01
|
| H A D | ex9_1.out | 5 errors on 9 x 9 grid: av |u-uexact| = 3.076e-03, |u-uexact|_inf = 1.334e-02
|
| H A D | ex9_2.out | 8 errors on 9 x 9 grid: av |u-uexact| = 3.076e-03, |u-uexact|_inf = 1.334e-02
|
| H A D | ex9_mg.out | 5 errors on 33 x 33 grid: av |u-uexact| = 8.182e-04, |u-uexact|_inf = 5.747e-03
|
| H A D | ex9_6.out | 18 errors on 17 x 17 grid: av |u-uexact| = 2.707e-03, |u-uexact|_inf = 1.428e-02
|
| H A D | ex9_3.out | 10 errors on 17 x 17 grid: av |u-uexact| = 2.707e-03, |u-uexact|_inf = 1.428e-02
|
| H A D | ex9_8.out | 18 errors: av |u-uexact| = 1.287e-03 |u-uexact|_inf = 7.882e-03
|
| H A D | ex9_7.out | 18 errors: av |u-uexact| = 1.287e-03 |u-uexact|_inf = 7.882e-03
|
| H A D | ex9_9.out | 21 errors: av |u-uexact| = 1.287e-03 |u-uexact|_inf = 7.882e-03
|
| /petsc/src/mat/impls/aij/mpi/mkl_cpardiso/ ! |
| H A D | mkl_cpardiso.c | 148 const PetscScalar *av, *bv; in MatConvertToTriples_mpiaij_mpiaij_MKL_CPARDISO() local 161 av = aa->a; in MatConvertToTriples_mpiaij_mpiaij_MKL_CPARDISO() 201 val[nz++] = *av++; in MatConvertToTriples_mpiaij_mpiaij_MKL_CPARDISO() 219 const PetscScalar *av, *bv; in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO() local 232 av = aa->a; in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO() 273 PetscCall(PetscArraycpy(val, av, countA * bs2)); in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO() 275 av += countA * bs2; in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO() 292 const PetscScalar *av, *bv; in MatConvertToTriples_mpisbaij_mpisbaij_MKL_CPARDISO() local 304 av = aa->a; in MatConvertToTriples_mpisbaij_mpisbaij_MKL_CPARDISO() 332 PetscCall(PetscArraycpy(val, av, countA * bs2)); in MatConvertToTriples_mpisbaij_mpisbaij_MKL_CPARDISO() [all …]
|
| /petsc/src/mat/impls/dense/seq/ ! |
| H A D | dense.c | 154 const MatScalar *av; in MatConvert_SeqAIJ_SeqDense() local 172 PetscCall(MatSeqAIJGetArrayRead(A, &av)); in MatConvert_SeqAIJ_SeqDense() 176 b->v[*aj * b->lda + i] = *av; in MatConvert_SeqAIJ_SeqDense() 178 av++; in MatConvert_SeqAIJ_SeqDense() 182 PetscCall(MatSeqAIJRestoreArrayRead(A, &av)); in MatConvert_SeqAIJ_SeqDense() 373 const PetscScalar *av; in MatDuplicateNoCreate_SeqDense() local 376 PetscCall(MatDenseGetArrayRead(A, &av)); in MatDuplicateNoCreate_SeqDense() 381 …etscArraycpy(PetscSafePointerPlusOffset(v, j * nlda), PetscSafePointerPlusOffset(av, j * lda), m)); in MatDuplicateNoCreate_SeqDense() 383 PetscCall(PetscArraycpy(v, av, A->rmap->n * A->cmap->n)); in MatDuplicateNoCreate_SeqDense() 386 PetscCall(MatDenseRestoreArrayRead(A, &av)); in MatDuplicateNoCreate_SeqDense() [all …]
|
| /petsc/src/mat/impls/aij/mpi/mumps/ ! |
| H A D | mumps.c | 789 const PetscScalar *av; in MatConvertToTriples_seqaij_seqaij() local 796 PetscCall(MatSeqAIJGetArrayRead(A, &av)); in MatConvertToTriples_seqaij_seqaij() 811 mumps->val = (PetscScalar *)av; in MatConvertToTriples_seqaij_seqaij() 815 …} else if (mumps->nest_vals) PetscCall(PetscArraycpy(mumps->val, av, aa->nz)); /* MatConvertToTrip… in MatConvertToTriples_seqaij_seqaij() 816 …else mumps->val = (PetscScalar *)av; /* in the default c… in MatConvertToTriples_seqaij_seqaij() 817 PetscCall(MatSeqAIJRestoreArrayRead(A, &av)); in MatConvertToTriples_seqaij_seqaij() 961 const PetscScalar *av, *v1; in MatConvertToTriples_seqaij_seqsbaij() local 975 PetscCall(MatSeqAIJGetArrayRead(A, &av)); in MatConvertToTriples_seqaij_seqsbaij() 1011 val[nz] = av[j]; in MatConvertToTriples_seqaij_seqsbaij() 1017 v1 = av + adiag[i]; in MatConvertToTriples_seqaij_seqsbaij() [all …]
|
| /petsc/src/mat/impls/dense/mpi/ ! |
| H A D | mpidense.c | 334 const PetscScalar *av = v + lda * icol[i]; in MatCreateSubMatrix_MPIDense() local 335 for (j = 0; j < nrows; j++) *bv++ = av[irow[j] - rstart]; in MatCreateSubMatrix_MPIDense() 678 const PetscScalar *av; in MatGetDiagonal_MPIDense() local 687 PetscCall(MatDenseGetArrayRead(a->A, &av)); in MatGetDiagonal_MPIDense() 689 for (i = 0; i < len; i++) x[i] = av[radd + i * lda + i]; in MatGetDiagonal_MPIDense() 690 PetscCall(MatDenseRestoreArrayRead(a->A, &av)); in MatGetDiagonal_MPIDense() 1028 const PetscScalar *av, *v; in MatNorm_MPIDense() local 1031 PetscCall(MatDenseGetArrayRead(mdn->A, &av)); in MatNorm_MPIDense() 1032 v = av; in MatNorm_MPIDense() 1050 v = av; in MatNorm_MPIDense() [all …]
|
| /petsc/src/mat/impls/aij/mpi/strumpack/ ! |
| H A D | strumpack.c | 938 const PetscScalar *av; in MatLUFactorNumeric_STRUMPACK() local 955 PetscCall(MatSeqAIJGetArrayRead(Aloc, &av)); in MatLUFactorNumeric_STRUMPACK() 960 …CK_set_distributed_csr_matrix", STRUMPACK_set_distributed_csr_matrix(*S, &m, ai, aj, av, dist, 0)); in MatLUFactorNumeric_STRUMPACK() 962 …StackCallExternalVoid("STRUMPACK_set_csr_matrix", STRUMPACK_set_csr_matrix(*S, &M, ai, aj, av, 0)); in MatLUFactorNumeric_STRUMPACK() 967 PetscCall(MatSeqAIJRestoreArrayRead(Aloc, &av)); in MatLUFactorNumeric_STRUMPACK()
|
| /petsc/src/mat/impls/aij/seq/klu/ ! |
| H A D | klu.c | 145 PetscScalar *av = a->a; in MatLUFactorNumeric_KLU() local 151 lu->Numeric = klu_K_factor(ai, aj, (PetscReal *)av, lu->Symbolic, &lu->Common); in MatLUFactorNumeric_KLU()
|
| /petsc/src/mat/impls/aij/mpi/superlu_dist/ ! |
| H A D | superlu_dist.c | 446 const PetscScalar *av; in MatLUFactorNumeric_SuperLU_DIST() local 469 PetscCall(MatSeqAIJGetArrayRead(Aloc, &av)); in MatLUFactorNumeric_SuperLU_DIST() 525 …for (PetscInt i = 0; i < nz; i++) lu->sval[i] = PetscRealPart(av[i]); /* PetscRealPart() is a no-o… in MatLUFactorNumeric_SuperLU_DIST() 528 PetscCall(PetscArraycpy(lu->val, av, nz)); in MatLUFactorNumeric_SuperLU_DIST() 531 PetscCall(MatSeqAIJRestoreArrayRead(Aloc, &av)); in MatLUFactorNumeric_SuperLU_DIST()
|
| /petsc/src/mat/impls/aij/seq/mkl_pardiso/ ! |
| H A D | mkl_pardiso.c | 195 PetscScalar *av = aav + adiag[i]; in MatMKLPardiso_Convert_seqaij() local 202 *vv = *av; in MatMKLPardiso_Convert_seqaij() 204 av++; in MatMKLPardiso_Convert_seqaij() 221 PetscScalar *av = aav + adiag[i]; in MatMKLPardiso_Convert_seqaij() local 224 *vv = *av; in MatMKLPardiso_Convert_seqaij() 226 av++; in MatMKLPardiso_Convert_seqaij()
|