1 // Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
4 // SPDX-License-Identifier: BSD-2-Clause
8 #include <ceed-impl.h>
16 /// ----------------------------------------------------------------------------
18 /// ----------------------------------------------------------------------------
27   @return An error code: 0 - success, otherwise - failure
38 /// ----------------------------------------------------------------------------
40 /// ----------------------------------------------------------------------------
50   @return An error code: 0 - success, otherwise - failure
55   if (!ceed->TensorContractCreate) {  in CeedTensorContractCreate()
66   CeedCall(CeedObjectCreate(ceed, NULL, CeedTensorContractDestroy_Object, &(*contract)->obj));  in CeedTensorContractCreate()
67   CeedCall(ceed->TensorContractCreate(*contract));  in CeedTensorContractCreate()
81   @param[in]  B        Middle index of `u`, one index of `t`
83   @param[in]  J        Middle index of `v`, one index of `t`
84   @param[in]  t        Tensor array to contract against
85 …@param[in]  t_mode   Transpose mode for `t`, @ref CEED_NOTRANSPOSE for `t_jb` @ref CEED_TRANSPOSE …
90   @return An error code: 0 - success, otherwise - failure
94 …dTensorContract contract, CeedInt A, CeedInt B, CeedInt C, CeedInt J, const CeedScalar *restrict t,  in CeedTensorContractApply()  argument
96   CeedCall(contract->Apply(contract, A, B, C, J, t, t_mode, add, u, v));  in CeedTensorContractApply()
110   @param[in]  B        Middle index of `u`, one of last two indices of `t`
112   @param[in]  D        First index of `v`, first index of `t`
113   @param[in]  J        Third index of `v`, one of last two indices of `t`
114   @param[in]  t        Tensor array to contract against
115 …@param[in]  t_mode   Transpose mode for `t`, @ref CEED_NOTRANSPOSE for `t_djb` @ref CEED_TRANSPOSE…
120   @return An error code: 0 - success, otherwise - failure
124 …ract contract, CeedInt A, CeedInt B, CeedInt C, CeedInt D, CeedInt J, const CeedScalar *restrict t,  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()
144   @return An error code: 0 - success, otherwise - failure
170   @return An error code: 0 - success, otherwise - failure
175   *(void **)data = contract->data;  in CeedTensorContractGetData()
185   @return An error code: 0 - success, otherwise - failure
190   contract->data = data;  in CeedTensorContractSetData()
199   @return An error code: 0 - success, otherwise - failure
213 …Note: If the value of `*tensor_copy` passed to this function is non-`NULL`, then it is assumed tha…
219   @return An error code: 0 - success, otherwise - failure
235   @return An error code: 0 - success, otherwise - failure
244   if ((*contract)->Destroy) {  in CeedTensorContractDestroy()
245     CeedCall((*contract)->Destroy(*contract));  in CeedTensorContractDestroy()
247   CeedCall(CeedObjectDestroy_Private(&(*contract)->obj));  in CeedTensorContractDestroy()