Lines Matching refs:eps
45 computeI1(maxpow):=block([sintegrand,sintegrandexp,s,sigma,tau1,k2,eps],
48 (1-eps)*subst([k2=4*eps/(1-eps)^2],sintegrand),
49 eps,maxpow),
53 tau1:ataylor(s/A1,eps,maxpow),
57 A1:A1/(1-eps),
62 static inline CT evaluate_series_A1(CT eps) {
63 CT eps2 = math::sqr(eps);
67 q:horner(ataylor(subst([eps=sqrt(eps2)],A1*(1-eps)-1),eps2,n)),
73 return (t + eps) / (1 - eps);
77 computeI2(maxpow):=block([sintegrand,sintegrandexp,s,sigma,tau1,k2,eps],
80 (1+eps)*subst([k2=4*eps/(1-eps)^2],sintegrand),
81 eps,maxpow),
85 tau1:ataylor(s/A2,eps,maxpow),
89 A2:A2/(1+eps),
94 CT evaluate_series_A2(CT const& eps)
96 CT const eps2 = math::sqr(eps);
100 q:horner(ataylor(subst([eps=sqrt(eps2)],A2*(1+eps)-1),eps2,n)),
106 return (t - eps) / (1 + eps);
110 computeI3(maxpow):=block([int,intexp,dlam,eta,del,eps,nu,f,z,n],
112 int:subst([k2=4*eps/(1-eps)^2],
115 intexp:jtaylor(int,n,eps,maxpow),
119 eta:jtaylor(dlam/A3,n,eps,maxpow),
120 A3:jtaylor(A3,n,eps,maxpow),
133 jtaylor(subst([n=n],A3),n,eps,nn-1),
138 string(horner(coeff(q,eps,i))),";")),
146 static inline evaluate_coeffs_C1(CT eps, CT c[]) {
147 CT eps2 = math::sqr(eps);
148 CT d = eps;
153 subst([eps=sqrt(eps2)],ataylor(C1[m],eps,n)/eps^m)),
155 if m>1 then print(concat(tab3,"d *= eps;")),
162 revertI1(maxpow):=block([tau,eps,tauacc:1,sigacc:0],
166 eps,maxpow)),
175 static inline evaluate_coeffs_C1p(CT eps, CT c[])
177 CT const eps2 = math::sqr(eps);
178 CT d = eps;
183 subst([eps=sqrt(eps2)],ataylor(C1p[m],eps,n)/eps^m)),
185 if m>1 then print(concat(tab3,"d *= eps;")),
194 static inline void evaluate_coeffs_C2(CT const& eps, CT c[])
196 CT const eps2 = math::sqr(eps);
197 CT d = eps;
202 subst([eps=sqrt(eps2)],ataylor(C2[m],eps,n)/eps^m)),
204 if m>1 then print(concat(tab3,"d *= eps;")),
222 jtaylor(subst([n=n],C3[m]),_n,eps,nn-1),
226 string(horner(coeff(q,eps,j))),";")),