From 39416ab47e61bf001e0a5c60e02e95616e94985e Mon Sep 17 00:00:00 2001
From: Giovanni La Mura <giovanni.lamura@inaf.it>
Date: Wed, 30 Apr 2025 17:17:59 +0200
Subject: [PATCH] Simplify LM recovery formula
---
src/libnptm/tfrfme.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/libnptm/tfrfme.cpp b/src/libnptm/tfrfme.cpp
index 37b872e..2b81215 100644
--- a/src/libnptm/tfrfme.cpp
+++ b/src/libnptm/tfrfme.cpp
@@ -81,7 +81,7 @@ Swap1* Swap1::from_hdf5(const std::string& file_name) {
if (status == 0) {
status = hdf_file->read("NLMMT", "INT32", &_nlmmt);
status = hdf_file->read("NKV", "INT32", &_nkv);
- lm = (int)((-2.0 + sqrt(4.0 + 2.0 * _nlmmt)) / 2.0);
+ lm = (int)(sqrt(4.0 + 2.0 * _nlmmt) / 2.0) - 1;
num_elements = 2 * _nlmmt * _nkv * _nkv;
instance = new Swap1(lm, _nkv);
elements = new double[num_elements]();
@@ -107,7 +107,7 @@ Swap1* Swap1::from_legacy(const std::string& file_name) {
input.open(file_name.c_str(), ios::in | ios::binary);
if (input.is_open()) {
input.read(reinterpret_cast<char *>(&_nlmmt), sizeof(int));
- lm = (int)((-2.0 + sqrt(4.0 + 2.0 * _nlmmt)) / 2.0);
+ lm = (int)(sqrt(4.0 + 2.0 * _nlmmt) / 2.0) - 1;
input.read(reinterpret_cast<char *>(&_nkv), sizeof(int));
instance = new Swap1(lm, _nkv);
int num_elements = _nlmmt * _nkv * _nkv;
--
GitLab