- 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 uintL* qsv = (uintL*) cl_alloca(N*sizeof(uintL));
- var uintL n;
- for (n = 0; n < N; n++) {
- init1(cl_I, av[n]) (A+n*B);
- if (n==0) {
- init1(cl_I, pv[n]) (1);
- init1(cl_I, qv[n]) (1);
- } else {
- init1(cl_I, pv[n]) (-((cl_I)(6*n-5)*(cl_I)(2*n-1)*(cl_I)(6*n-1)));
- init1(cl_I, qv[n]) ((cl_I)n*(cl_I)n*(cl_I)n*J1);
- }
- }
- var cl_pqa_series series;
- series.av = av;
- series.pv = pv; series.qv = qv;
- series.qsv = (len >= 35 ? qsv : 0); // 5% speedup for large len's
- var cl_LF fsum = eval_rational_series(N,series,actuallen);
- for (n = 0; n < N; n++) {
- av[n].~cl_I();
- pv[n].~cl_I();
- qv[n].~cl_I();
- }
+ var cl_LF fsum = eval_rational_series(N,series,actuallen,actuallen);