Параллельные вычисления, С++

[Anastasia.K]

Уважаемые форумчане!
У меня вопрос по параллельным вычислениям.

Код писать не надо, достаточно объяснить алгоритм словами.

Задача:
дана последовательность комплексных чисел длины N, нужно найти максимальное по модулю число (или числа), используя P потоков. Нужно выписать все элементы с максимальным по модулю значением и их позиции в векторе.
Сложность должна быть O( N/P+lg(P) )

Пишу в С++.

Мое решение:
последовательность пусть будет vector<double>.
я разбиваю этот вектор на части и посылаю каждому потоку по одной из частей (точнее позиции элементов), каждый поток находит максимальный по модулю элемент (или элементы) в своей части и запоминает их позиции в векторе. Далее используя эти же потоки, сравниваем данные, полученные каждым потоком. Например, первый поток обрабатывает результаты своей работы и второго потока ( например, если его максимальное значение больше, то про результаты второго забываем, если они равны, то объединяем результаты...), третий работает со своими результатами и результатами четвертого потока и.т.д, затем первый обрабатывает свои результаты и бывшего третьего, ну и т.д.

Что возвращают потоки? я бы написала структуру:

struct data{
  vector<int > index;   // тут сохраняем позиции элементов с макс. по модулю значением
  double max_val;  // а это само значение модуля
};

Вопрос: можно ли это все написать лучше?

Как бы вы решили эту задачу? Можно ли как-то оптимизировать решение, т.е. как-то лучше хранить результаты потоков, может быть лучше использовать другую структуру данных?

И ещё, какие есть хорошие книги по данной теме, по параллельным вычислениям и где можно найти хорошие примеры решения задач?

Заранее спасибо!

P.S. Я знаю про Google, мне хотелось бы узнать ваше мнение