root/library/tests/testResample.cpp @ 384

#include "estim/particles.h"
#include "math/square_mat.h"
#include "stat/exp_family.h"

using namespace bdm;

//These lines are needed for use of cout and endl
using std::cout;
using std::endl;

int main() {

        RV x("1");
        RV xm=x; xm.t(-1);const
        RV y("2");
        
        mat A = "1";
        vec vR = "1";
        ldmat R(vR);
                
        eEmp emp;
        euni eun;
        eun.set_parameters("0","1");
        emp.set_statistics(ones(10),&eun);
        vec &v=emp._w();
        Array<vec> &S=emp._samples();
        
        for (int i=0;i<10;i++){ v(i) = exp(-0.5*sum(pow(S(i)-1,2.0))*10);}
        v/=sum(v);
        
        cout << "p:" << S << endl;
        cout << "w:" << v << endl;
        
        ivec ind = emp.resample();
        
        cout << ind << endl;
        
        //Exit program:
        return 0;

}