// to more than double the floating-point precision because of the large
// extinction which takes place. But luckily we compute with integers.
var uintC actuallen = len+1; // 1 Schutz-Digit
// to more than double the floating-point precision because of the large
// extinction which takes place. But luckily we compute with integers.
var uintC actuallen = len+1; // 1 Schutz-Digit
CL_ALLOCA_STACK;
var cl_I* bv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* pv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* qv = (cl_I*) cl_alloca(N*sizeof(cl_I));
CL_ALLOCA_STACK;
var cl_I* bv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* pv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* qv = (cl_I*) cl_alloca(N*sizeof(cl_I));
for (n = 0; n < N; n++) {
init1(cl_I, bv[n]) (n+1);
init1(cl_I, pv[n]) (n==0 ? (cl_I)z : -(cl_I)z);
for (n = 0; n < N; n++) {
init1(cl_I, bv[n]) (n+1);
init1(cl_I, pv[n]) (n==0 ? (cl_I)z : -(cl_I)z);
// each.
// We compute f(x) classically and g(x) using the partial sums of f(x).
var uintC actuallen = len+2; // 2 Schutz-Digits
// each.
// We compute f(x) classically and g(x) using the partial sums of f(x).
var uintC actuallen = len+2; // 2 Schutz-Digits
var cl_LF one = cl_I_to_LF(1,actuallen);
var cl_LF fterm = one;
var cl_LF fsum = fterm;
var cl_LF gterm = cl_I_to_LF(0,actuallen);
var cl_LF gsum = gterm;
var cl_LF one = cl_I_to_LF(1,actuallen);
var cl_LF fterm = one;
var cl_LF fsum = fterm;
var cl_LF gterm = cl_I_to_LF(0,actuallen);
var cl_LF gsum = gterm;
// After n loops
// fterm = x^n/n!, fsum = 1 + x/1! + ... + x^n/n!,
// gterm = H_n*x^n/n!, gsum = H_1*x/1! + ... + H_n*x^n/n!.
// After n loops
// fterm = x^n/n!, fsum = 1 + x/1! + ... + x^n/n!,
// gterm = H_n*x^n/n!, gsum = H_1*x/1! + ... + H_n*x^n/n!.
const cl_LF compute_eulerconst_expintegral2 (uintC len)
{
var uintC actuallen = len+2; // 2 Schutz-Digits
const cl_LF compute_eulerconst_expintegral2 (uintC len)
{
var uintC actuallen = len+2; // 2 Schutz-Digits
for (n = 0; n < N; n++) {
init1(cl_I, args[n].p) (x);
init1(cl_I, args[n].q) (n+1);
for (n = 0; n < N; n++) {
init1(cl_I, args[n].p) (x);
init1(cl_I, args[n].q) (n+1);
{
// We compute f(x) classically and g(x) using the partial sums of f(x).
var uintC actuallen = len+1; // 1 Schutz-Digit
{
// We compute f(x) classically and g(x) using the partial sums of f(x).
var uintC actuallen = len+1; // 1 Schutz-Digit
- var uintL sx = (uintL)(0.25*0.693148*intDsize*actuallen)+1;
- var uintL N = (uintL)(3.591121477*sx);
+ var uintC sx = (uintC)(0.25*0.693148*intDsize*actuallen)+1;
+ var uintC N = (uintC)(3.591121477*sx);
var cl_LF fterm = cl_I_to_LF(1,actuallen);
var cl_LF fsum = fterm;
var cl_LF gterm = cl_I_to_LF(0,actuallen);
var cl_LF gsum = gterm;
var cl_LF fterm = cl_I_to_LF(1,actuallen);
var cl_LF fsum = fterm;
var cl_LF gterm = cl_I_to_LF(0,actuallen);
var cl_LF gsum = gterm;
// After n loops
// fterm = x^n/n!^2, fsum = 1 + x/1!^2 + ... + x^n/n!^2,
// gterm = H_n*x^n/n!^2, gsum = H_1*x/1!^2 + ... + H_n*x^n/n!^2.
// After n loops
// fterm = x^n/n!^2, fsum = 1 + x/1!^2 + ... + x^n/n!^2,
// gterm = H_n*x^n/n!^2, gsum = H_1*x/1!^2 + ... + H_n*x^n/n!^2.
// (Because HD_n grows like exp(n), hence HN_n grows like exp(n) as
// well, and we store all HN_n values in an array!)
var uintC actuallen = len+1; // 1 Schutz-Digit
// (Because HD_n grows like exp(n), hence HN_n grows like exp(n) as
// well, and we store all HN_n values in an array!)
var uintC actuallen = len+1; // 1 Schutz-Digit
- var uintL sx = (uintL)(0.25*0.693148*intDsize*actuallen)+1;
- var uintL N = (uintL)(3.591121477*sx);
+ var uintC sx = (uintC)(0.25*0.693148*intDsize*actuallen)+1;
+ var uintC N = (uintC)(3.591121477*sx);
var cl_I x = square((cl_I)sx);
CL_ALLOCA_STACK;
var cl_I* av = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* pv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* qv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I x = square((cl_I)sx);
CL_ALLOCA_STACK;
var cl_I* av = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* pv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* qv = (cl_I*) cl_alloca(N*sizeof(cl_I));
const cl_LF compute_eulerconst_besselintegral3 (uintC len)
{
var uintC actuallen = len+1; // 1 Schutz-Digit
const cl_LF compute_eulerconst_besselintegral3 (uintC len)
{
var uintC actuallen = len+1; // 1 Schutz-Digit
- var uintL sx = (uintL)(0.25*0.693148*intDsize*actuallen)+1;
- var uintL N = (uintL)(3.591121477*sx);
+ var uintC sx = (uintC)(0.25*0.693148*intDsize*actuallen)+1;
+ var uintC N = (uintC)(3.591121477*sx);
var cl_I x = square((cl_I)sx);
CL_ALLOCA_STACK;
var cl_I* pv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* qv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I x = square((cl_I)sx);
CL_ALLOCA_STACK;
var cl_I* pv = (cl_I*) cl_alloca(N*sizeof(cl_I));
var cl_I* qv = (cl_I*) cl_alloca(N*sizeof(cl_I));
const cl_LF compute_eulerconst_besselintegral4 (uintC len)
{
var uintC actuallen = len+2; // 2 Schutz-Digits
const cl_LF compute_eulerconst_besselintegral4 (uintC len)
{
var uintC actuallen = len+2; // 2 Schutz-Digits
- var uintL sx = (uintL)(0.25*0.693148*intDsize*actuallen)+1;
- var uintL N = (uintL)(3.591121477*sx);
+ var uintC sx = (uintC)(0.25*0.693148*intDsize*actuallen)+1;
+ var uintC N = (uintC)(3.591121477*sx);
var cl_I x = square((cl_I)sx);
CL_ALLOCA_STACK;
var cl_pqd_series_term* args = (cl_pqd_series_term*) cl_alloca(N*sizeof(cl_pqd_series_term));
var cl_I x = square((cl_I)sx);
CL_ALLOCA_STACK;
var cl_pqd_series_term* args = (cl_pqd_series_term*) cl_alloca(N*sizeof(cl_pqd_series_term));
for (n = 0; n < N; n++) {
init1(cl_I, args[n].p) (x);
init1(cl_I, args[n].q) (square((cl_I)(n+1)));
for (n = 0; n < N; n++) {
init1(cl_I, args[n].p) (x);
init1(cl_I, args[n].q) (square((cl_I)(n+1)));