Did the inner (that is, aij) summation in the generic RK method, and added stream fusion to the required libs