Changeset 1099

Timestamp:

Feb 18, 2008, 4:42:46 AM (15 years ago)

Author:

Jari Häkkinen

Message:

Readded resize() to vector.

Location:

trunk/yat/utility

Files:

• vector.cc (modified) (2 diffs)
• vector.h (modified) (1 diff)

trunk/yat/utility/vector.cc

@@ -195 +195 @@
   {
     delete_allocated_memory();
-    if (!other.size())
-      vec_=NULL;
-    else  {   
+    if (other.size()) {
       if (size()!=other.size())
         vec_ = gsl_vector_alloc(other.size());
@@ -203 +201 @@
         throw utility::GSL_error("vector failed to allocate memory");
       gsl_vector_memcpy(vec_,other.gsl_vector_p());
+      const_vec_ = vec_;
     }
-    const_vec_ = vec_;
     return *this;
   } 
+
+
+  void vector::resize(size_t n, double init_value)
+  {
+    delete_allocated_memory();
+    const_vec_ = vec_ = gsl_vector_alloc(n);
+    if (!vec_)
+      throw utility::GSL_error("vector::vector failed to allocate memory");
+    all(init_value);
+  }
 
 

trunk/yat/utility/vector.h

@@ -116 +116 @@
 
     /**
+      \brief Resize vector
+
+      All elements are set to \a init_value.
+
+      \note Underlying GSL vector is destroyed and a view into this
+      vector becomes invalid.
+    */
+    void resize(size_t, double init_value=0);
+
+    /**
        \brief The assignment operator.