- interpolate_gsl_line( x, 5, oY, aperiodic, periodic);
- interpolate_gsl_line((Y-1)<<YSHIFT | x, 1, oY, aperiodic, periodic);
+ interpolate_gsl_line( x, 5, Y, aperiodic, periodic);
+ interpolate_gsl_line((Y-1)<<YSHIFT | x, 1, Y, aperiodic, periodic);
+ }
+ for (y=0; y<(Y+1)/2; y++) {
+ interpolate_gsl_line(y<<YSHIFT | ((y>>1)&1), 0, X*2+1,
+ aperiodic, periodic);