Changeset 758

Timestamp:

Feb 19, 2007, 7:52:16 PM (17 years ago)

Author:

Jari Häkkinen

Message:

Fixes #2. Finalized adding GSL_error exceptions.

File:

• trunk/yat/utility/vector.cc (modified) (5 diffs)

trunk/yat/utility/vector.cc

@@ -92 +92 @@
                               gsl_matrix_row   (m.gsl_matrix_p(),i) :
                               gsl_matrix_column(m.gsl_matrix_p(),i)  );
+    if (!view_)
+      throw utility::GSL_error("vector::vector failed to setup view");
     proxy_v_ = v_ = &(view_->vector);
   }
@@ -102 +104 @@
                                    gsl_matrix_const_row   (m.gsl_matrix_p(),i) :
                                    gsl_matrix_const_column(m.gsl_matrix_p(),i) );
+    if (!view_const_)
+      throw utility::GSL_error("vector::vector failed to setup view");
     proxy_v_ = v_const_ = &(view_const_->vector);
   }
@@ -403 +407 @@
   double& vector::operator[](size_t i)
   {
-    return *gsl_vector_ptr(v_,i);
+    return this->operator()(i);
   }
 
@@ -409 +413 @@
   const double& vector::operator[](size_t i) const
   {
-    return *gsl_vector_const_ptr(proxy_v_,i);
+    return this->operator()(i);
   }
 
@@ -417 +421 @@
     if (size()!=a.size())
       return false;
+    // if gsl error handler disabled, out of bounds index will not
+    // abort the program.
     for (size_t i=0; i<size(); ++i)
       if (gsl_vector_get(proxy_v_,i)!=a(i))