src/float/transcendental/cl_LF_ratsumseries_pqcd.cc

   1 // eval_pqcd_series().
   2
   3 // General includes.
   4 #include "base/cl_sysdep.h"
   5
   6 // Specification.
   7 #include "float/transcendental/cl_LF_tran.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cln/lfloat.h"
  13 #include "cln/integer.h"
  14 #include "cln/real.h"
  15 #include "float/lfloat/cl_LF.h"
  16
  17 namespace cln {
  18
  19 const cl_LF eval_pqcd_series (uintC N, cl_pqcd_series_term* args, uintC len)
  20 {
  21         if (N==0)
  22                 return cl_I_to_LF(0,len);
  23         var cl_pqcd_series_result<cl_I> sums;
  24         eval_pqcd_series_aux(N,args,sums);
  25         // Instead of computing  fsum = T/Q  and  gsum = V/(D*Q)
  26         // and then dividing them, to compute  gsum/fsum, we save two
  27         // divisions by computing  V/(D*T).
  28         return
  29           cl_I_to_LF(sums.V,len) / The(cl_LF)(sums.D * cl_I_to_LF(sums.T,len));
  30 }
  31
  32 const cl_LF eval_pqcd_series (uintC N, cl_pqcd_series_stream& args, uintC len)
  33 {
  34         if (N==0)
  35                 return cl_I_to_LF(0,len);
  36         var cl_pqcd_series_result<cl_I> sums;
  37         eval_pqcd_series_aux(N,args,sums);
  38         // Instead of computing  fsum = T/Q  and  gsum = V/(D*Q)
  39         // and then dividing them, to compute  gsum/fsum, we save two
  40         // divisions by computing  V/(D*T).
  41         return
  42           cl_I_to_LF(sums.V,len) / The(cl_LF)(sums.D * cl_I_to_LF(sums.T,len));
  43 }
  44
  45 const cl_LF eval_pqcd_series (uintC N, cl_pqcd_series_stream& args, uintC len, uintC trunclen)
  46 {
  47         if (N==0)
  48                 return cl_I_to_LF(0,len);
  49         var cl_pqcd_series_result<cl_R> sums;
  50         eval_pqcd_series_aux(N,args,sums,trunclen);
  51         // Instead of computing  fsum = T/Q  and  gsum = V/(D*Q)
  52         // and then dividing them, to compute  gsum/fsum, we save two
  53         // divisions by computing  V/(D*T).
  54         return
  55           cl_R_to_LF(sums.V,len) / The(cl_LF)(sums.D * cl_R_to_LF(sums.T,len));
  56 }
  57
  58 }  // namespace cln