Changeset 755 for trunk/yat/utility/vector.cc

Timestamp:

Feb 18, 2007, 1:01:39 AM (16 years ago)

Author:

Jari Häkkinen

Message:

Addresses #2 and #65. Continued adding GSL_error exceptions, and added doxygen briefs.

File:

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

trunk/yat/utility/vector.cc

@@ -219 +219 @@
 
 
-  int vector::div(const vector& other)
-  {
-    return gsl_vector_div(v_,other.v_);
+  void vector::div(const vector& other)
+  {
+    int status=gsl_vector_div(v_,other.v_);
+    if (status)
+      throw utility::GSL_error(std::string("vector::div",status));
   }
 
@@ -289 +291 @@
 
 
-  int vector::mul(const vector& other)
-  {
-    return gsl_vector_mul(v_,other.v_);
-  }
-
-
-  int vector::reverse(void)
-  {
-    return gsl_vector_reverse(v_);
+  void vector::mul(const vector& other)
+  {
+    int status=gsl_vector_mul(v_,other.v_);
+    if (status)
+      throw utility::GSL_error(std::string("vector::mul",status));
+  }
+
+
+  void vector::reverse(void)
+  {
+    gsl_vector_reverse(v_);
   }
 
@@ -363 +367 @@
 
 
-  int vector::swap(vector& other)
-  {
-    return gsl_vector_swap(v_,other.v_);
-  }
-
-
-  int vector::swap_elements(size_t i, size_t j)
-  {
-    return gsl_vector_swap_elements(v_,i,j);
+  void vector::swap(vector& other)
+  {
+    int status=gsl_vector_swap(v_,other.v_);
+    if (status)
+      throw utility::GSL_error(std::string("vector::swap",status));
+  }
+
+
+  void vector::swap_elements(size_t i, size_t j)
+  {
+    int status=gsl_vector_swap_elements(v_, i, j);
+    if (status)
+      throw utility::GSL_error(std::string("vector::swap_elements",status));
   }
 
@@ -377 +385 @@
   double& vector::operator()(size_t i)
   {
-    return *gsl_vector_ptr(v_,i);
+    double* d=gsl_vector_ptr(v_, i);
+    if (!d)
+      throw utility::GSL_error("vector::operator()",GSL_EINVAL);
+    return *d;
   }
 
@@ -383 +394 @@
   const double& vector::operator()(size_t i) const
   {
-    return *gsl_vector_const_ptr(proxy_v_,i);
+    const double* d=gsl_vector_const_ptr(proxy_v_, i);
+    if (!d)
+      throw utility::GSL_error("vector::operator()",GSL_EINVAL);
+    return *d;
   }