Greetings! I'm trying to make the mentioned combo work together, EmbeddedLapack contains f2c-converted blas and lapack libs. Tests with dgesvd_ (6, 7) keep failing, test 2,3 with dgesvd are ok. VT and S matrixes differs from reference ones in 6,7. Correlation matrix is ok everywhere. So it seems that dgesvd breaks with bigger matrixes OR data is incorrectly handled between correlation and dgesvd. Could you, please, tell, why should we shuffle data in such way (is it transpose or smth else)? Sorry for stupid questions and thanks in advance!
for (i = 0; i < M; ++i) { for (j = 0; j < N; ++j) { A[j * M + i] = Y[i * yn + j]; } }
Greetings! I'm trying to make the mentioned combo work together, EmbeddedLapack contains f2c-converted blas and lapack libs. Tests with dgesvd_ (6, 7) keep failing, test 2,3 with dgesvd are ok. VT and S matrixes differs from reference ones in 6,7. Correlation matrix is ok everywhere. So it seems that dgesvd breaks with bigger matrixes OR data is incorrectly handled between correlation and dgesvd. Could you, please, tell, why should we shuffle data in such way (is it transpose or smth else)? Sorry for stupid questions and thanks in advance!
for (i = 0; i < M; ++i) { for (j = 0; j < N; ++j) { A[j * M + i] = Y[i * yn + j]; } }