add unit/integrted tests for dpc

numpy will happily do zero division (it  returns nan/inf).  The point of 5 is to set those values to 0 so the don't mess up the sums.  I think the two code blocks will yield identical results, the question is which is faster.