Распараллеливание кусочка кода на Си (MPI && openMP)

[Quew]

Есть примерно такой (в упрощенном варианте) кусок кода.

void function(...)
{
for (i1 = 0; i1 < N1; i1++){
for (i2 = 0; i2 < N2; i2++){
for (i3 = 0; i3 < N3; i3++){
for (i4 = 0; i4 < N4; i4++){
for (i5 = 0; i5 < N5; i5++){

tmp = func(a[i1], b[i2], c[i3], d[i4], e[i5]);
if (tmp < min) min = tmp;

}}}}}

printf("%f\n", min);
}

где i1..i5 - integer, N1..N5 - разные числа в определенном (не упорядочены), a..e - массивы с числами типа float.
func возвращает одно число, которое считается через кучу трёхэтажных уравнений. Необходимо просчитать все возможные варианты и найти минимальное.
Т.е. если объяснить на пальцах, то есть огромное дерево, и надо найти минимальный лист.

Необходимо распараллелить данный кусочек используя MPI и openMP (2 различные реализации, а не комбо) под линуксом, под Си. Компилятор gcc.

Для запуска вычислений есть кластер из 8машин по 8ядер на каждой. Вот желательно нагрузить его как-нибудь посильнее, чтобы сократить время расчётов.

Буду рад любым предложениям касаемо алгоритмов и примеров написания.
Заранее спасибо.

П.С. На openmp в голову приходить лишь разбиение первого цикла на 8 ветвей.