Я бы предпочел использовать C++ для выполнения кода, но я открыт для любых предложений по лучшему языку для данной ситуации. По сути, я хочу использовать алгоритм Штрассена для решения матриц, и я хочу знать, как я буду решать матрицы и измерять время их выполнения. # Версия 3.6
import numpy as np
def split(matrix):
"""
Splits a given matrix into quarters.
Input: nxn matrix
Output: tuple containing 4 n/2 x n/2 matrices corresponding to a, b, c, d
"""
row, col = matrix.shape
row2, col2 = row//2, col//2
return matrix[:row2, :col2], matrix[:row2, col2:], matrix[row2:, :col2],
matrix[row2:, col2:]
def strassen(x, y):
"""
Computes matrix product by divide and conquer approach, recursively.
Input: nxn matrices x and y
Output: nxn matrix, product of x and y
"""
# Base case when size of matrices is 1x1
if len(x) == 1:
return x * y
# Splitting the matrices into quadrants. This will be done recursively
# untill the base case is reached.
a, b, c, d = split(x)
e, f, g, h = split(y)
# Computing the 7 products, recursively (p1, p2...p7)
p1 = strassen(a, f - h)
p2 = strassen(a + b, h)
p3 = strassen(c + d, e)
p4 = strassen(d, g - e)
p5 = strassen(a + d, e + h)
p6 = strassen(b - d, g + h)
p7 = strassen(a - c, e + f)
# Computing the values of the 4 quadrants of the final matrix c
c11 = p5 + p4 - p2 + p6
c12 = p1 + p2
c21 = p3 + p4
c22 = p1 + p5 - p3 - p7
# Combining the 4 quadrants into a single matrix by stacking horizontally and vertically.
c = np.vstack((np.hstack((c11, c12)), np.hstack((c21, c22))))
return c
Я нашел код выше для алгоритма.
#include <time.h>
int main(void) {
clock_t tStart = clock();
/* Do your stuff here */
printf("Time taken: %.2fs\n", (double)(clock() - tStart)/CLOCKS_PER_SEC);
return 0;
}
Я нашел этот код для измерения времени выполнения кода. Однако я увидел, что могу использовать
/usr/bin/time ./MyProgram if I have cygwin installed.
Короче говоря, как бы я использовал свой код для решения реальных матриц с использованием алгоритма Штрассена и других алгоритмов решения матриц? И как мне запустить код? Спасибо за помощь, я новичок в кодировании, и я делаю это, чтобы проверить алгоритмическую эффективность различных алгоритмов решения матриц в разных сценариях.
