Copyright © 2007-2024 The OpenMx Project
The file omxSadmvnWrapper.cpp has the following hard-coded precision options for numerical integration:
double absEps = 1e-3; double relEps = 0; int MaxPts = 10000*cov->rows
The last of these deviates from classic Mx setting of
1000*ordCov->rows * ordCov->rows * ordCov->rows;
We have been experiencing persistent difficulties with optimization when the number of ordinal variables is large - more frequently than with classic Mx. Possibly the MaxPts setting could be changed to agree with classic Mx?
At the same time, providing mxOption() access to absEps and relEps would help to eliminate possible numerical precision related issues. Could that be added?