Currently, when sampled startpoints are outside the bounds, their value is set to the upper/lower bounds. This may put too much probability mass on the bounds. The problem is illustrated in https://petab--329.org.readthedocs.build/projects/libpetab-python/en/329/example/distributions.html. This should instead be handled by sampling from the respective truncated distributions (which we anyways want to support in PEtab2).