Failed assertion contiguous_solutions.size() == n_solutions in DDECal with DD intervals
When I run DDECal on a 8h obs, I get this error:
0%....10....20....30....40....50....60....70....80....90....100% Finishing processing ... DP3: /home/anoko/projects/DP3/ddecal/SolutionWriter.cc:107: void dp3::ddecal::SolutionWriter::Write(const std::vector<std::vector<std::vector<std::complex > > >&, const std::vector<std::vector<std::vectordp3::ddecal::Constraint::Result > >&, double, double, dp3::base::CalType, const std::vector<std::__cxx11::basic_string >&, const std::vectordp3::base::Direction&, const std::vector<std::vector<std::__cxx11::basic_string > >&, const std::vector&, const std::vector&, const string&): Assertion `contiguous_solutions.size() == n_solutions' failed. Aborted
This seems to be a result of not specifying all directions in the solutions_per_direction parset key. Dp3 is supposed to add values of 'one' for all unspecified direction.
I do not get this error when I provide a value for all 8 directions in ddecal.solutions_per_direction
. So in the parset below, if I change
ddecal.solutions_per_direction=[8,1,1,1,1]
to
ddecal.solutions_per_direction=[8,1,1,1,1,1,1,1]
all works fine (because my model file has 8 directions). This is the parset:
msin=L432696-NL.ms
msin.ntimes=51
msout=
steps=[ddecal]
ddecal.nchan=0
ddecal.solint=8
ddecal.solutions_per_direction=[8,1,1,1,1]
ddecal.h5parm=solutions.h5
ddecal.solveralgorithm=directioniterative
ddecal.sourcedb=3c196-field-bbs.sourcedb