(This was added to enable gain flagging in apply mode, IIRC.)
Unfortunately, the phase-only (and slope?) solver wants to have a model in precompute_attributes, so that it can precompute the constant parts of the Hessian. This causes it to crash in apply-only mode. SOP has been to use complex-2x2 in apply-only mode (since the gain tables are compatible), so we never noticed the problem.
The solution is to fix the phase-only solver to skip the Hessian precomputation if a model is not supplied.
In apply-only mode,
precompute_attributesis called here: https://github.com/ratt-ru/CubiCal/blob/master/cubical/solver.py#L615(This was added to enable gain flagging in apply mode, IIRC.)
Unfortunately, the phase-only (and slope?) solver wants to have a model in precompute_attributes, so that it can precompute the constant parts of the Hessian. This causes it to crash in apply-only mode. SOP has been to use complex-2x2 in apply-only mode (since the gain tables are compatible), so we never noticed the problem.
The solution is to fix the phase-only solver to skip the Hessian precomputation if a model is not supplied.