Home
last modified time | relevance | path

Searched refs:A (Results 1 – 25 of 69) sorted by relevance

123

/libCEED/include/ceed/jit-source/gallery/
H A Dceed-poisson3dbuild.h26 CeedScalar A[3][3]; in Poisson3DBuild() local
31 A[k][j] = J[(k + 1) % dim][(j + 1) % dim][i] * J[(k + 2) % dim][(j + 2) % dim][i] - in Poisson3DBuild()
35 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in Poisson3DBuild()
42 q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in Poisson3DBuild()
43 q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in Poisson3DBuild()
44 q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in Poisson3DBuild()
45 q_data[3][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in Poisson3DBuild()
46 q_data[4][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in Poisson3DBuild()
47 q_data[5][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in Poisson3DBuild()
/libCEED/examples/petsc/qfunctions/bps/
H A Dbp3.h49 CeedScalar A[3][3]; in SetupDiffGeo() local
54 A[k][j] = J[(k + 1) % dim][(j + 1) % dim][i] * J[(k + 2) % dim][(j + 2) % dim][i] - in SetupDiffGeo()
58 const CeedScalar detJ = J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]; in SetupDiffGeo()
62 qd[i + Q * 1] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in SetupDiffGeo()
63 qd[i + Q * 2] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in SetupDiffGeo()
64 qd[i + Q * 3] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in SetupDiffGeo()
65 qd[i + Q * 4] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in SetupDiffGeo()
66 qd[i + Q * 5] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in SetupDiffGeo()
67 qd[i + Q * 6] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in SetupDiffGeo()
H A Dbp1.h37 CeedScalar A[3][3]; in SetupMassGeo() local
42 A[k][j] = J[(k + 1) % dim][(j + 1) % dim][i] * J[(k + 2) % dim][(j + 2) % dim][i] - in SetupMassGeo()
46 const CeedScalar detJ = J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]; in SetupMassGeo()
/libCEED/examples/fluids/qfunctions/
H A Dutils_eigensolver_jacobi.h34 CEED_QFUNCTION_HELPER CeedInt MaxEntryRow(const CeedScalar *A, CeedInt N, CeedInt i) { in MaxEntryRow() argument
37 if (fabs(A[i * N + j]) > fabs(A[i * N + j_max])) j_max = j; in MaxEntryRow()
48 CEED_QFUNCTION_HELPER void MaxEntry(const CeedScalar *A, CeedInt N, CeedInt *max_idx_row, CeedInt *… in MaxEntry() argument
51 CeedScalar max_entry = fabs(A[*i_max * N + *j_max]); in MaxEntry()
54 if (fabs(A[i * N + j]) > max_entry) { in MaxEntry()
55 max_entry = fabs(A[i * N + j]); in MaxEntry()
71 CEED_QFUNCTION_HELPER void CalcRot(const CeedScalar *A, CeedInt N, CeedInt i, CeedInt j, CeedScalar… in CalcRot() argument
73 CeedScalar A_jj_ii = (A[j * N + j] - A[i * N + i]); in CalcRot()
78 CeedScalar A_ij = A[i * N + j]; in CalcRot()
159 CEED_QFUNCTION_HELPER void ApplyRot(CeedScalar *A, CeedInt N, CeedInt i, CeedInt j, CeedInt *max_id… in ApplyRot() argument
[all …]
H A Dutils.h44 …FUNCTION_HELPER void CopyMat3(const CeedScalar A[3][3], CeedScalar B[3][3]) { CopyN((const CeedSca… in CopyMat3()
71 CEED_QFUNCTION_HELPER void MatVecNM(const CeedScalar *A, const CeedScalar *x, const CeedInt N, cons… in MatVecNM() argument
75 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) b[i] += DotN(&A[i * M], x, M); in MatVecNM()
78 … i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) b[i] += A[j * M + i] * x[j]; } in MatVecNM()
84 CEED_QFUNCTION_HELPER void MatVec3(const CeedScalar A[3][3], const CeedScalar x[3], const CeedTrans… in MatVec3()
85 MatVecNM((const CeedScalar *)A, (const CeedScalar *)x, 3, 3, transpose_A, (CeedScalar *)b); in MatVec3()
90 CEED_QFUNCTION_HELPER void MatDiagNM(const CeedScalar *A, const CeedScalar *D, const CeedInt N, con… in MatDiagNM() argument
94 …i < N; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < M; j++) B[i * M + j] += D[i] * A[i * M + j]; } in MatDiagNM()
97 …i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) B[i * N + j] += D[i] * A[j * M + i]; } in MatDiagNM()
104 CEED_QFUNCTION_HELPER void MatDiag3(const CeedScalar A[3][3], const CeedScalar D[3], const CeedTran… in MatDiag3()
[all …]
/libCEED/examples/python/qfunctions/
H A Dex2-surface.h50 CeedScalar A[3][3]; in build_diff() local
56 A[k][j] = in build_diff()
62 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_diff()
69 q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_diff()
70 q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_diff()
71 q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in build_diff()
72 q_data[3][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in build_diff()
73 q_data[4][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in build_diff()
74 q_data[5][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in build_diff()
H A Dex3-volume.h60 CeedScalar A[3][3]; in build_mass_diff() local
63 A[k][j] = in build_mass_diff()
69 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
72 q_data[0][i] = w[i] * (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
79 q_data[1][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_mass_diff()
80 q_data[2][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_mass_diff()
81 q_data[3][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in build_mass_diff()
82 q_data[4][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in build_mass_diff()
83 q_data[5][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in build_mass_diff()
84 q_data[6][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in build_mass_diff()
/libCEED/examples/ceed/
H A Dex2-surface.h54 CeedScalar A[3][3]; in build_diff() local
60 A[k][j] = in build_diff()
66 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_diff()
73 q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_diff()
74 q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_diff()
75 q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in build_diff()
76 q_data[3][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in build_diff()
77 q_data[4][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in build_diff()
78 q_data[5][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in build_diff()
H A Dex3-volume.h64 CeedScalar A[3][3]; in build_mass_diff() local
70 A[k][j] = in build_mass_diff()
76 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
79 q_data[0][i] = w[i] * (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
86 q_data[1][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_mass_diff()
87 q_data[2][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_mass_diff()
88 q_data[3][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in build_mass_diff()
89 q_data[4][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in build_mass_diff()
90 q_data[5][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in build_mass_diff()
91 q_data[6][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in build_mass_diff()
/libCEED/tests/
H A Dt580-operator.h11 …NCTION_HELPER CeedScalar MatDet2x2(const CeedScalar A[2][2]) { return A[0][0] * A[1][1] - A[1][0] … in MatDet2x2()
14 CEED_QFUNCTION_HELPER int AlphaMatTransposeMatMult2x2(const CeedScalar alpha, const CeedScalar A[2]… in AlphaMatTransposeMatMult2x2()
20 C[j][k] += alpha * A[m][j] * B[m][k]; in AlphaMatTransposeMatMult2x2()
H A Dt301-basis.c13 CeedScalar A[12] = {1, -1, 4, 1, 4, -2, 1, 4, 2, 1, -1, 0}; in main() local
27 if (fabs(A_qr[i] - A[i]) > 100. * CEED_EPSILON) { in main()
29 … QR factorization A_qr[%" CeedInt_FMT "] = %f != A[%" CeedInt_FMT "] = %f\n", i, A_qr[i], i, A[i]); in main()
/libCEED/backends/opt/
H A Dceed-opt-tensor.c16 static inline int CeedTensorContractApply_Core_Opt(CeedTensorContract contract, CeedInt A, CeedInt … in CeedTensorContractApply_Core_Opt() argument
26 for (CeedInt a = 0; a < A; a++) { in CeedTensorContractApply_Core_Opt()
40 static int CeedTensorContractApply_Opt(CeedTensorContract contract, CeedInt A, CeedInt B, CeedInt C… in CeedTensorContractApply_Opt() argument
43 for (CeedInt q = 0; q < A * J * C; q++) v[q] = (CeedScalar)0.0; in CeedTensorContractApply_Opt()
46 if (C == 1) return CeedTensorContractApply_Core_Opt(contract, A, B, 1, J, t, t_mode, add, u, v); in CeedTensorContractApply_Opt()
47 else return CeedTensorContractApply_Core_Opt(contract, A, B, C, J, t, t_mode, add, u, v); in CeedTensorContractApply_Opt()
/libCEED/backends/avx/
H A Dceed-avx-tensor.c42 static inline int CeedTensorContract_Avx_Blocked(CeedTensorContract contract, CeedInt A, CeedInt B,… in CeedTensorContract_Avx_Blocked() argument
52 for (CeedInt a = 0; a < A; a++) { in CeedTensorContract_Avx_Blocked()
104 static inline int CeedTensorContract_Avx_Remainder(CeedTensorContract contract, CeedInt A, CeedInt … in CeedTensorContract_Avx_Remainder() argument
116 for (CeedInt a = 0; a < A; a++) { in CeedTensorContract_Avx_Remainder()
153 static inline int CeedTensorContract_Avx_Single(CeedTensorContract contract, CeedInt A, CeedInt B, … in CeedTensorContract_Avx_Single() argument
164 for (CeedInt a = 0; a < (A / AA) * AA; a += AA) { in CeedTensorContract_Avx_Single()
187 const CeedInt a = (A / AA) * AA; in CeedTensorContract_Avx_Single()
192 for (CeedInt aa = 0; aa < A - a; aa++) { in CeedTensorContract_Avx_Single()
200 for (CeedInt aa = 0; aa < A - a; aa++) { // unroll in CeedTensorContract_Avx_Single()
205 for (CeedInt aa = 0; aa < A - a; aa++) { in CeedTensorContract_Avx_Single()
[all …]
/libCEED/interface/
H A Dceed-tensor.c94 int CeedTensorContractApply(CeedTensorContract contract, CeedInt A, CeedInt B, CeedInt C, CeedInt J… in CeedTensorContractApply() argument
96 CeedCall(contract->Apply(contract, A, B, C, J, t, t_mode, add, u, v)); in CeedTensorContractApply()
124 int CeedTensorContractStridedApply(CeedTensorContract contract, CeedInt A, CeedInt B, CeedInt C, Ce… in CeedTensorContractStridedApply() argument
128 … CeedCall(contract->Apply(contract, A, J, C, B, t + d * B * J, t_mode, add, u + d * A * J * C, v)); in CeedTensorContractStridedApply()
132 … CeedCall(contract->Apply(contract, A, B, C, J, t + d * B * J, t_mode, add, u, v + d * A * J * C)); in CeedTensorContractStridedApply()
/libCEED/examples/solids/include/
H A Dmatops.h25 PetscErrorCode ApplyJacobian_Ceed(Mat A, Vec X, Vec Y);
28 PetscErrorCode Prolong_Ceed(Mat A, Vec X, Vec Y);
31 PetscErrorCode Restrict_Ceed(Mat A, Vec X, Vec Y);
34 PetscErrorCode GetDiag_Ceed(Mat A, Vec D);
/libCEED/examples/solids/qfunctions/
H A Dfinite-strain-neo-hookean.h79 CEED_QFUNCTION_HELPER int computeMatinvSym(const CeedScalar A[3][3], const CeedScalar detA, CeedSca… in computeMatinvSym()
82 A[1][1] * A[2][2] - A[1][2] * A[2][1], /* *NOPAD* */ in computeMatinvSym()
83 A[0][0] * A[2][2] - A[0][2] * A[2][0], /* *NOPAD* */ in computeMatinvSym()
84 A[0][0] * A[1][1] - A[0][1] * A[1][0], /* *NOPAD* */ in computeMatinvSym()
85 A[0][2] * A[1][0] - A[0][0] * A[1][2], /* *NOPAD* */ in computeMatinvSym()
86 A[0][1] * A[1][2] - A[0][2] * A[1][1], /* *NOPAD* */ in computeMatinvSym()
87 A[0][2] * A[2][1] - A[0][1] * A[2][2] /* *NOPAD* */ in computeMatinvSym()
H A Dfinite-strain-mooney-rivlin.h84 CEED_QFUNCTION_HELPER int computeMatinvSym(const CeedScalar A[3][3], const CeedScalar detA, CeedSca… in computeMatinvSym()
87 A[1][1] * A[2][2] - A[1][2] * A[2][1], /* *NOPAD* */ in computeMatinvSym()
88 A[0][0] * A[2][2] - A[0][2] * A[2][0], /* *NOPAD* */ in computeMatinvSym()
89 A[0][0] * A[1][1] - A[0][1] * A[1][0], /* *NOPAD* */ in computeMatinvSym()
90 A[0][2] * A[1][0] - A[0][0] * A[1][2], /* *NOPAD* */ in computeMatinvSym()
91 A[0][1] * A[1][2] - A[0][2] * A[1][1], /* *NOPAD* */ in computeMatinvSym()
92 A[0][2] * A[2][1] - A[0][1] * A[2][2] /* *NOPAD* */ in computeMatinvSym()
/libCEED/examples/petsc/include/
H A Dmatops.h20 PetscErrorCode MatGetDiag(Mat A, Vec D);
22 PetscErrorCode MatMult_Ceed(Mat A, Vec X, Vec Y);
24 PetscErrorCode MatMult_Prolong(Mat A, Vec X, Vec Y);
25 PetscErrorCode MatMult_Restrict(Mat A, Vec X, Vec Y);
/libCEED/examples/petsc/src/
H A Dmatops.c41 PetscErrorCode MatGetDiag(Mat A, Vec D) { in MatGetDiag() argument
45 PetscCall(MatShellGetContext(A, &op_apply_ctx)); in MatGetDiag()
93 PetscErrorCode MatMult_Ceed(Mat A, Vec X, Vec Y) { in MatMult_Ceed() argument
97 PetscCall(MatShellGetContext(A, &op_apply_ctx)); in MatMult_Ceed()
107 PetscErrorCode MatMult_Prolong(Mat A, Vec X, Vec Y) { in MatMult_Prolong() argument
112 PetscCall(MatShellGetContext(A, &pr_restr_ctx)); in MatMult_Prolong()
141 PetscErrorCode MatMult_Restrict(Mat A, Vec X, Vec Y) { in MatMult_Restrict() argument
146 PetscCall(MatShellGetContext(A, &pr_restr_ctx)); in MatMult_Restrict()
/libCEED/backends/xsmm/
H A Dceed-xsmm-tensor.c17 static int CeedTensorContractApply_Xsmm(CeedTensorContract contract, CeedInt A, CeedInt B, CeedInt … in CeedTensorContractApply_Xsmm() argument
25 …? libxsmm_create_gemm_shape(J, A, B, !t_mode ? B : J, B, J, LIBXSMM_DATATYPE_F64, LIBXSMM_DATATYPE… in CeedTensorContractApply_Xsmm()
27 …: libxsmm_create_gemm_shape(J, A, B, !t_mode ? B : J, B, J, LIBXSMM_DATATYPE_F32, LIBXSMM_DATATYPE… in CeedTensorContractApply_Xsmm()
56 for (CeedInt a = 0; a < A; a++) { in CeedTensorContractApply_Xsmm()
/libCEED/backends/ref/
H A Dceed-ref-tensor.c16 static int CeedTensorContractApply_Ref(CeedTensorContract contract, CeedInt A, CeedInt B, CeedInt C… in CeedTensorContractApply_Ref() argument
26 for (CeedInt q = 0; q < A * J * C; q++) v[q] = (CeedScalar)0.0; in CeedTensorContractApply_Ref()
29 for (CeedInt a = 0; a < A; a++) { in CeedTensorContractApply_Ref()
/libCEED/examples/fluids/include/
H A Dmat-ceed-impl.h42 PETSC_CEED_EXTERN PetscErrorCode MatGetDiagonal_Ceed(Mat A, Vec D);
43 PETSC_CEED_EXTERN PetscErrorCode MatMult_Ceed(Mat A, Vec X, Vec Y);
44 PETSC_CEED_EXTERN PetscErrorCode MatMultTranspose_Ceed(Mat A, Vec Y, Vec X);
/libCEED/examples/solids/src/
H A Dmatops.c92 PetscErrorCode ApplyJacobian_Ceed(Mat A, Vec X, Vec Y) { in ApplyJacobian_Ceed() argument
98 PetscCall(MatShellGetContext(A, &user)); in ApplyJacobian_Ceed()
108 PetscErrorCode Prolong_Ceed(Mat A, Vec X, Vec Y) { in Prolong_Ceed() argument
115 PetscCall(MatShellGetContext(A, &user)); in Prolong_Ceed()
145 PetscErrorCode Restrict_Ceed(Mat A, Vec X, Vec Y) { in Restrict_Ceed() argument
152 PetscCall(MatShellGetContext(A, &user)); in Restrict_Ceed()
182 PetscErrorCode GetDiag_Ceed(Mat A, Vec D) { in GetDiag_Ceed() argument
187 PetscCall(MatShellGetContext(A, &user)); in GetDiag_Ceed()
/libCEED/doc/sphinx/source/api/
H A DCeedOperator.rst6 A `CeedOperator` defines the finite/spectral element operator associated to a :ref:`CeedQFunction`.
7 A `CeedOperator` connects objects of the type :ref:`CeedElemRestriction`, :ref:`CeedBasis`, and :re…
/libCEED/doc/sphinx/source/
H A DCODE_OF_CONDUCT.md53 **Consequence**: A private, written warning from community leaders, providing clarity around the na…
54 A public apology may be requested.
58 **Community Impact**: A violation through a single incident or series of actions.
60 **Consequence**: A warning with consequences for continued behavior.
67 **Community Impact**: A serious violation of community standards, including sustained inappropriate…
69 **Consequence**: A temporary ban from any sort of interaction or public communication with the comm…
77 **Consequence**: A permanent ban from any sort of public interaction within the community.

123