Xtensor random::randn returning number number below lower bound and above upper bound

Question

I am trying to generate random numbers of shape {10} and range between -0.5 and 0.5

#include <iostream>
#include <xtensor/xio.hpp>
#include <xtensor/xrandom.hpp>

int main() {
    xt::xarray<double> a = xt::random::randn<double>({10}, -0.5, 0.5);
    std::cout << a;
    return 0;
}

but it's returning an array

{-0.432735, -0.573191, -0.269675, -1.435692, -0.418144, -0.607127,
 -0.350702, -0.913972, -0.494892,  0.027733}

where -0.573191 is smaller and -1.435692 is bigger than range

Ted Lyngmo · Accepted Answer · 2022-06-11T03:58:45.863

1

According to the documentation for randn:

template<class T, class S, class E = random::default_engine_type>
auto xt::random::randn(const S& shape, T mean = 0, T std_dev = 1,
                       E& engine = random::get_default_random_engine())

The function draws numbers from std::normal_distribution and you've set the mean value to be -0.5 and the standard deviation to be 0.5. These are not lower and upper bounds.

You may want to use the rand function instead:

template<class T, class S, class E = random::default_engine_type>
auto xt::random::rand(const S& shape, T lower = 0, T upper = 1,
                      E& engine = random::get_default_random_engine())

This draws numbers from std::uniform_real_distribution - so, it's not the same thing, but it has a lower and an upper bound.

edited Jun 11 '22 at 03:58

answered Jun 11 '22 at 03:45

Ted Lyngmo

93,841
5
60
108

yeah, that was fast. is there a way to set upper and lower bounds in randn or is there any other function – Jaisal Francis Jun 11 '22 at 03:49
@JaisalFrancis I added a function with `uniform_real_distribution` to the answer. – Ted Lyngmo Jun 11 '22 at 03:50
@JaisalFrancis It seems someone removed the comments. The answer did help you to solve the problem, right? It would be good if someone with a similar problem in the future is able to see that the answer helped in some way, if only in comments. – Ted Lyngmo Jun 11 '22 at 08:16
1

of course it helped me. I didn't remove the comments. It's strange. probably mods. – Jaisal Francis Jun 11 '22 at 09:44
@JaisalFrancis Yes, they do that sometimes if the comments doesn't seem to help the answer, but in this case, I thought they did :-) – Ted Lyngmo Jun 11 '22 at 10:05
Yep. Btw you DP look like Walter White :-) – Jaisal Francis Jun 11 '22 at 10:38
@JaisalFrancis Call me Heisenberg :-) – Ted Lyngmo Jun 11 '22 at 10:39

Xtensor random::randn returning number number below lower bound and above upper bound

1 Answers1