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