within_support() compares theta against lower/upper with sweep(),
which recycles a bound of the wrong length and warns instead of stopping.
The support test that comes back is then wrong, and it decides which
posterior draws are rejected as leakage and what log_prob() renormalizes
by. A length-1 bound is recycled here, once, so everything downstream sees
one bound per parameter.
