Re: [PPL-devel] [GIT] ppl/ppl(floating_point): Added a skeleton for the new linearization function.

Il 24/07/2010 12:42, Fabio Bossi ha scritto:

Module: ppl/ppl Branch: floating_point Commit: 1c4c91c248d74dc5abb8c5b45507dfe9665ff6b5 URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=1c4c91c248d74...

Author: Fabio Bossi bossi@cs.unipr.it Date: Sat Jul 24 12:39:51 2010 +0200

Added a skeleton for the new linearization function.

Also added an idea for the linearization of the unary negation expression.

---

+template <typename Target, typename FP_Interval_Type> +bool +linearize(const Concrete_Expression<Target>& expr,

•      const Box<FP_Interval_Type>& int_store,
  

•      const std::map<dimension_type, Linear_Form<FP_Interval_Type>>& lf_store,
  

•      Linear_Form<FP_Interval_Type>& result) {
  

• typedef typename FP_Interval_Type::boundary_type analyzer_format;
• typedef Linear_Form<FP_Interval_Type> FP_Linear_Form;
• typedef Box<FP_Interval_Type> FP_Interval_Abstract_Store;
• typedef std::map<dimension_type, FP_Linear_Form> FP_Linear_Form_Abstract_Store;
• // Check that analyzer_format is a floating point type.
• PPL_COMPILE_TIME_CHECK(!std::numeric_limits<analyzer_format>::is_exact,

•  "linearize<Target, FP_Interval_Type>:"
  

•  " FP_Interval_Type is not the type of an interval with floating point boundaries.");
  

• // Check that we are dealing with an expression of floating point type.
• PPL_ASSERT(expr.is_floating_point());
• /*
• Floating_Point_Format analyzed_format = expr.floating_point_format();
• FP_Interval_Type absolute_error =

•               compute_absolute_error<FP_Interval_Type>(analyzed_format);
  

• */
• switch(expr.kind()) {
• case INT_CON:

case Integer_Constant<Target>::KIND:

• // TODO.
• break;
• case FP_CON:
• Floating_Point_Constant<Target> fpc_expr =

•  static_cast<Floating_Point_Constant<Target>>(expr);
  

• result = FP_Linear_Form(FP_Interval(fpc_expr.get_value_as_string()));
• return true;
• break;
• case UOP:
• Unary_Operator<Target> uop_expr =

•  static_cast<Unary_Operator<Target>>(expr);
  

• switch (uop_expr.get_uop()) {
• case PLUS:

case Unary_Operator<Target>::PLUS:

•  return linearize(uop_expr.get_arg(), int_store, lf_store, result);
  

•  break;
  

• case MINUS:

•  if (!linearize(uop_expr.get_arg(), int_store, lf_store, result))
  

•    return false;
  

•  result.negate();
  

•  return true;
  

•  break;
  

• case BNOT:

•  return bnot_linearize(uop_expr, int_store, lf_store, result);
  

•  break;
  

• default:

•  throw std::runtime_error("PPL internal error");
  

• }
• break;
• case BOP:
• Binary_Operator<Target> bop_expr =

•  static_cast<Binary_Operator<Target>>(expr);
  

• switch (bop_expr.get_bop()) {
• case ADD:

•  return add_linearize(bop_expr, int_store, lf_store, result);
  

•  break;
  

• case SUB:

•  return sub_linearize(bop_expr, int_store, lf_store, result);
  

•  break;
  

• case MUL:

•  return mul_linearize(bop_expr, int_store, lf_store, result);
  

•  break;
  

• case DIV:

•  return div_linearize(bop_expr, int_store, lf_store, result);
  

•  break;
  

• case REM:
• case BAND:
• case BOR:
• case BXOR:
• case LSHIFT:
• case RSHIFT:
• default:

•  throw std::runtime_error("PPL internal error");
  

• }
• break;
• case CAST:
• // TODO.
• break;
• default:
• throw std::runtime_error("PPL internal error");
• }

+}

+} // namespace Parma_Polyhedra_Library

+#endif // !defined(PPL_linearize_hh)

... e similarmente per tutti gli altri.