Changeset 896 for library/bdm/estim

Timestamp:

04/09/10 09:31:37 (14 years ago)

Author:

mido

Message:

cleanup of MemDS and its descendants
bdmtoolbox/CMakeLists.txt slightly changed to avoid unnecessary MEX condition
"indeces" replaced by "indices"

Location:

library/bdm/estim

Files:

• arx.cpp (modified) (3 diffs)
• arx.h (modified) (1 diff)
• kalman.h (modified) (1 diff)

library/bdm/estim/arx.cpp

@@ -148 +148 @@
 @param Eg0 a copy of prior GiW density before estimation
 @param Egll likelihood of the current Eg
-@param indeces current indeces
+@param indices current indices
 \return best likelihood in the structure below the given one
 */
-double egiw_bestbelow ( egiw Eg, egiw Eg0, double Egll, ivec &indeces ) { //parameter Eg is a copy!
+double egiw_bestbelow ( egiw Eg, egiw Eg0, double Egll, ivec &indices ) { //parameter Eg is a copy!
         ldmat Vo = Eg._V(); //copy
         ldmat Vo0 = Eg._V(); //copy
@@ -164 +164 @@
         double belll = Egll;
 
-        ivec tmpindeces;
-        ivec maxindeces = indeces;
-
-
-        cout << "bb:(" << indeces << ") ll=" << Egll << endl;
+        ivec tmpindices;
+        ivec maxindices = indices;
+
+
+        cout << "bb:(" << indices << ") ll=" << Egll << endl;
 
         //try to remove only one rv
@@ -186 +186 @@
                 //
                 if ( tmpll > Egll ) { //increase of the likelihood
-                        tmpindeces = indeces;
-                        tmpindeces.del ( i );
+                        tmpindices = indices;
+                        tmpindices.del ( i );
                         //search for a better match in this substructure
-                        belll = egiw_bestbelow ( Eg, Eg0, tmpll, tmpindeces );
+                        belll = egiw_bestbelow ( Eg, Eg0, tmpll, tmpindices );
                         if ( belll > maxll ) { //better match found
                                 maxll = belll;
-                                maxindeces = tmpindeces;
+                                maxindices = tmpindices;
                         }
                 }
         }
-        indeces = maxindeces;
+        indices = maxindices;
         return maxll;
 }

library/bdm/estim/arx.h

@@ -97 +97 @@
         void ml_predictor_update ( mlnorm<sq_T> &pred ) const;
         mlstudent* predictor_student() const;
-        //! Brute force structure estimation.\return indeces of accepted regressors.
+        //! Brute force structure estimation.\return indices of accepted regressors.
         ivec structure_est ( egiw Eg0 );
-        //! Smarter structure estimation by Ludvik Tesar.\return indeces of accepted regressors.
+        //! Smarter structure estimation by Ludvik Tesar.\return indices of accepted regressors.
         ivec structure_est_LT ( egiw Eg0 );
         //!@}

library/bdm/estim/kalman.h

@@ -432 +432 @@
 Using Frobenius form, see [].
 
-For easier use in the future, indeces theta_in_A and theta_in_C are set. TODO - explain
+For easier use in the future, indices theta_in_A and theta_in_C are set. TODO - explain
 */
 //template<class sq_T>