Changeset 762 for trunk/yat/utility/matrix.h

Timestamp:

Feb 20, 2007, 8:18:48 PM (16 years ago)

Author:

Jari Häkkinen

Message:

Fixes #76. Creating/recreating BLAS support matrix only when needed.

File:

• trunk/yat/utility/matrix.h (modified) (3 diffs)

trunk/yat/utility/matrix.h

@@ -67 +67 @@
   /// the omitted functionality could be included.
   ///
-  /// @todo Maybe it would be smart to create temporary objects need
-  /// for BLAS in constructors. As it is now temporary objects are
-  /// called before BLAS functionality iss used, cf. const matrix&
-  /// matrix::operator*=(const matrix& other)
-  ///
   class matrix
   {
@@ -323 +318 @@
        \brief Swap row \a i and column \a j.
 
+       The matrix must be square.
+
        \throw GSL_error if either index is out of bounds, or if matrix
        is not square.
@@ -457 +454 @@
     gsl_matrix* create_gsl_matrix_copy(void) const;
 
+    // blas_result_ is used to temporarily store result in BLAS calls
+    // and when not NULL it should always have the same size as
+    // m_. Memory is not allocated for blas_result_ until it is
+    // needed.
+    gsl_matrix* blas_result_;
     gsl_matrix* m_;
     gsl_matrix_view* view_;