Hi I would like to propose a better implementation for 'test_indices':

We can remove the unneeded np.array casting:

Cleaner/New:
test_indices = list(set(range(len(texts))) - set(train_indices))

Old:
test_indices = np.array(list(set(range(len(texts))) - set(train_indices)))

We really should be testing each and every feature of Latexify and currently we are not.

Test generation is made easy by a macro that we supply:

using Latexify
@Latexify.generate_test latexify("x/y")

generates a test and puts it in your clipboard to be pasted:

@test latexify("x/y") == 
raw"$\frac{x}{y}$"

One just have to make sure that the test does ac

Currently, values for the llh, gradient, the computed trajectories of states and observables, or the sensitivities are checked in unit tests.
Unfortunately, some bugs, such as incorrect Jacobians or switched minus signs in the Newton solver, will not necessarily affect those quantities. However, they will substantially impact solver performance, by causing way too many steps to be taken. Hence, w