[GIT] ppl/ppl(floating_point): Added some tests.

16 Sep 2009

Module: ppl/ppl
Branch: floating_point
Commit: bd64da813c519c335d4aa434c3d275e766fda20d
URL:    http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=bd64da813c519...
Author: Roberto Amadini r.amadini@virgilio.it
Date:   Wed Sep 16 11:25:27 2009 +0200
Added some tests.
---
tests/Floating_Point_Expression/affineimage3.cc |  209 +++++++++++++++++++++--
 1 files changed, 196 insertions(+), 13 deletions(-)

diff --git a/tests/Floating_Point_Expression/affineimage3.cc b/tests/Floating_Point_Expression/affineimage3.cc
index d1ab707..a2bdc0f 100644
--- a/tests/Floating_Point_Expression/affineimage3.cc
+++ b/tests/Floating_Point_Expression/affineimage3.cc
@@ -24,17 +24,51 @@ site: http://www.cs.unipr.it/ppl/ . */
namespace {
-// tests affine_image(A, [-2, 1])
+// tests space_dimensions and trivial cases
 bool
 test01() {
   Variable A(0);
   Variable B(1);
+  Octagonal_Shape<double> oc1(0);
+  bool ok1 = false;
+  Linear_Form<db_r_oc> l(A);
+
+  try {
+      oc1.affine_image(A, l);
+  }
+  catch(std::invalid_argument e) {
+    nout << "oc1_space_dim < lf_space_dim" << endl;
+    ok1 = true;
+  }
+
+  bool ok2 = false;
+  Octagonal_Shape<double> oc2(1);
+
+  try {
+    oc2.affine_image(B, l);
+  }
+  catch(std::invalid_argument e) {
+    nout << "space_dim < var_id + 1" << endl;
+    oc2.affine_image(A, l);
+    Constraint_System cs(A < A);
+    oc2.add_constraints(cs);
+    oc2.affine_image(A, l);
+    ok2 = true;
+  }
+
+  return ok1 && ok2;
+}
+
+// tests affine_image(A, [-2, 1])
+bool
+test02() {
+  Variable A(0);
+  Variable B(1);
Octagonal_Shape<float> oc1(3);
   oc1.add_constraint(A <= 2);
   oc1.add_constraint(A - B <= 3);
   oc1.add_constraint(B <= 2);
-  print_constraints(oc1, "*** oc1  ***");
   fl_r_oc free_term(-2);
   free_term.join_assign(1);
   Linear_Form<fl_r_oc> l(free_term);
@@ -53,8 +87,36 @@ test01() {
 }
-// tests affine_image(A, [-1, -1]*A)
-bool test02() {
+// tests affine_image(A, [-1, -1]*A + [0.5, 2])
+bool test03() {
+  Variable A(0);
+  Variable B(1);
+
+  Octagonal_Shape<double> oc1(3);
+  oc1.add_constraint(A <= 2);
+  oc1.add_constraint(A - B <= 3);
+  oc1.add_constraint(B <= 2);
+  db_r_oc free_term(0.5);
+  free_term.join_assign(2);
+  Linear_Form<db_r_oc> l(-A);
+  l += free_term;
+  oc1.affine_image(A, l);
+  print_constraints(oc1, "*** oc1.affine_image(A, -A + [0.5, 2]) ***");
+
+  Octagonal_Shape<double> known_result(3);
+  known_result.add_constraint(-2 * A <= 3);
+  known_result.add_constraint(B <= 2);
+  known_result.add_constraint(2*B - 2*A <= 7);
+  print_constraints(known_result, "*** known_result ***");
+  bool ok = (oc1 == known_result);
+
+  return ok;
+
+}
+
+
+// tests affine_image(B, [1, 1]*B + [-1, 3])
+bool test04() {
   Variable A(0);
   Variable B(1);
@@ -62,12 +124,102 @@ bool test02() {
   oc1.add_constraint(A <= 2);
   oc1.add_constraint(A - B <= 3);
   oc1.add_constraint(B <= 2);
-  Linear_Form<fl_r_oc> l(-A);
+  fl_r_oc free_term(-1);
+  free_term.join_assign(3);
+  Linear_Form<fl_r_oc> l(B);
+  l += free_term;
+  oc1.affine_image(B, l);
+  print_constraints(oc1, "*** oc1.affine_image(B, B + [-1, 3]) ***");
+  oc1.ascii_dump();
+  Octagonal_Shape<float> known_result(3);
+ /* known_result.add_constraint(A >= -2);
+  known_result.add_constraint(B <= 2);
+  known_result.add_constraint(B - A <= 4);
+  print_constraints(known_result, "*** known_result ***");
+
+  bool ok = (oc1 == known_result);
+
+  return ok;
+  */
+  return true;
+}
+
+// tests affine_image(A, [1, 1]*B + [-2, 0.5])
+bool test05() {
+  Variable A(0);
+  Variable B(1);
+
+  Octagonal_Shape<double> oc1(3);
+  oc1.add_constraint(A <= 2);
+  oc1.add_constraint(A - B <= 3);
+  oc1.add_constraint(B <= 2);
+  db_r_oc free_term(-2);
+  free_term.join_assign(0.5);
+  Linear_Form<db_r_oc> l(B);
+  l += free_term;
   oc1.affine_image(A, l);
-  print_constraints(oc1, "*** oc1.affine_image(A, [-1, -1]*A) ***");
+  print_constraints(oc1, "*** oc1.affine_image(A, B + [-2, 0.5]) ***");
+  oc1.ascii_dump();
+  Octagonal_Shape<double> known_result(3);
+ /* known_result.add_constraint(A >= -2);
+  known_result.add_constraint(B <= 2);
+  known_result.add_constraint(B - A <= 4);
+  print_constraints(known_result, "*** known_result ***");
+  bool ok = (oc1 == known_result);
+
+  return ok;
+  */
+  return true;
+}
+
+// tests affine_image(B, [1, 1]*A + [-3, 1])
+bool test06() {
+  Variable A(0);
+  Variable B(1);
+
+  Octagonal_Shape<double> oc1(3);
+  oc1.add_constraint(A <= 2);
+  oc1.add_constraint(A - B <= 3);
+  oc1.add_constraint(B <= 2);
+  db_r_oc free_term(-3);
+  free_term.join_assign(1);
+  Linear_Form<db_r_oc> l(A);
+  l += free_term;
+  oc1.affine_image(B, l);
+  print_constraints(oc1, "*** oc1.affine_image(B, A + [-3, 1]) ***");
+  oc1.ascii_dump();
+  Octagonal_Shape<double> known_result(3);
+ /* known_result.add_constraint(A >= -2);
+  known_result.add_constraint(B <= 2);
+  known_result.add_constraint(B - A <= 4);
+  print_constraints(known_result, "*** known_result ***");
+
+  bool ok = (oc1 == known_result);
+
+  return ok;
+  */
+  return true;
+}
+
+// tests affine_image(B, [-1, -1]*A + [0, 4])
+bool test07() {
+  Variable A(0);
+  Variable B(1);
+
+  Octagonal_Shape<float> oc1(3);
+  oc1.add_constraint(A <= 2);
+  oc1.add_constraint(A - B <= 3);
+  oc1.add_constraint(B <= 2);
+  fl_r_oc free_term(0);
+  free_term.join_assign(4);
+  Linear_Form<fl_r_oc> l(-A);
+  l += free_term;
+  oc1.affine_image(B, l);
+  print_constraints(oc1, "*** oc1.affine_image(B, -A + [0, 4]) ***");
+  oc1.ascii_dump();
   Octagonal_Shape<float> known_result(3);
-  known_result.add_constraint(A >= -2);
+ /* known_result.add_constraint(A >= -2);
   known_result.add_constraint(B <= 2);
   known_result.add_constraint(B - A <= 4);
   print_constraints(known_result, "*** known_result ***");
@@ -75,23 +227,54 @@ bool test02() {
   bool ok = (oc1 == known_result);
return ok;
+  */
+  return true;
 }
-/*
-// tests affine_image(A, [-1, 1]*B)
-bool test03() {
+// tests affine_image(A, [-1, -1]*B + [0, 2])
+bool test08() {
+  Variable A(0);
+  Variable B(1);
+
+  Octagonal_Shape<float> oc1(3);
+  oc1.add_constraint(A <= 2);
+  oc1.add_constraint(A - B <= 3);
+  oc1.add_constraint(B <= 2);
+  fl_r_oc free_term(0);
+  free_term.join_assign(2);
+  Linear_Form<fl_r_oc> l(-B);
+  l += free_term;
+  oc1.affine_image(A, l);
+  print_constraints(oc1, "*** oc1.affine_image(A, -B + [0, 2]) ***");
+  oc1.ascii_dump();
+  Octagonal_Shape<float> known_result(3);
+ /* known_result.add_constraint(A >= -2);
+  known_result.add_constraint(B <= 2);
+  known_result.add_constraint(B - A <= 4);
+  print_constraints(known_result, "*** known_result ***");
+
+  bool ok = (oc1 == known_result);
+
+  return ok;
+  */
   return true;
 }
-// tests affine_image(A, i + i1*A + i2*B)
-bool test04() {
+// tests affine_image(A, i + i1 * A + i2 * B)
+bool test09() {
   return true;
 }
-*/
} //namespace
BEGIN_MAIN
   DO_TEST(test01);
   DO_TEST(test02);
+  DO_TEST(test03); /* FIXME: Not yet tested!
+  DO_TEST(test04);
+  DO_TEST(test05);
+  DO_TEST(test06);
+  DO_TEST(test07);
+  DO_TEST(test08);
+  DO_TEST(test09); */
 END_MAIN

    

Roberto Amadini

tags (0)

participants (1)