Changeset 1068 for library/bdm/stat/exp_family.cpp

Timestamp:

06/09/10 18:35:14 (14 years ago)

Author:

mido

Message:

patch of documentation - all conditional pdfs revised

Files:

• library/bdm/stat/exp_family.cpp (modified) (6 diffs)

library/bdm/stat/exp_family.cpp

@@ -134 +134 @@
 
     double m = nu - nPsi - dimx - 1;
-#define log2  0.693147180559945286226763983
-#define logpi 1.144729885849400163877476189
+#define    log2  0.693147180559945286226763983
+#define    logpi 1.144729885849400163877476189
 #define log2pi 1.83787706640935
 #define Inf std::numeric_limits<double>::infinity()
@@ -149 +149 @@
                  - 0.5 * dimx * ( m * log2 + 0.5 * ( dimx - 1 ) * log2pi )  - lg;
 
-//      bdm_assert_debug ( ( ( -nkG - nkW ) > -Inf ) && ( ( -nkG - nkW ) < Inf ), "ARX improper" );
+//    bdm_assert_debug ( ( ( -nkG - nkW ) > -Inf ) && ( ( -nkG - nkW ) < Inf ), "ARX improper" );
     if ( -nkG - nkW == Inf ) {
         cout << "??" << endl;
@@ -186 +186 @@
         M = iLsub * Lpsi;
     }
-    /*  if ( 0 ) { // test with Peterka
-                mat VF = V.to_mat();
-                mat Vf = VF ( 0, dimx - 1, 0, dimx - 1 );
-                mat Vzf = VF ( dimx, end, 0, dimx - 1 );
-                mat VZ = VF ( dimx, end, dimx, end );
-
-                mat Lam2 = Vf - Vzf.T() * inv ( VZ ) * Vzf;
-        }*/
+    /*    if ( 0 ) { // test with Peterka
+            mat VF = V.to_mat();
+            mat Vf = VF ( 0, dimx - 1, 0, dimx - 1 );
+            mat Vzf = VF ( dimx, end, 0, dimx - 1 );
+            mat VZ = VF ( dimx, end, dimx, end );
+
+            mat Lam2 = Vf - Vzf.T() * inv ( VZ ) * Vzf;
+        }*/
 }
 
@@ -203 +203 @@
 
         mat Lsub = L ( 1, end, 1, end );
-//              mat Dsub = diag ( D ( 1, end ) );
+//        mat Dsub = diag ( D ( 1, end ) );
 
         ldmat LD ( inv ( Lsub ).T(), 1.0 / D ( 1, end ) );
@@ -715 +715 @@
 }
 
-void    euni::to_setting  (Setting  &set) const {
+void     euni::to_setting  (Setting  &set) const {
     epdf::to_setting ( set );
     UI::save ( high, set, "high" );
@@ -731 +731 @@
     pdf::from_setting(set);
     g=UI::build<fnc>(set,"g",UI::compulsory);
-    validate();
-}
+}
+
 void mgdirac::to_setting(Setting &set) const {
     pdf::to_setting(set);
     UI::save(g.get(), set, "g");
 }
+
 void mgdirac::validate() {
     pdf::validate();