PLaSK library
Loading...
Searching...
No Matches
mradf5.c
Go to the documentation of this file.
1/* mradf5.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 mradf5_(integer *m, integer *ido, integer *l1,
46 doublereal *cc, integer *im1, integer *in1, doublereal *ch, integer *
47 im2, integer *in2, doublereal *wa1, doublereal *wa2, doublereal *wa3,
48 doublereal *wa4)
49{
50 /* System generated locals */
51 integer cc_dim1, cc_dim2, cc_dim3, cc_offset, ch_dim1, ch_dim2, ch_offset,
52 i__1, i__2, i__3, i__4;
53
54 /* Builtin functions */
55 double atan(doublereal), cos(doublereal), sin(doublereal);
56
57 /* Local variables */
58 integer i__, k, m1, m2, ic, m1d, m2s;
59 doublereal arg, ti11, ti12, tr11, tr12;
60 integer idp2;
61
62
63 /* Parameter adjustments */
64 --wa4;
65 --wa3;
66 --wa2;
67 --wa1;
68 cc_dim1 = *in1;
69 cc_dim2 = *ido;
70 cc_dim3 = *l1;
71 cc_offset = 1 + cc_dim1 * (1 + cc_dim2 * (1 + cc_dim3));
72 cc -= cc_offset;
73 ch_dim1 = *in2;
74 ch_dim2 = *ido;
75 ch_offset = 1 + ch_dim1 * (1 + ch_dim2 * 6);
76 ch -= ch_offset;
77
78 /* Function Body */
79 m1d = (*m - 1) * *im1 + 1;
80 m2s = 1 - *im2;
81 arg = atan(1.) * 8. / 5.;
82 tr11 = cos(arg);
83 ti11 = sin(arg);
84 tr12 = cos(arg * 2.);
85 ti12 = sin(arg * 2.);
86 i__1 = *l1;
87 for (k = 1; k <= i__1; ++k) {
88 m2 = m2s;
89 i__2 = m1d;
90 i__3 = *im1;
91 for (m1 = 1; i__3 < 0 ? m1 >= i__2 : m1 <= i__2; m1 += i__3) {
92 m2 += *im2;
93 ch[m2 + ((k * 5 + 1) * ch_dim2 + 1) * ch_dim1] = cc[m1 + ((k +
94 cc_dim3) * cc_dim2 + 1) * cc_dim1] + (cc[m1 + ((k +
95 cc_dim3 * 5) * cc_dim2 + 1) * cc_dim1] + cc[m1 + ((k + (
96 cc_dim3 << 1)) * cc_dim2 + 1) * cc_dim1]) + (cc[m1 + ((k
97 + (cc_dim3 << 2)) * cc_dim2 + 1) * cc_dim1] + cc[m1 + ((k
98 + cc_dim3 * 3) * cc_dim2 + 1) * cc_dim1]);
99 ch[m2 + (*ido + (k * 5 + 2) * ch_dim2) * ch_dim1] = cc[m1 + ((k +
100 cc_dim3) * cc_dim2 + 1) * cc_dim1] + tr11 * (cc[m1 + ((k
101 + cc_dim3 * 5) * cc_dim2 + 1) * cc_dim1] + cc[m1 + ((k + (
102 cc_dim3 << 1)) * cc_dim2 + 1) * cc_dim1]) + tr12 * (cc[m1
103 + ((k + (cc_dim3 << 2)) * cc_dim2 + 1) * cc_dim1] + cc[m1
104 + ((k + cc_dim3 * 3) * cc_dim2 + 1) * cc_dim1]);
105 ch[m2 + ((k * 5 + 3) * ch_dim2 + 1) * ch_dim1] = ti11 * (cc[m1 + (
106 (k + cc_dim3 * 5) * cc_dim2 + 1) * cc_dim1] - cc[m1 + ((k
107 + (cc_dim3 << 1)) * cc_dim2 + 1) * cc_dim1]) + ti12 * (cc[
108 m1 + ((k + (cc_dim3 << 2)) * cc_dim2 + 1) * cc_dim1] - cc[
109 m1 + ((k + cc_dim3 * 3) * cc_dim2 + 1) * cc_dim1]);
110 ch[m2 + (*ido + (k * 5 + 4) * ch_dim2) * ch_dim1] = cc[m1 + ((k +
111 cc_dim3) * cc_dim2 + 1) * cc_dim1] + tr12 * (cc[m1 + ((k
112 + cc_dim3 * 5) * cc_dim2 + 1) * cc_dim1] + cc[m1 + ((k + (
113 cc_dim3 << 1)) * cc_dim2 + 1) * cc_dim1]) + tr11 * (cc[m1
114 + ((k + (cc_dim3 << 2)) * cc_dim2 + 1) * cc_dim1] + cc[m1
115 + ((k + cc_dim3 * 3) * cc_dim2 + 1) * cc_dim1]);
116 ch[m2 + ((k * 5 + 5) * ch_dim2 + 1) * ch_dim1] = ti12 * (cc[m1 + (
117 (k + cc_dim3 * 5) * cc_dim2 + 1) * cc_dim1] - cc[m1 + ((k
118 + (cc_dim3 << 1)) * cc_dim2 + 1) * cc_dim1]) - ti11 * (cc[
119 m1 + ((k + (cc_dim3 << 2)) * cc_dim2 + 1) * cc_dim1] - cc[
120 m1 + ((k + cc_dim3 * 3) * cc_dim2 + 1) * cc_dim1]);
121/* L1001: */
122 }
123/* L101: */
124 }
125 if (*ido == 1) {
126 return 0;
127 }
128 idp2 = *ido + 2;
129 i__1 = *l1;
130 for (k = 1; k <= i__1; ++k) {
131 i__3 = *ido;
132 for (i__ = 3; i__ <= i__3; i__ += 2) {
133 ic = idp2 - i__;
134 m2 = m2s;
135 i__2 = m1d;
136 i__4 = *im1;
137 for (m1 = 1; i__4 < 0 ? m1 >= i__2 : m1 <= i__2; m1 += i__4) {
138 m2 += *im2;
139 ch[m2 + (i__ - 1 + (k * 5 + 1) * ch_dim2) * ch_dim1] = cc[m1
140 + (i__ - 1 + (k + cc_dim3) * cc_dim2) * cc_dim1] + (
141 wa1[i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 1))
142 * cc_dim2) * cc_dim1] + wa1[i__ - 1] * cc[m1 + (i__
143 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] + (wa4[
144 i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 5) *
145 cc_dim2) * cc_dim1] + wa4[i__ - 1] * cc[m1 + (i__ + (
146 k + cc_dim3 * 5) * cc_dim2) * cc_dim1])) + (wa2[i__ -
147 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 3) * cc_dim2) *
148 cc_dim1] + wa2[i__ - 1] * cc[m1 + (i__ + (k +
149 cc_dim3 * 3) * cc_dim2) * cc_dim1] + (wa3[i__ - 2] *
150 cc[m1 + (i__ - 1 + (k + (cc_dim3 << 2)) * cc_dim2) *
151 cc_dim1] + wa3[i__ - 1] * cc[m1 + (i__ + (k + (
152 cc_dim3 << 2)) * cc_dim2) * cc_dim1]));
153 ch[m2 + (i__ + (k * 5 + 1) * ch_dim2) * ch_dim1] = cc[m1 + (
154 i__ + (k + cc_dim3) * cc_dim2) * cc_dim1] + (wa1[i__
155 - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 1)) * cc_dim2)
156 * cc_dim1] - wa1[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
157 cc_dim3 << 1)) * cc_dim2) * cc_dim1] + (wa4[i__ - 2] *
158 cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2) *
159 cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1 + (k +
160 cc_dim3 * 5) * cc_dim2) * cc_dim1])) + (wa2[i__ - 2] *
161 cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2) *
162 cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1 + (k +
163 cc_dim3 * 3) * cc_dim2) * cc_dim1] + (wa3[i__ - 2] *
164 cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2) *
165 cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
166 cc_dim3 << 2)) * cc_dim2) * cc_dim1]));
167 ch[m2 + (i__ - 1 + (k * 5 + 3) * ch_dim2) * ch_dim1] = cc[m1
168 + (i__ - 1 + (k + cc_dim3) * cc_dim2) * cc_dim1] +
169 tr11 * (wa1[i__ - 2] * cc[m1 + (i__ - 1 + (k + (
170 cc_dim3 << 1)) * cc_dim2) * cc_dim1] + wa1[i__ - 1] *
171 cc[m1 + (i__ + (k + (cc_dim3 << 1)) * cc_dim2) *
172 cc_dim1] + wa4[i__ - 2] * cc[m1 + (i__ - 1 + (k +
173 cc_dim3 * 5) * cc_dim2) * cc_dim1] + wa4[i__ - 1] *
174 cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2) * cc_dim1]
175 ) + tr12 * (wa2[i__ - 2] * cc[m1 + (i__ - 1 + (k +
176 cc_dim3 * 3) * cc_dim2) * cc_dim1] + wa2[i__ - 1] *
177 cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2) * cc_dim1]
178 + wa3[i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 <<
179 2)) * cc_dim2) * cc_dim1] + wa3[i__ - 1] * cc[m1 + (
180 i__ + (k + (cc_dim3 << 2)) * cc_dim2) * cc_dim1]) +
181 ti11 * (wa1[i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 <<
182 1)) * cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (
183 i__ - 1 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1]
184 - (wa4[i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) *
185 cc_dim2) * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1
186 + (k + cc_dim3 * 5) * cc_dim2) * cc_dim1])) + ti12 * (
187 wa2[i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) *
188 cc_dim2) * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1
189 + (k + cc_dim3 * 3) * cc_dim2) * cc_dim1] - (wa3[i__
190 - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2)
191 * cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
192 cc_dim3 << 2)) * cc_dim2) * cc_dim1]));
193 ch[m2 + (ic - 1 + (k * 5 + 2) * ch_dim2) * ch_dim1] = cc[m1 +
194 (i__ - 1 + (k + cc_dim3) * cc_dim2) * cc_dim1] + tr11
195 * (wa1[i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 <<
196 1)) * cc_dim2) * cc_dim1] + wa1[i__ - 1] * cc[m1 + (
197 i__ + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] +
198 wa4[i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 5) *
199 cc_dim2) * cc_dim1] + wa4[i__ - 1] * cc[m1 + (i__ + (
200 k + cc_dim3 * 5) * cc_dim2) * cc_dim1]) + tr12 * (wa2[
201 i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 3) *
202 cc_dim2) * cc_dim1] + wa2[i__ - 1] * cc[m1 + (i__ + (
203 k + cc_dim3 * 3) * cc_dim2) * cc_dim1] + wa3[i__ - 2]
204 * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 2)) * cc_dim2) *
205 cc_dim1] + wa3[i__ - 1] * cc[m1 + (i__ + (k + (
206 cc_dim3 << 2)) * cc_dim2) * cc_dim1]) - (ti11 * (wa1[
207 i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 1)) *
208 cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (i__ - 1
209 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] - (wa4[
210 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2)
211 * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1 + (k +
212 cc_dim3 * 5) * cc_dim2) * cc_dim1])) + ti12 * (wa2[
213 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2)
214 * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1 + (k +
215 cc_dim3 * 3) * cc_dim2) * cc_dim1] - (wa3[i__ - 2] *
216 cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2) *
217 cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
218 cc_dim3 << 2)) * cc_dim2) * cc_dim1])));
219 ch[m2 + (i__ + (k * 5 + 3) * ch_dim2) * ch_dim1] = cc[m1 + (
220 i__ + (k + cc_dim3) * cc_dim2) * cc_dim1] + tr11 * (
221 wa1[i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 1)) *
222 cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (i__ - 1
223 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] + (wa4[
224 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2)
225 * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1 + (k +
226 cc_dim3 * 5) * cc_dim2) * cc_dim1])) + tr12 * (wa2[
227 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2)
228 * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1 + (k +
229 cc_dim3 * 3) * cc_dim2) * cc_dim1] + (wa3[i__ - 2] *
230 cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2) *
231 cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
232 cc_dim3 << 2)) * cc_dim2) * cc_dim1])) + (ti11 * (wa4[
233 i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 5) *
234 cc_dim2) * cc_dim1] + wa4[i__ - 1] * cc[m1 + (i__ + (
235 k + cc_dim3 * 5) * cc_dim2) * cc_dim1] - (wa1[i__ - 2]
236 * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 1)) * cc_dim2)
237 * cc_dim1] + wa1[i__ - 1] * cc[m1 + (i__ + (k + (
238 cc_dim3 << 1)) * cc_dim2) * cc_dim1])) + ti12 * (wa3[
239 i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 2)) *
240 cc_dim2) * cc_dim1] + wa3[i__ - 1] * cc[m1 + (i__ + (
241 k + (cc_dim3 << 2)) * cc_dim2) * cc_dim1] - (wa2[i__
242 - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 3) * cc_dim2)
243 * cc_dim1] + wa2[i__ - 1] * cc[m1 + (i__ + (k +
244 cc_dim3 * 3) * cc_dim2) * cc_dim1])));
245 ch[m2 + (ic + (k * 5 + 2) * ch_dim2) * ch_dim1] = ti11 * (wa4[
246 i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 5) *
247 cc_dim2) * cc_dim1] + wa4[i__ - 1] * cc[m1 + (i__ + (
248 k + cc_dim3 * 5) * cc_dim2) * cc_dim1] - (wa1[i__ - 2]
249 * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 1)) * cc_dim2)
250 * cc_dim1] + wa1[i__ - 1] * cc[m1 + (i__ + (k + (
251 cc_dim3 << 1)) * cc_dim2) * cc_dim1])) + ti12 * (wa3[
252 i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 2)) *
253 cc_dim2) * cc_dim1] + wa3[i__ - 1] * cc[m1 + (i__ + (
254 k + (cc_dim3 << 2)) * cc_dim2) * cc_dim1] - (wa2[i__
255 - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 3) * cc_dim2)
256 * cc_dim1] + wa2[i__ - 1] * cc[m1 + (i__ + (k +
257 cc_dim3 * 3) * cc_dim2) * cc_dim1])) - (cc[m1 + (i__
258 + (k + cc_dim3) * cc_dim2) * cc_dim1] + tr11 * (wa1[
259 i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 1)) *
260 cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (i__ - 1
261 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] + (wa4[
262 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2)
263 * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1 + (k +
264 cc_dim3 * 5) * cc_dim2) * cc_dim1])) + tr12 * (wa2[
265 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2)
266 * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1 + (k +
267 cc_dim3 * 3) * cc_dim2) * cc_dim1] + (wa3[i__ - 2] *
268 cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2) *
269 cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
270 cc_dim3 << 2)) * cc_dim2) * cc_dim1])));
271 ch[m2 + (i__ - 1 + (k * 5 + 5) * ch_dim2) * ch_dim1] = cc[m1
272 + (i__ - 1 + (k + cc_dim3) * cc_dim2) * cc_dim1] +
273 tr12 * (wa1[i__ - 2] * cc[m1 + (i__ - 1 + (k + (
274 cc_dim3 << 1)) * cc_dim2) * cc_dim1] + wa1[i__ - 1] *
275 cc[m1 + (i__ + (k + (cc_dim3 << 1)) * cc_dim2) *
276 cc_dim1] + (wa4[i__ - 2] * cc[m1 + (i__ - 1 + (k +
277 cc_dim3 * 5) * cc_dim2) * cc_dim1] + wa4[i__ - 1] *
278 cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2) * cc_dim1]
279 )) + tr11 * (wa2[i__ - 2] * cc[m1 + (i__ - 1 + (k +
280 cc_dim3 * 3) * cc_dim2) * cc_dim1] + wa2[i__ - 1] *
281 cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2) * cc_dim1]
282 + (wa3[i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 <<
283 2)) * cc_dim2) * cc_dim1] + wa3[i__ - 1] * cc[m1 + (
284 i__ + (k + (cc_dim3 << 2)) * cc_dim2) * cc_dim1])) + (
285 ti12 * (wa1[i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 <<
286 1)) * cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (
287 i__ - 1 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1]
288 - (wa4[i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) *
289 cc_dim2) * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1
290 + (k + cc_dim3 * 5) * cc_dim2) * cc_dim1])) - ti11 * (
291 wa2[i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) *
292 cc_dim2) * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1
293 + (k + cc_dim3 * 3) * cc_dim2) * cc_dim1] - (wa3[i__
294 - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2)
295 * cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
296 cc_dim3 << 2)) * cc_dim2) * cc_dim1])));
297 ch[m2 + (ic - 1 + (k * 5 + 4) * ch_dim2) * ch_dim1] = cc[m1 +
298 (i__ - 1 + (k + cc_dim3) * cc_dim2) * cc_dim1] + tr12
299 * (wa1[i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 <<
300 1)) * cc_dim2) * cc_dim1] + wa1[i__ - 1] * cc[m1 + (
301 i__ + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] + (
302 wa4[i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 5) *
303 cc_dim2) * cc_dim1] + wa4[i__ - 1] * cc[m1 + (i__ + (
304 k + cc_dim3 * 5) * cc_dim2) * cc_dim1])) + tr11 * (
305 wa2[i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 3) *
306 cc_dim2) * cc_dim1] + wa2[i__ - 1] * cc[m1 + (i__ + (
307 k + cc_dim3 * 3) * cc_dim2) * cc_dim1] + (wa3[i__ - 2]
308 * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 2)) * cc_dim2)
309 * cc_dim1] + wa3[i__ - 1] * cc[m1 + (i__ + (k + (
310 cc_dim3 << 2)) * cc_dim2) * cc_dim1])) - (ti12 * (wa1[
311 i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 1)) *
312 cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (i__ - 1
313 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] - (wa4[
314 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2)
315 * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1 + (k +
316 cc_dim3 * 5) * cc_dim2) * cc_dim1])) - ti11 * (wa2[
317 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2)
318 * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1 + (k +
319 cc_dim3 * 3) * cc_dim2) * cc_dim1] - (wa3[i__ - 2] *
320 cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2) *
321 cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
322 cc_dim3 << 2)) * cc_dim2) * cc_dim1])));
323 ch[m2 + (i__ + (k * 5 + 5) * ch_dim2) * ch_dim1] = cc[m1 + (
324 i__ + (k + cc_dim3) * cc_dim2) * cc_dim1] + tr12 * (
325 wa1[i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 1)) *
326 cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (i__ - 1
327 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] + (wa4[
328 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2)
329 * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1 + (k +
330 cc_dim3 * 5) * cc_dim2) * cc_dim1])) + tr11 * (wa2[
331 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2)
332 * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1 + (k +
333 cc_dim3 * 3) * cc_dim2) * cc_dim1] + (wa3[i__ - 2] *
334 cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2) *
335 cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
336 cc_dim3 << 2)) * cc_dim2) * cc_dim1])) + (ti12 * (wa4[
337 i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 5) *
338 cc_dim2) * cc_dim1] + wa4[i__ - 1] * cc[m1 + (i__ + (
339 k + cc_dim3 * 5) * cc_dim2) * cc_dim1] - (wa1[i__ - 2]
340 * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 1)) * cc_dim2)
341 * cc_dim1] + wa1[i__ - 1] * cc[m1 + (i__ + (k + (
342 cc_dim3 << 1)) * cc_dim2) * cc_dim1])) - ti11 * (wa3[
343 i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 2)) *
344 cc_dim2) * cc_dim1] + wa3[i__ - 1] * cc[m1 + (i__ + (
345 k + (cc_dim3 << 2)) * cc_dim2) * cc_dim1] - (wa2[i__
346 - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 3) * cc_dim2)
347 * cc_dim1] + wa2[i__ - 1] * cc[m1 + (i__ + (k +
348 cc_dim3 * 3) * cc_dim2) * cc_dim1])));
349 ch[m2 + (ic + (k * 5 + 4) * ch_dim2) * ch_dim1] = ti12 * (wa4[
350 i__ - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 5) *
351 cc_dim2) * cc_dim1] + wa4[i__ - 1] * cc[m1 + (i__ + (
352 k + cc_dim3 * 5) * cc_dim2) * cc_dim1] - (wa1[i__ - 2]
353 * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 1)) * cc_dim2)
354 * cc_dim1] + wa1[i__ - 1] * cc[m1 + (i__ + (k + (
355 cc_dim3 << 1)) * cc_dim2) * cc_dim1])) - ti11 * (wa3[
356 i__ - 2] * cc[m1 + (i__ - 1 + (k + (cc_dim3 << 2)) *
357 cc_dim2) * cc_dim1] + wa3[i__ - 1] * cc[m1 + (i__ + (
358 k + (cc_dim3 << 2)) * cc_dim2) * cc_dim1] - (wa2[i__
359 - 2] * cc[m1 + (i__ - 1 + (k + cc_dim3 * 3) * cc_dim2)
360 * cc_dim1] + wa2[i__ - 1] * cc[m1 + (i__ + (k +
361 cc_dim3 * 3) * cc_dim2) * cc_dim1])) - (cc[m1 + (i__
362 + (k + cc_dim3) * cc_dim2) * cc_dim1] + tr12 * (wa1[
363 i__ - 2] * cc[m1 + (i__ + (k + (cc_dim3 << 1)) *
364 cc_dim2) * cc_dim1] - wa1[i__ - 1] * cc[m1 + (i__ - 1
365 + (k + (cc_dim3 << 1)) * cc_dim2) * cc_dim1] + (wa4[
366 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 5) * cc_dim2)
367 * cc_dim1] - wa4[i__ - 1] * cc[m1 + (i__ - 1 + (k +
368 cc_dim3 * 5) * cc_dim2) * cc_dim1])) + tr11 * (wa2[
369 i__ - 2] * cc[m1 + (i__ + (k + cc_dim3 * 3) * cc_dim2)
370 * cc_dim1] - wa2[i__ - 1] * cc[m1 + (i__ - 1 + (k +
371 cc_dim3 * 3) * cc_dim2) * cc_dim1] + (wa3[i__ - 2] *
372 cc[m1 + (i__ + (k + (cc_dim3 << 2)) * cc_dim2) *
373 cc_dim1] - wa3[i__ - 1] * cc[m1 + (i__ - 1 + (k + (
374 cc_dim3 << 2)) * cc_dim2) * cc_dim1])));
375/* L1002: */
376 }
377/* L102: */
378 }
379/* L103: */
380 }
381 return 0;
382} /* mradf5_ */
383