I’ve written up a numeric function that computes the log-likelihood of a discrete-time hidden Markov model, which involves Gauss quadrature and an optimization step. The function using Julia 1.6.2 gives the same exact results as does another implementation written in MATLAB.
I noticed the function outputted different values once I upgraded from version 1.6.2 to version 1.7.3; the log-likelihood was different by around 100. When I went into the VS Code debugger everything worked correctly, except when I called the function from the debugging terminal it once again gave the wrong results.
Before I spend time deducing a MWE, does anyone have an idea what may be causing the problem? Is this the kind of issue that should be posted on the Github?