Index: src/unittests/Makefile.am
===================================================================
--- src/unittests/Makefile.am	(revision cca9efb7bf18583f35d0380acb97e9034c42bbe7)
+++ src/unittests/Makefile.am	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
@@ -22,6 +22,4 @@
   CountBondsUnitTest \
   FormulaUnittest \
-  GSLMatrixSymmetricUnitTest \
-  GSLMatrixUnitTest \
   GSLVectorUnitTest \
   InfoUnitTest \
@@ -32,4 +30,6 @@
   LogUnitTest \
   manipulateAtomsTest \
+  MatrixContentSymmetricUnitTest \
+  MatrixContentUnitTest \
   MatrixUnittest \
   MenuDescriptionUnitTest \
@@ -82,6 +82,4 @@
   CountBondsUnitTest.cpp \
   FormulaUnittest.cpp \
-  gslmatrixsymmetricunittest.cpp \
-  gslmatrixunittest.cpp \
   gslvectorunittest.cpp \
   infounittest.cpp \
@@ -91,4 +89,6 @@
   listofbondsunittest.cpp \
   logunittest.cpp \
+  MatrixContentSymmetricUnittest.cpp \
+  MatrixContentUnittest.cpp \
   MatrixUnittest.cpp \
   manipulateAtomsTest.cpp \
@@ -122,6 +122,4 @@
   CountBondsUnitTest.hpp \
   FormulaUnittest.hpp \
-  gslmatrixsymmetricunittest.hpp \
-  gslmatrixunittest.hpp \
   gslvectorunittest.hpp \
   infounittest.hpp \
@@ -132,4 +130,6 @@
   logunittest.hpp \
   manipulateAtomsTest.hpp \
+  MatrixContentSymmetricUnittest.hpp \
+  MatrixContentUnittest.hpp \
   MatrixUnittest.hpp \
   MenuDescriptionUnitTest.hpp \
@@ -187,9 +187,9 @@
 FormulaUnittest_LDADD = ${ALLLIBS}
 
-GSLMatrixSymmetricUnitTest_SOURCES = UnitTestMain.cpp gslmatrixsymmetricunittest.cpp gslmatrixsymmetricunittest.hpp
-GSLMatrixSymmetricUnitTest_LDADD = ${ALLLIBS}
-
-GSLMatrixUnitTest_SOURCES = UnitTestMain.cpp gslmatrixunittest.cpp gslmatrixunittest.hpp
-GSLMatrixUnitTest_LDADD = ${ALLLIBS}
+MatrixContentSymmetricUnitTest_SOURCES = UnitTestMain.cpp MatrixContentSymmetricUnittest.cpp MatrixContentSymmetricUnittest
+MatrixContentSymmetricUnitTest_LDADD = ${ALLLIBS}
+
+MatrixContentUnitTest_SOURCES = UnitTestMain.cpp MatrixContentUnittest.cpp MatrixContentUnittest.hpp
+MatrixContentUnitTest_LDADD = ${ALLLIBS}
 
 GSLVectorUnitTest_SOURCES = UnitTestMain.cpp gslvectorunittest.cpp gslvectorunittest.hpp
Index: src/unittests/MatrixContentSymmetricUnittest.cpp
===================================================================
--- src/unittests/MatrixContentSymmetricUnittest.cpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
+++ src/unittests/MatrixContentSymmetricUnittest.cpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
@@ -0,0 +1,284 @@
+/*
+ * Project: MoleCuilder
+ * Description: creates and alters molecular systems
+ * Copyright (C)  2010 University of Bonn. All rights reserved.
+ * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
+ */
+
+/*
+ * MatrixContentSymmetricUnittest.cpp
+ *
+ *  Created on: Jan 8, 2010
+ *      Author: heber
+ */
+
+// include config.h
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
+using namespace std;
+
+#include <cppunit/CompilerOutputter.h>
+#include <cppunit/extensions/TestFactoryRegistry.h>
+#include <cppunit/ui/text/TestRunner.h>
+
+#include "MatrixContentSymmetricUnittest.hpp"
+
+#include "LinearAlgebra/MatrixContent.hpp"
+
+#ifdef HAVE_TESTRUNNER
+#include "UnitTestMain.hpp"
+#endif /*HAVE_TESTRUNNER*/
+
+/********************************************** Test classes **************************************/
+
+// Registers the fixture into the 'registry'
+CPPUNIT_TEST_SUITE_REGISTRATION( MatrixContentSymmetricTest );
+
+
+void MatrixContentSymmetricTest::setUp()
+{
+  m = new MatrixContent(3,3);
+};
+
+void MatrixContentSymmetricTest::tearDown()
+{
+  delete(m);
+};
+
+/** Unit Test for accessing matrix elements.
+ *
+ */
+void MatrixContentSymmetricTest::AccessTest()
+{
+  // check whether all elements are initially zero
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( 0., m->at(i,j) );
+
+  // set
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, i*3+j );
+
+  // and check
+  double *ptr = NULL;
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++) {
+      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), m->at(i,j) );
+      ptr = m->Pointer(i,j);
+      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), *ptr );
+    }
+
+  // assignment
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, i*3+j );
+  MatrixContent *dest = new MatrixContent(3,3);
+  *dest = *m;
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( dest->at(i,j), m->at(i,j) );
+  delete(dest);
+
+  // out of bounds
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(4,2) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(2,17) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(1024,140040) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(-1,0) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(0,-1) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(-1,-1) );
+};
+
+/** Unit Test for initializating matrices.
+ *
+ */
+void MatrixContentSymmetricTest::InitializationTest()
+{
+  // set zero
+  m->setZero();
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( 0., m->at(i,j) );
+
+  // set all
+  m->setValue(1.5);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( 1.5, m->at(i,j) );
+
+  // set basis
+  m->setIdentity();
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->at(i,j) );
+
+  // set from array
+  double array[] = { 1., 0., 0.,
+                     0., 1., 0.,
+                     0., 0., 1. };
+  m->setFromDoubleArray(array);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->at(i,j) );
+
+};
+
+/** Unit Test for copying matrices.
+ *
+ */
+void MatrixContentSymmetricTest::CopyTest()
+{
+  // set basis
+  MatrixContent *dest = NULL;
+  for (int i=0;i<3;i++) {
+    m->setValue(i);
+    dest = new MatrixContent(m);
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( m->at(i,j) , dest->at(i,j) );
+
+    delete(dest);
+  }
+};
+
+/** Unit Test for exchanging rows and columns.
+ *
+ */
+void MatrixContentSymmetricTest::ExchangeTest()
+{
+  // set to 1,1,1,2, ...
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, i+1 );
+
+  // swap such that nothing happens
+  m->SwapColumns(1,2);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->at(i,j) );
+
+  // swap two rows
+  m->SwapRows(1,2);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      switch (j) {
+        case 0:
+          CPPUNIT_ASSERT_EQUAL( 1., m->at(j,i) );
+          break;
+        case 1:
+          CPPUNIT_ASSERT_EQUAL( 3., m->at(j,i) );
+          break;
+        case 2:
+          CPPUNIT_ASSERT_EQUAL( 2., m->at(j,i) );
+          break;
+        default:
+          CPPUNIT_ASSERT_EQUAL( -1., m->at(i,j) );
+      }
+  // check that op is reversable
+  m->SwapRows(1,2);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->at(i,j) );
+
+  // set to 1,2,3,1, ...
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, j+1. );
+
+  // swap such that nothing happens
+  m->SwapRows(0,2);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->at(i,j) );
+
+  // swap two columns
+  m->SwapColumns(0,2);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      switch (j) {
+        case 0:
+          CPPUNIT_ASSERT_EQUAL( 3., m->at(i,j) );
+          break;
+        case 1:
+          CPPUNIT_ASSERT_EQUAL( 2., m->at(i,j) );
+          break;
+        case 2:
+          CPPUNIT_ASSERT_EQUAL( 1., m->at(i,j) );
+          break;
+        default:
+          CPPUNIT_ASSERT_EQUAL( -1., m->at(i,j) );
+      }
+  // check that op is reversable
+  m->SwapColumns(0,2);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->at(i,j) );
+
+
+  // set to 1,2,3, ...
+  MatrixContent *n =  new MatrixContent(3,3);
+  for (int i=0;i<3;i++)
+    for (int j=0;j<3;j++) {
+      m->set(i,j, 3*i+j+1 );
+      n->set(i,j, 3*j+i+1 );
+    }
+  // transpose
+  MatrixContent res = ((const MatrixContent)*m).transpose();
+  CPPUNIT_ASSERT( *n == res );
+  // second transpose
+  res.transpose();
+  CPPUNIT_ASSERT( *m == res );
+};
+
+/** Unit Test for matrix properties.
+ *
+ */
+void MatrixContentSymmetricTest::PropertiesTest()
+{
+  // is zero
+  m->setZero();
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNull() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
+
+  // is positive
+  m->setValue(0.5);
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsPositive() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
+
+  // is negative
+  m->setValue(-0.1);
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNegative() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNonNegative() );
+
+  // is positive definite
+  double array[] = { 1., 0., 0.,
+                     0., 1., 1.,
+                     0., 0., 1. };
+  m->setFromDoubleArray(array);
+  CPPUNIT_ASSERT_EQUAL( true, m->IsPositiveDefinite() );
+
+  //determinant
+  m->setIdentity();
+  CPPUNIT_ASSERT_EQUAL( 1., m->Determinant() );
+
+  m->setZero();
+  CPPUNIT_ASSERT_EQUAL( 0., m->Determinant() );
+
+  m->set( 0, 0, 1.);
+  m->set( 1, 1, 1.);
+  m->set( 2, 1, 1.);
+  CPPUNIT_ASSERT_EQUAL( 0., m->Determinant() );
+  
+  double array2[] = { 2., 0., 1.,
+                     -3., 1., 1.,
+                     1., 5.5, 1. };
+  m->setFromDoubleArray(array2);
+  CPPUNIT_ASSERT_EQUAL( -26.5, m->Determinant() );
+};
Index: src/unittests/MatrixContentSymmetricUnittest.hpp
===================================================================
--- src/unittests/MatrixContentSymmetricUnittest.hpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
+++ src/unittests/MatrixContentSymmetricUnittest.hpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
@@ -0,0 +1,42 @@
+/*
+ * MatrixContentSymmetricUnittest.hpp
+ *
+ *  Created on: Jan 8, 2010
+ *      Author: heber
+ */
+
+#ifndef MATRIXCONTENTSYMMETRICUNITTEST_HPP_
+#define MATRIXCONTENTSYMMETRICUNITTEST_HPP_
+
+#include <cppunit/extensions/HelperMacros.h>
+
+class MatrixContent;
+
+/********************************************** Test classes **************************************/
+
+class MatrixContentSymmetricTest : public CppUnit::TestFixture
+{
+    CPPUNIT_TEST_SUITE( MatrixContentSymmetricTest) ;
+    CPPUNIT_TEST (AccessTest );
+    CPPUNIT_TEST (InitializationTest );
+    CPPUNIT_TEST (CopyTest );
+    CPPUNIT_TEST (ExchangeTest );
+    CPPUNIT_TEST (PropertiesTest );
+    CPPUNIT_TEST_SUITE_END();
+
+public:
+    void setUp();
+    void tearDown();
+
+    void AccessTest();
+    void InitializationTest();
+    void CopyTest();
+    void ExchangeTest();
+    void PropertiesTest();
+
+private:
+    MatrixContent *m;
+};
+
+
+#endif /* MATRIXCONTENTSYMMETRICUNITTEST_HPP_ */
Index: src/unittests/MatrixContentUnittest.cpp
===================================================================
--- src/unittests/MatrixContentUnittest.cpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
+++ src/unittests/MatrixContentUnittest.cpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
@@ -0,0 +1,266 @@
+/*
+ * Project: MoleCuilder
+ * Description: creates and alters molecular systems
+ * Copyright (C)  2010 University of Bonn. All rights reserved.
+ * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
+ */
+
+/*
+ * MatrixContentUnittest.cpp
+ *
+ *  Created on: Jan 8, 2010
+ *      Author: heber
+ */
+
+// include config.h
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
+using namespace std;
+
+#include <cppunit/CompilerOutputter.h>
+#include <cppunit/extensions/TestFactoryRegistry.h>
+#include <cppunit/ui/text/TestRunner.h>
+
+#include "MatrixContentUnittest.hpp"
+
+#include "LinearAlgebra/MatrixContent.hpp"
+
+#ifdef HAVE_TESTRUNNER
+#include "UnitTestMain.hpp"
+#endif /*HAVE_TESTRUNNER*/
+
+/********************************************** Test classes **************************************/
+
+// Registers the fixture into the 'registry'
+CPPUNIT_TEST_SUITE_REGISTRATION( MatrixContentTest );
+
+
+void MatrixContentTest::setUp()
+{
+  m = new MatrixContent(4,3);
+};
+
+void MatrixContentTest::tearDown()
+{
+  delete(m);
+};
+
+/** Unit Test for accessing matrix elements.
+ *
+ */
+void MatrixContentTest::AccessTest()
+{
+  // check whether all elements are initially zero
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( 0., m->at(i,j) );
+
+  // set
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, i*3+j );
+
+  // and check
+  double *ptr = NULL;
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++) {
+      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), m->at(i,j) );
+      ptr = m->Pointer(i,j);
+      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), *ptr );
+    }
+
+  // assignment
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, i*3+j );
+  MatrixContent *dest = new MatrixContent(4,3);
+  *dest = *m;
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( dest->at(i,j), m->at(i,j) );
+  delete(dest);
+
+  // out of bounds
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(5,2) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(2,17) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(1024,140040) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(-1,0) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(0,-1) );
+  //CPPUNIT_ASSERT_EQUAL(0., v->at(-1,-1) );
+};
+
+/** Unit Test for initializating matrices.
+ *
+ */
+void MatrixContentTest::InitializationTest()
+{
+  // set zero
+  m->setZero();
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( 0., m->at(i,j) );
+
+  // set all
+  m->setValue(1.5);
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( 1.5, m->at(i,j) );
+
+  // set basis
+  m->setIdentity();
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->at(i,j) );
+
+  // set from array
+  double array[] = { 1., 0., 0.,
+                     0., 1., 0.,
+                     0., 0., 1.,
+                     0., 0., 0. };
+  m->setFromDoubleArray(array);
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->at(i,j) );
+
+};
+
+/** Unit Test for copying matrices.
+ *
+ */
+void MatrixContentTest::CopyTest()
+{
+  // set basis
+  MatrixContent *dest = NULL;
+  for (int i=0;i<4;i++) {
+    m->setValue(i);
+    dest = new MatrixContent(m);
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( m->at(i,j) , dest->at(i,j) );
+
+    delete(dest);
+  }
+};
+
+/** Unit Test for exchanging rows and columns.
+ *
+ */
+void MatrixContentTest::ExchangeTest()
+{
+  // set to 1,1,1,2, ...
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, i+1 );
+
+  // swap such that nothing happens
+  CPPUNIT_ASSERT_EQUAL( true, m->SwapColumns(1,2) );
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->at(i,j) );
+
+  // swap two rows
+  CPPUNIT_ASSERT_EQUAL( true, m->SwapRows(1,2) );
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      switch (i) {
+        case 0:
+          CPPUNIT_ASSERT_EQUAL( 1., m->at(i,j) );
+          break;
+        case 1:
+          CPPUNIT_ASSERT_EQUAL( 3., m->at(i,j) );
+          break;
+        case 2:
+          CPPUNIT_ASSERT_EQUAL( 2., m->at(i,j) );
+          break;
+        case 3:
+          CPPUNIT_ASSERT_EQUAL( 4., m->at(i,j) );
+          break;
+        default:
+          CPPUNIT_ASSERT_EQUAL( -1., m->at(i,j) );
+      }
+  // check that op is reversable
+  CPPUNIT_ASSERT_EQUAL( true, m->SwapRows(1,2) );
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->at(i,j) );
+
+  // set to 1,2,3,1, ...
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      m->set(i,j, j+1. );
+
+  // swap such that nothing happens
+  CPPUNIT_ASSERT_EQUAL( true, m->SwapRows(0,2) );
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->at(i,j) );
+
+  // swap two columns
+  CPPUNIT_ASSERT_EQUAL( true, m->SwapColumns(0,2) );
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      switch (j) {
+        case 0:
+          CPPUNIT_ASSERT_EQUAL( 3., m->at(i,j) );
+          break;
+        case 1:
+          CPPUNIT_ASSERT_EQUAL( 2., m->at(i,j) );
+          break;
+        case 2:
+          CPPUNIT_ASSERT_EQUAL( 1., m->at(i,j) );
+          break;
+        default:
+          CPPUNIT_ASSERT_EQUAL( -1., m->at(i,j) );
+      }
+  // check that op is reversable
+  CPPUNIT_ASSERT_EQUAL( true, m->SwapColumns(0,2) );
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++)
+      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->at(i,j) );
+
+
+  // set to 1,2,3,4, ...
+  MatrixContent *n =  new MatrixContent(3,4);
+  for (int i=0;i<4;i++)
+    for (int j=0;j<3;j++) {
+      m->set(i,j, 3*i+j+1 );
+      n->set(j,i, 3*i+j+1 );
+    }
+  // transpose
+  MatrixContent res = ((const MatrixContent)(*m)).transpose();
+  CPPUNIT_ASSERT( *n == res );
+  // second transpose
+  MatrixContent res2 = ((const MatrixContent)res).transpose();
+  CPPUNIT_ASSERT( *m == res2 );
+};
+
+/** Unit Test for matrix properties.
+ *
+ */
+void MatrixContentTest::PropertiesTest()
+{
+  // is zero
+  m->setZero();
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNull() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
+
+  // is positive
+  m->setValue(0.5);
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsPositive() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
+
+  // is negative
+  m->setValue(-0.1);
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
+  CPPUNIT_ASSERT_EQUAL( true, m->IsNegative() );
+  CPPUNIT_ASSERT_EQUAL( false, m->IsNonNegative() );
+
+  // is positive definite
+  CPPUNIT_ASSERT_EQUAL( false, m->IsPositiveDefinite() );
+};
Index: src/unittests/MatrixContentUnittest.hpp
===================================================================
--- src/unittests/MatrixContentUnittest.hpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
+++ src/unittests/MatrixContentUnittest.hpp	(revision 0d4424de796ea0eb6de8098859070b5acdb6edae)
@@ -0,0 +1,42 @@
+/*
+ * MatrixContentUnittest.hpp
+ *
+ *  Created on: Jan 8, 2010
+ *      Author: heber
+ */
+
+#ifndef MATRIXCONTENTUNITTEST_HPP_
+#define MATRIXCONTENTUNITTEST_HPP_
+
+#include <cppunit/extensions/HelperMacros.h>
+
+class MatrixContent;
+
+/********************************************** Test classes **************************************/
+
+class MatrixContentTest : public CppUnit::TestFixture
+{
+    CPPUNIT_TEST_SUITE( MatrixContentTest) ;
+    CPPUNIT_TEST (AccessTest );
+    CPPUNIT_TEST (InitializationTest );
+    CPPUNIT_TEST (CopyTest );
+    CPPUNIT_TEST (ExchangeTest );
+    CPPUNIT_TEST (PropertiesTest );
+    CPPUNIT_TEST_SUITE_END();
+
+public:
+    void setUp();
+    void tearDown();
+
+    void AccessTest();
+    void InitializationTest();
+    void CopyTest();
+    void ExchangeTest();
+    void PropertiesTest();
+
+private:
+    MatrixContent *m;
+};
+
+
+#endif /* MATRIXCONTENTUNITTEST_HPP_ */
Index: src/unittests/gslmatrixsymmetricunittest.cpp
===================================================================
--- src/unittests/gslmatrixsymmetricunittest.cpp	(revision cca9efb7bf18583f35d0380acb97e9034c42bbe7)
+++ 	(revision )
@@ -1,283 +1,0 @@
-/*
- * Project: MoleCuilder
- * Description: creates and alters molecular systems
- * Copyright (C)  2010 University of Bonn. All rights reserved.
- * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
- */
-
-/*
- * gslmatrixunittest.cpp
- *
- *  Created on: Jan 8, 2010
- *      Author: heber
- */
-
-// include config.h
-#ifdef HAVE_CONFIG_H
-#include <config.h>
-#endif
-
-using namespace std;
-
-#include <cppunit/CompilerOutputter.h>
-#include <cppunit/extensions/TestFactoryRegistry.h>
-#include <cppunit/ui/text/TestRunner.h>
-
-#include "gslmatrixsymmetricunittest.hpp"
-
-#ifdef HAVE_TESTRUNNER
-#include "UnitTestMain.hpp"
-#endif /*HAVE_TESTRUNNER*/
-
-/********************************************** Test classes **************************************/
-
-// Registers the fixture into the 'registry'
-CPPUNIT_TEST_SUITE_REGISTRATION( GSLMatrixSymmetricTest );
-
-
-void GSLMatrixSymmetricTest::setUp()
-{
-  m = new GSLMatrix(3,3);
-};
-
-void GSLMatrixSymmetricTest::tearDown()
-{
-  delete(m);
-};
-
-/** Unit Test for accessing matrix elements.
- *
- */
-void GSLMatrixSymmetricTest::AccessTest()
-{
-  // check whether all elements are initially zero
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( 0., m->Get(i,j) );
-
-  // set
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, i*3+j );
-
-  // and check
-  double *ptr = NULL;
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++) {
-      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), m->Get(i,j) );
-      ptr = m->Pointer(i,j);
-      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), *ptr );
-    }
-
-  // assignment
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, i*3+j );
-  GSLMatrix *dest = new GSLMatrix(3,3);
-  *dest = *m;
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( dest->Get(i,j), m->Get(i,j) );
-  delete(dest);
-
-  // out of bounds
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(4,2) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(2,17) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(1024,140040) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(-1,0) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(0,-1) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(-1,-1) );
-};
-
-/** Unit Test for initializating matrices.
- *
- */
-void GSLMatrixSymmetricTest::InitializationTest()
-{
-  // set zero
-  m->SetZero();
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( 0., m->Get(i,j) );
-
-  // set all
-  m->SetAll(1.5);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( 1.5, m->Get(i,j) );
-
-  // set basis
-  m->SetIdentity();
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->Get(i,j) );
-
-  // set from array
-  double array[] = { 1., 0., 0.,
-                     0., 1., 0.,
-                     0., 0., 1. };
-  m->SetFromDoubleArray(array);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->Get(i,j) );
-
-};
-
-/** Unit Test for copying matrices.
- *
- */
-void GSLMatrixSymmetricTest::CopyTest()
-{
-  // set basis
-  GSLMatrix *dest = NULL;
-  for (int i=0;i<3;i++) {
-    m->SetAll(i);
-    dest = new GSLMatrix(m);
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( m->Get(i,j) , dest->Get(i,j) );
-
-    delete(dest);
-  }
-};
-
-/** Unit Test for exchanging rows and columns.
- *
- */
-void GSLMatrixSymmetricTest::ExchangeTest()
-{
-  // set to 1,1,1,2, ...
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, i+1 );
-
-  // swap such that nothing happens
-  m->SwapColumns(1,2);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->Get(i,j) );
-
-  // swap two rows
-  m->SwapRows(1,2);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      switch (j) {
-        case 0:
-          CPPUNIT_ASSERT_EQUAL( 1., m->Get(j,i) );
-          break;
-        case 1:
-          CPPUNIT_ASSERT_EQUAL( 3., m->Get(j,i) );
-          break;
-        case 2:
-          CPPUNIT_ASSERT_EQUAL( 2., m->Get(j,i) );
-          break;
-        default:
-          CPPUNIT_ASSERT_EQUAL( -1., m->Get(i,j) );
-      }
-  // check that op is reversable
-  m->SwapRows(1,2);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->Get(i,j) );
-
-  // set to 1,2,3,1, ...
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, j+1. );
-
-  // swap such that nothing happens
-  m->SwapRows(0,2);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->Get(i,j) );
-
-  // swap two columns
-  m->SwapColumns(0,2);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      switch (j) {
-        case 0:
-          CPPUNIT_ASSERT_EQUAL( 3., m->Get(i,j) );
-          break;
-        case 1:
-          CPPUNIT_ASSERT_EQUAL( 2., m->Get(i,j) );
-          break;
-        case 2:
-          CPPUNIT_ASSERT_EQUAL( 1., m->Get(i,j) );
-          break;
-        default:
-          CPPUNIT_ASSERT_EQUAL( -1., m->Get(i,j) );
-      }
-  // check that op is reversable
-  m->SwapColumns(0,2);
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->Get(i,j) );
-
-
-  // set to 1,2,3,4, ...
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, 3*i+j+1 );
-  // transpose
-  m->Transpose();
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)(3*i+j+1), m->Get(j,i) );
-  // second transpose
-  m->Transpose();
-  for (int i=0;i<3;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)(3*i+j+1), m->Get(i,j) );
-};
-
-/** Unit Test for matrix properties.
- *
- */
-void GSLMatrixSymmetricTest::PropertiesTest()
-{
-  // is zero
-  m->SetZero();
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNull() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
-
-  // is positive
-  m->SetAll(0.5);
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsPositive() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
-
-  // is negative
-  m->SetAll(-0.1);
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNegative() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNonNegative() );
-
-  // is positive definite
-  double array[] = { 1., 0., 0.,
-                     0., 1., 1.,
-                     0., 0., 1. };
-  m->SetFromDoubleArray(array);
-  CPPUNIT_ASSERT_EQUAL( true, m->IsPositiveDefinite() );
-
-  //determinant
-  m->SetIdentity();
-  CPPUNIT_ASSERT_EQUAL( 1., m->Determinant() );
-
-  m->SetZero();
-  CPPUNIT_ASSERT_EQUAL( 0., m->Determinant() );
-
-  m->Set( 0, 0, 1.);
-  m->Set( 1, 1, 1.);
-  m->Set( 2, 1, 1.);
-  CPPUNIT_ASSERT_EQUAL( 0., m->Determinant() );
-  
-  double array2[] = { 2., 0., 1.,
-                     -3., 1., 1.,
-                     1., 5.5, 1. };
-  m->SetFromDoubleArray(array2);
-  CPPUNIT_ASSERT_EQUAL( -26.5, m->Determinant() );
-};
Index: src/unittests/gslmatrixsymmetricunittest.hpp
===================================================================
--- src/unittests/gslmatrixsymmetricunittest.hpp	(revision cca9efb7bf18583f35d0380acb97e9034c42bbe7)
+++ 	(revision )
@@ -1,42 +1,0 @@
-/*
- * gslmatrixunittest.hpp
- *
- *  Created on: Jan 8, 2010
- *      Author: heber
- */
-
-#ifndef GSLMATRIXUNITTEST_HPP_
-#define GSLMATRIXUNITTEST_HPP_
-
-#include <cppunit/extensions/HelperMacros.h>
-
-#include "LinearAlgebra/gslmatrix.hpp"
-
-/********************************************** Test classes **************************************/
-
-class GSLMatrixSymmetricTest : public CppUnit::TestFixture
-{
-    CPPUNIT_TEST_SUITE( GSLMatrixSymmetricTest) ;
-    CPPUNIT_TEST (AccessTest );
-    CPPUNIT_TEST (InitializationTest );
-    CPPUNIT_TEST (CopyTest );
-    CPPUNIT_TEST (ExchangeTest );
-    CPPUNIT_TEST (PropertiesTest );
-    CPPUNIT_TEST_SUITE_END();
-
-public:
-    void setUp();
-    void tearDown();
-
-    void AccessTest();
-    void InitializationTest();
-    void CopyTest();
-    void ExchangeTest();
-    void PropertiesTest();
-
-private:
-    GSLMatrix *m;
-};
-
-
-#endif /* GSLMATRIXUNITTEST_HPP_ */
Index: src/unittests/gslmatrixunittest.cpp
===================================================================
--- src/unittests/gslmatrixunittest.cpp	(revision cca9efb7bf18583f35d0380acb97e9034c42bbe7)
+++ 	(revision )
@@ -1,265 +1,0 @@
-/*
- * Project: MoleCuilder
- * Description: creates and alters molecular systems
- * Copyright (C)  2010 University of Bonn. All rights reserved.
- * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
- */
-
-/*
- * gslmatrixunittest.cpp
- *
- *  Created on: Jan 8, 2010
- *      Author: heber
- */
-
-// include config.h
-#ifdef HAVE_CONFIG_H
-#include <config.h>
-#endif
-
-using namespace std;
-
-#include <cppunit/CompilerOutputter.h>
-#include <cppunit/extensions/TestFactoryRegistry.h>
-#include <cppunit/ui/text/TestRunner.h>
-
-#include "gslmatrixunittest.hpp"
-
-#ifdef HAVE_TESTRUNNER
-#include "UnitTestMain.hpp"
-#endif /*HAVE_TESTRUNNER*/
-
-/********************************************** Test classes **************************************/
-
-// Registers the fixture into the 'registry'
-CPPUNIT_TEST_SUITE_REGISTRATION( GSLMatrixTest );
-
-
-void GSLMatrixTest::setUp()
-{
-  m = new GSLMatrix(4,3);
-};
-
-void GSLMatrixTest::tearDown()
-{
-  delete(m);
-};
-
-/** Unit Test for accessing matrix elements.
- *
- */
-void GSLMatrixTest::AccessTest()
-{
-  // check whether all elements are initially zero
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( 0., m->Get(i,j) );
-
-  // set
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, i*3+j );
-
-  // and check
-  double *ptr = NULL;
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++) {
-      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), m->Get(i,j) );
-      ptr = m->Pointer(i,j);
-      CPPUNIT_ASSERT_EQUAL( (double)(i*3+j), *ptr );
-    }
-
-  // assignment
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, i*3+j );
-  GSLMatrix *dest = new GSLMatrix(4,3);
-  *dest = *m;
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( dest->Get(i,j), m->Get(i,j) );
-  delete(dest);
-
-  // out of bounds
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(5,2) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(2,17) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(1024,140040) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(-1,0) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(0,-1) );
-  //CPPUNIT_ASSERT_EQUAL(0., v->Get(-1,-1) );
-};
-
-/** Unit Test for initializating matrices.
- *
- */
-void GSLMatrixTest::InitializationTest()
-{
-  // set zero
-  m->SetZero();
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( 0., m->Get(i,j) );
-
-  // set all
-  m->SetAll(1.5);
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( 1.5, m->Get(i,j) );
-
-  // set basis
-  m->SetIdentity();
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->Get(i,j) );
-
-  // set from array
-  double array[] = { 1., 0., 0.,
-                     0., 1., 0.,
-                     0., 0., 1.,
-                     0., 0., 0. };
-  m->SetFromDoubleArray(array);
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( i == j ? 1. : 0. , m->Get(i,j) );
-
-};
-
-/** Unit Test for copying matrices.
- *
- */
-void GSLMatrixTest::CopyTest()
-{
-  // set basis
-  GSLMatrix *dest = NULL;
-  for (int i=0;i<4;i++) {
-    m->SetAll(i);
-    dest = new GSLMatrix(m);
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( m->Get(i,j) , dest->Get(i,j) );
-
-    delete(dest);
-  }
-};
-
-/** Unit Test for exchanging rows and columns.
- *
- */
-void GSLMatrixTest::ExchangeTest()
-{
-  // set to 1,1,1,2, ...
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, i+1 );
-
-  // swap such that nothing happens
-  CPPUNIT_ASSERT_EQUAL( true, m->SwapColumns(1,2) );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->Get(i,j) );
-
-  // swap two rows
-  CPPUNIT_ASSERT_EQUAL( true, m->SwapRows(1,2) );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      switch (i) {
-        case 0:
-          CPPUNIT_ASSERT_EQUAL( 1., m->Get(i,j) );
-          break;
-        case 1:
-          CPPUNIT_ASSERT_EQUAL( 3., m->Get(i,j) );
-          break;
-        case 2:
-          CPPUNIT_ASSERT_EQUAL( 2., m->Get(i,j) );
-          break;
-        case 3:
-          CPPUNIT_ASSERT_EQUAL( 4., m->Get(i,j) );
-          break;
-        default:
-          CPPUNIT_ASSERT_EQUAL( -1., m->Get(i,j) );
-      }
-  // check that op is reversable
-  CPPUNIT_ASSERT_EQUAL( true, m->SwapRows(1,2) );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)i+1., m->Get(i,j) );
-
-  // set to 1,2,3,1, ...
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, j+1. );
-
-  // swap such that nothing happens
-  CPPUNIT_ASSERT_EQUAL( true, m->SwapRows(0,2) );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->Get(i,j) );
-
-  // swap two columns
-  CPPUNIT_ASSERT_EQUAL( true, m->SwapColumns(0,2) );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      switch (j) {
-        case 0:
-          CPPUNIT_ASSERT_EQUAL( 3., m->Get(i,j) );
-          break;
-        case 1:
-          CPPUNIT_ASSERT_EQUAL( 2., m->Get(i,j) );
-          break;
-        case 2:
-          CPPUNIT_ASSERT_EQUAL( 1., m->Get(i,j) );
-          break;
-        default:
-          CPPUNIT_ASSERT_EQUAL( -1., m->Get(i,j) );
-      }
-  // check that op is reversable
-  CPPUNIT_ASSERT_EQUAL( true, m->SwapColumns(0,2) );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)j+1., m->Get(i,j) );
-
-
-  // set to 1,2,3,4, ...
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      m->Set(i,j, 3*i+j+1 );
-  // transpose
-  CPPUNIT_ASSERT_EQUAL( true, m->Transpose() );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)(3*i+j+1), m->Get(j,i) );
-  // second transpose
-  CPPUNIT_ASSERT_EQUAL( true, m->Transpose() );
-  for (int i=0;i<4;i++)
-    for (int j=0;j<3;j++)
-      CPPUNIT_ASSERT_EQUAL( (double)(3*i+j+1), m->Get(i,j) );
-};
-
-/** Unit Test for matrix properties.
- *
- */
-void GSLMatrixTest::PropertiesTest()
-{
-  // is zero
-  m->SetZero();
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNull() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
-
-  // is positive
-  m->SetAll(0.5);
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsPositive() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNegative() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNonNegative() );
-
-  // is negative
-  m->SetAll(-0.1);
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNull() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsPositive() );
-  CPPUNIT_ASSERT_EQUAL( true, m->IsNegative() );
-  CPPUNIT_ASSERT_EQUAL( false, m->IsNonNegative() );
-
-  // is positive definite
-  CPPUNIT_ASSERT_EQUAL( false, m->IsPositiveDefinite() );
-};
Index: src/unittests/gslmatrixunittest.hpp
===================================================================
--- src/unittests/gslmatrixunittest.hpp	(revision cca9efb7bf18583f35d0380acb97e9034c42bbe7)
+++ 	(revision )
@@ -1,42 +1,0 @@
-/*
- * gslmatrixunittest.hpp
- *
- *  Created on: Jan 8, 2010
- *      Author: heber
- */
-
-#ifndef GSLMATRIXUNITTEST_HPP_
-#define GSLMATRIXUNITTEST_HPP_
-
-#include <cppunit/extensions/HelperMacros.h>
-
-#include "LinearAlgebra/gslmatrix.hpp"
-
-/********************************************** Test classes **************************************/
-
-class GSLMatrixTest : public CppUnit::TestFixture
-{
-    CPPUNIT_TEST_SUITE( GSLMatrixTest) ;
-    CPPUNIT_TEST (AccessTest );
-    CPPUNIT_TEST (InitializationTest );
-    CPPUNIT_TEST (CopyTest );
-    CPPUNIT_TEST (ExchangeTest );
-    CPPUNIT_TEST (PropertiesTest );
-    CPPUNIT_TEST_SUITE_END();
-
-public:
-    void setUp();
-    void tearDown();
-
-    void AccessTest();
-    void InitializationTest();
-    void CopyTest();
-    void ExchangeTest();
-    void PropertiesTest();
-
-private:
-    GSLMatrix *m;
-};
-
-
-#endif /* GSLMATRIXUNITTEST_HPP_ */
