- omp_set_num_threads(nc >= 1 ? nc : omp_get_num_procs());
- #pragma omp parallel for
+ // TODO: stabilize this (for now, random result)
+// omp_set_num_threads(nc >= 1 ? nc : omp_get_num_procs());
+// #pragma omp parallel for
- W[mi(j,k,dim,dim)] = W[mi(k,j,dim,dim)];
+ W[mi(k,j,dim,dim)] = W[mi(j,k,dim,dim)];