SymEngine::NumerDenomVisitor Class Reference

Inheritance diagram for SymEngine::NumerDenomVisitor:

Collaboration diagram for SymEngine::NumerDenomVisitor:

Public Member Functions
	NumerDenomVisitor (const Ptr< RCP< const Basic > > &numer, const Ptr< RCP< const Basic > > &denom)
void	apply (const Basic &b)
void	bvisit (const Mul &x)
void	bvisit (const Add &x)
void	bvisit (const Pow &x)
void	bvisit (const Complex &x)
void	bvisit (const Rational &x)
void	bvisit (const Basic &x)

Private Attributes
Ptr< RCP< const Basic > >	numer_
Ptr< RCP< const Basic > >	denom_

Detailed Description

Definition at line 6 of file numer_denom.cpp.

Constructor & Destructor Documentation

NumerDenomVisitor()

SymEngine::NumerDenomVisitor::NumerDenomVisitor ( const Ptr< RCP< const Basic > > &  numer, const Ptr< RCP< const Basic > > &  denom  )

inline

Definition at line 12 of file numer_denom.cpp.

        : numer_{numer}, denom_{denom}
    {
    }

Member Function Documentation

apply()

void SymEngine::NumerDenomVisitor::apply ( const Basic &  b )

inline

Definition at line 18 of file numer_denom.cpp.

    {
        b.accept(*this);
    }

bvisit() [1/6]

void SymEngine::NumerDenomVisitor::bvisit ( const Add &  x )

inline

Definition at line 49 of file numer_denom.cpp.

    {
 
        RCP<const Basic> curr_num = zero;
        RCP<const Basic> curr_den = one;
        RCP<const Basic> arg_num, arg_den, den_mul, divx;
        RCP<const Basic> divx_num, divx_den;
 
        for (const auto &arg : x.get_args()) {
            // TODO: This is naive and slow. Fix it
            as_numer_denom(arg, outArg(arg_num), outArg(arg_den));
 
            divx = div(arg_den, curr_den);
            as_numer_denom(divx, outArg(divx_num), outArg(divx_den));
            if (eq(*divx_den, *one)) {
                // the curr_den completely divides the arg_den
                curr_den = arg_den;
                curr_num = add(mul(curr_num, divx), arg_num);
                continue;
            }
 
            divx = div(curr_den, arg_den);
            as_numer_denom(divx, outArg(divx_num), outArg(divx_den));
            // the below two lines, cover the general case, as well as the case
            // where arg_den completely divides curr_den
            curr_den = mul(curr_den, divx_den);
            curr_num = add(mul(curr_num, divx_den), mul(arg_num, divx_num));
        }
 
        *numer_ = curr_num;
        *denom_ = curr_den;
    }

bvisit() [2/6]

void SymEngine::NumerDenomVisitor::bvisit ( const Basic &  x )

inline

Definition at line 126 of file numer_denom.cpp.

    {
        *numer_ = x.rcp_from_this();
        *denom_ = one;
    }

bvisit() [3/6]

void SymEngine::NumerDenomVisitor::bvisit ( const Complex &  x )

inline

Definition at line 100 of file numer_denom.cpp.

    {
 
        RCP<const Integer> den, den1, den2;
        RCP<const Integer> num1, num2;
 
        num1 = integer(get_num(x.real_));
        num2 = integer(get_num(x.imaginary_));
 
        den1 = integer(get_den(x.real_));
        den2 = integer(get_den(x.imaginary_));
        den = lcm(*den1, *den2);
 
        num1 = rcp_static_cast<const Integer>(mul(num1, div(den, den1)));
        num2 = rcp_static_cast<const Integer>(mul(num2, div(den, den2)));
 
        *numer_ = Complex::from_two_nums(*num1, *num2);
        *denom_ = den;
    }

bvisit() [4/6]

void SymEngine::NumerDenomVisitor::bvisit ( const Mul &  x )

inline

Definition at line 23 of file numer_denom.cpp.

    {
        RCP<const Basic> curr = one;
        RCP<const Basic> arg_num, arg_den, t;
 
        for (const auto &arg : x.get_args()) {
            as_numer_denom(arg, outArg(arg_num), outArg(arg_den));
            curr = div(mul(curr, arg_num), arg_den);
        }
 
        if (not is_a<Mul>(*curr)) {
            apply(*curr);
            return;
        }
 
        RCP<const Basic> curr_num = one, curr_den = one;
        for (const auto &arg : curr->get_args()) {
            as_numer_denom(arg, outArg(arg_num), outArg(arg_den));
            curr_num = mul(curr_num, arg_num);
            curr_den = mul(curr_den, arg_den);
        }
 
        *numer_ = curr_num;
        *denom_ = curr_den;
    }

bvisit() [5/6]

void SymEngine::NumerDenomVisitor::bvisit ( const Pow &  x )

inline

Definition at line 82 of file numer_denom.cpp.

    {
 
        RCP<const Basic> base_, exp_, num, den;
        base_ = x.get_base();
        exp_ = x.get_exp();
        as_numer_denom(base_, outArg(num), outArg(den));
 
        // if the exp is a negative numer, or is intuitively 'negative'
        if (handle_minus(exp_, outArg(exp_))) {
            *numer_ = pow(den, exp_);
            *denom_ = pow(num, exp_);
        } else {
            *numer_ = pow(num, exp_);
            *denom_ = pow(den, exp_);
        }
    }

bvisit() [6/6]

void SymEngine::NumerDenomVisitor::bvisit ( const Rational &  x )

inline

Definition at line 120 of file numer_denom.cpp.

    {
        *numer_ = x.get_num();
        *denom_ = x.get_den();
    }

Field Documentation

denom_

Ptr<RCP<const Basic> > SymEngine::NumerDenomVisitor::denom_

private

Definition at line 9 of file numer_denom.cpp.

numer_

Ptr<RCP<const Basic> > SymEngine::NumerDenomVisitor::numer_

private

Definition at line 9 of file numer_denom.cpp.

---

The documentation for this class was generated from the following file:

• /home/runner/work/symengine/symengine/symengine/numer_denom.cpp