Renamed the functions and tried to document. Basically density() decides whether to call get_compressible_density() or get_corrected_density().