still not faster than fortran :/