Comment Re:Ehhh Meh (Score 1) 127
Those numeric and simulation libraries were optimized in conjunnction with the Fortran compiler to take advantage of the hardware. The most obvious example; having fixed sized multi-dimensional arrays as global variables. For regular grids, the compiler can then decide which way to slice that data block up so that every processing node gets assigned a chunk of data. Since each function is not more than a few loop counters reading previous and current state for each grid cell, those get optimized into parallelised for-each calls. Once the scientists and engineers started writing and sharing numeric and simulation libraries and building on top of them, there isn't any need to change.