Merge branch 'master' into doe

1 files changed, 81 insertions, 0 deletions

diff --git a/omega/examples/old_test/lu_ijk.oc-rt b/omega/examples/old_test/lu_ijk.oc-rt
new file mode 100644
index 0000000..d067dff
--- /dev/null
+++ b/omega/examples/old_test/lu_ijk.oc-rt
@@ -0,0 +1,81 @@
+>>> symbolic n,lb,ub;
+>>> I5 := {[k,j] : 1 <= k < j <= n};
+>>> I7 := {[k,j,i] : 1 <= k < j,i <= n};
+>>> T5 := {[k,j] -> [k,j,1,0]};
+>>> T7 := {[k,j,i] -> [i,j,0,k]};
+>>> D57 := {[k,j] -> [k,j,i] : 1 <= k < j,i <= n};
+>>> D75 := {[k,j,k+1] -> [k+1,j] : 1 <= k <= j-2 && j <= n} union
+>>> 	{[k,k+1,k+1] -> [k+1,j'] : 1 <= k <= j'-2 && j' <= n};
+>>> D77 := {[k,j,i] -> [k+1,j,i] : 1 <= k <= i-2, j-2 && j <= n && i <= n}
+>>> 	union
+>>> 	{[k,k+1,i] -> [k+1,j',i] : 1 <= k <= i-2, j'-2 && j' <= n && i <= n};
+>>> difference ( (inverse T5) join D57 join T7 );
+{[In_1,0,-1,In_4]: 1 <= In_4 && In_1+In_4 <= n && 1 <= In_1}
+>>> difference ( (inverse T7) join D75 join T5 );
+{[0,0,1,In_4]: -n+2 <= In_4 <= -1} union
+ {[0,In_2,1,In_4]: In_4 <= -1 && 1 <= In_2 && In_2 < n+In_4}
+>>> difference ( (inverse T7) join D77 join T7 );
+{[0,0,0,1]: 3 <= n} union
+ {[0,In_2,0,1]: 1 <= In_2 <= n-2}
+>>> codegen T5:I5,T7:I7;
+if (n >= 2) {
+  for(t1 = 1; t1 <= n; t1++) {
+    for(t2 = 2; t2 <= n; t2++) {
+      for(t4 = 1; t4 <= min(t2-1,t1-1); t4++) {
+        s2(t4,t2,t1);
+      }
+      if (t1 <= t2-1) {
+        s1(t1,t2);
+      }
+    }
+  }
+}
+
+>>> codegen 2 T5:I5,T7:I7;
+if (n >= 2) {
+  for(t1 = 1; t1 <= n; t1++) {
+    if (t1 <= 1) {
+      for(t2 = 2; t2 <= n; t2++) {
+        s1(1,t2);
+      }
+    }
+    for(t2 = 2; t2 <= t1; t2++) {
+      for(t4 = 1; t4 <= t2-1; t4++) {
+        s2(t4,t2,t1);
+      }
+    }
+    if (t1 >= 2) {
+      for(t2 = t1+1; t2 <= n; t2++) {
+        for(t4 = 1; t4 <= t1-1; t4++) {
+          s2(t4,t2,t1);
+        }
+        s1(t1,t2);
+      }
+    }
+  }
+}
+
+>>> codegen 3 T5:I5,T7:I7;
+for(t2 = 2; t2 <= n; t2++) {
+  s1(1,t2);
+}
+for(t1 = 2; t1 <= n-1; t1++) {
+  for(t2 = 2; t2 <= t1; t2++) {
+    for(t4 = 1; t4 <= t2-1; t4++) {
+      s2(t4,t2,t1);
+    }
+  }
+  for(t2 = t1+1; t2 <= n; t2++) {
+    for(t4 = 1; t4 <= t1-1; t4++) {
+      s2(t4,t2,t1);
+    }
+    s1(t1,t2);
+  }
+}
+for(t2 = 2; t2 <= n; t2++) {
+  for(t4 = 1; t4 <= t2-1; t4++) {
+    s2(t4,t2,n);
+  }
+}
+
+