Vladif5c.cin

// vladif5c.cin is the test routine in C++ for evaluation of the natural tetration through the Cauchi integral, using the displaced contour. // z_type Zo=z_type(.31813150520476413, 1.3372357014306895); z_type Zc=z_type(.31813150520476413,-1.3372357014306895);
 * 1) include "GLxw2048.inc"

z_type f4(z_type z){ //NOT SHIFTED FOR x1 !!!! #include"f2048natu.inc" //K and A are defined there int j,k,m,n; DB x,y, u, t; z_type c,d, cu,cd; //     z_type E[K],G[K]; z_type E[2048],G[2048]; DO(k,K){c=F[k];E[k]=log(c);G[k]=exp(c);} //     the initioalization abouve should run at the compillation c=0.; DO(k,K){t=A*GLx[k]; c+= GLw[k]*( G[k]/(z_type( 1.,t)-z) - E[k]/(z_type(-1.,t)-z) );} cu=.5-I/(2.*M_PI)*log( (z_type(1.,-A)+z)/(z_type(1., A)-z) ); cd=.5-I/(2.*M_PI)*log( (z_type(1.,-A)-z)/(z_type(1., A)+z) ); return c*(A/(2.*M_PI)) +Zo*cu+Zc*cd; }

//#include "x1.inc" //DB x1= -0.02454385682582; 1024;20; DB x1= 0.00743143611046; //2048;24;

z_type F4(z_type z){   DB x=Re(z); if(x<-.5) return log(F4(z+1.)); if(x> .5) return exp(F4(z-1.)); return f4(z+x1); }

z_type f5(z_type z){ //NOT SHIFTED FOR x1 !!!! //printf("f5 is called at, %9.4f %9.6f\n",Re(z), Im(z)); #include"f2048minus05.inc" //K and A are defined there int j,k,m,n; DB x,y, u, t; z_type c,d, cu,cd; //     z_type E[K],G[K]; z_type E[2048],G[2048]; z+=.5; DO(k,K){c=F[k];E[k]=log(c);G[k]=exp(c);} //     the initioalization abouve should run at the compillation c=0.; DO(k,K){t=A*GLx[k]; c+= GLw[k]*( G[k]/(z_type( 1.,t)-z) - E[k]/(z_type(-1.,t)-z) );} cu=.5-I/(2.*M_PI)*log( (z_type(1.,-A)+z)/(z_type(1., A)-z) ); cd=.5-I/(2.*M_PI)*log( (z_type(1.,-A)-z)/(z_type(1., A)+z) ); c= c*(A/(2.*M_PI)) +Zo*cu+Zc*cd; //printf("f5 is about to retunv value, %9.4f %9.6f\n",Re(c), Im(c)); return c;       }

//#include "x1.inc" //DB x1= 0.00743143611046; //2048;24; //#include "test20.inc" DB x1minus05= 0.00776084382757;

z_type F5(z_type z){   DB x=Re(z); //             printf("F5 is called at, %9.4f %9.6f\n",Re(z), Im(z)); if(x<-1.) return log(F5(z+1.)); if(x> 0.) return exp(F5(z-1.)); //             printf("F5 sends to f5\n"); return f5(z+x1minus05); } //#include"f4d1.cin" // //