Changeset 1172

Timestamp:

Feb 27, 2008, 3:32:15 PM (16 years ago)

Author:

Peter

Message:

allowing resize to zero size. In Matrix if resize is called with only one of columns or rows set to zero a zero by zero matrix is created. I think we can live with that. Though an assertion is called because it is strange.

Location:

trunk/yat/utility

Files:

• Matrix.cc (modified) (1 diff)
• Vector.cc (modified) (1 diff)

trunk/yat/utility/Matrix.cc

@@ -320 +320 @@
     delete_allocated_memory();
 
-    m_ = gsl_matrix_alloc(r,c);
-    if (!m_)
-      throw utility::GSL_error("Matrix::Matrix failed to allocate memory");
-    all(init_value);
+    assert( (r&&c) || (!r&&!c) );
+    if (r && c){
+      m_ = gsl_matrix_alloc(r,c);
+      if (!m_)
+        throw utility::GSL_error("Matrix::Matrix failed to allocate memory");
+      all(init_value);
+    }
 
     // no need to delete blas_result_ if the number of rows fit, it

trunk/yat/utility/Vector.cc

@@ -220 +220 @@
   {
     delete_allocated_memory();
+    if (!n)
+      return;
     const_vec_ = vec_ = gsl_vector_alloc(n);
     if (!vec_)