Problem
The function is doing too much at this point. We will only accept sample mass density moving forward.
If users want to use packing fraction for computing sample mass density, then there should be a separate function for that.
Proposed solution