Sqrt2f23l.cin

// Sqrt2f23l.cin suggests routine F23L for evaluation of real–holomorphic abelexponential to base $b\!=\!\sqrt{2}$.

//In order to evaluate $\mathrm{tet}_{\sqrt{2}}(z)$, the routine should be called as F23L(z)

z_type f23L(z_type z){ int n; z_type e,s,k;
DB TcL[23]={1., //coeff. of expansion of exp(-q(z+1.2 ...) by powers of (2-F).
 -.56472283831773236365, 0.29964618138408807683,
 -.15593239048925425850, 0.8035187974815443609e-1,
 -0.411584960662439279e-1, 0.2099852095441203541e-1,
 -0.1068258032026355653e-1, 0.542288102231591005e-2,
 -0.2748252661868267e-2, 0.13909151872677962e-2,
 -0.703181586212482131e-3, 0.35517006776480e-3,
 -0.1792537427481520668e-3, 0.9040887657183e-4,
 -0.45572543028501136e-4, 0.2296022632181e-4,
 -0.1156277075032e-4, 0.5820169657e-5,
        -0.291e-5, 0.144e-5, -.71e-6 };
z=2.-z; s=TcL[22]; for(n=21; n>=0; n--){ s*=z; s+=TcL[n]; }
// return -log(s*z)/0.36651292058166432701 -1.251551478822190;};
        return -log(-s*z)/0.36651292058166432701+2.131917787095039;};
                //.32663425997828098238;

//z_type TQ2L3(z_type z){ DB b=sqrt(2.); if(abs(z-2.)>9999.) return 9999.;
z_type F23L(z_type z){ DB b=sqrt(2.); if(abs(z-2.)>9999.) return 9999.;
                        if(abs(z-2.)>.4) return F23L(exp(z*log(b)))-1. ;
                         return f23L(z); }

/**/