[MaL] copyValues resizes target vector

Update and unify documentation.

MathLib/LinAlg/Eigen/EigenVector.h

@@ -19,6 +19,8 @@
 #include <Eigen/Eigen>
 #include <Eigen/Sparse>
 
+#include "EigenMapTools.h"
+
 namespace MathLib
 {
 /// Global vector based on Eigen vector
@@ -90,11 +92,13 @@ public:
         }
     }
 
-    /// Copy vector values.
+    /// Copy local entries to a vector.
+    /// \param u a vector for the values of local entries. It will be resized to
+    /// hold the current vector data.
     void copyValues(std::vector<double>& u) const
     {
-        assert(u.size() == (std::size_t)vec_.size());
-        copy_n(vec_.data(), vec_.size(), u.begin());
+        u.resize(size());
+        toVector(u) = vec_;
     }
 
 #ifndef NDEBUG

MathLib/LinAlg/Lis/LisVector.h

@@ -99,10 +99,12 @@ public:
         }
     }
 
-    /// Copy vector values.
+    /// Copy local entries to a vector.
+    /// \param u a vector for the values of local entries. It will be resized to
+    /// hold the current vector data.
     void copyValues(std::vector<double>& u) const
     {
-        assert(u.size() == size());
+        u.resize(size());
         lis_vector_get_values(vec_, 0, size(), u.data());
     }
 

MathLib/LinAlg/PETSc/PETScVector.cpp

@@ -196,7 +196,7 @@ void PETScVector::setLocalAccessibleVector() const
 
 void PETScVector::copyValues(std::vector<PetscScalar>& u) const
 {
-    assert(u.size() == (std::size_t)(getLocalSize() + getGhostSize()));
+    u.resize(getLocalSize() + getGhostSize());
 
     PetscScalar* loc_x = getLocalVector();
     std::copy_n(loc_x, getLocalSize() + getGhostSize(), u.begin());

MathLib/LinAlg/PETSc/PETScVector.h

@@ -177,8 +177,9 @@ public:
      */
     PETSc_Vec const& getRawVector() const { return v_; }
     /*!
-       Copy local entries including ghost ones to an array
-       \param u Preallocated vector for the values of local entries.
+       Copy local entries including ghost ones to a vector.
+       \param u a vector for the values of local entries. It will be resized to
+       hold the current vector data.
     */
     void copyValues(std::vector<PetscScalar>& u) const;