Changeset 1478 for trunk/yat/utility/SVD.cc

Timestamp:

Sep 6, 2008, 8:55:46 PM (13 years ago)

Author:

Peter

Message:

reverting [1474] as original docs were correct. Clarified that SVD only works when rows>=columns for all algorithms and added an assert checking this.

File:

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

trunk/yat/utility/SVD.cc

@@ -31 +31 @@
 #include "VectorBase.h"
 
+#include <cassert>
 #include <sstream>
 
@@ -51 +52 @@
   void SVD::decompose(SVDalgorithm algo)
   {
+    assert(U_.rows()>=U_.columns());
     int status=0;
     switch (algo) {