Changeset 175 for bdm/stat/libBM.cpp

Timestamp:

10/09/08 11:32:11 (16 years ago)

Author:

smidl

Message:

mprod is now based on compositepdf

Files:

• bdm/stat/libBM.cpp (modified) (1 diff)

bdm/stat/libBM.cpp

@@ -183 +183 @@
 void RV::newids(){ids=linspace ( RVcounter+1, RVcounter+len ),RVcounter+=len;}
 
+RV compositepdf::getrv(bool checkoverlap){
+        RV rv; //empty rv
+        bool rvaddok;
+        for (int i = 0;i < n;i++ ) {
+                rvaddok=rv.add ( mpdfs ( i )->_rv() ); //add rv to common rvs.
+                        // If rvaddok==false, mpdfs overlap => assert error.
+                it_assert_debug(rvaddok||(!checkoverlap),"mprod::mprod() input mpdfs overlap in rv!");
+        };
+        return rv;
+}
+
+void compositepdf::setrvc(const RV &rv, RV &rvc){
+        for (int i = 0;i < n;i++ ) {
+                rvc.add ( mpdfs ( i )->_rvc().subt ( rv ) ); //add rv to common rvc
+        };
+}
+
+void compositepdf::setindices(const RV &rv){
+        for (int i = 0;i < n;i++ ) {
+                        // find ith rv in common rv
+                rvsinrv ( i ) = mpdfs ( i )->_rv().dataind ( rv );
+                rvcsinrv ( i ) = mpdfs ( i )->_rvc().dataind ( rv );
+        }
+}
+
+void compositepdf::setrvcinrv(const RV &rvc, Array<ivec> &rvinds){
+        for (int i = 0;i < n;i++ ) {
+                        // find ith rv in common rv
+                rvinds ( i ) = mpdfs ( i )->_rvc().dataind ( rvc );
+        }
+}
+
+
 void BM::bayesB(const mat &Data){
         for(int t=0;t<Data.cols();t++){bayes(Data.get_col(t));}