lu_spmd-1.c

Go to the documentation of this file.

if (ub >= lb)
  for (int c0 = 1; c0 <= ub; c0 += 1)
    for (int c1 = c0; c1 <= n; c1 += 1) {
      if (c0 >= lb && c1 >= c0 + 1) {
        s0(c0, c1);
        if (n >= ub + 1)
          s2(c0, c1);
      } else if (lb >= c0 + 1) {
        s3(c0, c1, lb, c0, c1);
      }
      for (int c3 = max(lb, c0); c3 <= ub; c3 += 1)
        s1(c0, c1, c3);
    }