replace doNormal by doNormalMulti

src/doNormalMulti.h

@@ -0,0 +1,44 @@
+#ifndef _doNormalMulti_h
+#define _doNormalMulti_h
+
+#include <boost/numeric/ublas/symmetric.hpp>
+#include <boost/numeric/ublas/lu.hpp>
+
+#include "doDistrib.h"
+
+namespace ublas = boost::numeric::ublas;
+
+template < typename EOT >
+class doNormalMulti : public doDistrib< EOT >
+{
+public:
+    typedef typename EOT::AtomType AtomType;
+
+    doNormalMulti
+    (
+     const ublas::vector< AtomType >& mean,
+     const ublas::symmetric_matrix< AtomType, ublas::lower >& varcovar
+     )
+	: _mean(mean), _varcovar(varcovar)
+    {
+	assert(_mean.size() > 0);
+	assert(_mean.size() == _varcovar.size1());
+	assert(_mean.size() == _varcovar.size2());
+    }
+
+    unsigned int size()
+    {
+	assert(_mean.size() == _varcovar.size1());
+	assert(_mean.size() == _varcovar.size2());
+	return _mean.size();
+    }
+
+    ublas::vector< AtomType >& mean(){return _mean;}
+    ublas::symmetric_matrix< AtomType, ublas::lower >& varcovar(){return _varcovar;}
+
+private:
+    ublas::vector< AtomType > _mean;
+    ublas::symmetric_matrix< AtomType, ublas::lower > _varcovar;
+};
+
+#endif // !_doNormalMulti_h