add basic sampling test
(#601)
very mild tweak to rho array definition
(#601)