Changeset 1190 for trunk/yat/utility/VectorMutable.cc

Timestamp:

Feb 29, 2008, 7:17:25 PM (16 years ago)

Author:

Peter

Message:

adding missing implementation of operator += and added some throwing

File:

• trunk/yat/utility/VectorMutable.cc (modified) (2 diffs)

trunk/yat/utility/VectorMutable.cc

@@ -122 +122 @@
 
 
+  const VectorMutable& VectorMutable::operator+=(const VectorBase& other)
+  {
+    if (size()!=other.size())
+      throw std::runtime_error("VectorMutable::operator+= size must be same.");
+    assert(vec_);
+    int status=gsl_vector_add(vec_, other.gsl_vector_p());
+    if (status)
+      throw utility::GSL_error(std::string("VectorMutable::sub", status));
+    return *this;
+  }
+
+
   const VectorMutable& VectorMutable::operator+=(double d)
   {
@@ -132 +144 @@
   const VectorMutable& VectorMutable::operator-=(const VectorBase& other)
   {
+    if (size()!=other.size())
+      throw std::runtime_error("VectorMutable::operator-= size must be same.");
     assert(vec_);
     int status=gsl_vector_sub(vec_, other.gsl_vector_p());