PLaSK library
Loading...
Searching...
No Matches
r1f5kb.c
Go to the documentation of this file.
1/* r1f5kb.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 r1f5kb_(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_offset, ch_dim1, ch_dim2, ch_dim3, 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_offset = 1 + cc_dim1 * (1 + cc_dim2 * 6);
70 cc -= cc_offset;
71 ch_dim1 = *in2;
72 ch_dim2 = *ido;
73 ch_dim3 = *l1;
74 ch_offset = 1 + ch_dim1 * (1 + ch_dim2 * (1 + ch_dim3));
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 + ch_dim3) * ch_dim2 + 1) * ch_dim1 + 1] = cc[((k * 5 + 1) *
86 cc_dim2 + 1) * cc_dim1 + 1] + cc[(*ido + (k * 5 + 2) *
87 cc_dim2) * cc_dim1 + 1] * 2. + cc[(*ido + (k * 5 + 4) *
88 cc_dim2) * cc_dim1 + 1] * 2.;
89 ch[((k + (ch_dim3 << 1)) * ch_dim2 + 1) * ch_dim1 + 1] = cc[((k * 5 +
90 1) * cc_dim2 + 1) * cc_dim1 + 1] + tr11 * 2. * cc[(*ido + (k *
91 5 + 2) * cc_dim2) * cc_dim1 + 1] + tr12 * 2. * cc[(*ido + (k
92 * 5 + 4) * cc_dim2) * cc_dim1 + 1] - (ti11 * 2. * cc[((k * 5
93 + 3) * cc_dim2 + 1) * cc_dim1 + 1] + ti12 * 2. * cc[((k * 5 +
94 5) * cc_dim2 + 1) * cc_dim1 + 1]);
95 ch[((k + ch_dim3 * 3) * ch_dim2 + 1) * ch_dim1 + 1] = cc[((k * 5 + 1)
96 * cc_dim2 + 1) * cc_dim1 + 1] + tr12 * 2. * cc[(*ido + (k * 5
97 + 2) * cc_dim2) * cc_dim1 + 1] + tr11 * 2. * cc[(*ido + (k *
98 5 + 4) * cc_dim2) * cc_dim1 + 1] - (ti12 * 2. * cc[((k * 5 +
99 3) * cc_dim2 + 1) * cc_dim1 + 1] - ti11 * 2. * cc[((k * 5 + 5)
100 * cc_dim2 + 1) * cc_dim1 + 1]);
101 ch[((k + (ch_dim3 << 2)) * ch_dim2 + 1) * ch_dim1 + 1] = cc[((k * 5 +
102 1) * cc_dim2 + 1) * cc_dim1 + 1] + tr12 * 2. * cc[(*ido + (k *
103 5 + 2) * cc_dim2) * cc_dim1 + 1] + tr11 * 2. * cc[(*ido + (k
104 * 5 + 4) * cc_dim2) * cc_dim1 + 1] + (ti12 * 2. * cc[((k * 5
105 + 3) * cc_dim2 + 1) * cc_dim1 + 1] - ti11 * 2. * cc[((k * 5 +
106 5) * cc_dim2 + 1) * cc_dim1 + 1]);
107 ch[((k + ch_dim3 * 5) * ch_dim2 + 1) * ch_dim1 + 1] = cc[((k * 5 + 1)
108 * cc_dim2 + 1) * cc_dim1 + 1] + tr11 * 2. * cc[(*ido + (k * 5
109 + 2) * cc_dim2) * cc_dim1 + 1] + tr12 * 2. * cc[(*ido + (k *
110 5 + 4) * cc_dim2) * cc_dim1 + 1] + (ti11 * 2. * cc[((k * 5 +
111 3) * cc_dim2 + 1) * cc_dim1 + 1] + ti12 * 2. * cc[((k * 5 + 5)
112 * 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 + ch_dim3) * ch_dim2) * ch_dim1 + 1] = cc[(i__ -
125 1 + (k * 5 + 1) * cc_dim2) * cc_dim1 + 1] + (cc[(i__ - 1
126 + (k * 5 + 3) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k
127 * 5 + 2) * cc_dim2) * cc_dim1 + 1]) + (cc[(i__ - 1 + (k *
128 5 + 5) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 +
129 4) * cc_dim2) * cc_dim1 + 1]);
130 ch[(i__ + (k + ch_dim3) * ch_dim2) * ch_dim1 + 1] = cc[(i__ + (k *
131 5 + 1) * cc_dim2) * cc_dim1 + 1] + (cc[(i__ + (k * 5 + 3)
132 * cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 2) *
133 cc_dim2) * cc_dim1 + 1]) + (cc[(i__ + (k * 5 + 5) *
134 cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 4) * cc_dim2)
135 * cc_dim1 + 1]);
136 ch[(i__ - 1 + (k + (ch_dim3 << 1)) * ch_dim2) * ch_dim1 + 1] =
137 wa1[i__ - 2] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) *
138 cc_dim1 + 1] + tr11 * (cc[(i__ - 1 + (k * 5 + 3) *
139 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) *
140 cc_dim2) * cc_dim1 + 1]) + tr12 * (cc[(i__ - 1 + (k * 5 +
141 5) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
142 cc_dim2) * cc_dim1 + 1]) - (ti11 * (cc[(i__ + (k * 5 + 3)
143 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
144 cc_dim2) * cc_dim1 + 1]) + ti12 * (cc[(i__ + (k * 5 + 5) *
145 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
146 * cc_dim1 + 1]))) - wa1[i__ - 1] * (cc[(i__ + (k * 5 + 1)
147 * cc_dim2) * cc_dim1 + 1] + tr11 * (cc[(i__ + (k * 5 + 3)
148 * cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 2) *
149 cc_dim2) * cc_dim1 + 1]) + tr12 * (cc[(i__ + (k * 5 + 5) *
150 cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 4) * cc_dim2)
151 * cc_dim1 + 1]) + (ti11 * (cc[(i__ - 1 + (k * 5 + 3) *
152 cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 2) *
153 cc_dim2) * cc_dim1 + 1]) + ti12 * (cc[(i__ - 1 + (k * 5 +
154 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 4) *
155 cc_dim2) * cc_dim1 + 1])));
156 ch[(i__ + (k + (ch_dim3 << 1)) * ch_dim2) * ch_dim1 + 1] = wa1[
157 i__ - 2] * (cc[(i__ + (k * 5 + 1) * cc_dim2) * cc_dim1 +
158 1] + tr11 * (cc[(i__ + (k * 5 + 3) * cc_dim2) * cc_dim1 +
159 1] - cc[(ic + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) +
160 tr12 * (cc[(i__ + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] -
161 cc[(ic + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]) + (ti11 *
162 (cc[(i__ - 1 + (k * 5 + 3) * cc_dim2) * cc_dim1 + 1] - cc[
163 (ic - 1 + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) + ti12 *
164 (cc[(i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] - cc[
165 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]))) + wa1[
166 i__ - 1] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) *
167 cc_dim1 + 1] + tr11 * (cc[(i__ - 1 + (k * 5 + 3) *
168 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) *
169 cc_dim2) * cc_dim1 + 1]) + tr12 * (cc[(i__ - 1 + (k * 5 +
170 5) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
171 cc_dim2) * cc_dim1 + 1]) - (ti11 * (cc[(i__ + (k * 5 + 3)
172 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
173 cc_dim2) * cc_dim1 + 1]) + ti12 * (cc[(i__ + (k * 5 + 5) *
174 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
175 * cc_dim1 + 1])));
176 ch[(i__ - 1 + (k + ch_dim3 * 3) * ch_dim2) * ch_dim1 + 1] = wa2[
177 i__ - 2] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) *
178 cc_dim1 + 1] + tr12 * (cc[(i__ - 1 + (k * 5 + 3) *
179 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) *
180 cc_dim2) * cc_dim1 + 1]) + tr11 * (cc[(i__ - 1 + (k * 5 +
181 5) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
182 cc_dim2) * cc_dim1 + 1]) - (ti12 * (cc[(i__ + (k * 5 + 3)
183 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
184 cc_dim2) * cc_dim1 + 1]) - ti11 * (cc[(i__ + (k * 5 + 5) *
185 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
186 * cc_dim1 + 1]))) - wa2[i__ - 1] * (cc[(i__ + (k * 5 + 1)
187 * cc_dim2) * cc_dim1 + 1] + tr12 * (cc[(i__ + (k * 5 + 3)
188 * cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 2) *
189 cc_dim2) * cc_dim1 + 1]) + tr11 * (cc[(i__ + (k * 5 + 5) *
190 cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 4) * cc_dim2)
191 * cc_dim1 + 1]) + (ti12 * (cc[(i__ - 1 + (k * 5 + 3) *
192 cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 2) *
193 cc_dim2) * cc_dim1 + 1]) - ti11 * (cc[(i__ - 1 + (k * 5 +
194 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 4) *
195 cc_dim2) * cc_dim1 + 1])));
196 ch[(i__ + (k + ch_dim3 * 3) * ch_dim2) * ch_dim1 + 1] = wa2[i__ -
197 2] * (cc[(i__ + (k * 5 + 1) * cc_dim2) * cc_dim1 + 1] +
198 tr12 * (cc[(i__ + (k * 5 + 3) * cc_dim2) * cc_dim1 + 1] -
199 cc[(ic + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) + tr11 * (
200 cc[(i__ + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic
201 + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]) + (ti12 * (cc[(
202 i__ - 1 + (k * 5 + 3) * cc_dim2) * cc_dim1 + 1] - cc[(ic
203 - 1 + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) - ti11 * (cc[
204 (i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic
205 - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]))) + wa2[i__
206 - 1] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) * cc_dim1 +
207 1] + tr12 * (cc[(i__ - 1 + (k * 5 + 3) * cc_dim2) *
208 cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) * cc_dim2) *
209 cc_dim1 + 1]) + tr11 * (cc[(i__ - 1 + (k * 5 + 5) *
210 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
211 cc_dim2) * cc_dim1 + 1]) - (ti12 * (cc[(i__ + (k * 5 + 3)
212 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
213 cc_dim2) * cc_dim1 + 1]) - ti11 * (cc[(i__ + (k * 5 + 5) *
214 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
215 * cc_dim1 + 1])));
216 ch[(i__ - 1 + (k + (ch_dim3 << 2)) * ch_dim2) * ch_dim1 + 1] =
217 wa3[i__ - 2] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) *
218 cc_dim1 + 1] + tr12 * (cc[(i__ - 1 + (k * 5 + 3) *
219 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) *
220 cc_dim2) * cc_dim1 + 1]) + tr11 * (cc[(i__ - 1 + (k * 5 +
221 5) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
222 cc_dim2) * cc_dim1 + 1]) + (ti12 * (cc[(i__ + (k * 5 + 3)
223 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
224 cc_dim2) * cc_dim1 + 1]) - ti11 * (cc[(i__ + (k * 5 + 5) *
225 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
226 * cc_dim1 + 1]))) - wa3[i__ - 1] * (cc[(i__ + (k * 5 + 1)
227 * cc_dim2) * cc_dim1 + 1] + tr12 * (cc[(i__ + (k * 5 + 3)
228 * cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 2) *
229 cc_dim2) * cc_dim1 + 1]) + tr11 * (cc[(i__ + (k * 5 + 5) *
230 cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 4) * cc_dim2)
231 * cc_dim1 + 1]) - (ti12 * (cc[(i__ - 1 + (k * 5 + 3) *
232 cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 2) *
233 cc_dim2) * cc_dim1 + 1]) - ti11 * (cc[(i__ - 1 + (k * 5 +
234 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 4) *
235 cc_dim2) * cc_dim1 + 1])));
236 ch[(i__ + (k + (ch_dim3 << 2)) * ch_dim2) * ch_dim1 + 1] = wa3[
237 i__ - 2] * (cc[(i__ + (k * 5 + 1) * cc_dim2) * cc_dim1 +
238 1] + tr12 * (cc[(i__ + (k * 5 + 3) * cc_dim2) * cc_dim1 +
239 1] - cc[(ic + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) +
240 tr11 * (cc[(i__ + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] -
241 cc[(ic + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]) - (ti12 *
242 (cc[(i__ - 1 + (k * 5 + 3) * cc_dim2) * cc_dim1 + 1] - cc[
243 (ic - 1 + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) - ti11 *
244 (cc[(i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] - cc[
245 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]))) + wa3[
246 i__ - 1] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) *
247 cc_dim1 + 1] + tr12 * (cc[(i__ - 1 + (k * 5 + 3) *
248 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) *
249 cc_dim2) * cc_dim1 + 1]) + tr11 * (cc[(i__ - 1 + (k * 5 +
250 5) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
251 cc_dim2) * cc_dim1 + 1]) + (ti12 * (cc[(i__ + (k * 5 + 3)
252 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
253 cc_dim2) * cc_dim1 + 1]) - ti11 * (cc[(i__ + (k * 5 + 5) *
254 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
255 * cc_dim1 + 1])));
256 ch[(i__ - 1 + (k + ch_dim3 * 5) * ch_dim2) * ch_dim1 + 1] = wa4[
257 i__ - 2] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) *
258 cc_dim1 + 1] + tr11 * (cc[(i__ - 1 + (k * 5 + 3) *
259 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) *
260 cc_dim2) * cc_dim1 + 1]) + tr12 * (cc[(i__ - 1 + (k * 5 +
261 5) * cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
262 cc_dim2) * cc_dim1 + 1]) + (ti11 * (cc[(i__ + (k * 5 + 3)
263 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
264 cc_dim2) * cc_dim1 + 1]) + ti12 * (cc[(i__ + (k * 5 + 5) *
265 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
266 * cc_dim1 + 1]))) - wa4[i__ - 1] * (cc[(i__ + (k * 5 + 1)
267 * cc_dim2) * cc_dim1 + 1] + tr11 * (cc[(i__ + (k * 5 + 3)
268 * cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 2) *
269 cc_dim2) * cc_dim1 + 1]) + tr12 * (cc[(i__ + (k * 5 + 5) *
270 cc_dim2) * cc_dim1 + 1] - cc[(ic + (k * 5 + 4) * cc_dim2)
271 * cc_dim1 + 1]) - (ti11 * (cc[(i__ - 1 + (k * 5 + 3) *
272 cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 2) *
273 cc_dim2) * cc_dim1 + 1]) + ti12 * (cc[(i__ - 1 + (k * 5 +
274 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic - 1 + (k * 5 + 4) *
275 cc_dim2) * cc_dim1 + 1])));
276 ch[(i__ + (k + ch_dim3 * 5) * ch_dim2) * ch_dim1 + 1] = wa4[i__ -
277 2] * (cc[(i__ + (k * 5 + 1) * cc_dim2) * cc_dim1 + 1] +
278 tr11 * (cc[(i__ + (k * 5 + 3) * cc_dim2) * cc_dim1 + 1] -
279 cc[(ic + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) + tr12 * (
280 cc[(i__ + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic
281 + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]) - (ti11 * (cc[(
282 i__ - 1 + (k * 5 + 3) * cc_dim2) * cc_dim1 + 1] - cc[(ic
283 - 1 + (k * 5 + 2) * cc_dim2) * cc_dim1 + 1]) + ti12 * (cc[
284 (i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1 + 1] - cc[(ic
285 - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1 + 1]))) + wa4[i__
286 - 1] * (cc[(i__ - 1 + (k * 5 + 1) * cc_dim2) * cc_dim1 +
287 1] + tr11 * (cc[(i__ - 1 + (k * 5 + 3) * cc_dim2) *
288 cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 2) * cc_dim2) *
289 cc_dim1 + 1]) + tr12 * (cc[(i__ - 1 + (k * 5 + 5) *
290 cc_dim2) * cc_dim1 + 1] + cc[(ic - 1 + (k * 5 + 4) *
291 cc_dim2) * cc_dim1 + 1]) + (ti11 * (cc[(i__ + (k * 5 + 3)
292 * cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 2) *
293 cc_dim2) * cc_dim1 + 1]) + ti12 * (cc[(i__ + (k * 5 + 5) *
294 cc_dim2) * cc_dim1 + 1] + cc[(ic + (k * 5 + 4) * cc_dim2)
295 * cc_dim1 + 1])));
296/* L102: */
297 }
298/* L103: */
299 }
300 return 0;
301} /* r1f5kb_ */
302