You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
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]; } }
The text was updated successfully, but these errors were encountered:
Thank you very much for your comments @vaalberith Both lmtest2.c and lmtest6.c use dgesvd_ to compute job parameters and job result. And yes: for (i = 0; i < M; ++i) { for (j = 0; j < N; ++j) { A[j * M + i] = Y[i * yn + j]; } }
is to transpose Y[M, N] into A[N, M] cause this is what dgesvd_ is expecting.
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]; } }
The text was updated successfully, but these errors were encountered: