For testing purposes let the vacuum step/matrix functions be used.

Constant density is still used if cons:0 is requested.  For an density <0 the vacuum versions are used.