VPF::Алгоритм циклической редукции

peretc

Дата 14.1.2012, 19:05 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 1
Регистрация: 14.1.2012

Репутация: нет
Всего: нет

Доброго времени суток!

Необходимо реализовать решение системы линейных уравнений методом циклической редукции, причем нужны последовательная и параллельная реализации.

С последовательным вариантом я разобрался, а вот с параллельным никак не догоняю. Информации на этот счет достаточно мало, удается найти только то, что отсутствует фаза замещения.

Если кто-то может мне пояснить, как реализовать параллельный вариант или как модифицировать последовательный - буду очень благодарен.

PS Последовательный вариант реализую так:

Код


int tridag_cpu(double *a, double *b, double *c, double *r, double *x, int n, int stride) {
// 1. прямой ход
for(int nn = n, low = 2; nn > 1; nn /= 2, low *= 2, stride *= 2) {
for(int i = low - 1; i < n; i += stride * 2) {
    float alpha = -a[i] / b[i - stride];
    if(alpha!=alpha) { return 1;  }
    float gamma = -c[i] / b[i + stride];
    if(gamma!=gamma) { return 1; }
    a[i] = alpha * a[i - stride];
    b[i] = alpha * c[i - stride] + b[i] + gamma * a[i + stride];
    c[i] = gamma * c[i + stride];
    r[i] = alpha * r[i - stride] + r[i] + gamma * r[i + stride];

 } 
 } 

 // 2. обратный ход
 x[n/2] = r[n/2] / b[n/2];
 for(stride /= 2; stride >= 1; stride /= 2) {
 for(int i = stride - 1; i < n; i += stride * 2) {
 x[i] = (r[i]
 - (i - stride > 0 ? a[i] * x[i - stride] : 0.0)
 - (i + stride < n ? c[i] * x[i + stride] : 0.0)
 ) / b[i];
 }
 }
 return 0;
 }

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Алгоритмы \| Следующая тема »