cc_dim1 = *ido;    cc_offset = 1 + (cc_dim1 << 2);    cc -= cc_offset;    --wa1;    --wa2;    /* Function Body */    if (*ido != 2) {	goto L102;    }    i__1 = *l1;    for (k = 1; k <= i__1; ++k) {	tr2 = cc[(k * 3 + 2) * cc_dim1 + 1] + cc[(k * 3 + 3) * cc_dim1 + 1];	cr2 = cc[(k * 3 + 1) * cc_dim1 + 1] + taur * tr2;	ch[(k + ch_dim2) * ch_dim1 + 1] = cc[(k * 3 + 1) * cc_dim1 + 1] + tr2;	ti2 = cc[(k * 3 + 2) * cc_dim1 + 2] + cc[(k * 3 + 3) * cc_dim1 + 2];	ci2 = cc[(k * 3 + 1) * cc_dim1 + 2] + taur * ti2;	ch[(k + ch_dim2) * ch_dim1 + 2] = cc[(k * 3 + 1) * cc_dim1 + 2] + ti2;	cr3 = taui * (cc[(k * 3 + 2) * cc_dim1 + 1] - cc[(k * 3 + 3) * 		cc_dim1 + 1]);	ci3 = taui * (cc[(k * 3 + 2) * cc_dim1 + 2] - cc[(k * 3 + 3) * 		cc_dim1 + 2]);	ch[(k + (ch_dim2 << 1)) * ch_dim1 + 1] = cr2 - ci3;	ch[(k + ch_dim2 * 3) * ch_dim1 + 1] = cr2 + ci3;	ch[(k + (ch_dim2 << 1)) * ch_dim1 + 2] = ci2 + cr3;	ch[(k + ch_dim2 * 3) * ch_dim1 + 2] = ci2 - cr3;/* L101: */    }    return;L102:    i__1 = *l1;    for (k = 1; k <= i__1; ++k) {	i__2 = *ido;	for (i__ = 2; i__ <= i__2; i__ += 2) {	    tr2 = cc[i__ - 1 + (k * 3 + 2) * cc_dim1] + cc[i__ - 1 + (k * 3 + 		    3) * cc_dim1];	    cr2 = cc[i__ - 1 + (k * 3 + 1) * cc_dim1] + taur * tr2;	    ch[i__ - 1 + (k + ch_dim2) * ch_dim1] = cc[i__ - 1 + (k * 3 + 1) *		     cc_dim1] + tr2;	    ti2 = cc[i__ + (k * 3 + 2) * cc_dim1] + cc[i__ + (k * 3 + 3) * 		    cc_dim1];	    ci2 = cc[i__ + (k * 3 + 1) * cc_dim1] + taur * ti2;	    ch[i__ + (k + ch_dim2) * ch_dim1] = cc[i__ + (k * 3 + 1) * 		    cc_dim1] + ti2;	    cr3 = taui * (cc[i__ - 1 + (k * 3 + 2) * cc_dim1] - cc[i__ - 1 + (		    k * 3 + 3) * cc_dim1]);	    ci3 = taui * (cc[i__ + (k * 3 + 2) * cc_dim1] - cc[i__ + (k * 3 + 		    3) * cc_dim1]);	    dr2 = cr2 - ci3;	    dr3 = cr2 + ci3;	    di2 = ci2 + cr3;	    di3 = ci2 - cr3;	    ch[i__ + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i__ - 1] * di2 + 		    wa1[i__] * dr2;	    ch[i__ - 1 + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i__ - 1] * dr2 		    - wa1[i__] * di2;	    ch[i__ + (k + ch_dim2 * 3) * ch_dim1] = wa2[i__ - 1] * di3 + wa2[		    i__] * dr3;	    ch[i__ - 1 + (k + ch_dim2 * 3) * ch_dim1] = wa2[i__ - 1] * dr3 - 		    wa2[i__] * di3;/* L103: */	}/* L104: */    }    return;} /* passb3_ *//* Subroutine */ static void passb4(integer_t *ido, integer_t *l1, real_t *cc, 	real_t *ch, real_t *wa1, real_t *wa2, real_t *wa3){    /* System generated locals */    integer_t cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;    /* Local variables */    integer_t i__, k;    real_t ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2, 	    tr3, tr4;    /* Parameter adjustments */    ch_dim1 = *ido;    ch_dim2 = *l1;    ch_offset = 1 + ch_dim1 * (1 + ch_dim2);    ch -= ch_offset;    cc_dim1 = *ido;    cc_offset = 1 + cc_dim1 * 5;    cc -= cc_offset;    --wa1;    --wa2;    --wa3;    /* Function Body */    if (*ido != 2) {	goto L102;    }    i__1 = *l1;    for (k = 1; k <= i__1; ++k) {	ti1 = cc[((k << 2) + 1) * cc_dim1 + 2] - cc[((k << 2) + 3) * cc_dim1 		+ 2];	ti2 = cc[((k << 2) + 1) * cc_dim1 + 2] + cc[((k << 2) + 3) * cc_dim1 		+ 2];	tr4 = cc[((k << 2) + 4) * cc_dim1 + 2] - cc[((k << 2) + 2) * cc_dim1 		+ 2];	ti3 = cc[((k << 2) + 2) * cc_dim1 + 2] + cc[((k << 2) + 4) * cc_dim1 		+ 2];	tr1 = cc[((k << 2) + 1) * cc_dim1 + 1] - cc[((k << 2) + 3) * cc_dim1 		+ 1];	tr2 = cc[((k << 2) + 1) * cc_dim1 + 1] + cc[((k << 2) + 3) * cc_dim1 		+ 1];	ti4 = cc[((k << 2) + 2) * cc_dim1 + 1] - cc[((k << 2) + 4) * cc_dim1 		+ 1];	tr3 = cc[((k << 2) + 2) * cc_dim1 + 1] + cc[((k << 2) + 4) * cc_dim1 		+ 1];	ch[(k + ch_dim2) * ch_dim1 + 1] = tr2 + tr3;	ch[(k + ch_dim2 * 3) * ch_dim1 + 1] = tr2 - tr3;	ch[(k + ch_dim2) * ch_dim1 + 2] = ti2 + ti3;	ch[(k + ch_dim2 * 3) * ch_dim1 + 2] = ti2 - ti3;	ch[(k + (ch_dim2 << 1)) * ch_dim1 + 1] = tr1 + tr4;	ch[(k + (ch_dim2 << 2)) * ch_dim1 + 1] = tr1 - tr4;	ch[(k + (ch_dim2 << 1)) * ch_dim1 + 2] = ti1 + ti4;	ch[(k + (ch_dim2 << 2)) * ch_dim1 + 2] = ti1 - ti4;/* L101: */    }    return;L102:    i__1 = *l1;    for (k = 1; k <= i__1; ++k) {	i__2 = *ido;	for (i__ = 2; i__ <= i__2; i__ += 2) {	    ti1 = cc[i__ + ((k << 2) + 1) * cc_dim1] - cc[i__ + ((k << 2) + 3)		     * cc_dim1];	    ti2 = cc[i__ + ((k << 2) + 1) * cc_dim1] + cc[i__ + ((k << 2) + 3)		     * cc_dim1];	    ti3 = cc[i__ + ((k << 2) + 2) * cc_dim1] + cc[i__ + ((k << 2) + 4)		     * cc_dim1];	    tr4 = cc[i__ + ((k << 2) + 4) * cc_dim1] - cc[i__ + ((k << 2) + 2)		     * cc_dim1];	    tr1 = cc[i__ - 1 + ((k << 2) + 1) * cc_dim1] - cc[i__ - 1 + ((k <<		     2) + 3) * cc_dim1];	    tr2 = cc[i__ - 1 + ((k << 2) + 1) * cc_dim1] + cc[i__ - 1 + ((k <<		     2) + 3) * cc_dim1];	    ti4 = cc[i__ - 1 + ((k << 2) + 2) * cc_dim1] - cc[i__ - 1 + ((k <<		     2) + 4) * cc_dim1];	    tr3 = cc[i__ - 1 + ((k << 2) + 2) * cc_dim1] + cc[i__ - 1 + ((k <<		     2) + 4) * cc_dim1];	    ch[i__ - 1 + (k + ch_dim2) * ch_dim1] = tr2 + tr3;	    cr3 = tr2 - tr3;	    ch[i__ + (k + ch_dim2) * ch_dim1] = ti2 + ti3;	    ci3 = ti2 - ti3;	    cr2 = tr1 + tr4;	    cr4 = tr1 - tr4;	    ci2 = ti1 + ti4;	    ci4 = ti1 - ti4;	    ch[i__ - 1 + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i__ - 1] * cr2 		    - wa1[i__] * ci2;	    ch[i__ + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i__ - 1] * ci2 + 		    wa1[i__] * cr2;	    ch[i__ - 1 + (k + ch_dim2 * 3) * ch_dim1] = wa2[i__ - 1] * cr3 - 		    wa2[i__] * ci3;	    ch[i__ + (k + ch_dim2 * 3) * ch_dim1] = wa2[i__ - 1] * ci3 + wa2[		    i__] * cr3;	    ch[i__ - 1 + (k + (ch_dim2 << 2)) * ch_dim1] = wa3[i__ - 1] * cr4 		    - wa3[i__] * ci4;	    ch[i__ + (k + (ch_dim2 << 2)) * ch_dim1] = wa3[i__ - 1] * ci4 + 		    wa3[i__] * cr4;/* L103: */	}/* L104: */    }    return;} /* passb4_ *//* Subroutine */ static void passb5(integer_t *ido, integer_t *l1, real_t *cc, 	real_t *ch, real_t *wa1, real_t *wa2, real_t *wa3, 	real_t *wa4){    /* Initialized data */    static real_t tr11 = 	    REAL_CONSTANT(0.3090169943749474241022934171828195886015458990288143106772431137);    static real_t ti11 = 	    REAL_CONSTANT(0.9510565162951535721164393337938214340569863412575022244730564442);    static real_t tr12 = 	    REAL_CONSTANT(-0.8090169943749474241022934171828190588601545899028814310677431135);    static real_t ti12 = 	    REAL_CONSTANT(0.5877852522924731291687059546390727685976524376431459107227248076);    /* System generated locals */    integer_t cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;    /* Local variables */    integer_t i__, k;    real_t ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4, 	    ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;    /* Parameter adjustments */    ch_dim1 = *ido;    ch_dim2 = *l1;    ch_offset = 1 + ch_dim1 * (1 + ch_dim2);    ch -= ch_offset;    cc_dim1 = *ido;    cc_offset = 1 + cc_dim1 * 6;    cc -= cc_offset;    --wa1;    --wa2;    --wa3;    --wa4;    /* Function Body */    if (*ido != 2) {	goto L102;    }    i__1 = *l1;    for (k = 1; k <= i__1; ++k) {	ti5 = cc[(k * 5 + 2) * cc_dim1 + 2] - cc[(k * 5 + 5) * cc_dim1 + 2];	ti2 = cc[(k * 5 + 2) * cc_dim1 + 2] + cc[(k * 5 + 5) * cc_dim1 + 2];	ti4 = cc[(k * 5 + 3) * cc_dim1 + 2] - cc[(k * 5 + 4) * cc_dim1 + 2];	ti3 = cc[(k * 5 + 3) * cc_dim1 + 2] + cc[(k * 5 + 4) * cc_dim1 + 2];	tr5 = cc[(k * 5 + 2) * cc_dim1 + 1] - cc[(k * 5 + 5) * cc_dim1 + 1];	tr2 = cc[(k * 5 + 2) * cc_dim1 + 1] + cc[(k * 5 + 5) * cc_dim1 + 1];	tr4 = cc[(k * 5 + 3) * cc_dim1 + 1] - cc[(k * 5 + 4) * cc_dim1 + 1];	tr3 = cc[(k * 5 + 3) * cc_dim1 + 1] + cc[(k * 5 + 4) * cc_dim1 + 1];	ch[(k + ch_dim2) * ch_dim1 + 1] = cc[(k * 5 + 1) * cc_dim1 + 1] + tr2 		+ tr3;	ch[(k + ch_dim2) * ch_dim1 + 2] = cc[(k * 5 + 1) * cc_dim1 + 2] + ti2 		+ ti3;	cr2 = cc[(k * 5 + 1) * cc_dim1 + 1] + tr11 * tr2 + tr12 * tr3;	ci2 = cc[(k * 5 + 1) * cc_dim1 + 2] + tr11 * ti2 + tr12 * ti3;	cr3 = cc[(k * 5 + 1) * cc_dim1 + 1] + tr12 * tr2 + tr11 * tr3;	ci3 = cc[(k * 5 + 1) * cc_dim1 + 2] + tr12 * ti2 + tr11 * ti3;	cr5 = ti11 * tr5 + ti12 * tr4;	ci5 = ti11 * ti5 + ti12 * ti4;	cr4 = ti12 * tr5 - ti11 * tr4;	ci4 = ti12 * ti5 - ti11 * ti4;	ch[(k + (ch_dim2 << 1)) * ch_dim1 + 1] = cr2 - ci5;	ch[(k + ch_dim2 * 5) * ch_dim1 + 1] = cr2 + ci5;	ch[(k + (ch_dim2 << 1)) * ch_dim1 + 2] = ci2 + cr5;	ch[(k + ch_dim2 * 3) * ch_dim1 + 2] = ci3 + cr4;	ch[(k + ch_dim2 * 3) * ch_dim1 + 1] = cr3 - ci4;	ch[(k + (ch_dim2 << 2)) * ch_dim1 + 1] = cr3 + ci4;	ch[(k + (ch_dim2 << 2)) * ch_dim1 + 2] = ci3 - cr4;	ch[(k + ch_dim2 * 5) * ch_dim1 + 2] = ci2 - cr5;/* L101: */    }    return;L102:    i__1 = *l1;    for (k = 1; k <= i__1; ++k) {	i__2 = *ido;	for (i__ = 2; i__ <= i__2; i__ += 2) {	    ti5 = cc[i__ + (k * 5 + 2) * cc_dim1] - cc[i__ + (k * 5 + 5) * 		    cc_dim1];	    ti2 = cc[i__ + (k * 5 + 2) * cc_dim1] + cc[i__ + (k * 5 + 5) * 		    cc_dim1];	    ti4 = cc[i__ + (k * 5 + 3) * cc_dim1] - cc[i__ + (k * 5 + 4) * 		    cc_dim1];	    ti3 = cc[i__ + (k * 5 + 3) * cc_dim1] + cc[i__ + (k * 5 + 4) * 		    cc_dim1];	    tr5 = cc[i__ - 1 + (k * 5 + 2) * cc_dim1] - cc[i__ - 1 + (k * 5 + 		    5) * cc_dim1];	    tr2 = cc[i__ - 1 + (k * 5 + 2) * cc_dim1] + cc[i__ - 1 + (k * 5 + 		    5) * cc_dim1];	    tr4 = cc[i__ - 1 + (k * 5 + 3) * cc_dim1] - cc[i__ - 1 + (k * 5 + 		    4) * cc_dim1];	    tr3 = cc[i__ - 1 + (k * 5 + 3) * cc_dim1] + cc[i__ - 1 + (k * 5 + 		    4) * cc_dim1];	    ch[i__ - 1 + (k + ch_dim2) * ch_dim1] = cc[i__ - 1 + (k * 5 + 1) *		     cc_dim1] + tr2 + tr3;	    ch[i__ + (k + ch_dim2) * ch_dim1] = cc[i__ + (k * 5 + 1) * 		    cc_dim1] + ti2 + ti3;	    cr2 = cc[i__ - 1 + (k * 5 + 1) * cc_dim1] + tr11 * tr2 + tr12 * 		    tr3;	    ci2 = cc[i__ + (k * 5 + 1) * cc_dim1] + tr11 * ti2 + tr12 * ti3;	    cr3 = cc[i__ - 1 + (k * 5 + 1) * cc_dim1] + tr12 * tr2 + tr11 * 		    tr3;	    ci3 = cc[i__ + (k * 5 + 1) * cc_dim1] + tr12 * ti2 + tr11 * ti3;	    cr5 = ti11 * tr5 + ti12 * tr4;	    ci5 = ti11 * ti5 + ti12 * ti4;	    cr4 = ti12 * tr5 - ti11 * tr4;	    ci4 = ti12 * ti5 - ti11 * ti4;	    dr3 = cr3 - ci4;	    dr4 = cr3 + ci4;	    di3 = ci3 + cr4;	    di4 = ci3 - cr4;	    dr5 = cr2 + ci5;	    dr2 = cr2 - ci5;	    di5 = ci2 - cr5;	    di2 = ci2 + cr5;	    ch[i__ - 1 + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i__ - 1] * dr2 		    - wa1[i__] * di2;	    ch[i__ + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i__ - 1] * di2 + 		    wa1[i__] * dr2;	    ch[i__ - 1 + (k + ch_dim2 * 3) * ch_dim1] = wa2[i__ - 1] * dr3 - 		    wa2[i__] * di3;	    ch[i__ + (k + ch_dim2 * 3) * ch_dim1] = wa2[i__ - 1] * di3 + wa2[		    i__] * dr3;	    ch[i__ - 1 + (k + (ch_dim2 << 2)) * ch_dim1] = wa3[i__ - 1] * dr4 		    - wa3[i__] * di4;	    ch[i__ + (k + (ch_dim2 << 2)) * ch_dim1] = wa3[i__ - 1] * di4 + 		    wa3[i__] * dr4;	    ch[i__ - 1 + (k + ch_dim2 * 5) * ch_dim1] = wa4[i__ - 1] * dr5 - 		    wa4[i__] * di5;	    ch[i__ + (k + ch_dim2 * 5) * ch_dim1] = wa4[i__ - 1] * di5 + wa4[		    i__] * dr5;/* L103: */	}/* L104: */    }    return;} /* passb5_ *//* Subroutine */ static void passf(integer_t *nac, integer_t *ido, integer_t *ip, integer_t *	l1, integer_t *idl1, real_t *cc, real_t *c1, real_t *c2, 	real_t *ch, real_t *ch2, real_t *wa){    /* System generated locals */    integer_t ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,	     c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset, 	    i__1, i__2, i__3;    /* Local variables */    integer_t i__, j, k, l, jc, lc, ik, nt, idj, idl, inc, idp;    real_t wai, war;    integer_t ipp2, idij, idlj, idot, ipph;    /* Parameter adjustments */    ch_dim1 = *ido;    ch_dim2 = *l1;    ch_offset = 1 + ch_dim1 * (1 + ch_dim2);    ch -= ch_offset;    c1_dim1 = *ido;    c1_dim2 = *l1;    c1_offset = 1 + c1_dim1 * (1 + c1_dim2);    c1 -= c1_offset;    cc_dim1 = *ido;    cc_dim2 = *ip;    cc_offset = 1 + cc_dim1 * (1 + cc_dim2);    cc -= cc_offset;    ch2_dim1 = *idl1;    ch2_offset = 1 + ch2_dim1;    ch2 -= ch2_offset;    c2_dim1 = *idl1;    c2_offset = 1 + c2_dim1;    c2 -= c2_offset;    --wa;    /* Function Body */    idot = *ido / 2;    nt = *ip * *idl1;    ipp2 = *ip + 2;    ipph = (*ip + 1) / 2;    idp = *ip * *ido;    if (*ido < *l1) {	goto L106;    }    i__1 = ipph;    for (j = 2; j <= i__1; ++j) {	jc = ipp2 - j;	i__2 = *l1;	for (k = 1; k <= i__2; ++k) {	    i__3 = *ido;	    for (i__ = 1; i__ <= i__3; ++i__) {		ch[i__ + (k + j * ch_dim2) * ch_dim1] = cc[i__ + (j + k * 			cc_dim2) * cc_dim1] + cc[i__ + (jc + k * cc_dim2) * 			cc_dim1];		ch[i__ + (k + jc * ch_dim2) * ch_dim1] = cc[i__ + (j + k * 			cc_dim2) * cc_dim1] - cc[i__ + (jc + k * cc_dim2) * 			cc_dim1];/* L101: */	    }/* L102: */	}/* L103: */    }    i__1 = *l1;    for (k = 1; k <= i__1; ++k) {	i__2 = *ido;	for (i__ = 1; i__ <= i__2; ++i__) {	    ch[i__ + (k + ch_dim2) * ch_dim1] = cc[i__ + (k * cc_dim2 + 1) * 		    cc_dim1];/* L104: */	}/* L105: */    }    goto L112;L106:    i__1 = ipph;    for (j = 2; j <= i__1; ++j) {	jc = ipp2 - j;	i__2 = *ido;	for (i__ = 1; i__ <= i__2; ++i__) {	    i__3 = *l1;	    for (k = 1; k <= i__3; ++k) {		ch[i__ + (k + j * ch_dim2) * ch_dim1] = cc[i__ + (j + k * 			cc_dim2) * cc_dim1] + cc[i__ + (jc + k * cc_dim2) * 			cc_dim1];		ch[i__ + (k + jc * ch_dim2) * ch_dim1] = cc[i__ + (j + k * 			cc_dim2) * cc_dim1] - cc[i__ + (jc + k * cc_dim2) * 			cc_dim1];/* L107: */	    }/* L108: */	}/* L109: */    }    i__1 = *ido;    for (i__ = 1; i__ <= i__1; ++i__) {	i__2 = *l1;	for (k = 1; k <= i__2; ++k) {	    ch[i__ + (k + ch_dim2) * ch_dim1] = cc[i__ + (k * cc_dim2 + 1) * 		    cc_dim1];/* L110: */	}/* L111: */    }L112:    idl = 2 - *ido;    inc = 0;    i__1 = ipph;    for (l = 2; l <= i__1; ++l) {	lc = ipp2 - l;	idl += *ido;	i__2 = *idl1;	for (ik = 1; ik <= i__2; ++ik) {	    c2[ik + l * c2_dim1] = ch2[ik + ch2_dim1] + wa[idl - 1] * ch2[ik 		    + (ch2_dim1 << 1)];	    c2[ik + lc * c2_dim1] = -wa[idl] * ch2[ik + *ip * ch2_dim1];/* L113: */	}	idlj = idl;	inc += *ido;	i__2 = ipph;	for (j = 3; j <= i__2; ++j) {	    jc = ipp2 - j;	    idlj += inc;	    if (idlj > idp) {		idlj -= idp;	    }	    war = wa[idlj - 1];	    wai = wa[idlj];	    i__3 = *idl1;	    for (ik = 1; ik <= i__3; ++ik) {		c2[ik + l * c2_dim1] += war * ch2[ik + j * ch2_dim1];		c2[ik + lc * c2_dim1] -= wai * ch2[ik + jc * ch2_dim1];/* L114: */	    }/* L115: */	}/* L116: */    }    i__1 = ipph;    for (j = 2; j <= i__1; ++j) {	i__2 = *idl1;	for (ik = 1; ik <= i__2; ++ik) {	    ch2[ik + ch2_dim1] += ch2[ik + j * ch2_dim1];/* L117: */	}/* L118: */    }    i__1 = ipph;    for (j = 2; j <= i__1; ++j) {	jc = ipp2 - j;	i__2 = *idl1;	for (ik = 2; ik <= i__2; ik += 2) {	    ch2[ik - 1 + j * ch2_dim1] = c2[ik - 1 + j * c2_dim1] - c2[ik + 		    jc * c2_dim1];	    ch2[ik - 1 + jc * ch2_dim1] = c2[ik - 1 + j * c2_dim1] + c2[ik + 		    jc * c2_dim1];	    ch2[ik + j * ch2_dim1] = c2[ik + j * c2_dim1] + c2[ik - 1 + jc * 		    c2_dim1];	    ch2[ik + jc * ch2_dim1] = c2[ik + j * c2_dim1] - c2[ik - 1 + jc * 		    c2_dim1];/* L119: */	}/* L120: */    }    *nac = 1;    if (*ido == 2) {	return;