source: to-imperative/trunk/library/Arithm/mult.cc @ 1010

// $Source$
// $Revision: 1010 $
// $Date: 2003-07-15 12:12:14 +0000 (Tue, 15 Jul 2003) $
// $Author: pooh $

#include <rf_core.hh>
#include <rf_short_int.hh>
#include <rf_integer.ih>
#include <inttypes.h>
#include <gmp.h>

namespace refal
{

using namespace rfrt;
using namespace rftype;

namespace Arithm
{

RF_FUNC (_a_, (RF_ARG s_Int1, s_Int2), (RF_RES s_Int))

  Term* p_Int1 = s_Int1.get_first();
  Term* p_Int2 = s_Int2.get_first();

  intptr_t int1;
  intptr_t int2;
  Expr result;

  if (p_Int1->get_type() == type_short_int) {
    int1 = ((ShortInt*)p_Int1)->to_int();
    if (p_Int2->get_type() == type_short_int) {
      int2 = ((ShortInt*)p_Int2)->to_int();
      if (imaxabs(int1) < INTPTR_MAX / imaxabs(int2)) {
        result = ShortInt::create_expr(int1 * int2);
      } else {
        result = Integer::create_expr(int1);
        Integer* res = static_cast<Integer*>(result.get_first());
        mpz_t* resp = res->get_mpz_ptr();
        mpz_mul_si(*resp, *resp, int2);
      }
    } else if (p_Int2->get_type() == type_int) {
      result = Integer::create_expr(int1);
      Integer* res = static_cast<Integer*>(result.get_first());
      mpz_t* resp = res->get_mpz_ptr();
      Integer* i2 = static_cast<Integer*>(p_Int2);
      mpz_t* i2p = i2->get_mpz_ptr();
      mpz_mul(*resp, *resp, *i2p);
    } else {
      RF_LIB_ERROR("Invalid argument");
    }
  } else if (p_Int1->get_type() == type_int) {
    Integer* i1 = static_cast<Integer*>(p_Int1);
    mpz_t* i1p = i1->get_mpz_ptr();
    if (p_Int2->get_type() == type_short_int) {
      int2 = ((ShortInt*)p_Int2)->to_int();
      result = Integer::create_expr(int2);
      Integer* res = static_cast<Integer*>(result.get_first());
      mpz_t* resp = res->get_mpz_ptr();
      mpz_mul(*resp, *resp, *i1p);
    } else if (p_Int2->get_type() == type_int) {
      result = Integer::create_expr(0);
      Integer* res = static_cast<Integer*>(result.get_first());
      mpz_t* resp = res->get_mpz_ptr();
      Integer* i2 = static_cast<Integer*>(p_Int2);
      mpz_t* i2p = i2->get_mpz_ptr();
      mpz_mul(*resp, *i1p, *i2p);
    } else {
      RF_LIB_ERROR("Invalid argument");
    }
  } else {
    RF_LIB_ERROR("Invalid argument");
  }
  s_Int = result;
RF_END
}

}