convolution.c

su 的源代码库

		}
	} else if (lx==11) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==12) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==13) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==14) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==15) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==16) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==17) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		x16 = x[ifx+16];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;  z1 += x16*ya;
			yb = y[i-ifx-16];  z0 += x16*yb;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==18) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		x16 = x[ifx+16];
		x17 = x[ifx+17];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;  z1 += x16*ya;
			yb = y[i-ifx-16];  z0 += x16*yb;  z1 += x17*yb;
			ya = y[i-ifx-17];  z0 += x17*ya;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==19) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		x16 = x[ifx+16];
		x17 = x[ifx+17];
		x18 = x[ifx+18];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;  z1 += x16*ya;
			yb = y[i-ifx-16];  z0 += x16*yb;  z1 += x17*yb;
			ya = y[i-ifx-17];  z0 += x17*ya;  z1 += x18*ya;
			yb = y[i-ifx-18];  z0 += x18*yb;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==20) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		x16 = x[ifx+16];
		x17 = x[ifx+17];
		x18 = x[ifx+18];
		x19 = x[ifx+19];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;  z1 += x16*ya;
			yb = y[i-ifx-16];  z0 += x16*yb;  z1 += x17*yb;
			ya = y[i-ifx-17];  z0 += x17*ya;  z1 += x18*ya;
			yb = y[i-ifx-18];  z0 += x18*yb;  z1 += x19*yb;
			ya = y[i-ifx-19];  z0 += x19*ya;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==21) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		x16 = x[ifx+16];
		x17 = x[ifx+17];
		x18 = x[ifx+18];
		x19 = x[ifx+19];
		x20 = x[ifx+20];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;  z1 += x16*ya;
			yb = y[i-ifx-16];  z0 += x16*yb;  z1 += x17*yb;
			ya = y[i-ifx-17];  z0 += x17*ya;  z1 += x18*ya;
			yb = y[i-ifx-18];  z0 += x18*yb;  z1 += x19*yb;
			ya = y[i-ifx-19];  z0 += x19*ya;  z1 += x20*ya;
			yb = y[i-ifx-20];  z0 += x20*yb;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==22) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		x16 = x[ifx+16];
		x17 = x[ifx+17];
		x18 = x[ifx+18];
		x19 = x[ifx+19];
		x20 = x[ifx+20];
		x21 = x[ifx+21];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;  z1 += x16*ya;
			yb = y[i-ifx-16];  z0 += x16*yb;  z1 += x17*yb;
			ya = y[i-ifx-17];  z0 += x17*ya;  z1 += x18*ya;
			yb = y[i-ifx-18];  z0 += x18*yb;  z1 += x19*yb;
			ya = y[i-ifx-19];  z0 += x19*ya;  z1 += x20*ya;
			yb = y[i-ifx-20];  z0 += x20*yb;  z1 += x21*yb;
			ya = y[i-ifx-21];  z0 += x21*ya;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==23) {
		x0 = x[ifx];
		x1 = x[ifx+1];
		x2 = x[ifx+2];
		x3 = x[ifx+3];
		x4 = x[ifx+4];
		x5 = x[ifx+5];
		x6 = x[ifx+6];
		x7 = x[ifx+7];
		x8 = x[ifx+8];
		x9 = x[ifx+9];
		x10 = x[ifx+10];
		x11 = x[ifx+11];
		x12 = x[ifx+12];
		x13 = x[ifx+13];
		x14 = x[ifx+14];
		x15 = x[ifx+15];
		x16 = x[ifx+16];
		x17 = x[ifx+17];
		x18 = x[ifx+18];
		x19 = x[ifx+19];
		x20 = x[ifx+20];
		x21 = x[ifx+21];
		x22 = x[ifx+22];
		for (i=ilow; i<=ihigh-1; i+=2) {
			ya = y[i+1-ifx];  z1 = x0*ya;
			yb = y[i-ifx];  z0 = x0*yb;  z1 += x1*yb;
			ya = y[i-ifx-1];  z0 += x1*ya;  z1 += x2*ya;
			yb = y[i-ifx-2];  z0 += x2*yb;  z1 += x3*yb;
			ya = y[i-ifx-3];  z0 += x3*ya;  z1 += x4*ya;
			yb = y[i-ifx-4];  z0 += x4*yb;  z1 += x5*yb;
			ya = y[i-ifx-5];  z0 += x5*ya;  z1 += x6*ya;
			yb = y[i-ifx-6];  z0 += x6*yb;  z1 += x7*yb;
			ya = y[i-ifx-7];  z0 += x7*ya;  z1 += x8*ya;
			yb = y[i-ifx-8];  z0 += x8*yb;  z1 += x9*yb;
			ya = y[i-ifx-9];  z0 += x9*ya;  z1 += x10*ya;
			yb = y[i-ifx-10];  z0 += x10*yb;  z1 += x11*yb;
			ya = y[i-ifx-11];  z0 += x11*ya;  z1 += x12*ya;
			yb = y[i-ifx-12];  z0 += x12*yb;  z1 += x13*yb;
			ya = y[i-ifx-13];  z0 += x13*ya;  z1 += x14*ya;
			yb = y[i-ifx-14];  z0 += x14*yb;  z1 += x15*yb;
			ya = y[i-ifx-15];  z0 += x15*ya;  z1 += x16*ya;
			yb = y[i-ifx-16];  z0 += x16*yb;  z1 += x17*yb;
			ya = y[i-ifx-17];  z0 += x17*ya;  z1 += x18*ya;
			yb = y[i-ifx-18];  z0 += x18*yb;  z1 += x19*yb;
			ya = y[i-ifx-19];  z0 += x19*ya;  z1 += x20*ya;
			yb = y[i-ifx-20];  z0 += x20*yb;  z1 += x21*yb;
			ya = y[i-ifx-21];  z0 += x21*ya;  z1 += x22*ya;
			yb = y[i-ifx-22];  z0 += x22*yb;
			z[i+1] = z1;
			z[i] = z0;
		}
	} else if (lx==24) {