add unit/integrted tests for dpc

You can also do this as 

```python
non_zero_arr = (div_v != 0)
c[non_zero_arr] /= div_v[non_zero_arr]
```

You should bench mark those two formulations (I suggest `%timeit` in IPython) to sort out which is faster.