При разработке параллельных алгоритмов решения сложных научно-технических задач принципиальным моментом является анализ эффективности использования параллелизма, состоящий обычно в оценке получаемого ускорения процесса вычислений (сокращения времени решения задачи). Формирование подобных оценок ускорения может осуществляться применительно к выбранному вычислительному алгоритму (оценка эффективности распараллеливания конкретного алгоритма). Другой важный подход состоит в построении оценок максимально возможного ускорения процесса решения задачи конкретного типа (оценка эффективности параллельного способа решения задачи).
В данной лекции рассматривается
Для описания существующих информационных зависимостей в
выбираемых алгоритмах решения задач может быть использована
модель в виде
Представим множество операций, выполняемых в исследуемом алгоритме решения вычислительной задачи, и существующие между операциями информационные зависимости в виде ациклического ориентированного графа
G = (V, R),
(рис 2.1) Пример вычислительной модели алгоритма в виде графа "операции – операнды"где V = {1,...,|V|} есть множество вершин графа,
представляющих выполняемые операции алгоритма, а R
есть множество дуг графа (при этом дуга r = (i, j)
принадлежит графу только в том случае, если операция j
использует результат выполнения операции i ). Для
примера на рис. 2.1 показан
граф алгоритма вычисления площади прямоугольника, заданного
координатами двух противолежащих углов. Как можно заметить по
приведенному примеру, для выполнения выбранного алгоритма решения
задачи могут быть использованы разные схемы вычислений и построены
соответственно разные вычислительные модели. Как будет показано
далее, разные схемы вычислений обладают различными возможностями
для распараллеливания и, тем самым, при построении модели
вычислений может быть поставлена задача выбора наиболее подходящей
для параллельного исполнения вычислительной схемы алгоритма.
В рассматриваемой вычислительной модели алгоритма вершины без
входных дуг могут использоваться для задания операций ввода, а
вершины без выходных дуг – для операций вывода. Обозначим через V множество вершин графа без вершин ввода, а через d(G) — диаметр (длину максимального пути) графа.
Операции алгоритма, между которыми нет пути в рамках выбранной схемы вычислений, могут быть выполнены параллельно (для вычислительной схемы на рис. 2.1, например, параллельно могут быть реализованы сначала все операции умножения, а затем первые две операции вычитания). Возможный способ описания параллельного выполнения алгоритма может состоять в следующем (см., например, [2, 22]).
Пусть p есть количество процессоров, используемых
для выполнения алгоритма. Тогда для параллельного выполнения
вычислений необходимо задать множество (
в котором для каждой операции $$i\in V$$ указывается
номер используемого для выполнения операции процессора Pi и время начала выполнения операции ti. Для того чтобы Hp:
Вычислительная схема алгоритма G совместно с Hp может рассматриваться как
модель параллельного алгоритма Ap(G,Hp),
исполняемого с использованием p процессоров. Время
выполнения параллельного алгоритма определяется максимальным
значением времени, применяемым в расписании$$T_p(G,H_p)=\max_{i \in V}(t_i + 1).$$
Для выбранной схемы вычислений желательно использование расписания, обеспечивающего минимальное время исполнения алгоритма$$T_p(G)=\min_{H_p} T_p(G,H_p).$$
Уменьшение времени выполнения может быть обеспечено и путем подбора наилучшей вычислительной схемы$$T_p = \min_{G} T_p(G).$$
Оценки Tp(G,Hp), Tp(G)
и Tp могут быть применены в качестве
показателей времени выполнения параллельного алгоритма. Кроме
того, для анализа максимально возможного параллелизма можно
определить оценку наиболее быстрого исполнения алгоритма$$T_{\infty} = \min_{p \ge 1} T_p.$$
Оценку $$T_{\infty }$$ можно рассматривать
как минимально возможное
Оценка T1 определяет
Приведем без доказательства теоретические положения, характеризующие свойства оценок времени выполнения параллельного алгоритма (см. [22]).
Теорема 1. Минимально возможное
Теорема 2. Пусть для некоторой вершины вывода в вычислительной схеме алгоритма существует путь из каждой вершины ввода. Кроме того, пусть входная степень вершин схемы (количество входящих дуг) не превышает 2. Тогда минимально возможное время выполнения параллельного алгоритма ограничено снизу значением
$$T_{\infty }(G)=log_{2}n,$$где n есть количество вершин ввода в схеме алгоритма.
Теорема 3. При уменьшении числа используемых процессоров
Теорема 4. Для любого количества используемых процессоров справедлива следующая верхняя оценка для времени выполнения параллельного алгоритма
$$\forall p\Rightarrow T_{p}<T_{\infty }+T_{1}/p.$$Теорема 5. Времени выполнения алгоритма, которое сопоставимо с минимально возможным временем $$T_{\infty }$$, можно достичь при количестве процессоров порядка $$p \sim T_{1}/T_{\infty }$$, а именно,
$$p\ge T_{1}/T_{\infty }\Rightarrow T_{p}\le 2T_{\infty }.$$При меньшем количестве процессоров
Приведенные утверждения позволяют дать следующие рекомендации по правилам формирования параллельных алгоритмов:
Для вывода рекомендаций по формированию расписания по параллельному выполнению алгоритма приведем доказательство теоремы 4.
Доказательство теоремы 4. Пусть $$H_{\infty }$$
есть p процессоров может быть построено
следующим образом. Выполнение алгоритма разделим на $$T_{\infty }$$ шагов; на каждом шаге $$\tau$$ следует выполнить все $$n_{\tau }$$
операций, которые выполнялись на итерации $$\tau$$ расписания $$H_{\infty }$$. Эти операции могут быть выполнены не
более чем за $$\lceil n_{\tau }/p\rceil$$ итераций
при использовании p процессоров. Как результат, время
выполнения алгоритма Tp может быть оценено
следующим образом$$T_p = \sum_{\tau=1}^{T_{\infty}}
\left{\lceil}
\frac{n_{tau}}{p}
\right{\rceil}
<
\sum_{\tau=1}^{T_{\infty}}
\left(
\frac{n_{tau}}{p} + 1
\right)
= \frac{T_1}{p} + T_{\infty} .$$
Доказательство теоремы дает практический способ построения
расписания параллельного алгоритма. Первоначально может быть
построено
p процессоров, по сравнению
с последовательным вариантом выполнения вычислений определяется
величиной
Sp(n)=T1(n)/Tp(n),
т.е. как отношение времени решения задач на скалярной ЭВМ к
времени выполнения параллельного алгоритма (величина n применяется для параметризации вычислительной
сложности решаемой задачи и может пониматься, например, как
количество входных данных задачи).
Ep(n)=T1(n)/(pTp(n))=Sp(n)/p
(величина эффективности определяет среднюю долю времени выполнения алгоритма, в течение которой процессоры реально задействованы для решения задачи).
Из приведенных соотношений можно показать, что в наилучшем
случае Sp(n)=p и Ep(n)=1.
При практическом применении данных показателей для оценки эффективности
параллельных вычислений следует учитывать два важных момента:
Sp(n)>p
- в этом случае говорят о существовании сверхлинейного
(superlinear) ускорения. Несмотря на парадоксальность таких
ситуаций ( Ep(n)=1 легко обеспечивается при использовании
одного процессора). Как результат, разработка методов параллельных
вычислений часто предполагает выбор некоторого компромиссного варианта
с учетом желаемых показателей ускорения и эффективности.При выборе надлежащего параллельного способа решения задачи может оказаться полезной оценка стоимости (cost) вычислений, определяемой как произведение времени параллельного решения задачи и числа используемых процессоров
Cp=pTp.
В связи с этим можно определить понятие стоимостно-оптимального
(cost-optimal) параллельного алгоритма как метода,
Далее для иллюстрации введенных понятий в следующем пункте будет рассмотрен учебный пример решения задачи вычисления частных сумм для последовательности числовых значений. Кроме того, данные показатели будут использоваться для характеристики эффективности всех рассматриваемых в лекциях 6 – 11 параллельных алгоритмов при решении ряда типовых задач вычислительной математики.
Рассмотрим для демонстрации ряда проблем, возникающих при
разработке параллельных методов вычислений, сравнительно простую задачу нахождения частных сумм
последовательности числовых значений$$S_k = \sum_{i=1}^k x_i, \quad 1 \le k \le n,$$
где n есть количество суммируемых значений (данная
задача известна также под названием prefix sum
problem ).
Изучение возможных параллельных методов решения данной задачи начнем с еще более простого варианта ее постановки – с задачи вычисления общей суммы имеющегося набора значений (в таком виде задача суммирования является частным случаем общей задачи редукции )$$S = \sum_{i=1}^n x_i.$$
Традиционный алгоритм для решения этой задачи состоит в последовательном суммировании элементов числового набора
S=0, S=S+x1,...
Вычислительная схема данного алгоритма может быть представлена следующим образом (см. рис. 2.2):
G1=(V1,R1),
где V1={v01,...,v0n, v11,...,v1n}
есть множество операций (вершины v01,...,v0n
обозначают операции ввода, каждая вершина v1i, 1<=i<=n, соответствует прибавлению
значения xi к накапливаемой сумме S ), а
R1={(v0i,v1i),(v1i,v1i+1), 1<=i<=n–1}
есть множество дуг, определяющих информационные зависимости операций.
(рис 2.2) Последовательная вычислительная схема алгоритма суммированияКак можно заметить, данный "стандартный" алгоритм суммирования допускает только строго последовательное исполнение и не может быть распараллелен.
Параллелизм алгоритма суммирования становится возможным только
при ином способе построения процесса вычислений, основанном на
использовании ассоциативности операции сложения. Получаемый новый
вариант суммирования (известный в литературе как
Данная вычислительная схема может быть определена как граф
(пусть n=2k )
G2(V2,R2),
(рис 2.3) Каскадная схема алгоритма суммированиягде V2={(vi1,...,vli), 0<=i<=k, 1<=li<=2-1n}
есть вершины графа ( (v01,...v0n)
- операции ввода, (v1l,...,v1n/2) -
операции суммирования первой итерации и т.д.), а множество дуг
графа определяется соотношениями:
R2={(vi-1,2j-1vij),(vi-1,2jvij), 1<=i<=k, 1<=j<=2-in}.
Как нетрудно оценить, количество итераций
k=log2n,
а общее количество операций суммирования
Kпосл=n/2+n/4+...+1=n–1
совпадает с количеством операций последовательного варианта
алгоритма суммирования. При параллельном исполнении отдельных
итераций
Kпар=log2n.
Поскольку считается, что T1=Kпосл, Tp=Kпар,
поэтому показатели ускорения и эффективности
Sp=T1/Tp=(n–1)/log2n, Ep=T1/pTp=(n–1)/(plog2n)=(n–1)/((n/2)log2n),
где p=n/2 есть необходимое для выполнения
Анализируя полученные характеристики, можно отметить, что
время параллельного выполнения
Получение асимптотически ненулевой эффективности может быть
обеспечено, например, при использовании модифицированной n=2k, k=2s. Тогда
в новом варианте
(n/log2n ) групп, в
каждой из которых содержится log2n
элементов; далее для каждой группы вычисляется сумма значений
при помощи последовательного алгоритма суммирования; вычисления
в каждой группе могут выполняться независимо друг от друга
(т.е. параллельно – для этого необходимо наличие не менее (n/log2n ) процессоров);(n/log2n)
сумм отдельных групп применяется обычная
(рис 2.4) Модифицированная каскадная схема суммированияТогда для выполнения первого этапа требуется log2n
параллельных операций при использовании p1=(n/log2n)
процессоров. Для выполнения второго этапа необходимо
log2(n/log2n)<=log2n
параллельных операций для p2=(n/log2n)/2
процессоров. Как результат, данный способ суммирования характеризуется
следующими показателями:
Tp=2log2n, p=(n/log2n).
С учетом полученных оценок показатели ускорения и эффективности
модифицированной
Sp=T1/Tp=(n–1)/2log2n, Ep=T1/pTp=(n–1)/(2(n/log2n)log2n)=(n–1)/2n.
Сравнивая данные оценки с показателями обычной
Можно отметить также, что данные значения показателей достигаются
при количестве процессоров, определенном в теореме 5. Кроме того,
необходимо подчеркнуть, что, в отличие от обычной
Cp=pTp=(n/log2n)(2log2n)
является пропорциональной времени выполнения последовательного алгоритма.
Вернемся к исходной задаче вычисления всех частных сумм последовательности значений и проведем анализ возможных способов последовательной и параллельной организации вычислений. Вычисление всех частных сумм на скалярном компьютере может быть получено при помощи обычного последовательного алгоритма суммирования при том же количестве операций (!)
T1=n.
При параллельном исполнении применение log2n параллельных операций (как и в случае
вычисления общей суммы), может состоять в следующем (см.
рис. 2.5, а также [22]):
S вектора
суммируемых значений ( S=x );i, 1<=i<=log2n,
формируется вспомогательный вектор Q путем сдвига вправо
вектора S на 2i-1 позиций
(освобождающиеся при сдвиге позиции слева устанавливаются в нулевые
значения); итерация алгоритма завершается параллельной операцией
суммирования векторов S и Q.
(рис 2.5) Схема параллельного алгоритма вычисления всех частных сумм(величины Si-j означают суммы значений
от i до j элементов числовой
последовательности)
Всего параллельный алгоритм выполняется за log2n параллельных операций сложения.
На каждой итерации алгоритма параллельно выполняются n
скалярных операций сложения и, таким образом, общее количество
скалярных операций определяется величиной
Kпар=nlog2n
(параллельный алгоритм содержит большее (!) количество операций
по сравнению с последовательным способом суммирования). Необходимое
количество процессоров определяется количеством суммируемых значений
( p=n ).
С учетом полученных соотношений показатели ускорения и эффективности параллельного алгоритма вычисления всех частных сумм оцениваются следующим образом:
Sp=T1/Tp=n/log2n, Ep=T1/pTp=n/(plog2n)=n/(nlog2n)=1/log2n.
Как следует из построенных оценок,
Оценка качества параллельных вычислений предполагает знание
наилучших (максимально достижимых) значений показателей ускорения
и эффективности, однако получение идеальных величин Sp=p для ускорения и Ep=1 для эффективности может быть
обеспечено не для всех вычислительно трудоемких задач. Так, для
рассматриваемого учебного примера в предыдущем пункте минимально
достижимое время параллельного вычисления суммы числовых значений
составляет log2n. Определенное содействие
в решении этой проблемы могут оказать теоретические утверждения,
приведенные в начале данной лекции. В дополнение к ним рассмотрим
еще ряд закономерностей, которые могут быть чрезвычайно полезны
при построении оценок максимально достижимого
1. f есть доля последовательных
вычислений в применяемом алгоритме обработки данных, тогда
в соответствии с p процессоров ограничивается величиной$$S_p \le \frac{1}{f + (1-f)/p} \le S^* = \frac{1}{f}.$$
Так, например, при наличии всего 10% последовательных команд в
выполняемых вычислениях эффект использования параллелизма не
может превышать 10-кратного ускорения обработки данных. В
рассмотренном учебном примере вычисления суммы значений для f=log2n/n и, как результат, величина
возможного ускорения ограничена оценкой S*=n/log2n.
Следует отметить также, что рассмотрение f
является постоянной величиной и не зависит от параметра n,
определяющего вычислительную сложность решаемой задачи. Однако для
большого ряда задач доля f=f(n) является убывающей
функцией от n, и в этом случае Sp=Sp(n) является возрастающей
функцией от параметра n (данное утверждение часто
именуется эффект Амдаля ). Так, например, для
учебного примера вычисления суммы значений при использовании
фиксированного числа процессоров p суммируемый набор
данных может быть разделен на блоки размера n/p, для
которых сначала параллельно могут быть вычислены частные суммы, а
далее эти суммы можно сложить при помощи
Tp=(n/p)+log2p,
что приводит к оценке доли последовательных расчетов как величины
f=(1/p)+log2p/n.
Как следует из полученного выражения, доля последовательных
расчетов f убывает с ростом n и в предельном
случае мы получаемом идеальную оценку максимально возможного ускорения S*=p.
2.
С учетом введенной величины g можно получить
что позволяет построить оценку для ускорения$$S_p = \frac{T_1}{T_p} = \frac{\tau (n) + \pi (n)}{\tau (n) + \pi (n) / p} =
\frac{(\tau (n) + \pi (n) / p)(g + (1-g)p)}{\tau (n) + \pi (n) / p} ,$$
которая после упрощения приводится к виду
Sp = g+(1–g)p = p+(1–p)g.
Применительно к учебному примеру суммирования значений при
использовании p процессоров время параллельного
выполнения, как уже отмечалось выше, составляет
Tp = (n/p)+log2p,
что соответствует последовательной доле$$g = \frac{\log_2 p}{(n/p) + \log_2 p} .$$
За счет увеличения числа суммируемых значений величина g может быть пренебрежимо малой, обеспечивая
получение идеального возможного ускорения Sp=p.
При рассмотрении
Целью применения параллельных вычислений во многих случаях
является не только уменьшение времени выполнения расчетов, но
и обеспечение возможности решения более сложных вариантов задач
(таких постановок, решение которых не представляется возможным
при использовании однопроцессорных вычислительных систем).
Способность параллельного алгоритма эффективно использовать
процессоры при повышении сложности вычислений является важной
характеристикой выполняемых расчетов. Поэтому параллельный
алгоритм называют масштабируемым
( scalable ), если при росте числа
процессоров он обеспечивает увеличение ускорения при сохранении
постоянного уровня эффективности использования процессоров.
Возможный способ характеристики свойств
Оценим накладные расходы ( total overhead ), которые имеют место при выполнении параллельного алгоритма
T0=pTp–T1.
Накладные расходы появляются за счет необходимости организации взаимодействия процессоров, выполнения некоторых дополнительных действий, синхронизации параллельных вычислений и т.п. Используя введенное обозначение, можно получить новые выражения для времени параллельного решения задачи и соответствующего ускорения:$$T_p = \frac{T_1 + T_0}{p}, \quad S_p = \frac{T_1}{T_p} = \frac{pT_1}{T_1 + T_0} .$$
Применяя полученные соотношения,
Последнее выражение показывает, что если сложность решаемой
задачи является фиксированной ( T1=const ),
то при росте числа процессоров T0.
При фиксации числа процессоров T1 (предполагается, что при росте
параметра сложности n накладные расходы T0
увеличиваются медленнее, чем объем вычислений T1 ).
Как результат, при увеличении числа процессоров в большинстве
случаев можно обеспечить определенный уровень эффективности при
помощи соответствующего повышения сложности решаемых задач.
Поэтому важной характеристикой параллельных вычислений становится
соотношение необходимых темпов роста сложности расчетов и числа
используемых процессоров.
Пусть E=const есть желаемый уровень эффективности
выполняемых вычислений. Из выражения для эффективности можно получить$$\frac{T_0}{T_1} = \frac{1-E}{E} \; \text{или } T_1 = KT_0, K=E/(1-E).$$
Порождаемую последним соотношением зависимость n=F(p)
между сложностью решаемой задачи и числом процессоров обычно
называют
Покажем в качестве иллюстрации вывод
Как результат, например, при числе процессоров p=16
для обеспечения уровня эффективности E=0,5 (т.е. K=1 ) количество суммируемых значений должно быть не
менее n=64. Или же, при увеличении числа процессоров с p до q (q>p) для обеспечения
пропорционального роста ускорения (Sq/Sp)=(q/p) необходимо увеличить
число суммируемых значений n в (qlog2q)/(plog2p) раз.
В лекции рассматривается
В основу данной модели положен ациклический ориентированный
граф, в котором вершины представляют операции, а дуги
соответствуют зависимостям операций по данным. При наличии
такого графа для определения параллельного алгоритма достаточно
задать
Представление вычислений при помощи моделей подобного вида
позволяет получить аналитически ряд характеристик
разрабатываемых параллельных алгоритмов, среди которых время
выполнения, схема оптимального расписания, оценки максимально
возможного быстродействия методов решения поставленных задач.
Для более простого построения теоретических оценок в лекции
рассматривается понятие
Для оценки оптимальности разрабатываемых методов параллельных
вычислений в лекции приводятся широко используемые в теории и
практике параллельного программирования основные показатели
качества -
Для демонстрации применимости рассмотренных моделей и методов
анализа параллельных алгоритмов в лекции рассматривается задача
нахождения частных сумм последовательности числовых значений. На
данном примере отмечается проблема сложности распараллеливания
последовательных алгоритмов, которые изначально не были
ориентированы на возможность организации параллельных вычислений.
Для выделения "скрытого" параллелизма показывается
возможность преобразования исходной последовательной схемы
вычислений и приводится получаемая в результате таких
преобразований
В завершение лекции рассматривается вопрос построения оценок
максимально достижимых значений показателей эффективности. Для
получения таких оценок может быть использован
Дополнительная информация по моделированию и анализу параллельных вычислений может быть получена, например, в [2, 22]), полезная информация содержится также в [51, 63].
Рассмотрение учебной задачи суммирования последовательности числовых значений было выполнено в [22].
Впервые
Систематическое изложение (на момент издания работы) вопросов моделирования и анализа параллельных вычислений приводится в [77].
1. Разработайте модель и выполните оценку показателей ускорения и эффективности параллельных вычислений:
2. Выполните в соответствии с
3. Выполните оценку ускорения масштабирования для задач п.1.
4. Выполните построение
5. Разработайте модель и выполните полный анализ эффективности
параллельных вычислений (
При разработке параллельных алгоритмов решения сложных научно-технических задач принципиальным моментом является анализ эффективности использования параллелизма, состоящий обычно в оценке получаемого ускорения процесса вычислений (сокращения времени решения задачи). Формирование подобных оценок ускорения может осуществляться применительно к выбранному вычислительному алгоритму (оценка эффективности распараллеливания конкретного алгоритма). Другой важный подход состоит в построении оценок максимально возможного ускорения процесса решения задачи конкретного типа (оценка эффективности параллельного способа решения задачи).
В данной лекции рассматривается
Для описания существующих информационных зависимостей в
выбираемых алгоритмах решения задач может быть использована
модель в виде
Представим множество операций, выполняемых в исследуемом алгоритме решения вычислительной задачи, и существующие между операциями информационные зависимости в виде ациклического ориентированного графа
G = (V, R),
(рис 2.1) Пример вычислительной модели алгоритма в виде графа "операции – операнды"где V = {1,...,|V|} есть множество вершин графа,
представляющих выполняемые операции алгоритма, а R
есть множество дуг графа (при этом дуга r = (i, j)
принадлежит графу только в том случае, если операция j
использует результат выполнения операции i ). Для
примера на рис. 2.1 показан
граф алгоритма вычисления площади прямоугольника, заданного
координатами двух противолежащих углов. Как можно заметить по
приведенному примеру, для выполнения выбранного алгоритма решения
задачи могут быть использованы разные схемы вычислений и построены
соответственно разные вычислительные модели. Как будет показано
далее, разные схемы вычислений обладают различными возможностями
для распараллеливания и, тем самым, при построении модели
вычислений может быть поставлена задача выбора наиболее подходящей
для параллельного исполнения вычислительной схемы алгоритма.
В рассматриваемой вычислительной модели алгоритма вершины без
входных дуг могут использоваться для задания операций ввода, а
вершины без выходных дуг – для операций вывода. Обозначим через V множество вершин графа без вершин ввода, а через d(G) — диаметр (длину максимального пути) графа.
Операции алгоритма, между которыми нет пути в рамках выбранной схемы вычислений, могут быть выполнены параллельно (для вычислительной схемы на рис. 2.1, например, параллельно могут быть реализованы сначала все операции умножения, а затем первые две операции вычитания). Возможный способ описания параллельного выполнения алгоритма может состоять в следующем (см., например, [2, 22]).
Пусть p есть количество процессоров, используемых
для выполнения алгоритма. Тогда для параллельного выполнения
вычислений необходимо задать множество (
в котором для каждой операции $$i\in V$$ указывается
номер используемого для выполнения операции процессора Pi и время начала выполнения операции ti. Для того чтобы Hp:
Вычислительная схема алгоритма G совместно с Hp может рассматриваться как
модель параллельного алгоритма Ap(G,Hp),
исполняемого с использованием p процессоров. Время
выполнения параллельного алгоритма определяется максимальным
значением времени, применяемым в расписании$$T_p(G,H_p)=\max_{i \in V}(t_i + 1).$$
Для выбранной схемы вычислений желательно использование расписания, обеспечивающего минимальное время исполнения алгоритма$$T_p(G)=\min_{H_p} T_p(G,H_p).$$
Уменьшение времени выполнения может быть обеспечено и путем подбора наилучшей вычислительной схемы$$T_p = \min_{G} T_p(G).$$
Оценки Tp(G,Hp), Tp(G)
и Tp могут быть применены в качестве
показателей времени выполнения параллельного алгоритма. Кроме
того, для анализа максимально возможного параллелизма можно
определить оценку наиболее быстрого исполнения алгоритма$$T_{\infty} = \min_{p \ge 1} T_p.$$
Оценку $$T_{\infty }$$ можно рассматривать
как минимально возможное
Оценка T1 определяет
Приведем без доказательства теоретические положения, характеризующие свойства оценок времени выполнения параллельного алгоритма (см. [22]).
Теорема 1. Минимально возможное
Теорема 2. Пусть для некоторой вершины вывода в вычислительной схеме алгоритма существует путь из каждой вершины ввода. Кроме того, пусть входная степень вершин схемы (количество входящих дуг) не превышает 2. Тогда минимально возможное время выполнения параллельного алгоритма ограничено снизу значением
$$T_{\infty }(G)=log_{2}n,$$где n есть количество вершин ввода в схеме алгоритма.
Теорема 3. При уменьшении числа используемых процессоров
Теорема 4. Для любого количества используемых процессоров справедлива следующая верхняя оценка для времени выполнения параллельного алгоритма
$$\forall p\Rightarrow T_{p}<T_{\infty }+T_{1}/p.$$Теорема 5. Времени выполнения алгоритма, которое сопоставимо с минимально возможным временем $$T_{\infty }$$, можно достичь при количестве процессоров порядка $$p \sim T_{1}/T_{\infty }$$, а именно,
$$p\ge T_{1}/T_{\infty }\Rightarrow T_{p}\le 2T_{\infty }.$$При меньшем количестве процессоров
Приведенные утверждения позволяют дать следующие рекомендации по правилам формирования параллельных алгоритмов:
Для вывода рекомендаций по формированию расписания по параллельному выполнению алгоритма приведем доказательство теоремы 4.
Доказательство теоремы 4. Пусть $$H_{\infty }$$
есть p процессоров может быть построено
следующим образом. Выполнение алгоритма разделим на $$T_{\infty }$$ шагов; на каждом шаге $$\tau$$ следует выполнить все $$n_{\tau }$$
операций, которые выполнялись на итерации $$\tau$$ расписания $$H_{\infty }$$. Эти операции могут быть выполнены не
более чем за $$\lceil n_{\tau }/p\rceil$$ итераций
при использовании p процессоров. Как результат, время
выполнения алгоритма Tp может быть оценено
следующим образом$$T_p = \sum_{\tau=1}^{T_{\infty}}
\left{\lceil}
\frac{n_{tau}}{p}
\right{\rceil}
<
\sum_{\tau=1}^{T_{\infty}}
\left(
\frac{n_{tau}}{p} + 1
\right)
= \frac{T_1}{p} + T_{\infty} .$$
Доказательство теоремы дает практический способ построения
расписания параллельного алгоритма. Первоначально может быть
построено
p процессоров, по сравнению
с последовательным вариантом выполнения вычислений определяется
величиной
Sp(n)=T1(n)/Tp(n),
т.е. как отношение времени решения задач на скалярной ЭВМ к
времени выполнения параллельного алгоритма (величина n применяется для параметризации вычислительной
сложности решаемой задачи и может пониматься, например, как
количество входных данных задачи).
Ep(n)=T1(n)/(pTp(n))=Sp(n)/p
(величина эффективности определяет среднюю долю времени выполнения алгоритма, в течение которой процессоры реально задействованы для решения задачи).
Из приведенных соотношений можно показать, что в наилучшем
случае Sp(n)=p и Ep(n)=1.
При практическом применении данных показателей для оценки эффективности
параллельных вычислений следует учитывать два важных момента:
Sp(n)>p
- в этом случае говорят о существовании сверхлинейного
(superlinear) ускорения. Несмотря на парадоксальность таких
ситуаций ( Ep(n)=1 легко обеспечивается при использовании
одного процессора). Как результат, разработка методов параллельных
вычислений часто предполагает выбор некоторого компромиссного варианта
с учетом желаемых показателей ускорения и эффективности.При выборе надлежащего параллельного способа решения задачи может оказаться полезной оценка стоимости (cost) вычислений, определяемой как произведение времени параллельного решения задачи и числа используемых процессоров
Cp=pTp.
В связи с этим можно определить понятие стоимостно-оптимального
(cost-optimal) параллельного алгоритма как метода,
Далее для иллюстрации введенных понятий в следующем пункте будет рассмотрен учебный пример решения задачи вычисления частных сумм для последовательности числовых значений. Кроме того, данные показатели будут использоваться для характеристики эффективности всех рассматриваемых в лекциях 6 – 11 параллельных алгоритмов при решении ряда типовых задач вычислительной математики.
Рассмотрим для демонстрации ряда проблем, возникающих при
разработке параллельных методов вычислений, сравнительно простую задачу нахождения частных сумм
последовательности числовых значений$$S_k = \sum_{i=1}^k x_i, \quad 1 \le k \le n,$$
где n есть количество суммируемых значений (данная
задача известна также под названием prefix sum
problem ).
Изучение возможных параллельных методов решения данной задачи начнем с еще более простого варианта ее постановки – с задачи вычисления общей суммы имеющегося набора значений (в таком виде задача суммирования является частным случаем общей задачи редукции )$$S = \sum_{i=1}^n x_i.$$
Традиционный алгоритм для решения этой задачи состоит в последовательном суммировании элементов числового набора
S=0, S=S+x1,...
Вычислительная схема данного алгоритма может быть представлена следующим образом (см. рис. 2.2):
G1=(V1,R1),
где V1={v01,...,v0n, v11,...,v1n}
есть множество операций (вершины v01,...,v0n
обозначают операции ввода, каждая вершина v1i, 1<=i<=n, соответствует прибавлению
значения xi к накапливаемой сумме S ), а
R1={(v0i,v1i),(v1i,v1i+1), 1<=i<=n–1}
есть множество дуг, определяющих информационные зависимости операций.
(рис 2.2) Последовательная вычислительная схема алгоритма суммированияКак можно заметить, данный "стандартный" алгоритм суммирования допускает только строго последовательное исполнение и не может быть распараллелен.
Параллелизм алгоритма суммирования становится возможным только
при ином способе построения процесса вычислений, основанном на
использовании ассоциативности операции сложения. Получаемый новый
вариант суммирования (известный в литературе как
Данная вычислительная схема может быть определена как граф
(пусть n=2k )
G2(V2,R2),
(рис 2.3) Каскадная схема алгоритма суммированиягде V2={(vi1,...,vli), 0<=i<=k, 1<=li<=2-1n}
есть вершины графа ( (v01,...v0n)
- операции ввода, (v1l,...,v1n/2) -
операции суммирования первой итерации и т.д.), а множество дуг
графа определяется соотношениями:
R2={(vi-1,2j-1vij),(vi-1,2jvij), 1<=i<=k, 1<=j<=2-in}.
Как нетрудно оценить, количество итераций
k=log2n,
а общее количество операций суммирования
Kпосл=n/2+n/4+...+1=n–1
совпадает с количеством операций последовательного варианта
алгоритма суммирования. При параллельном исполнении отдельных
итераций
Kпар=log2n.
Поскольку считается, что T1=Kпосл, Tp=Kпар,
поэтому показатели ускорения и эффективности
Sp=T1/Tp=(n–1)/log2n, Ep=T1/pTp=(n–1)/(plog2n)=(n–1)/((n/2)log2n),
где p=n/2 есть необходимое для выполнения
Анализируя полученные характеристики, можно отметить, что
время параллельного выполнения
Получение асимптотически ненулевой эффективности может быть
обеспечено, например, при использовании модифицированной n=2k, k=2s. Тогда
в новом варианте
(n/log2n ) групп, в
каждой из которых содержится log2n
элементов; далее для каждой группы вычисляется сумма значений
при помощи последовательного алгоритма суммирования; вычисления
в каждой группе могут выполняться независимо друг от друга
(т.е. параллельно – для этого необходимо наличие не менее (n/log2n ) процессоров);(n/log2n)
сумм отдельных групп применяется обычная
(рис 2.4) Модифицированная каскадная схема суммированияТогда для выполнения первого этапа требуется log2n
параллельных операций при использовании p1=(n/log2n)
процессоров. Для выполнения второго этапа необходимо
log2(n/log2n)<=log2n
параллельных операций для p2=(n/log2n)/2
процессоров. Как результат, данный способ суммирования характеризуется
следующими показателями:
Tp=2log2n, p=(n/log2n).
С учетом полученных оценок показатели ускорения и эффективности
модифицированной
Sp=T1/Tp=(n–1)/2log2n, Ep=T1/pTp=(n–1)/(2(n/log2n)log2n)=(n–1)/2n.
Сравнивая данные оценки с показателями обычной
Можно отметить также, что данные значения показателей достигаются
при количестве процессоров, определенном в теореме 5. Кроме того,
необходимо подчеркнуть, что, в отличие от обычной
Cp=pTp=(n/log2n)(2log2n)
является пропорциональной времени выполнения последовательного алгоритма.
Вернемся к исходной задаче вычисления всех частных сумм последовательности значений и проведем анализ возможных способов последовательной и параллельной организации вычислений. Вычисление всех частных сумм на скалярном компьютере может быть получено при помощи обычного последовательного алгоритма суммирования при том же количестве операций (!)
T1=n.
При параллельном исполнении применение log2n параллельных операций (как и в случае
вычисления общей суммы), может состоять в следующем (см.
рис. 2.5, а также [22]):
S вектора
суммируемых значений ( S=x );i, 1<=i<=log2n,
формируется вспомогательный вектор Q путем сдвига вправо
вектора S на 2i-1 позиций
(освобождающиеся при сдвиге позиции слева устанавливаются в нулевые
значения); итерация алгоритма завершается параллельной операцией
суммирования векторов S и Q.
(рис 2.5) Схема параллельного алгоритма вычисления всех частных сумм(величины Si-j означают суммы значений
от i до j элементов числовой
последовательности)
Всего параллельный алгоритм выполняется за log2n параллельных операций сложения.
На каждой итерации алгоритма параллельно выполняются n
скалярных операций сложения и, таким образом, общее количество
скалярных операций определяется величиной
Kпар=nlog2n
(параллельный алгоритм содержит большее (!) количество операций
по сравнению с последовательным способом суммирования). Необходимое
количество процессоров определяется количеством суммируемых значений
( p=n ).
С учетом полученных соотношений показатели ускорения и эффективности параллельного алгоритма вычисления всех частных сумм оцениваются следующим образом:
Sp=T1/Tp=n/log2n, Ep=T1/pTp=n/(plog2n)=n/(nlog2n)=1/log2n.
Как следует из построенных оценок,
Оценка качества параллельных вычислений предполагает знание
наилучших (максимально достижимых) значений показателей ускорения
и эффективности, однако получение идеальных величин Sp=p для ускорения и Ep=1 для эффективности может быть
обеспечено не для всех вычислительно трудоемких задач. Так, для
рассматриваемого учебного примера в предыдущем пункте минимально
достижимое время параллельного вычисления суммы числовых значений
составляет log2n. Определенное содействие
в решении этой проблемы могут оказать теоретические утверждения,
приведенные в начале данной лекции. В дополнение к ним рассмотрим
еще ряд закономерностей, которые могут быть чрезвычайно полезны
при построении оценок максимально достижимого
1. f есть доля последовательных
вычислений в применяемом алгоритме обработки данных, тогда
в соответствии с p процессоров ограничивается величиной$$S_p \le \frac{1}{f + (1-f)/p} \le S^* = \frac{1}{f}.$$
Так, например, при наличии всего 10% последовательных команд в
выполняемых вычислениях эффект использования параллелизма не
может превышать 10-кратного ускорения обработки данных. В
рассмотренном учебном примере вычисления суммы значений для f=log2n/n и, как результат, величина
возможного ускорения ограничена оценкой S*=n/log2n.
Следует отметить также, что рассмотрение f
является постоянной величиной и не зависит от параметра n,
определяющего вычислительную сложность решаемой задачи. Однако для
большого ряда задач доля f=f(n) является убывающей
функцией от n, и в этом случае Sp=Sp(n) является возрастающей
функцией от параметра n (данное утверждение часто
именуется эффект Амдаля ). Так, например, для
учебного примера вычисления суммы значений при использовании
фиксированного числа процессоров p суммируемый набор
данных может быть разделен на блоки размера n/p, для
которых сначала параллельно могут быть вычислены частные суммы, а
далее эти суммы можно сложить при помощи
Tp=(n/p)+log2p,
что приводит к оценке доли последовательных расчетов как величины
f=(1/p)+log2p/n.
Как следует из полученного выражения, доля последовательных
расчетов f убывает с ростом n и в предельном
случае мы получаемом идеальную оценку максимально возможного ускорения S*=p.
2.
С учетом введенной величины g можно получить
что позволяет построить оценку для ускорения$$S_p = \frac{T_1}{T_p} = \frac{\tau (n) + \pi (n)}{\tau (n) + \pi (n) / p} =
\frac{(\tau (n) + \pi (n) / p)(g + (1-g)p)}{\tau (n) + \pi (n) / p} ,$$
которая после упрощения приводится к виду
Sp = g+(1–g)p = p+(1–p)g.
Применительно к учебному примеру суммирования значений при
использовании p процессоров время параллельного
выполнения, как уже отмечалось выше, составляет
Tp = (n/p)+log2p,
что соответствует последовательной доле$$g = \frac{\log_2 p}{(n/p) + \log_2 p} .$$
За счет увеличения числа суммируемых значений величина g может быть пренебрежимо малой, обеспечивая
получение идеального возможного ускорения Sp=p.
При рассмотрении
Целью применения параллельных вычислений во многих случаях
является не только уменьшение времени выполнения расчетов, но
и обеспечение возможности решения более сложных вариантов задач
(таких постановок, решение которых не представляется возможным
при использовании однопроцессорных вычислительных систем).
Способность параллельного алгоритма эффективно использовать
процессоры при повышении сложности вычислений является важной
характеристикой выполняемых расчетов. Поэтому параллельный
алгоритм называют масштабируемым
( scalable ), если при росте числа
процессоров он обеспечивает увеличение ускорения при сохранении
постоянного уровня эффективности использования процессоров.
Возможный способ характеристики свойств
Оценим накладные расходы ( total overhead ), которые имеют место при выполнении параллельного алгоритма
T0=pTp–T1.
Накладные расходы появляются за счет необходимости организации взаимодействия процессоров, выполнения некоторых дополнительных действий, синхронизации параллельных вычислений и т.п. Используя введенное обозначение, можно получить новые выражения для времени параллельного решения задачи и соответствующего ускорения:$$T_p = \frac{T_1 + T_0}{p}, \quad S_p = \frac{T_1}{T_p} = \frac{pT_1}{T_1 + T_0} .$$
Применяя полученные соотношения,
Последнее выражение показывает, что если сложность решаемой
задачи является фиксированной ( T1=const ),
то при росте числа процессоров T0.
При фиксации числа процессоров T1 (предполагается, что при росте
параметра сложности n накладные расходы T0
увеличиваются медленнее, чем объем вычислений T1 ).
Как результат, при увеличении числа процессоров в большинстве
случаев можно обеспечить определенный уровень эффективности при
помощи соответствующего повышения сложности решаемых задач.
Поэтому важной характеристикой параллельных вычислений становится
соотношение необходимых темпов роста сложности расчетов и числа
используемых процессоров.
Пусть E=const есть желаемый уровень эффективности
выполняемых вычислений. Из выражения для эффективности можно получить$$\frac{T_0}{T_1} = \frac{1-E}{E} \; \text{или } T_1 = KT_0, K=E/(1-E).$$
Порождаемую последним соотношением зависимость n=F(p)
между сложностью решаемой задачи и числом процессоров обычно
называют
Покажем в качестве иллюстрации вывод
Как результат, например, при числе процессоров p=16
для обеспечения уровня эффективности E=0,5 (т.е. K=1 ) количество суммируемых значений должно быть не
менее n=64. Или же, при увеличении числа процессоров с p до q (q>p) для обеспечения
пропорционального роста ускорения (Sq/Sp)=(q/p) необходимо увеличить
число суммируемых значений n в (qlog2q)/(plog2p) раз.
В лекции рассматривается
В основу данной модели положен ациклический ориентированный
граф, в котором вершины представляют операции, а дуги
соответствуют зависимостям операций по данным. При наличии
такого графа для определения параллельного алгоритма достаточно
задать
Представление вычислений при помощи моделей подобного вида
позволяет получить аналитически ряд характеристик
разрабатываемых параллельных алгоритмов, среди которых время
выполнения, схема оптимального расписания, оценки максимально
возможного быстродействия методов решения поставленных задач.
Для более простого построения теоретических оценок в лекции
рассматривается понятие
Для оценки оптимальности разрабатываемых методов параллельных
вычислений в лекции приводятся широко используемые в теории и
практике параллельного программирования основные показатели
качества -
Для демонстрации применимости рассмотренных моделей и методов
анализа параллельных алгоритмов в лекции рассматривается задача
нахождения частных сумм последовательности числовых значений. На
данном примере отмечается проблема сложности распараллеливания
последовательных алгоритмов, которые изначально не были
ориентированы на возможность организации параллельных вычислений.
Для выделения "скрытого" параллелизма показывается
возможность преобразования исходной последовательной схемы
вычислений и приводится получаемая в результате таких
преобразований
В завершение лекции рассматривается вопрос построения оценок
максимально достижимых значений показателей эффективности. Для
получения таких оценок может быть использован
Дополнительная информация по моделированию и анализу параллельных вычислений может быть получена, например, в [2, 22]), полезная информация содержится также в [51, 63].
Рассмотрение учебной задачи суммирования последовательности числовых значений было выполнено в [22].
Впервые
Систематическое изложение (на момент издания работы) вопросов моделирования и анализа параллельных вычислений приводится в [77].
1. Разработайте модель и выполните оценку показателей ускорения и эффективности параллельных вычислений:
2. Выполните в соответствии с
3. Выполните оценку ускорения масштабирования для задач п.1.
4. Выполните построение
5. Разработайте модель и выполните полный анализ эффективности
параллельных вычислений (
Для получения официальных документов о завершении программы дополнительного профессионального образования (удостоверения о повышении квалификации, дипломов о профессиональной переподготовке и MBA) необходимо предоставить:
Внимание! Вы можете не заказывать доставку бумажной версии официального документы, а скачать его в электронном виде и распечатать самостоятельно. Информация о выданном документе в течение 1 месяца загружается в Федеральную информационную систему «Федеральный реестр сведений о документах об образовании и (или) о квалификации, документах об обучении» - ФИС ФРДО.
Доступ на новый сайт осуществляется с использованием адреса электронной почты, который был указан вами при регистрации на "старом". Мы постарались перенести все ваши данные с прежнего ресурса, однако не исключена вероятность потери части информации.
При возникновении проблемы со входом, воспользуйтесь функцией сброса пароля
Если вы обнаружите несоответствия, пожалуйста, сообщите нам.