int main(){
CRNG rng;
scalar_t m1[2][1] = {0, 0};
Mat mean_matrix(2, 1, CV_64FC1, m1);
scalar_t m2[2][2] = {small_deviation, 0, 0, small_deviation};
Mat sigma_matrix(2, 2, CV_64FC1, m2);
CCSV file1(".\\Matlab\\test1.csv");
CCSV file2(".\\Matlab\\test2.csv");
for(int i = 0; i < 1000; ++i){
vector<scalar_t> mean(2); mean[0] = 0; mean[1] = 0;
vector<scalar_t> d(2); d[0] = 0; d[1] = 0;
vector<vector<scalar_t> > body(100);
for(int j = 0; j < 100; ++j){
vector<scalar_t> stem = get(mean_matrix, sigma_matrix, rng);
mean[0] += stem[0]; mean[1] += stem[1];
body[j] = stem;
}
mean[0] /= 100; mean[1] /= 100;
file1.put(mean);
for(int j = 0; j < 100; ++j){
d[0] += (body[j][0]-mean[0])*(body[j][0]-mean[0]); d[1] += (body[j][1]-mean[1])*(body[j][1]-mean[1]);
}
d[0] /= 100; d[1] /= 100;
file2.put(d);
}
}
shared_ptr<arma::mat> mean_matrix() {
return mean_matrix(all_examples);
}
