Changeset 180 for bdm/math

Timestamp:

10/15/08 19:08:06 (16 years ago)

Author:

smidl

Message:

Modifications of BDM to reflect changes in basics

Files:

• bdm/math/libDC.h (modified) (7 diffs)

bdm/math/libDC.h

@@ -15 +15 @@
 
 #include <itpp/itbase.h>
+#include "../itpp_ext.h"
 
 using namespace itpp;
@@ -172 +173 @@
 };
 
-/*! \brief Matrix stored in LD form, (typically known as UD) 
+/*! \brief Matrix stored in LD form, (commonly known as UD) 
 
 Matrix is decomposed as follows: \f[M = L'DL\f] where only \f$L\f$ and \f$D\f$ matrices are stored.
@@ -180 +181 @@
 {
         public:
-
                 //! Construct by copy of L and D.
                 ldmat ( const mat &L, const vec &D );
                 //! Construct by decomposition of full matrix V.
                 ldmat (const mat &V );
+                //! Construct by restructuring of V0 accordint to permutation vector perm.
+                ldmat (const ldmat &V0, const ivec &perm):sqmat(V0.rows()){     ldform(V0.L.get_cols(perm), V0.D);};
                 //! Construct diagonal matrix with diagonal D0
                 ldmat ( vec D0 );
@@ -191 +193 @@
                 //! Default initialization with proper size
                 ldmat(const int dim0);
-
+                
                 //! Destructor for future use;
                 virtual ~ldmat(){};
@@ -206 +208 @@
                 double qform (const vec &v ) const;
                 double invqform (const vec &v ) const;
-//      sqmat& operator -= ( const sqmat & ld2 );
                 void clear();
                 int cols() const;
@@ -212 +213 @@
                 vec sqrt_mult ( const vec &v ) const;
 
+                
                 /*! \brief Matrix inversion preserving the chosen form.
                 @param Inv a space where the inverse is stored.
@@ -267 +269 @@
 
 };
-
 
 //////// Operations: