PLaSK library
Loading...
Searching...
No Matches
cmf5kf.c
Go to the documentation of this file.
1/* cmf5kf.f -- translated by f2c (version 20100827).
2 You must link the resulting object file with libf2c:
3 on Microsoft Windows system, link with libf2c.lib;
4 on Linux or Unix systems, link with .../path/to/libf2c.a -lm
5 or, if you install libf2c.a in a standard place, with -lf2c -lm
6 -- in that order, at the end of the command line, as in
7 cc *.o -lf2c -lm
8 Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
9
10 http://www.netlib.org/f2c/libf2c.zip
11*/
12
13#include "f2c.h"
14
15/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
16/* * * */
17/* * copyright (c) 2011 by UCAR * */
18/* * * */
19/* * University Corporation for Atmospheric Research * */
20/* * * */
21/* * all rights reserved * */
22/* * * */
23/* * FFTPACK version 5.1 * */
24/* * * */
25/* * A Fortran Package of Fast Fourier * */
26/* * * */
27/* * Subroutines and Example Programs * */
28/* * * */
29/* * by * */
30/* * * */
31/* * Paul Swarztrauber and Dick Valent * */
32/* * * */
33/* * of * */
34/* * * */
35/* * the National Center for Atmospheric Research * */
36/* * * */
37/* * Boulder, Colorado (80307) U.S.A. * */
38/* * * */
39/* * which is sponsored by * */
40/* * * */
41/* * the National Science Foundation * */
42/* * * */
43/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
44
45/* Subroutine */ int cmf5kf_(integer *lot, integer *ido, integer *l1, integer
46 *na, doublereal *cc, integer *im1, integer *in1, doublereal *ch,
47 integer *im2, integer *in2, doublereal *wa)
48{
49 /* Initialized data */
50
51 static doublereal tr11 = .3090169943749474;
52 static doublereal ti11 = -.9510565162951536;
53 static doublereal tr12 = -.8090169943749474;
54 static doublereal ti12 = -.5877852522924731;
55
56 /* System generated locals */
57 integer cc_dim2, cc_dim3, cc_dim4, cc_offset, ch_dim2, ch_dim3, ch_offset,
58 wa_dim1, wa_offset, i__1, i__2, i__3, i__4;
59
60 /* Local variables */
61 integer i__, k, m1, m2;
62 doublereal sn, ci2, ci3, ci4, ci5;
63 integer m1d;
64 doublereal di3, di4, di5, di2, cr2, cr3, cr5, cr4, ti2, ti3, ti4;
65 integer m2s;
66 doublereal ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5, chold1, chold2;
67
68 /* Parameter adjustments */
69 wa_dim1 = *ido;
70 wa_offset = 1 + wa_dim1 * 5;
71 wa -= wa_offset;
72 cc_dim2 = *in1;
73 cc_dim3 = *l1;
74 cc_dim4 = *ido;
75 cc_offset = 1 + 2 * (1 + cc_dim2 * (1 + cc_dim3 * (1 + cc_dim4)));
76 cc -= cc_offset;
77 ch_dim2 = *in2;
78 ch_dim3 = *l1;
79 ch_offset = 1 + 2 * (1 + ch_dim2 * (1 + ch_dim3 * 6));
80 ch -= ch_offset;
81
82 /* Function Body */
83
84/* FFTPACK 5.0 auxiliary routine */
85
86 m1d = (*lot - 1) * *im1 + 1;
87 m2s = 1 - *im2;
88 if (*ido > 1) {
89 goto L102;
90 }
91 sn = 1. / (doublereal) (*l1 * 5);
92 if (*na == 1) {
93 goto L106;
94 }
95 i__1 = *l1;
96 for (k = 1; k <= i__1; ++k) {
97 i__2 = m1d;
98 i__3 = *im1;
99 for (m1 = 1; i__3 < 0 ? m1 >= i__2 : m1 <= i__2; m1 += i__3) {
100 ti5 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
101 1) + 2] - cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
102 cc_dim2 << 1) + 2];
103 ti2 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
104 1) + 2] + cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
105 cc_dim2 << 1) + 2];
106 ti4 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
107 + 2] - cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
108 cc_dim2 << 1) + 2];
109 ti3 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
110 + 2] + cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
111 cc_dim2 << 1) + 2];
112 tr5 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
113 1) + 1] - cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
114 cc_dim2 << 1) + 1];
115 tr2 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
116 1) + 1] + cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
117 cc_dim2 << 1) + 1];
118 tr4 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
119 + 1] - cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
120 cc_dim2 << 1) + 1];
121 tr3 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
122 + 1] + cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
123 cc_dim2 << 1) + 1];
124 chold1 = sn * (cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 <<
125 1) + 1] + tr2 + tr3);
126 chold2 = sn * (cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 <<
127 1) + 2] + ti2 + ti3);
128 cr2 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1]
129 + tr11 * tr2 + tr12 * tr3;
130 ci2 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2]
131 + tr11 * ti2 + tr12 * ti3;
132 cr3 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1]
133 + tr12 * tr2 + tr11 * tr3;
134 ci3 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2]
135 + tr12 * ti2 + tr11 * ti3;
136 cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1] =
137 chold1;
138 cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2] =
139 chold2;
140 cr5 = ti11 * tr5 + ti12 * tr4;
141 ci5 = ti11 * ti5 + ti12 * ti4;
142 cr4 = ti12 * tr5 - ti11 * tr4;
143 ci4 = ti12 * ti5 - ti11 * ti4;
144 cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 << 1) + 1]
145 = sn * (cr2 - ci5);
146 cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) * cc_dim2 << 1) + 1] =
147 sn * (cr2 + ci5);
148 cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 << 1) + 2]
149 = sn * (ci2 + cr5);
150 cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1) + 2] =
151 sn * (ci3 + cr4);
152 cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1) + 1] =
153 sn * (cr3 - ci4);
154 cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) * cc_dim2 << 1) + 1]
155 = sn * (cr3 + ci4);
156 cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) * cc_dim2 << 1) + 2]
157 = sn * (ci3 - cr4);
158 cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) * cc_dim2 << 1) + 2] =
159 sn * (ci2 - cr5);
160/* L101: */
161 }
162 }
163 return 0;
164L106:
165 i__3 = *l1;
166 for (k = 1; k <= i__3; ++k) {
167 m2 = m2s;
168 i__2 = m1d;
169 i__1 = *im1;
170 for (m1 = 1; i__1 < 0 ? m1 >= i__2 : m1 <= i__2; m1 += i__1) {
171 m2 += *im2;
172 ti5 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
173 1) + 2] - cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
174 cc_dim2 << 1) + 2];
175 ti2 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
176 1) + 2] + cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
177 cc_dim2 << 1) + 2];
178 ti4 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
179 + 2] - cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
180 cc_dim2 << 1) + 2];
181 ti3 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
182 + 2] + cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
183 cc_dim2 << 1) + 2];
184 tr5 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
185 1) + 1] - cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
186 cc_dim2 << 1) + 1];
187 tr2 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
188 1) + 1] + cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
189 cc_dim2 << 1) + 1];
190 tr4 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
191 + 1] - cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
192 cc_dim2 << 1) + 1];
193 tr3 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
194 + 1] + cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
195 cc_dim2 << 1) + 1];
196 ch[(m2 + (k + ch_dim3 * 6) * ch_dim2 << 1) + 1] = sn * (cc[(m1 + (
197 k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1] + tr2 +
198 tr3);
199 ch[(m2 + (k + ch_dim3 * 6) * ch_dim2 << 1) + 2] = sn * (cc[(m1 + (
200 k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2] + ti2 +
201 ti3);
202 cr2 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1]
203 + tr11 * tr2 + tr12 * tr3;
204 ci2 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2]
205 + tr11 * ti2 + tr12 * ti3;
206 cr3 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1]
207 + tr12 * tr2 + tr11 * tr3;
208 ci3 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2]
209 + tr12 * ti2 + tr11 * ti3;
210 cr5 = ti11 * tr5 + ti12 * tr4;
211 ci5 = ti11 * ti5 + ti12 * ti4;
212 cr4 = ti12 * tr5 - ti11 * tr4;
213 ci4 = ti12 * ti5 - ti11 * ti4;
214 ch[(m2 + (k + ch_dim3 * 7) * ch_dim2 << 1) + 1] = sn * (cr2 - ci5)
215 ;
216 ch[(m2 + (k + ch_dim3 * 10) * ch_dim2 << 1) + 1] = sn * (cr2 +
217 ci5);
218 ch[(m2 + (k + ch_dim3 * 7) * ch_dim2 << 1) + 2] = sn * (ci2 + cr5)
219 ;
220 ch[(m2 + (k + (ch_dim3 << 3)) * ch_dim2 << 1) + 2] = sn * (ci3 +
221 cr4);
222 ch[(m2 + (k + (ch_dim3 << 3)) * ch_dim2 << 1) + 1] = sn * (cr3 -
223 ci4);
224 ch[(m2 + (k + ch_dim3 * 9) * ch_dim2 << 1) + 1] = sn * (cr3 + ci4)
225 ;
226 ch[(m2 + (k + ch_dim3 * 9) * ch_dim2 << 1) + 2] = sn * (ci3 - cr4)
227 ;
228 ch[(m2 + (k + ch_dim3 * 10) * ch_dim2 << 1) + 2] = sn * (ci2 -
229 cr5);
230/* L107: */
231 }
232 }
233 return 0;
234L102:
235 i__1 = *l1;
236 for (k = 1; k <= i__1; ++k) {
237 m2 = m2s;
238 i__2 = m1d;
239 i__3 = *im1;
240 for (m1 = 1; i__3 < 0 ? m1 >= i__2 : m1 <= i__2; m1 += i__3) {
241 m2 += *im2;
242 ti5 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
243 1) + 2] - cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
244 cc_dim2 << 1) + 2];
245 ti2 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
246 1) + 2] + cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
247 cc_dim2 << 1) + 2];
248 ti4 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
249 + 2] - cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
250 cc_dim2 << 1) + 2];
251 ti3 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
252 + 2] + cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
253 cc_dim2 << 1) + 2];
254 tr5 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
255 1) + 1] - cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
256 cc_dim2 << 1) + 1];
257 tr2 = cc[(m1 + (k + ((cc_dim4 << 1) + 1) * cc_dim3) * cc_dim2 <<
258 1) + 1] + cc[(m1 + (k + (cc_dim4 * 5 + 1) * cc_dim3) *
259 cc_dim2 << 1) + 1];
260 tr4 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
261 + 1] - cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
262 cc_dim2 << 1) + 1];
263 tr3 = cc[(m1 + (k + (cc_dim4 * 3 + 1) * cc_dim3) * cc_dim2 << 1)
264 + 1] + cc[(m1 + (k + ((cc_dim4 << 2) + 1) * cc_dim3) *
265 cc_dim2 << 1) + 1];
266 ch[(m2 + (k + ch_dim3 * 6) * ch_dim2 << 1) + 1] = cc[(m1 + (k + (
267 cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1] + tr2 + tr3;
268 ch[(m2 + (k + ch_dim3 * 6) * ch_dim2 << 1) + 2] = cc[(m1 + (k + (
269 cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2] + ti2 + ti3;
270 cr2 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1]
271 + tr11 * tr2 + tr12 * tr3;
272 ci2 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2]
273 + tr11 * ti2 + tr12 * ti3;
274 cr3 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 1]
275 + tr12 * tr2 + tr11 * tr3;
276 ci3 = cc[(m1 + (k + (cc_dim4 + 1) * cc_dim3) * cc_dim2 << 1) + 2]
277 + tr12 * ti2 + tr11 * ti3;
278 cr5 = ti11 * tr5 + ti12 * tr4;
279 ci5 = ti11 * ti5 + ti12 * ti4;
280 cr4 = ti12 * tr5 - ti11 * tr4;
281 ci4 = ti12 * ti5 - ti11 * ti4;
282 ch[(m2 + (k + ch_dim3 * 7) * ch_dim2 << 1) + 1] = cr2 - ci5;
283 ch[(m2 + (k + ch_dim3 * 10) * ch_dim2 << 1) + 1] = cr2 + ci5;
284 ch[(m2 + (k + ch_dim3 * 7) * ch_dim2 << 1) + 2] = ci2 + cr5;
285 ch[(m2 + (k + (ch_dim3 << 3)) * ch_dim2 << 1) + 2] = ci3 + cr4;
286 ch[(m2 + (k + (ch_dim3 << 3)) * ch_dim2 << 1) + 1] = cr3 - ci4;
287 ch[(m2 + (k + ch_dim3 * 9) * ch_dim2 << 1) + 1] = cr3 + ci4;
288 ch[(m2 + (k + ch_dim3 * 9) * ch_dim2 << 1) + 2] = ci3 - cr4;
289 ch[(m2 + (k + ch_dim3 * 10) * ch_dim2 << 1) + 2] = ci2 - cr5;
290/* L103: */
291 }
292 }
293 i__3 = *ido;
294 for (i__ = 2; i__ <= i__3; ++i__) {
295 i__2 = *l1;
296 for (k = 1; k <= i__2; ++k) {
297 m2 = m2s;
298 i__1 = m1d;
299 i__4 = *im1;
300 for (m1 = 1; i__4 < 0 ? m1 >= i__1 : m1 <= i__1; m1 += i__4) {
301 m2 += *im2;
302 ti5 = cc[(m1 + (k + (i__ + (cc_dim4 << 1)) * cc_dim3) *
303 cc_dim2 << 1) + 2] - cc[(m1 + (k + (i__ + cc_dim4 * 5)
304 * cc_dim3) * cc_dim2 << 1) + 2];
305 ti2 = cc[(m1 + (k + (i__ + (cc_dim4 << 1)) * cc_dim3) *
306 cc_dim2 << 1) + 2] + cc[(m1 + (k + (i__ + cc_dim4 * 5)
307 * cc_dim3) * cc_dim2 << 1) + 2];
308 ti4 = cc[(m1 + (k + (i__ + cc_dim4 * 3) * cc_dim3) * cc_dim2
309 << 1) + 2] - cc[(m1 + (k + (i__ + (cc_dim4 << 2)) *
310 cc_dim3) * cc_dim2 << 1) + 2];
311 ti3 = cc[(m1 + (k + (i__ + cc_dim4 * 3) * cc_dim3) * cc_dim2
312 << 1) + 2] + cc[(m1 + (k + (i__ + (cc_dim4 << 2)) *
313 cc_dim3) * cc_dim2 << 1) + 2];
314 tr5 = cc[(m1 + (k + (i__ + (cc_dim4 << 1)) * cc_dim3) *
315 cc_dim2 << 1) + 1] - cc[(m1 + (k + (i__ + cc_dim4 * 5)
316 * cc_dim3) * cc_dim2 << 1) + 1];
317 tr2 = cc[(m1 + (k + (i__ + (cc_dim4 << 1)) * cc_dim3) *
318 cc_dim2 << 1) + 1] + cc[(m1 + (k + (i__ + cc_dim4 * 5)
319 * cc_dim3) * cc_dim2 << 1) + 1];
320 tr4 = cc[(m1 + (k + (i__ + cc_dim4 * 3) * cc_dim3) * cc_dim2
321 << 1) + 1] - cc[(m1 + (k + (i__ + (cc_dim4 << 2)) *
322 cc_dim3) * cc_dim2 << 1) + 1];
323 tr3 = cc[(m1 + (k + (i__ + cc_dim4 * 3) * cc_dim3) * cc_dim2
324 << 1) + 1] + cc[(m1 + (k + (i__ + (cc_dim4 << 2)) *
325 cc_dim3) * cc_dim2 << 1) + 1];
326 ch[(m2 + (k + (i__ * 5 + 1) * ch_dim3) * ch_dim2 << 1) + 1] =
327 cc[(m1 + (k + (i__ + cc_dim4) * cc_dim3) * cc_dim2 <<
328 1) + 1] + tr2 + tr3;
329 ch[(m2 + (k + (i__ * 5 + 1) * ch_dim3) * ch_dim2 << 1) + 2] =
330 cc[(m1 + (k + (i__ + cc_dim4) * cc_dim3) * cc_dim2 <<
331 1) + 2] + ti2 + ti3;
332 cr2 = cc[(m1 + (k + (i__ + cc_dim4) * cc_dim3) * cc_dim2 << 1)
333 + 1] + tr11 * tr2 + tr12 * tr3;
334 ci2 = cc[(m1 + (k + (i__ + cc_dim4) * cc_dim3) * cc_dim2 << 1)
335 + 2] + tr11 * ti2 + tr12 * ti3;
336 cr3 = cc[(m1 + (k + (i__ + cc_dim4) * cc_dim3) * cc_dim2 << 1)
337 + 1] + tr12 * tr2 + tr11 * tr3;
338 ci3 = cc[(m1 + (k + (i__ + cc_dim4) * cc_dim3) * cc_dim2 << 1)
339 + 2] + tr12 * ti2 + tr11 * ti3;
340 cr5 = ti11 * tr5 + ti12 * tr4;
341 ci5 = ti11 * ti5 + ti12 * ti4;
342 cr4 = ti12 * tr5 - ti11 * tr4;
343 ci4 = ti12 * ti5 - ti11 * ti4;
344 dr3 = cr3 - ci4;
345 dr4 = cr3 + ci4;
346 di3 = ci3 + cr4;
347 di4 = ci3 - cr4;
348 dr5 = cr2 + ci5;
349 dr2 = cr2 - ci5;
350 di5 = ci2 - cr5;
351 di2 = ci2 + cr5;
352 ch[(m2 + (k + (i__ * 5 + 2) * ch_dim3) * ch_dim2 << 1) + 1] =
353 wa[i__ + wa_dim1 * 5] * dr2 + wa[i__ + wa_dim1 * 9] *
354 di2;
355 ch[(m2 + (k + (i__ * 5 + 2) * ch_dim3) * ch_dim2 << 1) + 2] =
356 wa[i__ + wa_dim1 * 5] * di2 - wa[i__ + wa_dim1 * 9] *
357 dr2;
358 ch[(m2 + (k + (i__ * 5 + 3) * ch_dim3) * ch_dim2 << 1) + 1] =
359 wa[i__ + wa_dim1 * 6] * dr3 + wa[i__ + wa_dim1 * 10] *
360 di3;
361 ch[(m2 + (k + (i__ * 5 + 3) * ch_dim3) * ch_dim2 << 1) + 2] =
362 wa[i__ + wa_dim1 * 6] * di3 - wa[i__ + wa_dim1 * 10] *
363 dr3;
364 ch[(m2 + (k + (i__ * 5 + 4) * ch_dim3) * ch_dim2 << 1) + 1] =
365 wa[i__ + wa_dim1 * 7] * dr4 + wa[i__ + wa_dim1 * 11] *
366 di4;
367 ch[(m2 + (k + (i__ * 5 + 4) * ch_dim3) * ch_dim2 << 1) + 2] =
368 wa[i__ + wa_dim1 * 7] * di4 - wa[i__ + wa_dim1 * 11] *
369 dr4;
370 ch[(m2 + (k + (i__ * 5 + 5) * ch_dim3) * ch_dim2 << 1) + 1] =
371 wa[i__ + (wa_dim1 << 3)] * dr5 + wa[i__ + wa_dim1 *
372 12] * di5;
373 ch[(m2 + (k + (i__ * 5 + 5) * ch_dim3) * ch_dim2 << 1) + 2] =
374 wa[i__ + (wa_dim1 << 3)] * di5 - wa[i__ + wa_dim1 *
375 12] * dr5;
376/* L104: */
377 }
378 }
379/* L105: */
380 }
381 return 0;
382} /* cmf5kf_ */
383