Otros ¿Sabéis de algún código para Ising 2D multiparalelo?

[idontknow]

Estoy haciendo una práctica en la uni en la que nos han dado un código en fortran para calcular un modelo de N espines de Ising 2D de red cuadrada y lado L. El problema es que cerca del punto crítico como hay muchas fluctuaciones se necesita poner un gran número de iteraciones y esto es lentísimo. He mirado por google códigos y bueno he encontrado alguno programado en c++ para cuda. En especial, he encontrado esta tesis (ir a la página 37) pero este código para mi es algo ilegible. Como sólo quiero ejecutarlo, la verdad que en estos momentos me da igual entender ese código pero claro, no entiendo que narices significan los parámetros EMD, ranD y por qué le pasa el int spins si después no lo utiliza y usa block_size... Indiferentemente de este código, ¿Os sabéis algún código los más mayorcitos en esto?

Opero con windows y con linux via máquina virtual. Tengo IDE de c++, java (a parte de fortran que no necesita).