DAS_2023_1/antonov_dmitry_lab_6/parallel_determ.py

import random
from multiprocessing import Pool
import time


def submatrix(matrix, row, col):
    return [[matrix[i][j] for j in range(len(matrix[i])) if j != col] for i in range(len(matrix)) if i != row]


def determinant(matrix):
    size = len(matrix)

    # Простой случай: детерминант матрицы 1x1
    if size == 1:
        return matrix[0][0]

    # Простой случай: детерминант матрицы 2x2
    if size == 2:
        return matrix[0][0] * matrix[1][1] - matrix[0][1] * matrix[1][0]

    det = 0
    for col in range(size):
        det += ((-1) ** col) * matrix[0][col] * determinant(submatrix(matrix, 0, col))

    return det


def generate_random_matrix(size, lower_limit, upper_limit):
    return [[random.uniform(lower_limit, upper_limit) for _ in range(size)] for _ in range(size)]


def sequential_determinant_calculation(matrix_size, lower_limit, upper_limit):
    random_matrix = generate_random_matrix(matrix_size, lower_limit, upper_limit)

    start_time = time.time()
    result = determinant(random_matrix)
    end_time = time.time()

    print(f"Последовательный детерминант: {result}")
    print(f"Последовательное время: {end_time - start_time} секунд")


def parallel_determinant_calculation(matrix_size, lower_limit, upper_limit, num_processes):
    random_matrix = generate_random_matrix(matrix_size, lower_limit, upper_limit)

    matrices_to_process = [submatrix(random_matrix, 0, col) for col in range(matrix_size)]

    start_time = time.time()
    with Pool(processes=num_processes) as pool:
        determinants = pool.map(determinant, matrices_to_process)

    result = sum(((-1) ** col) * random_matrix[0][col] * det for col, det in enumerate(determinants))
    end_time = time.time()

    print(f"Параллельный детерминант: {result}")
    print(f"Параллельное время: {end_time - start_time} секунд")


if __name__ == "__main__":
    matrix_size = 20  # размер матрицы
    lower_limit = 10  # числа в матрице от
    upper_limit = 20  # и до
    processes = [1, 2, 4, 8, 16, 32]
    #
    # # последовательное вычисление
    # sequential_determinant_calculation(matrix_size, lower_limit, upper_limit)

    # параллельное вычисление
    for i in processes:
        print("Потоков " + str(i))
        parallel_determinant_calculation(matrix_size, lower_limit, upper_limit, i)