1145 static double tr11 = .309016994374947;
1146 static double ti11 = -.951056516295154;
1147 static double tr12 = -.809016994374947;
1148 static double ti12 = -.587785252292473;
1151 int cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i_1, i_2;
1155 static double ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5,
1156 cr4, ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;
1160 cc_offset = cc_dim1 * 6 + 1;
1164 ch_offset = ch_dim1 * (ch_dim2 + 1) + 1;
1176 for (k = 1; k <= i_1; ++k) {
1177 ti5 = cc[(k * 5 + 2) * cc_dim1 + 2] - cc[(k * 5 + 5) * cc_dim1 + 2];
1178 ti2 = cc[(k * 5 + 2) * cc_dim1 + 2] + cc[(k * 5 + 5) * cc_dim1 + 2];
1179 ti4 = cc[(k * 5 + 3) * cc_dim1 + 2] - cc[(k * 5 + 4) * cc_dim1 + 2];
1180 ti3 = cc[(k * 5 + 3) * cc_dim1 + 2] + cc[(k * 5 + 4) * cc_dim1 + 2];
1181 tr5 = cc[(k * 5 + 2) * cc_dim1 + 1] - cc[(k * 5 + 5) * cc_dim1 + 1];
1182 tr2 = cc[(k * 5 + 2) * cc_dim1 + 1] + cc[(k * 5 + 5) * cc_dim1 + 1];
1183 tr4 = cc[(k * 5 + 3) * cc_dim1 + 1] - cc[(k * 5 + 4) * cc_dim1 + 1];
1184 tr3 = cc[(k * 5 + 3) * cc_dim1 + 1] + cc[(k * 5 + 4) * cc_dim1 + 1];
1185 ch[(k + ch_dim2) * ch_dim1 + 1] = cc[(k * 5 + 1) * cc_dim1 + 1] + tr2
1187 ch[(k + ch_dim2) * ch_dim1 + 2] = cc[(k * 5 + 1) * cc_dim1 + 2] + ti2
1189 cr2 = cc[(k * 5 + 1) * cc_dim1 + 1] + tr11 * tr2 + tr12 * tr3;
1190 ci2 = cc[(k * 5 + 1) * cc_dim1 + 2] + tr11 * ti2 + tr12 * ti3;
1191 cr3 = cc[(k * 5 + 1) * cc_dim1 + 1] + tr12 * tr2 + tr11 * tr3;
1192 ci3 = cc[(k * 5 + 1) * cc_dim1 + 2] + tr12 * ti2 + tr11 * ti3;
1193 cr5 = ti11 * tr5 + ti12 * tr4;
1194 ci5 = ti11 * ti5 + ti12 * ti4;
1195 cr4 = ti12 * tr5 - ti11 * tr4;
1196 ci4 = ti12 * ti5 - ti11 * ti4;
1197 ch[(k + (ch_dim2 << 1)) * ch_dim1 + 1] = cr2 - ci5;
1198 ch[(k + ch_dim2 * 5) * ch_dim1 + 1] = cr2 + ci5;
1199 ch[(k + (ch_dim2 << 1)) * ch_dim1 + 2] = ci2 + cr5;
1200 ch[(k + ch_dim2 * 3) * ch_dim1 + 2] = ci3 + cr4;
1201 ch[(k + ch_dim2 * 3) * ch_dim1 + 1] = cr3 - ci4;
1202 ch[(k + (ch_dim2 << 2)) * ch_dim1 + 1] = cr3 + ci4;
1203 ch[(k + (ch_dim2 << 2)) * ch_dim1 + 2] = ci3 - cr4;
1204 ch[(k + ch_dim2 * 5) * ch_dim1 + 2] = ci2 - cr5;
1210 for (k = 1; k <= i_1; ++k) {
1212 for (i = 2; i <= i_2; i += 2) {
1213 ti5 = cc[i + (k * 5 + 2) * cc_dim1] - cc[i + (k * 5 + 5) *
1215 ti2 = cc[i + (k * 5 + 2) * cc_dim1] + cc[i + (k * 5 + 5) *
1217 ti4 = cc[i + (k * 5 + 3) * cc_dim1] - cc[i + (k * 5 + 4) *
1219 ti3 = cc[i + (k * 5 + 3) * cc_dim1] + cc[i + (k * 5 + 4) *
1221 tr5 = cc[i - 1 + (k * 5 + 2) * cc_dim1] - cc[i - 1 + (k * 5 + 5) *
1223 tr2 = cc[i - 1 + (k * 5 + 2) * cc_dim1] + cc[i - 1 + (k * 5 + 5) *
1225 tr4 = cc[i - 1 + (k * 5 + 3) * cc_dim1] - cc[i - 1 + (k * 5 + 4) *
1227 tr3 = cc[i - 1 + (k * 5 + 3) * cc_dim1] + cc[i - 1 + (k * 5 + 4) *
1229 ch[i - 1 + (k + ch_dim2) * ch_dim1] = cc[i - 1 + (k * 5 + 1) *
1230 cc_dim1] + tr2 + tr3;
1231 ch[i + (k + ch_dim2) * ch_dim1] = cc[i + (k * 5 + 1) * cc_dim1] +
1233 cr2 = cc[i - 1 + (k * 5 + 1) * cc_dim1] + tr11 * tr2 + tr12 * tr3;
1235 ci2 = cc[i + (k * 5 + 1) * cc_dim1] + tr11 * ti2 + tr12 * ti3;
1236 cr3 = cc[i - 1 + (k * 5 + 1) * cc_dim1] + tr12 * tr2 + tr11 * tr3;
1238 ci3 = cc[i + (k * 5 + 1) * cc_dim1] + tr12 * ti2 + tr11 * ti3;
1239 cr5 = ti11 * tr5 + ti12 * tr4;
1240 ci5 = ti11 * ti5 + ti12 * ti4;
1241 cr4 = ti12 * tr5 - ti11 * tr4;
1242 ci4 = ti12 * ti5 - ti11 * ti4;
1251 ch[i - 1 + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i - 1] * dr2 +
1253 ch[i + (k + (ch_dim2 << 1)) * ch_dim1] = wa1[i - 1] * di2 - wa1[i]
1255 ch[i - 1 + (k + ch_dim2 * 3) * ch_dim1] = wa2[i - 1] * dr3 + wa2[
1257 ch[i + (k + ch_dim2 * 3) * ch_dim1] = wa2[i - 1] * di3 - wa2[i] *
1259 ch[i - 1 + (k + (ch_dim2 << 2)) * ch_dim1] = wa3[i - 1] * dr4 +
1261 ch[i + (k + (ch_dim2 << 2)) * ch_dim1] = wa3[i - 1] * di4 - wa3[i]
1263 ch[i - 1 + (k + ch_dim2 * 5) * ch_dim1] = wa4[i - 1] * dr5 + wa4[
1265 ch[i + (k + ch_dim2 * 5) * ch_dim1] = wa4[i - 1] * di5 - wa4[i] *