hadamard_test.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25

#include "hadamard_pt.hpp"
#include <iostream>

int main() {
  double β0 = 1.0;
  double β1 = 10.0;
  unsigned N = 16;
  unsigned n = 4;

  ParallelMeasurement m1;
  std::vector<Measurement> m2(N);
  std::vector<Measurement*> m3;
  for (Measurement& m : m2) {
    m3.push_back(&m);
  }

  PT a(β0, β1, N, n, m1, m3);

  a.tune(1000, 10, 0.1, 0.1);

  for (MCMC& M : a.Ms) {
    std::cout << M.β << " ";
  }
  std::cout << "\n";
}