From 7ad971abd32d767c02abdbb26e599ac90e4b9c50 Mon Sep 17 00:00:00 2001
From: Giovanni La Mura <giovanni.lamura@inaf.it>
Date: Mon, 13 Nov 2023 10:57:32 +0100
Subject: [PATCH] Adapt common blocks structure to cluster case
---
src/include/Commons.h | 2 ++
src/libnptm/Commons.cpp | 1 +
2 files changed, 3 insertions(+)
diff --git a/src/include/Commons.h b/src/include/Commons.h
index db9a134d..4aeb5ef0 100644
--- a/src/include/Commons.h
+++ b/src/include/Commons.h
@@ -214,6 +214,8 @@ public:
//! \brief QUESTION: definition?
std::complex<double> *vj0;
//! \brief QUESTION: definition?
+ std::complex<double> *vj;
+ //! \brief QUESTION: definition?
std::complex<double> *vyhj;
//! \brief QUESTION: definition?
std::complex<double> *vyj0;
diff --git a/src/libnptm/Commons.cpp b/src/libnptm/Commons.cpp
index 004ee59e..342381f8 100644
--- a/src/libnptm/Commons.cpp
+++ b/src/libnptm/Commons.cpp
@@ -99,6 +99,7 @@ C1_AddOns::C1_AddOns(C4 *c4) {
nlemt = 2 * c4->nlem;
vh = new complex<double>[(nsph * nsph - 1) * c4->litpo];
vj0 = new complex<double>[nsph * c4->lmtpo];
+ vj = new complex<double>[1]; // QUESTION: is 1 really enough for a general case?
vyhj = new complex<double>[nsph * (nsph - 1) * c4->litpo];
vyj0 = new complex<double>[nsph * c4->lmtpo];
am0m = new complex<double>*[nlemt];
--
GitLab