curious lack of warning/error

The following do check out ok against classic Mx. I suspect that a) OpenMx only pays attention to the triangle below the diagonal of the covariance matrix (second example); and b) that it does something to curtail out-of-range values for the correlation so that optimization doesn't come to a grinding halt. I've logged a bug report: http://openmx.psyc.virginia.edu/issue/2013/02/omxmnor-problem
and also fixed a couple of minor errors in the example script (reading & mxFactoring the data)


> omxMnor(matrix(c(1,.5,.5,1),2,2), c(0, 0), c(-Inf, -Inf), c(0,0))
[,1]
[1,] 0.3333333
> omxMnor(matrix(c(1,.3,.5,1),2,2), c(0, 0), c(-Inf, -Inf), c(1, 1))
[,1]
[1,] 0.7281473
> omxMnor(matrix(c(1,.3,.3,1),2,2), c(0, 0), c(-Inf, -Inf), c(1, 1))
[,1]
[1,] 0.7281473


I don't see anything wrong with the standardization code (in omxAlgebraFunctions.c) called before the call to SADMVN to get the integral:

void omxStandardizeCovMatrix(omxMatrix* cov, double* corList, double* weights) {
// Maybe coerce this into an algebra or sequence of algebras?

if(OMX_DEBUG) { Rprintf("Standardizing matrix."); }

int rows = cov->rows;

for(int i = 0; i < rows; i++) {
weights[i] = sqrt(omxMatrixElement(cov, i, i));
}

for(int i = 0; i < rows; i++) {
for(int j = 0; j < i; j++) {
corList[((i*(i-1))/2) + j] = omxMatrixElement(cov, i, j) / (weights[i] * weights[j]);
}
}
}