def prepare_ivi(N=100, dtype=np.double): A = np.random.rand(3260, 3260) B = np.random.rand(3260, 3000) return (A, B) #return 1.0*toc/(trials), 2*N*N*N*1e-9, times