Commit 5dac1d6a authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Reconstruct angle vectors when building InclusionOutputInfo from HDF5

parent 4908ff55

src/libnptm/outputs.cpp

+4 −0
Original line number Diff line number Diff line
@@ -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;