Better handling of empty vector constants
Fix estimation of variance for small sample sizes
Add tests for scalar outputs (size=())
Add tolerance for std, different formula for mean tolerance
Merge pull request #1608 from lamblin/fix_rng_empty_shape

Fix rng empty shape