The distribution bounds does not seem to actually scale the distribution from its default [0,1]. Example code below produces pce.samples that are [0,1].
dimension = 3
alpha = 2.#a = b = 2 makes a normal dist
beta = 2.
bounds = np.zeros([2,3])
bounds[0,0],bounds[1,0] = -3, 3
bounds[0,1],bounds[1,1] = -3, 3
bounds[0,2],bounds[1,2] = -3, 3
dist = BetaDistribution(alpha=alpha, beta=beta, dim=dimension,domain = bounds)
The distribution bounds does not seem to actually scale the distribution from its default [0,1]. Example code below produces pce.samples that are [0,1].