From 5dac1d6a2232c7c615ce37c53bb3df446d2586f2 Mon Sep 17 00:00:00 2001
From: Giovanni La Mura <giovanni.lamura@inaf.it>
Date: Wed, 19 Feb 2025 12:31:33 +0100
Subject: [PATCH] Reconstruct angle vectors when building InclusionOutputInfo
 from HDF5

---
 src/libnptm/outputs.cpp | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/libnptm/outputs.cpp b/src/libnptm/outputs.cpp
index 58f40cd5..68e0e62b 100644
--- a/src/libnptm/outputs.cpp
+++ b/src/libnptm/outputs.cpp
@@ -3065,6 +3065,10 @@ InclusionOutputInfo::InclusionOutputInfo(const std::string &hdf5_name) {
     vec_fsas12 = new dcomplex[xi_block_size];
     status = hdf_file->read("VEC_FSAS12", str_type, vec_fsas12);
     // Initialize directions (they are scale-independent)
+    vec_dir_tidg = new double[_num_theta];
+    vec_dir_tsdg = new double[_num_thetas];
+    vec_dir_pidg = new double[_num_phi];
+    vec_dir_psdg = new double[_num_phis];
     double cti = th, cpi = ph, cts = ths, cps = phs;
     for (int di = 0; di < _num_theta; di++) {
       vec_dir_tidg[di] = cti;
-- 
GitLab