summaryrefslogtreecommitdiff
path: root/omegalib/examples/old_test/lu_ijk.oc-rt
diff options
context:
space:
mode:
Diffstat (limited to 'omegalib/examples/old_test/lu_ijk.oc-rt')
-rw-r--r--omegalib/examples/old_test/lu_ijk.oc-rt81
1 files changed, 0 insertions, 81 deletions
diff --git a/omegalib/examples/old_test/lu_ijk.oc-rt b/omegalib/examples/old_test/lu_ijk.oc-rt
deleted file mode 100644
index d067dff..0000000
--- a/omegalib/examples/old_test/lu_ijk.oc-rt
+++ /dev/null
@@ -1,81 +0,0 @@
->>> 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);
- }
-}
-
-