0

Здравствуйте! Есть функция, которая умножает матрицы, но с разделением на потоки. Вот код:

unsigned Mult()
{
    int id = (int)param;
    int start = n / NumberOfThreads * id;
    int end = n / NumberOfThreads * (id + 1); 
    cout << "id=" << id;
    cout << "start=" << start;
    cout << "end=" << end;

    for (int i = start; i < end; i++)
            {
                for (int j = 0; j < k; j++)
                {
                    matr_3[i][j] = 0;
                    for (int z = 0; z < m; z++)
                        matr_3[i][j] += matr_1[i][z] * matr_2[z][j];
                }
            }
    return 0;
}

Проблема есть в этом месте:

int start = n / NumberOfThreads * id;
int end = n / NumberOfThreads * (id + 1);

Т.е. где потокам объясняется, как делится. При одном потоке для всех размерностей матриц вроде все хорошо. Но когда количество потоков четное, а хотя бы одна размерность матриц нечетная и наоборот, возникает проблема (последние строки или столбцы считаются неправильно). На примере можно понять о чем я: пусть кол-во потоков = 2, а n = 3. Тогда по алгоритму для первого потока start = 0, end = 1,5. А int только целое, и он округляет до 1. И в конце алгоритма end второга потока будет равно 2, а нужно 3. Т.е. нужно каким-то образом оптимизировать алгоритм деления на потоки, чтобы и для четных, и для нечетных всё было верно, т.е. универсальный. Каким образом это можно сделать?

1
  • Кстати, вы уверены, что int id = (int)param; делает то, что надо? Вроде бы param у вас указатель, а вы его приводите к целому числу. По идее, должно быть int id = *((int*)(param));. Посмотрите в отладчике, точно ли в id у вас получаются числа типа 1,2,3, а не что-то явно большее.
    – ASten
    Commented 9 мар. 2015 в 13:22

2 ответа 2

0

"Лишние операции" последнему потоку можно скинуть, добавив проверку при определении конца интервала строк, с которым работает функция:

int end = n / NumberOfThreads * (id + 1); // ваша строка
// дополнительная проверка и увеличение интервала для последнего потока
if( (n-1) - end < n/NumberOfThreads ) { end = n-1; }
2
  • Спасибо. Немного исправил, и вроде правильно.
    – compl
    Commented 9 мар. 2015 в 10:49
  • Если вы планируете развивать дальше вашу программу, я бы советовал вам выносить логику по определению входных данных за пределы функции потока. То есть, рассчитывать интервалы начала и конца в вызывающем коде, и передавать эти границы интервала в функцию потока как входные данные. Сейчас вы передаёте логический номер потока (int id = (int)param;) и по нему считаете границы. Советую вам определять start и end из param: int start = ((std::pair<int,int>*)(param))->first; int end = ((std::pair<int,int>*)(param))->second;. param должен быть указателем на std::pair с началом и концом.
    – ASten
    Commented 9 мар. 2015 в 13:17
0

Попробуйте в последний поток кидать "лишние" операции, т.е. если потоков 2 - делите количество операций на 2 и лишние кидаете в последний. Тоже самое и с тремя потоками.

1
  • Спасибо. Но не совсем понимаю, как кидать эти лишние операции.
    – compl
    Commented 8 мар. 2015 в 21:48

Ваш ответ

Нажимая «Отправить ответ», вы соглашаетесь с условиями пользования и подтверждаете, что прочитали политику конфиденциальности.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.