Home
last modified time | relevance | path

Searched refs:ay (Results 1 – 17 of 17) sorted by relevance

/petsc/src/mat/impls/baij/seq/
H A Ddgefa2.c18 MatScalar *aa, *ax, *ay, work[4], stmp; in PetscKernel_A_gets_inverse_A_2() local
82 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_2()
83 for (ll = 0; ll < i__3; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_2()
110 ay = &a[j3 + 1]; in PetscKernel_A_gets_inverse_A_2()
111 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_2()
127 ay = &a[k3 + 1]; in PetscKernel_A_gets_inverse_A_2()
128 ay[0] += stmp * ax[0]; in PetscKernel_A_gets_inverse_A_2()
129 ay[1] += stmp * ax[1]; in PetscKernel_A_gets_inverse_A_2()
134 ay = &a[2 * l + 1]; in PetscKernel_A_gets_inverse_A_2()
136 ax[0] = ay[0]; in PetscKernel_A_gets_inverse_A_2()
[all …]
H A Ddgefa7.c18 MatScalar *aa, *ax, *ay, work[49], stmp; in PetscKernel_A_gets_inverse_A_7() local
83 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_7()
84 for (ll = 0; ll < i__3; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_7()
111 ay = &a[j3 + 1]; in PetscKernel_A_gets_inverse_A_7()
112 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_7()
129 ay = &a[k3 + 1]; in PetscKernel_A_gets_inverse_A_7()
130 ay[0] += stmp * ax[0]; in PetscKernel_A_gets_inverse_A_7()
131 ay[1] += stmp * ax[1]; in PetscKernel_A_gets_inverse_A_7()
132 ay[2] += stmp * ax[2]; in PetscKernel_A_gets_inverse_A_7()
133 ay[3] += stmp * ax[3]; in PetscKernel_A_gets_inverse_A_7()
[all …]
H A Ddgefa6.c18 MatScalar *aa, *ax, *ay, work[36], stmp; in PetscKernel_A_gets_inverse_A_6() local
83 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_6()
84 for (ll = 0; ll < i__3; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_6()
111 ay = &a[j3 + 1]; in PetscKernel_A_gets_inverse_A_6()
112 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_6()
129 ay = &a[k3 + 1]; in PetscKernel_A_gets_inverse_A_6()
130 ay[0] += stmp * ax[0]; in PetscKernel_A_gets_inverse_A_6()
131 ay[1] += stmp * ax[1]; in PetscKernel_A_gets_inverse_A_6()
132 ay[2] += stmp * ax[2]; in PetscKernel_A_gets_inverse_A_6()
133 ay[3] += stmp * ax[3]; in PetscKernel_A_gets_inverse_A_6()
[all …]
H A Ddgefa5.c18 MatScalar *aa, *ax, *ay, stmp; in PetscKernel_A_gets_inverse_A_5() local
83 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_5()
84 for (ll = 0; ll < i__3; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_5()
111 ay = &a[j3 + 1]; in PetscKernel_A_gets_inverse_A_5()
112 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_5()
129 ay = &a[k3 + 1]; in PetscKernel_A_gets_inverse_A_5()
130 ay[0] += stmp * ax[0]; in PetscKernel_A_gets_inverse_A_5()
131 ay[1] += stmp * ax[1]; in PetscKernel_A_gets_inverse_A_5()
132 ay[2] += stmp * ax[2]; in PetscKernel_A_gets_inverse_A_5()
133 ay[3] += stmp * ax[3]; in PetscKernel_A_gets_inverse_A_5()
[all …]
H A Ddgefa4.c18 MatScalar *aa, *ax, *ay, work[16], stmp; in PetscKernel_A_gets_inverse_A_4() local
83 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_4()
84 for (ll = 0; ll < i__3; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_4()
111 ay = &a[j3 + 1]; in PetscKernel_A_gets_inverse_A_4()
112 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_4()
129 ay = &a[k3 + 1]; in PetscKernel_A_gets_inverse_A_4()
130 ay[0] += stmp * ax[0]; in PetscKernel_A_gets_inverse_A_4()
131 ay[1] += stmp * ax[1]; in PetscKernel_A_gets_inverse_A_4()
132 ay[2] += stmp * ax[2]; in PetscKernel_A_gets_inverse_A_4()
133 ay[3] += stmp * ax[3]; in PetscKernel_A_gets_inverse_A_4()
[all …]
H A Ddgefa3.c18 MatScalar *aa, *ax, *ay, work[9], stmp; in PetscKernel_A_gets_inverse_A_3() local
83 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_3()
84 for (ll = 0; ll < i__3; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_3()
111 ay = &a[j3 + 1]; in PetscKernel_A_gets_inverse_A_3()
112 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_3()
129 ay = &a[k3 + 1]; in PetscKernel_A_gets_inverse_A_3()
130 ay[0] += stmp * ax[0]; in PetscKernel_A_gets_inverse_A_3()
131 ay[1] += stmp * ax[1]; in PetscKernel_A_gets_inverse_A_3()
132 ay[2] += stmp * ax[2]; in PetscKernel_A_gets_inverse_A_3()
137 ay = &a[3 * l + 1]; in PetscKernel_A_gets_inverse_A_3()
[all …]
H A Ddgedi.c20 MatScalar *aa, *ax, *ay, tmp; in PetscLINPACKgedi() local
45 ay = &a[1 + jn1]; in PetscLINPACKgedi()
46 for (ll = 0; ll < k; ll++) ay[ll] += t * ax[ll]; in PetscLINPACKgedi()
66 ay = &a[kn + 1]; in PetscLINPACKgedi()
67 for (ll = 0; ll < n; ll++) ay[ll] += t * ax[ll]; in PetscLINPACKgedi()
72 ay = &a[l * n + 1]; in PetscLINPACKgedi()
75 ax[ll] = ay[ll]; in PetscLINPACKgedi()
76 ay[ll] = tmp; in PetscLINPACKgedi()
H A Ddgefa.c18 MatScalar t, *ax, *ay, *aa; in PetscLINPACKgefa() local
81 ay = &a[1 + k + jn1]; in PetscLINPACKgefa()
82 for (ll = 0; ll < i__3; ll++) ay[ll] += t * ax[ll]; in PetscLINPACKgefa()
/petsc/src/mat/tests/
H A Dex229.c8 PetscScalar *ay; in myF() local
13 PetscCall(VecGetArray(y, &ay)); in myF()
24 ay[i] = yy; in myF()
26 PetscCall(VecRestoreArray(y, &ay)); in myF()
/petsc/src/dm/impls/stag/
H A Dstagmulti.c94 static PetscErrorCode SetInterpolationWeight2d_Private(PetscScalar ax, PetscScalar ay, PetscScalar … in SetInterpolationWeight2d_Private() argument
97 weight[0] = (1.0 - ax) * (1.0 - ay); in SetInterpolationWeight2d_Private()
98 weight[1] = ax * (1.0 - ay); in SetInterpolationWeight2d_Private()
99 weight[2] = (1.0 - ax) * ay; in SetInterpolationWeight2d_Private()
100 weight[3] = ax * ay; in SetInterpolationWeight2d_Private()
140 static PetscErrorCode SetInterpolationWeight3d_Private(PetscScalar ax, PetscScalar ay, PetscScalar … in SetInterpolationWeight3d_Private() argument
143 weight[0] = (1.0 - ax) * (1.0 - ay) * (1.0 - az); in SetInterpolationWeight3d_Private()
144 weight[1] = ax * (1.0 - ay) * (1.0 - az); in SetInterpolationWeight3d_Private()
145 weight[2] = (1.0 - ax) * ay * (1.0 - az); in SetInterpolationWeight3d_Private()
146 weight[3] = ax * ay * (1.0 - az); in SetInterpolationWeight3d_Private()
[all …]
/petsc/include/
H A Dpetscdmstag.h147 …tatic inline PetscErrorCode DMStagGet1dCoordinateArraysDOFRead(DM dm, void *ax, void *ay, void *az) in DMStagGet1dCoordinateArraysDOFRead() argument
149 return DMStagGetProductCoordinateArraysRead(dm, ax, ay, az); in DMStagGet1dCoordinateArraysDOFRead()
159 …c inline PetscErrorCode DMStagRestore1dCoordinateArraysDOFRead(DM dm, void *ax, void *ay, void *az) in DMStagRestore1dCoordinateArraysDOFRead() argument
161 return DMStagRestoreProductCoordinateArraysRead(dm, ax, ay, az); in DMStagRestore1dCoordinateArraysDOFRead()
/petsc/src/mat/impls/dense/mpi/
H A Dmpidense.c477 PetscScalar *ay; in MatMultColumnRange_MPIDense() local
483 PetscCall(VecGetArrayWriteAndMemType(mdn->lvec, &ay, &aymtype)); in MatMultColumnRange_MPIDense()
484 …PetscCall(PetscSFBcastWithMemTypeBegin(mdn->Mvctx, MPIU_SCALAR, axmtype, ax, aymtype, ay, MPI_REPL… in MatMultColumnRange_MPIDense()
485 PetscCall(PetscSFBcastEnd(mdn->Mvctx, MPIU_SCALAR, ax, ay, MPI_REPLACE)); in MatMultColumnRange_MPIDense()
486 PetscCall(VecRestoreArrayWriteAndMemType(mdn->lvec, &ay)); in MatMultColumnRange_MPIDense()
496 PetscScalar *ay; in MatMult_MPIDense() local
502 PetscCall(VecGetArrayWriteAndMemType(mdn->lvec, &ay, &aymtype)); in MatMult_MPIDense()
503 …PetscCall(PetscSFBcastWithMemTypeBegin(mdn->Mvctx, MPIU_SCALAR, axmtype, ax, aymtype, ay, MPI_REPL… in MatMult_MPIDense()
504 PetscCall(PetscSFBcastEnd(mdn->Mvctx, MPIU_SCALAR, ax, ay, MPI_REPLACE)); in MatMult_MPIDense()
505 PetscCall(VecRestoreArrayWriteAndMemType(mdn->lvec, &ay)); in MatMult_MPIDense()
[all …]
/petsc/src/mat/impls/aij/seq/seqhipsparse/
H A Daijhipsparse.hip.cxx3264 PetscScalar *ay; in MatAXPY_SeqAIJHIPSPARSE() local
3299 PetscCall(MatSeqAIJHIPSPARSEGetArray(Y, &ay)); in MatAXPY_SeqAIJHIPSPARSE()
3302 …ta().get(), csrx->column_indices->data().get(), &b, cy->mat->descr, y->nz, ay, csry->row_offsets->… in MatAXPY_SeqAIJHIPSPARSE()
3303 …csry->column_indices->data().get(), cy->mat->descr, ay, csry->row_offsets->data().get(), csry->col… in MatAXPY_SeqAIJHIPSPARSE()
3306 …ta().get(), csrx->column_indices->data().get(), &b, cy->mat->descr, y->nz, ay, csry->row_offsets->… in MatAXPY_SeqAIJHIPSPARSE()
3307 …csry->column_indices->data().get(), cy->mat->descr, ay, csry->row_offsets->data().get(), csry->col… in MatAXPY_SeqAIJHIPSPARSE()
3313 …ta().get(), csrx->column_indices->data().get(), &b, cy->mat->descr, y->nz, ay, csry->row_offsets->… in MatAXPY_SeqAIJHIPSPARSE()
3314 …csry->column_indices->data().get(), cy->mat->descr, ay, csry->row_offsets->data().get(), csry->col… in MatAXPY_SeqAIJHIPSPARSE()
3320 PetscCall(MatSeqAIJHIPSPARSERestoreArray(Y, &ay)); in MatAXPY_SeqAIJHIPSPARSE()
3326 PetscCall(MatSeqAIJHIPSPARSEGetArray(Y, &ay)); in MatAXPY_SeqAIJHIPSPARSE()
[all …]
/petsc/src/vec/vec/impls/seq/cupm/
H A Dvecseqcupm_impl.hpp1724 const auto ay = PetscAbsScalar(y); in operator ()() local
1725 const auto skip = au < this->ignore_max_ || ay < this->ignore_max_; in operator ()()
1727 const auto tolr = skip ? 0.0 : PetscRealPart(thrust::get<3>(x)) * PetscMax(au, ay); in operator ()()
1749 const auto ay = PetscAbsScalar(thrust::get<1>(x)); in operator ()() local
1750 const auto skip = au < this->ignore_max_ || ay < this->ignore_max_; in operator ()()
1752 const auto tolr = skip ? 0.0 : PetscRealPart(thrust::get<4>(x)) * PetscMax(au, ay); in operator ()()
/petsc/src/mat/impls/aij/seq/seqcusparse/
H A Daijcusparse.cu3871 PetscScalar *ay; in MatAXPY_SeqAIJCUSPARSE() local
3907 PetscCall(MatSeqAIJCUSPARSEGetArray(Y, &ay)); in MatAXPY_SeqAIJCUSPARSE()
3910 …ta().get(), csrx->column_indices->data().get(), &b, cy->mat->descr, y->nz, ay, csry->row_offsets->… in MatAXPY_SeqAIJCUSPARSE()
3911 …csry->column_indices->data().get(), cy->mat->descr, ay, csry->row_offsets->data().get(), csry->col… in MatAXPY_SeqAIJCUSPARSE()
3914 …ta().get(), csrx->column_indices->data().get(), &b, cy->mat->descr, y->nz, ay, csry->row_offsets->… in MatAXPY_SeqAIJCUSPARSE()
3915 …csry->column_indices->data().get(), cy->mat->descr, ay, csry->row_offsets->data().get(), csry->col… in MatAXPY_SeqAIJCUSPARSE()
3921 …ta().get(), csrx->column_indices->data().get(), &b, cy->mat->descr, y->nz, ay, csry->row_offsets->… in MatAXPY_SeqAIJCUSPARSE()
3922 …csry->column_indices->data().get(), cy->mat->descr, ay, csry->row_offsets->data().get(), csry->col… in MatAXPY_SeqAIJCUSPARSE()
3928 PetscCall(MatSeqAIJCUSPARSERestoreArray(Y, &ay)); in MatAXPY_SeqAIJCUSPARSE()
3934 PetscCall(MatSeqAIJCUSPARSEGetArray(Y, &ay)); in MatAXPY_SeqAIJCUSPARSE()
[all …]
/petsc/src/ksp/pc/impls/bddc/
H A Dbddcprivate.c2218 PetscScalar *ay; in PCBDDCBenignMatMult_Private_Private() local
2221 PetscCall(VecGetArray(y, &ay)); in PCBDDCBenignMatMult_Private_Private()
2228 val = -ay[idxs[nz - 1]]; in PCBDDCBenignMatMult_Private_Private()
2232 sum += ay[idxs[j]]; in PCBDDCBenignMatMult_Private_Private()
2233 ay[idxs[j]] += val; in PCBDDCBenignMatMult_Private_Private()
2235 ay[idxs[nz - 1]] += sum; in PCBDDCBenignMatMult_Private_Private()
2237 for (j = 0; j < nz - 1; j++) ay[idxs[j]] += val; in PCBDDCBenignMatMult_Private_Private()
2238 ay[idxs[nz - 1]] = 0.; in PCBDDCBenignMatMult_Private_Private()
2242 PetscCall(VecRestoreArray(y, &ay)); in PCBDDCBenignMatMult_Private_Private()
/petsc/doc/
H A Dpetsc.bib27582 title = {Predicting Salinity in the {C}hesapeake {B}ay Using Backpropagation},