Changeset 787 for library/bdm/stat

Timestamp:

01/14/10 21:37:57 (14 years ago)

Author:

smidl

Message:

arena experiment + numerical fixes

Files:

• library/bdm/stat/merger.cpp (modified) (2 diffs)

library/bdm/stat/merger.cpp

@@ -163 +163 @@
         mat Tmp = zeros ( dim, dim );
         for ( int i = 0; i < Npoints; i++ ) {
-                Tmp += w ( i ) * outer_product ( S ( i ), S ( i ) );
-        }
-        return Tmp - outer_product ( mea, mea );
+                vec tmp=S ( i )-mea; //inefficient but numerically stable
+                Tmp += w ( i ) * outer_product (tmp , tmp );
+        }
+        return Tmp;
 }
 
@@ -174 +175 @@
 void merger_mix::merge ( ) {
         if(Npoints<1){
-                set_support(enorm<fsqmat>(zeros(dim), eye(dim)), 100);
+                set_support(enorm<fsqmat>(zeros(dim), eye(dim)), 1000);
         }