Changeset 772

Timestamp:

01/12/10 20:25:54 (14 years ago)

Author:

smidl

Message:

reordering in chmat

Files:

• library/bdm/math/chmat.h (modified) (1 diff)

library/bdm/math/chmat.h

@@ -76 +76 @@
           Some templates require this constructor to compile, but
           it shouldn't actually be called.
+          \note may be numerically unstable
         */
-        chmat ( const chmat &M, const ivec &perm ) NOT_IMPLEMENTED_VOID;
+        chmat ( const chmat &M, const ivec &perm ): sqmat(perm.length()) {
+                mat complete = M.to_mat ();
+                mat subs = complete.get_cols(perm);
+                Ch = chol(subs.get_rows(perm));
+        }
 
         //! Access function