Index: library/tests/stresssuite/square_mat_stress.cpp
===================================================================
--- library/tests/stresssuite/square_mat_stress.cpp (revision 766)
+++ library/tests/stresssuite/square_mat_stress.cpp (revision 1064)
@@ -37,11 +37,11 @@
                   const itpp::mat &actual, double tolerance,
                   TestDetails const& details ) {
-	if ( !AreClose ( expected, actual, tolerance ) ) {
-		MemoryOutStream stream;
-		stream << "failed at " << expected.rows()
-		<< " x " << expected.cols();
-
-		results.OnTestFailure ( details, stream.GetText() );
-	}
+    if ( !AreClose ( expected, actual, tolerance ) ) {
+        MemoryOutStream stream;
+        stream << "failed at " << expected.rows()
+               << " x " << expected.cols();
+
+        results.OnTestFailure ( details, stream.GetText() );
+    }
 }
 
@@ -52,201 +52,201 @@
 template<typename TMatrix>
 void test_matrix ( int index, square_mat_point *point ) {
-	Real_Timer tt;
-
-	cout << "agenda[" << index << "]:" << endl;
-	mat A = point->get_matrix();
-	int sz = A.rows();
-	CHECK_EQUAL ( A.cols(), sz );
-
-	tt.tic();
-	TMatrix sq_mat ( A );
-	double elapsed = tt.toc();
-	cout << "ctor(" << sz << " x " << sz << "): " << elapsed << " s" << endl;
-
-	tt.tic();
-	mat res = sq_mat.to_mat();
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		CHECK_CLOSE ( A, res, epsilon );
-	}
-
-	cout << "to_mat: " << elapsed << " s" << endl;
-
-	vec v = point->get_vector();
-	double w = point->get_scalar();
-	TMatrix sq_mat2 = sq_mat;
-
-	tt.tic();
-	sq_mat2.opupdt ( v, w );
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		mat expA = A + w * outer_product ( v, v );
-		CHECK_CLOSE ( expA, sq_mat2.to_mat(), epsilon );
-	}
-
-	cout << "opupdt: " << elapsed << " s" << endl;
-
-	TMatrix invmat ( sz );
-
-	tt.tic();
-	sq_mat.inv ( invmat );
-	elapsed = tt.toc();
-
-	mat invA;
-	if ( !fast ) {
-		invA = inv ( A );
-		CHECK_CLOSE ( invA, invmat.to_mat(), epsilon );
-	}
-
-	cout << "inv: " << elapsed << " s" << endl;
-
-	tt.tic();
-	double ld = sq_mat.logdet();
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		double d = det ( A );
-		CHECK_CLOSE ( log ( d ), ld, epsilon );
-	}
-
-	cout << "logdet: " << elapsed << " s" << endl;
-
-	tt.tic();
-	double q = sq_mat.qform ( ones ( sz ) );
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		CHECK_CLOSE ( sumsum ( A ), q, epsilon );
-	}
-
-	cout << "qform(1): " << elapsed << " s" << endl;
-
-	tt.tic();
-	q = sq_mat.qform ( v );
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		double r = ( A * v ) * v;
-		CHECK_CLOSE ( r, q, epsilon );
-	}
-
-	cout << "qform(v): " << elapsed << " s" << endl;
-
-	tt.tic();
-	q = sq_mat.invqform ( v );
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		double r = ( invA * v ) * v;
-		CHECK_CLOSE ( r, q, epsilon );
-	}
-
-	cout << "invqform: " << elapsed << " s" << endl;
-
-	TMatrix twice = sq_mat;
-
-	tt.tic();
-	twice += sq_mat;
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		res = 2 * A;
-		CHECK_CLOSE ( res, twice.to_mat(), epsilon );
-	}
-
-	cout << "+=: " << elapsed << " s" << endl;
-
-	sq_mat2 = sq_mat;
-
-	tt.tic();
-	sq_mat2.mult_sym ( A );
-	elapsed = tt.toc();
-
-	if ( !fast ) {
-		res = ( A * A ) * A.T();
-		CHECK_CLOSE ( res, sq_mat2.to_mat(), epsilon );
-	}
-
-	cout << "^2: " << elapsed << " s" << endl;
+    Real_Timer tt;
+
+    cout << "agenda[" << index << "]:" << endl;
+    mat A = point->get_matrix();
+    int sz = A.rows();
+    CHECK_EQUAL ( A.cols(), sz );
+
+    tt.tic();
+    TMatrix sq_mat ( A );
+    double elapsed = tt.toc();
+    cout << "ctor(" << sz << " x " << sz << "): " << elapsed << " s" << endl;
+
+    tt.tic();
+    mat res = sq_mat.to_mat();
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        CHECK_CLOSE ( A, res, epsilon );
+    }
+
+    cout << "to_mat: " << elapsed << " s" << endl;
+
+    vec v = point->get_vector();
+    double w = point->get_scalar();
+    TMatrix sq_mat2 = sq_mat;
+
+    tt.tic();
+    sq_mat2.opupdt ( v, w );
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        mat expA = A + w * outer_product ( v, v );
+        CHECK_CLOSE ( expA, sq_mat2.to_mat(), epsilon );
+    }
+
+    cout << "opupdt: " << elapsed << " s" << endl;
+
+    TMatrix invmat ( sz );
+
+    tt.tic();
+    sq_mat.inv ( invmat );
+    elapsed = tt.toc();
+
+    mat invA;
+    if ( !fast ) {
+        invA = inv ( A );
+        CHECK_CLOSE ( invA, invmat.to_mat(), epsilon );
+    }
+
+    cout << "inv: " << elapsed << " s" << endl;
+
+    tt.tic();
+    double ld = sq_mat.logdet();
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        double d = det ( A );
+        CHECK_CLOSE ( log ( d ), ld, epsilon );
+    }
+
+    cout << "logdet: " << elapsed << " s" << endl;
+
+    tt.tic();
+    double q = sq_mat.qform ( ones ( sz ) );
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        CHECK_CLOSE ( sumsum ( A ), q, epsilon );
+    }
+
+    cout << "qform(1): " << elapsed << " s" << endl;
+
+    tt.tic();
+    q = sq_mat.qform ( v );
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        double r = ( A * v ) * v;
+        CHECK_CLOSE ( r, q, epsilon );
+    }
+
+    cout << "qform(v): " << elapsed << " s" << endl;
+
+    tt.tic();
+    q = sq_mat.invqform ( v );
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        double r = ( invA * v ) * v;
+        CHECK_CLOSE ( r, q, epsilon );
+    }
+
+    cout << "invqform: " << elapsed << " s" << endl;
+
+    TMatrix twice = sq_mat;
+
+    tt.tic();
+    twice += sq_mat;
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        res = 2 * A;
+        CHECK_CLOSE ( res, twice.to_mat(), epsilon );
+    }
+
+    cout << "+=: " << elapsed << " s" << endl;
+
+    sq_mat2 = sq_mat;
+
+    tt.tic();
+    sq_mat2.mult_sym ( A );
+    elapsed = tt.toc();
+
+    if ( !fast ) {
+        res = ( A * A ) * A.T();
+        CHECK_CLOSE ( res, sq_mat2.to_mat(), epsilon );
+    }
+
+    cout << "^2: " << elapsed << " s" << endl;
 }
 
 void test_agenda ( FTestMatrix test ) {
-	UIFile fag ( agenda_file_name );
-	Array<shared_ptr<square_mat_point> > mag;
-	UI::get ( mag, fag, "agenda", UI::compulsory );
-	int sz = mag.size();
-	CHECK ( sz > 0 );
-	for ( int i = 0; i < sz; ++i ) {
-		test ( i, mag ( i ).get() );
-	}
+    UIFile fag ( agenda_file_name );
+    Array<shared_ptr<square_mat_point> > mag;
+    UI::get ( mag, fag, "agenda", UI::compulsory );
+    int sz = mag.size();
+    CHECK ( sz > 0 );
+    for ( int i = 0; i < sz; ++i ) {
+        test ( i, mag ( i ).get() );
+    }
 }
 
 SUITE ( ldmat ) {
-	TEST ( agenda ) {
-		test_agenda ( test_matrix<ldmat> );
-	}
+    TEST ( agenda ) {
+        test_agenda ( test_matrix<ldmat> );
+    }
 }
 
 SUITE ( fsqmat ) {
-	TEST ( agenda ) {
-		test_agenda ( test_matrix<fsqmat> );
-	}
+    TEST ( agenda ) {
+        test_agenda ( test_matrix<fsqmat> );
+    }
 }
 
 SUITE ( chmat ) {
-	TEST ( agenda ) {
-		test_agenda ( test_matrix<chmat> );
-	}
+    TEST ( agenda ) {
+        test_agenda ( test_matrix<chmat> );
+    }
 }
 
 int main ( int argc, char const *argv[] ) {
-	bool unknown = false;
-	int update_next = 0; // 1 suite, 2 epsilon, 3 agenda file
-	const char *suite = "ldmat";
-	const char **param = argv + 1;
-	while ( *param && !unknown ) {
-		if ( update_next ) {
-			if ( update_next == 1 ) {
-				suite = *param;
-			} else if ( update_next == 2 ) {
-				double eps = atof ( *param );
-				if ( eps > 0 ) {
-					epsilon = eps;
-				} else {
-					cerr << "invalid epsilon value ignored" << endl;
-				}
-			} else {
-				agenda_file_name = *param;
-			}
-
-			update_next = 0;
-		} else {
-			if ( !strcmp ( *param, "-a" ) ) {
-				update_next = 3;
-			} else if ( !strcmp ( *param, "-c" ) ) {
-				update_next = 1;
-			} else if ( !strcmp ( *param, "-e" ) ) {
-				update_next = 2;
-			} else if ( !strcmp ( *param, "-f" ) ) {
-				fast = true;
-			} else {
-				unknown = true;
-			}
-		}
-
-		++param;
-	}
-
-	if ( unknown || update_next ) {
-		cerr << "usage: " << argv[0] << " [ -f ] [ -e epsilon ] [ -a agenda_input.cfg ] [ -c class ]" << endl;
-	} else {
-		UnitTest::TestReporterStdout reporter;
-		UnitTest::TestRunner runner ( reporter );
-		return runner.RunTestsIf ( UnitTest::Test::GetTestList(),
-		                           suite,
-		                           UnitTest::True(),
-		                           0 );
-	}
-}
+    bool unknown = false;
+    int update_next = 0; // 1 suite, 2 epsilon, 3 agenda file
+    const char *suite = "ldmat";
+    const char **param = argv + 1;
+    while ( *param && !unknown ) {
+        if ( update_next ) {
+            if ( update_next == 1 ) {
+                suite = *param;
+            } else if ( update_next == 2 ) {
+                double eps = atof ( *param );
+                if ( eps > 0 ) {
+                    epsilon = eps;
+                } else {
+                    cerr << "invalid epsilon value ignored" << endl;
+                }
+            } else {
+                agenda_file_name = *param;
+            }
+
+            update_next = 0;
+        } else {
+            if ( !strcmp ( *param, "-a" ) ) {
+                update_next = 3;
+            } else if ( !strcmp ( *param, "-c" ) ) {
+                update_next = 1;
+            } else if ( !strcmp ( *param, "-e" ) ) {
+                update_next = 2;
+            } else if ( !strcmp ( *param, "-f" ) ) {
+                fast = true;
+            } else {
+                unknown = true;
+            }
+        }
+
+        ++param;
+    }
+
+    if ( unknown || update_next ) {
+        cerr << "usage: " << argv[0] << " [ -f ] [ -e epsilon ] [ -a agenda_input.cfg ] [ -c class ]" << endl;
+    } else {
+        UnitTest::TestReporterStdout reporter;
+        UnitTest::TestRunner runner ( reporter );
+        return runner.RunTestsIf ( UnitTest::Test::GetTestList(),
+                                   suite,
+                                   UnitTest::True(),
+                                   0 );
+    }
+}
