root/library/tests/stresssuite/resample_stress.cpp

#include "../mat_checks.h"
#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;

TEST ( resample_stress ) {

    RV x ( "1" );
    RV xm = x;
    xm.t_plus ( -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;
    resample ( v,ind );

    cout << ind << endl;
}