Test case that triggers the bug.

main.cc (text/plain), 914 bytes, created by Gert Wollny on 2011-06-07 16:07:42 UTC

#include <itpp/signal/fastica.h> 
const int rows = 5; const int cols = 10; const int nica = 3; double test_signal[rows][cols] = { {-2, 
2, 2, 1.6, -3.2, 1.6, -3.57628e-08, -3.57628e-08, 8.34465e-08, -2}, {-1, 1, 1, -0.4, 0.8, -0.4, 2.14577e-07, -2.38419e-08, 
-2.38419e-08, -1}, {2.38419e-08, 2.38419e-08, 2.38419e-08, -2.4, 4.8, -2.4, 2.38419e-08, 2.38419e-08, 2.38419e-08, 2.38419e-08}, {1, 
-1, -1, -0.4, 0.8, -0.4, 2.38419e-08, 2.38419e-08, 2.38419e-08, 1}, {2, -2, -2, 1.6, -3.2, 1.6, 2.38419e-08, 2.38419e-08, 2.38419e-08, 
2} }; int main(int argc, const char **args) { itpp::mat Signal(rows,cols); for (int r = 0; r < rows; ++r) { for (int c = 0; c < cols; 
++c) { Signal(r,c) = test_signal[c][r]; } } itpp::Fast_ICA fastICA(Signal); fastICA.set_nrof_independent_components(nica); 
fastICA.set_non_linearity( FICA_NONLIN_TANH ); fastICA.set_approach( FICA_APPROACH_DEFL ); fastICA.separate(); return 0; }