Changeset 980

Timestamp:

05/25/10 20:55:20 (15 years ago)

Author:

smidl

Message:

new way of copying data: vec_from_vec

Files:

• library/bdm/base/bdmbase.h (modified) (3 diffs)

library/bdm/base/bdmbase.h

@@ -1020 +1020 @@
 };
 
+class vec_from_vec: public vec {
+        protected:
+                datalink_part dl;
+        public:
+                void update(const vec &v1){
+                        bdm_assert_debug(length()==dl._downsize(),"vec_from_vec incompatible");
+                        bdm_assert_debug(v1.length()>=dl._upsize(),"vec_from_vec incompatible");//does not guarantee anything!
+                        dl.filldown(v1,*this);
+                };
+                void connect(const RV &vecrv, const RV & vec1){dl.set_connection(vecrv,vec1);set_length(vecrv._dsize());};
+};
+
+class vec_from_2vec: public vec {
+        protected:
+                datalink_part dl1;
+                datalink_part dl2;
+        public:
+                void update(const vec &v1, const vec &v2){
+                        bdm_assert_debug(length()==dl1._downsize()+dl2._downsize(),"vec_from_vec incompatible");
+                        bdm_assert_debug(v1.length()>=dl1._upsize(),"vec_from_vec incompatible");
+                        bdm_assert_debug(v2.length()>=dl2._upsize(),"vec_from_vec incompatible");
+                        dl1.filldown(v1,*this);
+                        dl2.filldown(v2,*this);
+                };
+                void connect(const RV &rv, const RV & rv1, const RV &rv2){
+                        dl1.set_connection(rv,rv1);
+                        dl2.set_connection(rv,rv2);
+                        set_length(rv._dsize());
+                };
+};
+
 //! buffered datalink from 2 vectors to 1
 class datalink_2to1_buffered {
@@ -1382 +1413 @@
                 UI::save( &rvc, set, "rvc" );           
                 UI::save( &posterior()._rv(), set, "rv" );
+                if (prior0) UI::save( prior0, set, "prior" );
                 UI::save( log_level, set );
         }
@@ -1392 +1424 @@
                         const_cast<epdf&> ( posterior() ).log_level[epdf::logmean] = true;;
                 }
-                if (prior0)
-                        set_prior(prior0.get());
+                set_prior(prior0.get());
         }
 };