Try with OpenMP - not yet finished!