Programming: Monte-Carlo simulation of spin-glass model

Write a complete Monte-Carlo simulation program for the Sherrington–Kirkpatrick model of spin glass.

The model is defined by the Hamiltonian, \[H=-\sum_{i<j}J_{i,j}s_i s_j,\] where $s_i=\pm 1$ with $i=0,..,N-1$ are the $N$ spins of the system.

Both structure ($h_i$,$J_{i,j}$) and state ($s_i$) data are stored in one dimensional arrays.

int sz = 64;
int * spin = new int[sz];
double * h = new double [sz];
double * j = new double [sz*(sz-1)/2];