StatErrEst.cpp

Go to the documentation of this file.

#include "StatErrEst.h"

/*
void StatErrEst::MakeSamples()
{
        boost::lagged_fibonacci607 generator(static_cast<unsigned int>(std::time(0)));
        boost::normal_distribution<> Zxxr_dist(Z.GetZxx().real(),Z.GetdZxx()), Zxxi_dist(Z.GetZxx().imag(),Z.GetdZxx());
        boost::normal_distribution<> Zxyr_dist(Z.GetZxy().real(),Z.GetdZxy()), Zxyi_dist(Z.GetZxy().imag(),Z.GetdZxy());
        boost::normal_distribution<> Zyxr_dist(Z.GetZyx().real(),Z.GetdZyx()), Zyxi_dist(Z.GetZyx().imag(),Z.GetdZyx());
        boost::normal_distribution<> Zyyr_dist(Z.GetZyy().real(),Z.GetdZyy()), Zyyi_dist(Z.GetZyy().imag(),Z.GetdZyy());
   
   
    boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<> > Zxxr(generator, Zxxr_dist);
        boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<> > Zxxi(generator, Zxxi_dist);
        boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<>  > Zxyr(generator, Zxyr_dist);
        boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<>  > Zxyi(generator, Zxyi_dist);
        boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<> > Zyxr(generator, Zyxr_dist);
        boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<> > Zyxi(generator, Zyxi_dist);
        boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<> > Zyyr(generator, Zyyr_dist);
        boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<> > Zyyi(generator, Zyyi_dist);
        
        Samples.reserve(nrealizations);
        for (int i = 0; i < nrealizations; ++i)
        {
                //MTTensor Point(std::complex<double>(Zxxr(),Zxxi()),std::complex<double>(Zxyr(),Zxyi()),
                //                      std::complex<double>(Zyxr(),Zyxi()),std::complex<double>(Zyyr(),Zyyi()),Z.GetFrequency());
                //std::cout <<"Point: " << Point.GetZxy() << " " << Point.GetRhoxy() << std::endl;
                Samples.push_back(Generator());
                //std::cout << value << std::endl;
        }
}*/