diff --git a/include/Directions.h b/include/Directions.h index 1534eb8527491a132dbba5d47c7c14ed8ef118f2..214f218146536caaa2da47330cd40908012465ac 100644 --- a/include/Directions.h +++ b/include/Directions.h @@ -184,9 +184,9 @@ void radec2lmn::operator()(Arch, const Direction &reference, const C &ra, double sin_dec = std::sin(dec[i]); double cos_dec = std::cos(dec[i]); - lmn(i, 0) = cos_dec * sin_delta_ra; - lmn(i, 1) = sin_dec * cos_dec0 - cos_dec * sin_dec0 * cos_delta_ra; - lmn(i, 2) = sin_dec * sin_dec0 + cos_dec * cos_dec0 * cos_delta_ra; + lmn(0, i) = cos_dec * sin_delta_ra; + lmn(1, i) = sin_dec * cos_dec0 - cos_dec * sin_dec0 * cos_delta_ra; + lmn(2, i) = sin_dec * sin_dec0 + cos_dec * cos_dec0 * cos_delta_ra; } } @@ -197,7 +197,7 @@ inline void Directions::radec2lmn<Directions::computation_strategy::MULTI_SIMD>( const Direction &reference, xt::xtensor<double, 2> &lmn) { xt::xtensor<double, 2> lmn_tmp({3, ra.size()}); xsimd::dispatch(xsimd::radec2lmn{})(reference, ra, dec, lmn_tmp, - xsimd::aligned_mode()); + xsimd::unaligned_mode()); lmn = xt::transpose(lmn_tmp); } #endif \ No newline at end of file