From 33424eb6446f8450ceaf7c281285405826bbe41b Mon Sep 17 00:00:00 2001
From: Yen-Chen Chen <yen-chen.chen@tum.de>
Date: Thu, 13 Mar 2025 17:19:38 +0100
Subject: [PATCH] Fix bug
---
src/horizontal/mo_lib_divrot.cpp | 56 ++++++++++++++++----------------
1 file changed, 28 insertions(+), 28 deletions(-)
diff --git a/src/horizontal/mo_lib_divrot.cpp b/src/horizontal/mo_lib_divrot.cpp
index 82884a5..b63415b 100644
--- a/src/horizontal/mo_lib_divrot.cpp
+++ b/src/horizontal/mo_lib_divrot.cpp
@@ -408,28 +408,28 @@ void recon_lsq_cell_q_svd(const T *p_cc, const int *lsq_idx_c,
auto lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 4, std::make_pair(0, 9), jb);
p_coeff_view(5, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_pseudoinv_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b_subview), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b_subview), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 3, std::make_pair(0, 9), jb);
p_coeff_view(4, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_pseudoinv_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b_subview), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b_subview), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 2, std::make_pair(0, 9), jb);
p_coeff_view(3, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_pseudoinv_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b_subview), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b_subview), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 1, std::make_pair(0, 9), jb);
p_coeff_view(2, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_pseudoinv_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b_subview), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b_subview), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 0, std::make_pair(0, 9), jb);
p_coeff_view(1, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_pseudoinv_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b_subview), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b_subview), 0.0);
p_coeff_view(0, jc, jk, jb) =
p_cc_view(jc, jk, jb) -
p_coeff_view(1, jc, jk, jb) * lsq_moments_view(jc, jb, 0) -
@@ -733,48 +733,48 @@ void recon_lsq_cell_c_svd(const T *p_cc, const int *lsq_idx_c,
auto lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 8, std::make_pair(0, 9), jb);
p_coeff_view(9, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 7, std::make_pair(0, 9), jb);
p_coeff_view(8, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 6, std::make_pair(0, 9), jb);
p_coeff_view(7, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 5, std::make_pair(0, 9), jb);
p_coeff_view(6, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 4, std::make_pair(0, 9), jb);
p_coeff_view(5, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 3, std::make_pair(0, 9), jb);
p_coeff_view(4, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 2, std::make_pair(0, 9), jb);
p_coeff_view(3, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 1, std::make_pair(0, 9), jb);
p_coeff_view(2, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
lsq_pseudoinv_subview =
subview(lsq_pseudoinv_view, jc, 0, std::make_pair(0, 9), jb);
p_coeff_view(1, jc, jk, jb) = std::inner_product(
- KE::cbegin(lsq_qtmat_c_subview), KE::cend(lsq_qtmat_c_subview),
- KE::cbegin(z_b), 0.0);
+ KE::cbegin(lsq_pseudoinv_subview),
+ KE::cend(lsq_pseudoinv_subview), KE::cbegin(z_b), 0.0);
p_coeff_view(0, jc, jk, jb) =
p_cc_view(jc, jk, jb) -
p_coeff_view(1, jc, jk, jb) * lsq_moments_view(jc, jb, 0) -
--
GitLab