Crash with additional constraints #341
Description
When adding an additional constraint (past https://github.com/brendanjmeade/wna/commit/c4c1131da5d4013cd9d671b088cd8fdb0963727b) in the WNA model, the MCMC solver fails immediately (see error below). At first, I figured this was because I was trying to do something kinematically implausible. To test this, I reverted to the previous running model and then tried adding what I would consider and easy to satisfy constraint that would have negligible kinematic importance. Interestingly, this led to a failed solve as well. A branch for this can be found at: https://github.com/brendanjmeade/wna/tree/extra_con
INFO Stations at full weight (area >= threshold): 1237
- celeri.solve_mcmc:_add_station_velocity_likelihood:446 - 2025-12-13 18:16:33
─────────────────────────────────────────────────────────────────────────────────────────────────────────────
0 1s 0s ERROR An error has been caught in function '<module>', process 'MainProcess' (33224), thread 'MainThread' (8716725312):
- __main__:<module>:10 - 2025-12-13 18:16:48
Traceback (most recent call last):
File "/Users/meade/Desktop/celeri-org/celeri/.pixi/envs/default/lib/python3.13/site-packages/loguru/_logger.py", line 1297, in catch_wrapper
return function(*args, **kwargs)
File "/Users/meade/Desktop/celeri-org/celeri/celeri/scripts/celeri_solve.py", line 40, in main
estimation = celeri.solve_mcmc(model)
File "/Users/meade/Desktop/celeri-org/celeri/celeri/solve_mcmc.py", line 680, in solve_mcmc
trace = nutpie.sample(compiled, **kwargs)
File "/Users/meade/Desktop/celeri-org/celeri/.pixi/envs/default/lib/python3.13/site-packages/nutpie/sample.py", line 858, in sample
result = sampler.wait()
File "/Users/meade/Desktop/celeri-org/celeri/.pixi/envs/default/lib/python3.13/site-packages/nutpie/sample.py", line 519, in wait
self._sampler.wait(timeout)
~~~~~~~~~~~~~~~~~~^^^^^^^^^
RuntimeError: All initialization points failed
Caused by:
Logp function returned error: ErrorCode(3)
Traceback (most recent call last):
File "/Users/meade/Desktop/celeri-org/celeri/.pixi/envs/default/bin/celeri-solve", line 10, in <module>
sys.exit(main())
~~~~^^
File "/Users/meade/Desktop/celeri-org/celeri/.pixi/envs/default/lib/python3.13/site-packages/loguru/_logger.py", line 1297, in catch_wrapper
return function(*args, **kwargs)
File "/Users/meade/Desktop/celeri-org/celeri/celeri/scripts/celeri_solve.py", line 40, in main
estimation = celeri.solve_mcmc(model)
File "/Users/meade/Desktop/celeri-org/celeri/celeri/solve_mcmc.py", line 680, in solve_mcmc
trace = nutpie.sample(compiled, **kwargs)
File "/Users/meade/Desktop/celeri-org/celeri/.pixi/envs/default/lib/python3.13/site-packages/nutpie/sample.py", line 858, in sample
result = sampler.wait()
File "/Users/meade/Desktop/celeri-org/celeri/.pixi/envs/default/lib/python3.13/site-packages/nutpie/sample.py", line 519, in wait
self._sampler.wait(timeout)
~~~~~~~~~~~~~~~~~~^^^^^^^^^
RuntimeError: All initialization points failed
Caused by:
Logp function returned error: ErrorCode(3)