From 4be85f376d4939f8e5534b7985917e4cfccedfaf Mon Sep 17 00:00:00 2001 From: Michael B Kuhn <31661049+mbkuhn@users.noreply.github.com> Date: Mon, 22 Jul 2024 16:23:16 -0600 Subject: [PATCH] fix tensor identity usage (#1150) --- amr-wind/overset/overset_ops_K.H | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/amr-wind/overset/overset_ops_K.H b/amr-wind/overset/overset_ops_K.H index 143731ad36..bd8780d3f0 100644 --- a/amr-wind/overset/overset_ops_K.H +++ b/amr-wind/overset/overset_ops_K.H @@ -168,7 +168,7 @@ vs::Tensor AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE gp_flux_tensor( fz(i, j - 1, k, 8) + fz(i - 1, j - 1, k, 8) + tiny; - auto f_otimes_gradp = vs::Tensor::I(); + auto f_otimes_gradp = vs::Tensor::identity(); // Average fluxes (faces) to nodes where pressure exists f_otimes_gradp.xx() = (fx(i, j, k, 5) + fx(i, j, k - 1, 5) + fx(i, j - 1, k, 5) + fx(i, j - 1, k - 1, 5)) / @@ -249,12 +249,12 @@ vs::Tensor AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE normal_reinit_tensor( n_zeta.normalize(); // To do outer product between vectors, set up tensors - auto n_tensor = vs::Tensor::I(); - auto n_tensor_T = vs::Tensor::I(); + auto n_tensor = vs::Tensor::identity(); + auto n_tensor_T = vs::Tensor::identity(); n_tensor.rows(n_zeta, n_zeta, n_zeta); n_tensor_T.cols(n_zeta, n_zeta, n_zeta); // Multiply tensors elementwise - auto n_otimes_n = vs::Tensor::I(); + auto n_otimes_n = vs::Tensor::identity(); n_otimes_n.rows( n_tensor.x() * n_tensor_T.x(), n_tensor.y() * n_tensor_T.y(), n_tensor.z() * n_tensor_T.z());