add unit/integrted tests for dpc

I preserved the original code block as it is almost 3 times faster. c[non_zero_arr] /= div_v[non_zero_arr] is time-consuming.