Changeset 793 for trunk/test/matrix_test.cc

Timestamp:

Mar 11, 2007, 4:40:13 PM (17 years ago)

Author:

Jari Häkkinen

Message:

Fixes #193. Last test for the new functionality added to matrix.

File:

• trunk/test/matrix_test.cc (modified) (3 diffs)

trunk/test/matrix_test.cc

@@ -119 +119 @@
   }
 
+  // Test of const view implementation (make sure no zero pointer is
+  // used). Here we use the bad style of not making the view const!
+  // If test fails with a null pointer exception it is an uncatchable
+  // core dump!
+  {
+    *error << "\tconst view implementation" << std::endl;
+    // Change the next statement to
+    // const utility::matrix m(10,10,3.0);
+    // when ticket:202 is fixes
+    utility::matrix m(10,10,3.0);
+    utility::matrix mview(m,3,3,3,3);
+    // const utility::vector vview(vv,0,5,1); // this is the proper line
+    utility::matrix m2(3,3,2.0);
+    m2.mul(mview); // should work even without const since const arg passing
+    m2.div(mview); // should work even without const since const arg passing
+  }
+
+  // checking that copy constructor creates an independent object when
+  // a non-view matrix is copied
+  {
+    *error << "\tcopy constructor" << std::endl;
+    utility::matrix m2(m5);
+    ok &= (m2.rows()==m5.rows());
+    ok &= (m2.columns()==m5.columns());
+    ok &= (m2==m5);
+    ok &= (&m2 != &m5);
+    ok &= !m2.isview();
+  }
+
+  // checking that copy constructor creates an independent object when
+  // a view matrix is copied
+  {
+    *error << "\tcopy contructor on view" << std::endl;
+    // Change the next statement to
+    // const utility::matrix m(10,10,3.0);
+    // when ticket:202 is fixes
+    utility::matrix m(10,10,3.0);
+    utility::matrix mview(m,3,3,3,3);
+    utility::matrix m3(mview);
+    ok &= (mview.rows()==m3.rows());
+    ok &= (mview.columns()==m3.columns());
+    ok &= (m3 == mview);
+    ok &= (&m3 != &mview);
+    ok &= !m3.isview();
+  }
+
+  // checking that assignment operator throws an exception if matrices
+  // differ in size
+  {
+    *error << "\tassignment operator" << std::endl;
+    // GSL will catch the error in this test there for the GSL error
+    // handling must be disabled until after the exception is
+    // catched. The GSL error handler is reinstated after the
+    // try-catch construct.
+    gsl_error_handler_t* err_handler=gsl_set_error_handler_off();
+    bool exception_happens=false;
+    try {
+      utility::matrix m(m5.rows()+1,3,0.0);
+      m=m5;
+    } catch (utility::GSL_error& err) {
+      exception_happens=true;
+    }
+    if (!exception_happens) {
+      *error << "Matrix assignment operator did not throw expected exception"
+             << std::endl;
+      ok=false;
+    }
+    gsl_set_error_handler(err_handler);
+  }
+
+  // checking that assignment operator changes the underlying object when
+  // a view is changed.
+  {
+    *error << "\tassignment operator on view" << std::endl;
+    bool this_ok(true);
+    utility::matrix mat_view(m5,3,3,3,3);
+    utility::matrix m2(3,3,12.0);
+    mat_view=m2;
+    for (u_int i=0; i<mat_view.rows(); ++i)
+      for (u_int j=0; j<mat_view.columns(); ++j)
+        if (m5(i+3,j+3)!=mat_view(i,j))
+          this_ok=false;
+    if (!this_ok) {
+      *error << "FAIL: assignemnt operator on view" << std::endl;
+      ok=false;
+    }
+  }
+
+  // checking clone functionality
+  {
+    *error << "\tclone functionality" << std::endl;
+    bool this_ok=true;
+    *error << "\t\tcloning normal matrix" << std::endl;
+    utility::matrix mat2;
+    mat2.clone(m5);
+    if (mat2.rows()!=m5.rows() || mat2.columns()!=m5.columns())
+      this_ok=false;
+    else
+      for (u_int i=0; i<m5.rows(); ++i)
+        for (u_int j=0; j<m5.columns(); ++j)
+          if (mat2(i,j)!=m5(i,j))
+            this_ok=false;
+    if (mat2.gsl_matrix_p()==m5.gsl_matrix_p())
+      this_ok=false;
+    *error << "\t\tcloning matrix view (sub matrix)" << std::endl;
+    utility::matrix* mat_view=new utility::matrix(m5,3,3,3,3);
+    utility::matrix mat_view2;
+    mat_view2.clone(*mat_view);
+    if (!mat_view2.isview())
+      this_ok=false;
+    if ( (mat_view->rows()!=mat_view2.rows()) ||
+         (mat_view->columns()!=mat_view2.columns()) )
+      this_ok=false;
+    else
+      for (u_int i=0; i<mat_view2.rows(); ++i)
+        for (u_int j=0; j<mat_view2.columns(); ++j)
+          if ((*mat_view)(i,j)!=mat_view2(i,j))
+            this_ok=false;
+    *error << "\t\tcloned matrix view independence" << std::endl;
+    delete mat_view;
+    for (u_int i=0; i<mat_view2.rows(); ++i)
+      for (u_int j=0; j<mat_view2.columns(); ++j)
+        if (m5(i+3,j+3)!=mat_view2(i,j))
+          this_ok=false;
+    if (!this_ok) {
+      *error << "FAIL: clone test" << std::endl;
+      ok=false;
+    }
+  }
+
   *error << "\tsub-(row)vector" << std::endl;
   // Checking that the row view works, i.e. mutation to the view are
   // reflected in the viewed object.
+  m5_sum2=0;
+  for (size_t i=0; i<m5.rows(); ++i)
+    for (size_t j=0; j<m5.columns(); ++j)
+      m5_sum2+=m5(i,j);
   utility::vector v5subrow(m5,3);
   double v5subrow_sum=0;
@@ -140 +274 @@
   // Checking that the column view works, i.e. mutation to the view
   // are reflected in the viewed object.
+  m5_sum3=0;
+  for (size_t i=0; i<m5.rows(); ++i)
+    for (size_t j=0; j<m5.columns(); ++j)
+      m5_sum3+=m5(i,j);
   utility::vector v5subcolumn(m5,0,false);
   double v5subcolumn_sum=0;
@@ -154 +292 @@
     *error << "error sub-vector test 2" << std::endl;
   }
+
   // Checking that the column view above mutates the values in the row
   // view.