To do list for closure_phase_bias

• remove large matrix in memory

• make nl selection automatic instead of giving a default value

• loosen the requirement for complete sequential network

• add more explanation in the code, mainly when eq.29 and eq.30 are used.