def timing_demo1(N=12): domain = 2**np.arange(1, N + 1) times = [] for n in domain: start = time() solutions.random_matrix(n) times.append(time() - start) plt.plot(domain, times, 'g.-', linewidth=2, markersize=15) plt.xlabel("n", fontsize=14) plt.ylabel("Seconds", fontsize=14) return domain, times
def timing_demo1(N=12): domain = 2**np.arange(1,N+1) times = [] for n in domain: start = time() solutions.random_matrix(n) times.append(time() - start) plt.plot(domain, times, 'g.-', linewidth=2, markersize=15) plt.xlabel("n", fontsize=14) plt.ylabel("Seconds", fontsize=14) return domain, times
def prob1A(N=9): domain = 2**np.arange(1,N+1) vector_times, matrix_times = [], [] for n in domain: A = solutions.random_matrix(n) x = solutions.random_vector(n) B = solutions.random_matrix(n) start = time() solutions.matrix_vector_product(A, x) vector_times.append(time() - start) start = time() solutions.matrix_matrix_product(A, B) matrix_times.append(time() - start) plt.plot(domain, vector_times, 'b.-', lw=2, ms=15) plt.xlabel("n", fontsize=14); plt.ylabel("Seconds", fontsize=14) plt.title("Matrix-Vector Multiplication") return domain, vector_times, matrix_times
def prob1A(N=9): domain = 2**np.arange(1, N + 1) vector_times, matrix_times = [], [] for n in domain: A = solutions.random_matrix(n) x = solutions.random_vector(n) B = solutions.random_matrix(n) start = time() solutions.matrix_vector_product(A, x) vector_times.append(time() - start) start = time() solutions.matrix_matrix_product(A, B) matrix_times.append(time() - start) plt.plot(domain, vector_times, 'b.-', lw=2, ms=15) plt.xlabel("n", fontsize=14) plt.ylabel("Seconds", fontsize=14) plt.title("Matrix-Vector Multiplication") return domain, vector_times, matrix_times