`model.initial_point()` returns inf if using `pm.distributions.transforms.Ordered()`

[pipme]

Thanks! pm.sample(1000, model=pymc_model, initvals = {"mu": np.array([0.0, 1.0])}) also works when sampling, but I think specifying initval in the way you suggested is better!