Home
last modified time | relevance | path

Searched refs:ap (Results 1 – 25 of 26) sorted by relevance

12

/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact3.c11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local
57 ap = aa + j * 9; /* ptr to the beginning of j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_3()
58 for (k = 0; k < 9; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_3()
60 for (k1 = 0; k1 < 3; k1++) *ap++ = dk[k + 3 * k1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
74 ap = aa + jmin * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
78 for (i = 0; i < 9; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
H A Dsbaijfact6.c11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_4() local
58 ap = aa + j * 16; /* ptr to the beginning of j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_4()
59 for (k = 0; k < 16; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_4()
61 for (k1 = 0; k1 < 4; k1++) *ap++ = dk[k + 4 * k1]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
75 ap = aa + jmin * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
79 for (i = 0; i < 16; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
H A Dsbaij.c516 MatScalar *ap, *aa = a->a; in MatGetValues_SeqSBAIJ() local
528 ap = aa + bs2 * ai[brow]; in MatGetValues_SeqSBAIJ()
550 *v++ = ap[bs2 * i + bs * cidx + ridx]; in MatGetValues_SeqSBAIJ()
583 MatScalar *ap, *aa = a->a, *bap; in MatSetValuesBlocked_SeqSBAIJ() local
593 ap = aa + bs2 * ai[row]; in MatSetValuesBlocked_SeqSBAIJ()
621 bap = ap + bs2 * i; in MatSetValuesBlocked_SeqSBAIJ()
648 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqSBAIJ()
653 PetscCall(PetscArraymove(ap + bs2 * (i + 1), ap + bs2 * i, bs2 * (N - i + 1))); in MatSetValuesBlocked_SeqSBAIJ()
654 PetscCall(PetscArrayzero(ap + bs2 * i, bs2)); in MatSetValuesBlocked_SeqSBAIJ()
656 bap = ap + bs2 * i; in MatSetValuesBlocked_SeqSBAIJ()
[all …]
H A Dsbaijfact4.c12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
36 ap = aa + jmin * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
40 for (i = 0; i < 9; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
H A Dsbaijfact.c566 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
615 ap = aa + j * bs2; /* ptr to the beginning of j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_N()
616 for (k = 0; k < bs2; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_N()
618 for (k1 = 0; k1 < bs; k1++) *ap++ = dk[k + bs * k1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
632 ap = aa + jmin * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
636 for (i = 0; i < bs2; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
739 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
764 ap = aa + jmin * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
768 for (i = 0; i < bs2; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
872 MatScalar *ba = b->a, *aa, *ap; in MatCholeskyFactorNumeric_SeqSBAIJ_2() local
[all …]
H A Dsbaijfact7.c11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
57 ap = aa + j * 25; /* ptr to the beginning of j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_5()
58 for (k = 0; k < 25; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_5()
60 for (k1 = 0; k1 < 5; k1++) *ap++ = dk[k + 5 * k1]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
74 ap = aa + jmin * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
78 for (i = 0; i < 25; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
H A Dsbaijfact5.c12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() local
37 ap = aa + jmin * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
41 for (i = 0; i < 16; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
H A Dsbaijfact9.c11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
59 ap = aa + j * 36; /* ptr to the beginning of j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_6()
60 for (k = 0; k < 36; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_6()
62 for (k1 = 0; k1 < 6; k1++) *ap++ = dk[k + 6 * k1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
76 ap = aa + jmin * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
80 for (i = 0; i < 36; i++) *wp++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
H A Dsbaijfact8.c12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
36 ap = aa + jmin * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
40 for (i = 0; i < 25; i++) *rtmp_ptr++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
H A Dsbaijfact11.c11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
60 ap = aa + j * 49; /* ptr to the beginning of j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_7()
61 for (k = 0; k < 49; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */ in MatCholeskyFactorNumeric_SeqSBAIJ_7()
63 for (k1 = 0; k1 < 7; k1++) *ap++ = dk[k + 7 * k1]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
77 ap = aa + jmin * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
81 for (i = 0; i < 49; i++) *wp++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
H A Dsbaijfact12.c12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
39 ap = aa + jmin * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
43 PetscCall(PetscArraycpy(wp, ap, 49)); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
44 ap += 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
H A Dsbaijfact10.c12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
44 ap = aa + jmin * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
48 for (i = 0; i < 36; i++) *wp++ = *ap++; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
/petsc/src/mat/graphops/coarsen/impls/hem/
H A Dhem.c430 MatScalar *ap; in MatCoarsenApply_HEM_private() local
560 ap = PetscSafePointerPlusOffset(matA->a, ii[lid]); in MatCoarsenApply_HEM_private()
564 if ((tt = PetscRealPart(ap[jj])) > threshold && lidj != lid) { in MatCoarsenApply_HEM_private()
572 ap = matB->a + ii[ix]; in MatCoarsenApply_HEM_private()
575 if ((tt = PetscRealPart(ap[jj])) > threshold) { in MatCoarsenApply_HEM_private()
632 ap = PetscSafePointerPlusOffset(matA->a, ii[lid]); in MatCoarsenApply_HEM_private()
636 if ((tt = PetscRealPart(ap[jj])) > threshold && lidj != lid) { in MatCoarsenApply_HEM_private()
649 ap = matB->a + ii[ix]; in MatCoarsenApply_HEM_private()
652 if ((tt = PetscRealPart(ap[jj])) > threshold) { in MatCoarsenApply_HEM_private()
741 ap = matB->a + ii[ix]; in MatCoarsenApply_HEM_private()
[all …]
/petsc/src/mat/impls/blockmat/seq/
H A Dblockmat.c223 Mat *ap, *aa = a->a; in MatSetValues_BlockMat() local
232 ap = aa + ai[brow]; in MatSetValues_BlockMat()
262 …MatSeqXAIJReallocateAIJ(A, a->mbs, 1, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValues_BlockMat()
268 ap[ii + 1] = ap[ii]; in MatSetValues_BlockMat()
270 if (N >= i) ap[i] = NULL; in MatSetValues_BlockMat()
275 if (!*(ap + i)) PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF, bs, bs, 0, NULL, ap + i)); in MatSetValues_BlockMat()
276 PetscCall(MatSetValues(ap[i], 1, &ridx, 1, &cidx, &value, is)); in MatSetValues_BlockMat()
601 Mat *aa = a->a, *ap; in MatAssemblyEnd_BlockMat() local
613 ap = aa + ai[i]; in MatAssemblyEnd_BlockMat()
617 ap[j - fshift] = ap[j]; in MatAssemblyEnd_BlockMat()
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c205 ap = aa + bs2 * ai[brow]; \
221 bap = ap + bs2 * _i + bs * cidx + ridx; \
229 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
233 PetscCall(PetscArraymove(ap + bs2 * (_i + 1), ap + bs2 * _i, bs2 * (N - _i + 1))); \
234 PetscCall(PetscArrayzero(ap + bs2 * _i, bs2)); \
236 ap[bs2 * _i + bs * cidx + ridx] = value; \
245 ap = ba + bs2 * bi[brow]; \
261 bap = ap + bs2 * _i + bs * cidx + ridx; \
269 …MatSeqXAIJReallocateAIJ(B, b->mbs, bs2, nrow, brow, bcol, rmax, ba, bi, bj, rp, ap, bimax, b->none…
273 PetscCall(PetscArraymove(ap + bs2 * (_i + 1), ap + bs2 * _i, bs2 * (N - _i + 1))); \
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaij.c1296 MatScalar *ap, *aa = a->a, *bap; in matsetvaluesblocked4_() local
1304 ap = aa + 16 * ai[row]; in matsetvaluesblocked4_()
1322 bap = ap + 16 * i; in matsetvaluesblocked4_()
1334 PetscCallVoid(PetscArraycpy(ap + 16 * (ii + 1), ap + 16 * (ii), 16)); in matsetvaluesblocked4_()
1336 if (N >= i) PetscCallVoid(PetscArrayzero(ap + 16 * i, 16)); in matsetvaluesblocked4_()
1338 bap = ap + 16 * i; in matsetvaluesblocked4_()
1364 MatScalar *ap, value, *aa = a->a, *bap; in matsetvalues4_() local
1371 ap = aa + 16 * ai[brow]; in matsetvalues4_()
1392 bap = ap + 16 * i + 4 * cidx + ridx; in matsetvalues4_()
1401 PetscCallVoid(PetscArraymove(ap + 16 * i + 16, ap + 16 * i, 16 * (N - i + 1))); in matsetvalues4_()
[all …]
/petsc/src/ksp/pc/impls/tfs/
H A Divec.c170 va_list ap; in PCTFS_ivec_non_uniform() local
173 va_start(ap, n); in PCTFS_ivec_non_uniform()
174 arg3 = va_arg(ap, PetscInt *); in PCTFS_ivec_non_uniform()
175 va_end(ap); in PCTFS_ivec_non_uniform()
/petsc/src/mat/impls/aij/seq/
H A Daij.c367 MatScalar *ap, value, *aa; in MatSeqAIJSetValuesLocalFast() local
374 ap = aa + ai[row]; in MatSeqAIJSetValuesLocalFast()
392 ap[i] += value; in MatSeqAIJSetValuesLocalFast()
408 MatScalar *ap = NULL, value = 0.0, *aa; in MatSetValues_SeqAIJ() local
419 if (!A->structure_only) ap = PetscSafePointerPlusOffset(aa, ai[row]); in MatSetValues_SeqAIJ()
444 ap[i] += value; in MatSetValues_SeqAIJ()
446 } else ap[i] = value; in MatSetValues_SeqAIJ()
458 …MatSeqXAIJReallocateAIJ(A, A->rmap->n, 1, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, M… in MatSetValues_SeqAIJ()
467 PetscCall(PetscArraymove(ap + i + 1, ap + i, N - i + 1)); in MatSetValues_SeqAIJ()
468 ap[i] = value; in MatSetValues_SeqAIJ()
[all …]
/petsc/src/mat/impls/baij/mpi/
H A Dmpibaij.c164 ap = PetscSafePointerPlusOffset(aa, bs2 * ai[brow]); \
180 bap = ap + bs2 * _i + bs * cidx + ridx; \
188 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
192 PetscCall(PetscArraymove(ap + bs2 * (_i + 1), ap + bs2 * _i, bs2 * (N - _i + 1))); \
193 PetscCall(PetscArrayzero(ap + bs2 * _i, bs2)); \
195 ap[bs2 * _i + bs * cidx + ridx] = value; \
204 ap = PetscSafePointerPlusOffset(ba, bs2 * bi[brow]); \
220 bap = ap + bs2 * _i + bs * cidx + ridx; \
228 …MatSeqXAIJReallocateAIJ(B, b->mbs, bs2, nrow, brow, bcol, rmax, ba, bi, bj, rp, ap, bimax, b->none…
232 PetscCall(PetscArraymove(ap + bs2 * (_i + 1), ap + bs2 * _i, bs2 * (N - _i + 1))); \
[all …]
/petsc/src/mat/impls/aij/mpi/
H A Dmpiptap.c56 Mat_SeqAIJ *ap = (Mat_SeqAIJ *)ptap->AP_loc->data; in MatProductCtxDestroy_MPIAIJ_PtAP() local
57 PetscCall(PetscFree(ap->i)); in MatProductCtxDestroy_MPIAIJ_PtAP()
58 PetscCall(PetscFree2(ap->j, ap->a)); in MatProductCtxDestroy_MPIAIJ_PtAP()
100 Mat_SeqAIJ *ap, *p_loc, *p_oth = NULL, *c_seq; in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable() local
124 ap = (Mat_SeqAIJ *)AP_loc->data; in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
138 api = ap->i; in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
139 apj = ap->j; in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
144 apa = ap->a + api[i]; in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
1861 Mat_SeqAIJ *ap, *p_loc, *p_oth = NULL, *c_seq; in MatPtAPNumeric_MPIAIJ_MPIAIJ() local
1885 ap = (Mat_SeqAIJ *)AP_loc->data; in MatPtAPNumeric_MPIAIJ_MPIAIJ()
[all …]
H A Dmpimatmatmult.c1594 MatProductCtx_APMPI *ap; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ() local
1614 ap = (MatProductCtx_APMPI *)C->product->data; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1615 …PetscCheck(ap, PetscObjectComm((PetscObject)C), PETSC_ERR_ARG_WRONGSTATE, "PtA cannot be computed.… in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1616 …PetscCheck(ap->A_loc, PetscObjectComm((PetscObject)C), PETSC_ERR_ARG_WRONGSTATE, "PtA cannot be re… in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1621 merge = ap->merge; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1634 A_loc = ap->A_loc; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1760 MatProductCtx_APMPI *ap; in MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ() local
1793 PetscCall(PetscNew(&ap)); in MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ()
1798 ap->A_loc = A_loc; in MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ()
2074 C->product->data = ap; in MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ()
[all …]
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DPETSc.pyx445 FILE *fd, const char fmt[], va_list ap, argument
452 CHKERR(PetscVSNPrintf(&cstring[0], stringlen, fmt, &final_pos, ap))
458 CHKERR(PetscVSNPrintf(&cstring[0], stringlen, fmt, &final_pos, ap))
464 CHKERR(PetscVFPrintfDefault(fd, fmt, ap))
/petsc/src/dm/impls/plex/
H A Dplexproject.c129 PetscInt ap; in DMProjectPoint_Func_Private() local
130 PetscCall(DMPlexGetActivePoint(dm, &ap)); in DMProjectPoint_Func_Private()
131 … PetscCall(PetscPrintf(PETSC_COMM_SELF, "Project point %" PetscInt_FMT ", analytic: ref (", ap)); in DMProjectPoint_Func_Private()
/petsc/share/petsc/datafiles/meshes/
H A Dmesh-3d-box-innersphere_bin.msh11 …�jrϐ�¿Z�B���ɿ4���ӮS$(�?.�70+��0����?5����{���?Gr�V�������'ο6������ST�?*�>�ap�? ��1�,�?7����+
/petsc/src/ksp/pc/impls/ml/
H A Dml.c125 …ic int PetscML_matvec(ML_Operator *ML_data, int in_length, double p[], int out_length, double ap[]) in PetscML_matvec() argument
142 PetscCall(VecPlaceArray(ml->y, ap)); in PetscML_matvec()

12