Multiple observed value for a single RV

[xidulu]

The code below runs without exceptions, is it desired?

w1 = normal_scale(1, 1)
x = bernoulli(w1)
y = 1
calc = Calculate("numpyro",niter=10000)
ys = calc.sample(w1, [x, x, x], [1, 0, 1])

[justindomke]

Hmmm, looks like the earlier pieces of evidence are simply being ignored now, so the last call is equivalent to:

ys = calc.sample(w1, x, 1)

That's surely bad—we should throw an exception rather than accepting this.