Changeset 1083 for library/bdm/stat/merger.h

Timestamp:

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

Author:

smidl

Message:

Mergers + fixes

Files:

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

library/bdm/stat/merger.h

@@ -38 +38 @@
                 //! 
                 //! check if all epdfs are on the same support
-                void validate_sources() {
-                        int n=sources.length();
+                virtual void set_sources(const Array<shared_ptr<pdf> > &A) {
+                        int n=A.length();
                         
                         bdm_assert(n>0,"merger has no sources to merge");
@@ -48 +48 @@
                         
                         // check compatibility of sources -- no types are needed
-                        int dim0 = sources(0)->dimension();
+                        int dim0 = A(0)->dimension();
                         for (int i=0; i<n; i++){
+                                const epdf *si=dynamic_cast<const epdf*>(A(i).get());
+                                if (si){
+                                        sources(i) = dynamic_cast<epdf*>(A(i)->_copy());
+                                }
                                 bdm_assert(sources(i)->dimension()==dim0, "Merger: Incompatible dimensions of sources");
                                 if (sources(i)->isnamed()){