TFracNum.h

Go to the documentation of this file.

#ifndef TFracNum_h
#define TFracNum_h

//
// Uncomment the following line
// if you want to work in CINT (root.cern.ch)
//
//#define __JCINT__
#ifndef __JCINT__
#define Int_t    long long
#define Double_t double
#define Bool_t   bool
#endif

#include "Primes.h"

#ifndef __JCINT__
const Int_t MAXPRIMSQUARED=PRIMES[NPRIMFIELD-1]*PRIMES[NPRIMFIELD-1]*1000;
const char IOUTSTRING[5]="%lld";
#else
const Int_t MAXPRIMSQUARED=PRIMES[NPRIMFIELD-1]*PRIMES[NPRIMFIELD-1];
const char IOUTSTRING[3]="%d";
#endif

//
// Fractional number representation by the prime number
// decomposition of the nominator and the denominator
//
class TFracNum {
  
 private:
  //
  // since Num is appearing as short form of "Number",
  // nom/NOM is taken when the numerator is meant
  //
  // maximum prime index of numerator.
  Int_t maxPrimNom;
  
  // maximum prime index of denominator.
  Int_t maxPrimDen;
  
  // Prime number decomposition of numerator. Field length is maxPrimNom,
  //  NOM[0] is the exponent of 2, NOM[1] of 3, and so on.
  Int_t *NOM;
  // Prime number decomposition of denominator, analogue to NOM 
  Int_t *DEN;

  // Prefactor, including sign
  // Negative fractional number have sign_prefrac=-1
  // Special cases: 
  // Division by zero      (<=> infinity)     => sign_prefac=-7777
  // Division zero by zero (<=> undetermined) => sign_prefac=-6666
  Int_t sign_prefac;

  // Integers of numerator and denominator
  Int_t NOM_INT;
  Int_t DEN_INT;
  Double_t dvalue;

 public:
  TFracNum(){
    maxPrimNom=0;
    maxPrimDen=0;
    NOM=0;
    DEN=0; 
    sign_prefac=1;
  };
  
  TFracNum(
           Int_t mN,
           Int_t mD, 
           Int_t* N, 
           Int_t* D, 
           Int_t s);

  TFracNum(
           Int_t inom,
           Int_t iden);
  

  TFracNum(
           Int_t N, 
           Int_t D,  
           const char* s);
  
  Int_t DenomCommonDivisor(const TFracNum &) const;
  

  TFracNum* SumSignedRoots(TFracNum*);
  
  const char* FracString();


  const char* FracStringSqrt();

  Double_t Print() const;

  Double_t PrintToErr() const;

  Double_t Dval(){return dvalue;};
  

  Bool_t Sqrt();

  Bool_t FlipSign();

  Bool_t Abs();

  Bool_t Invert();
  
  Int_t GetSign();
  
  Int_t GetNumerator(){return NOM_INT;};
  
  Int_t GetDenominator(){return DEN_INT;};
  
  Bool_t PrintDifference(const TFracNum &) const;

  char* HeaderString();
  
  Bool_t   operator== (const TFracNum &) const;
  Bool_t   operator>  (const TFracNum &) const;
  TFracNum operator*  (const TFracNum &) const;
  TFracNum operator+  (const TFracNum &) const;
  
  Bool_t SetINTs();
  //ClassDef(TFracNum,1);
  
};

const TFracNum TFracNum_Zero( 0,1);
const TFracNum TFracNum_One ( 1,1);
const TFracNum TFracNum_Two ( 2,1);
const TFracNum TFracNum_mTwo(-2,1);
const TFracNum TFracNum_Half( 1,2);
const TFracNum TFracNum_Quarter( 1,4);

TFracNum am0_to_J(Int_t J, Int_t m, Int_t m0);
TFracNum c_sub_ell(Int_t ell);
TFracNum cm0_sub_ell(Int_t ell, Int_t m0);
TFracNum cm0_sub_ell_2(Int_t ell, Int_t m0);

#endif