add unit/integrted tests for dpc

I will rewrite 414 as:  _gy = dpc_fit(ffy, ref_fy, fy, start_point, solver)[1] to save some memory I/O