Счетные циклы и OpenMP
Некоторые документы, связанные с OpenMP, заявляют, что для того, чтобы цикл обрабатывался OpenMP, он должен быть «счетным», обеспечивая различные определения для цикла, который является «счетным»:
• the number of iterations in the loop must be countable
with an integer and loop use a fixed increment.
• the loop count can be “determined” ( what does it mean “determined”? )
Действительно ли это требование OpenMP? Или это требование конкретной реализации компилятора OpenMP?
Кроме того, можно ли распараллелить следующий цикл - кажется, не подсчитываемый - OpenMP (обратите внимание, что вопрос в том, можно ли распараллелить код, а не в том, есть ли способ создать параллельный эквивалент этого кода)
for ( i = 0; i < cnt; )
{
x1 = 2.0 * x - 1.;
if ( x1 < 1.0 )
{
i = i + 3;
x = x*2.;
}
else // if ( x1 >= 1. )
{
i = i + 2;
x = x/2.;
}
}
Ответы 1
Стандарт OpenMP требует Каноническая петля
The canonical form allows the iteration count of all associated loops to be computed before executing the outermost loop. The computation is performed for each loop in an integer type.
Это точно не ограничивается целыми числами, вы можете использовать типы указателей в C и типы итераторов произвольного доступа в C++, которые логически преобразуются в целые числа.
Переменная цикла не должна изменяться внутри цикла for (за исключением выражения приращения). Поэтому ваш пример не является допустимым кодом OpenMP.
Дополнительные ограничения заключаются в том, что тест цикла должен быть сравнением отношения <, <=, > или >= - не != или чего-то более сложного. Приращение должно быть постоянным на протяжении всего цикла. Отсутствующее выражение приращения в вашем примере также делает его недействительным.
Вы можете прочитать об этом в разделе 2.6 Стандарт OpenMP.
Вот чего я не понимаю: что значит «вычислено»? (или иногда называемый: «решительный»).
Скажем, количество циклов похоже на (upper_limit - lower_limit) / increment. Вычисление означает, что вы знаете, что значение является определенным числом, например. 42 перед выполнением цикла.
Вы не можете распараллелить этот цикл как есть, потому что индекс зависит от результата предыдущей итерации, т.е. у вас есть зависимость. Этот счетчик относится только к тому, что
cntизвестен в тот момент, когда мы входим в цикл, аiувеличивается некоторым тривиальным образом.