PLaSK library
Loading...
Searching...
No Matches
mradb5.c
Go to the documentation of this file.
1/* mradb5.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 mradb5_(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_offset, ch_dim1, ch_dim2, ch_dim3, 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_offset = 1 + cc_dim1 * (1 + cc_dim2 * 6);
71 cc -= cc_offset;
72 ch_dim1 = *in2;
73 ch_dim2 = *ido;
74 ch_dim3 = *l1;
75 ch_offset = 1 + ch_dim1 * (1 + ch_dim2 * (1 + ch_dim3));
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 + ch_dim3) * ch_dim2 + 1) * ch_dim1] = cc[m1 + ((k *
94 5 + 1) * cc_dim2 + 1) * cc_dim1] + cc[m1 + (*ido + (k * 5
95 + 2) * cc_dim2) * cc_dim1] * 2. + cc[m1 + (*ido + (k * 5
96 + 4) * cc_dim2) * cc_dim1] * 2.;
97 ch[m2 + ((k + (ch_dim3 << 1)) * ch_dim2 + 1) * ch_dim1] = cc[m1 +
98 ((k * 5 + 1) * cc_dim2 + 1) * cc_dim1] + tr11 * 2. * cc[
99 m1 + (*ido + (k * 5 + 2) * cc_dim2) * cc_dim1] + tr12 *
100 2. * cc[m1 + (*ido + (k * 5 + 4) * cc_dim2) * cc_dim1] - (
101 ti11 * 2. * cc[m1 + ((k * 5 + 3) * cc_dim2 + 1) * cc_dim1]
102 + ti12 * 2. * cc[m1 + ((k * 5 + 5) * cc_dim2 + 1) *
103 cc_dim1]);
104 ch[m2 + ((k + ch_dim3 * 3) * ch_dim2 + 1) * ch_dim1] = cc[m1 + ((
105 k * 5 + 1) * cc_dim2 + 1) * cc_dim1] + tr12 * 2. * cc[m1
106 + (*ido + (k * 5 + 2) * cc_dim2) * cc_dim1] + tr11 * 2. *
107 cc[m1 + (*ido + (k * 5 + 4) * cc_dim2) * cc_dim1] - (ti12
108 * 2. * cc[m1 + ((k * 5 + 3) * cc_dim2 + 1) * cc_dim1] -
109 ti11 * 2. * cc[m1 + ((k * 5 + 5) * cc_dim2 + 1) * cc_dim1]
110 );
111 ch[m2 + ((k + (ch_dim3 << 2)) * ch_dim2 + 1) * ch_dim1] = cc[m1 +
112 ((k * 5 + 1) * cc_dim2 + 1) * cc_dim1] + tr12 * 2. * cc[
113 m1 + (*ido + (k * 5 + 2) * cc_dim2) * cc_dim1] + tr11 *
114 2. * cc[m1 + (*ido + (k * 5 + 4) * cc_dim2) * cc_dim1] + (
115 ti12 * 2. * cc[m1 + ((k * 5 + 3) * cc_dim2 + 1) * cc_dim1]
116 - ti11 * 2. * cc[m1 + ((k * 5 + 5) * cc_dim2 + 1) *
117 cc_dim1]);
118 ch[m2 + ((k + ch_dim3 * 5) * ch_dim2 + 1) * ch_dim1] = cc[m1 + ((
119 k * 5 + 1) * cc_dim2 + 1) * cc_dim1] + tr11 * 2. * cc[m1
120 + (*ido + (k * 5 + 2) * cc_dim2) * cc_dim1] + tr12 * 2. *
121 cc[m1 + (*ido + (k * 5 + 4) * cc_dim2) * cc_dim1] + (ti11
122 * 2. * cc[m1 + ((k * 5 + 3) * cc_dim2 + 1) * cc_dim1] +
123 ti12 * 2. * cc[m1 + ((k * 5 + 5) * cc_dim2 + 1) * cc_dim1]
124 );
125/* L1001: */
126 }
127/* L101: */
128 }
129 if (*ido == 1) {
130 return 0;
131 }
132 idp2 = *ido + 2;
133 i__1 = *l1;
134 for (k = 1; k <= i__1; ++k) {
135 i__3 = *ido;
136 for (i__ = 3; i__ <= i__3; i__ += 2) {
137 ic = idp2 - i__;
138 m2 = m2s;
139 i__2 = m1d;
140 i__4 = *im1;
141 for (m1 = 1; i__4 < 0 ? m1 >= i__2 : m1 <= i__2; m1 += i__4) {
142 m2 += *im2;
143 ch[m2 + (i__ - 1 + (k + ch_dim3) * ch_dim2) * ch_dim1] = cc[
144 m1 + (i__ - 1 + (k * 5 + 1) * cc_dim2) * cc_dim1] + (
145 cc[m1 + (i__ - 1 + (k * 5 + 3) * cc_dim2) * cc_dim1]
146 + cc[m1 + (ic - 1 + (k * 5 + 2) * cc_dim2) * cc_dim1])
147 + (cc[m1 + (i__ - 1 + (k * 5 + 5) * cc_dim2) *
148 cc_dim1] + cc[m1 + (ic - 1 + (k * 5 + 4) * cc_dim2) *
149 cc_dim1]);
150 ch[m2 + (i__ + (k + ch_dim3) * ch_dim2) * ch_dim1] = cc[m1 + (
151 i__ + (k * 5 + 1) * cc_dim2) * cc_dim1] + (cc[m1 + (
152 i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic
153 + (k * 5 + 2) * cc_dim2) * cc_dim1]) + (cc[m1 + (i__
154 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (
155 k * 5 + 4) * cc_dim2) * cc_dim1]);
156 ch[m2 + (i__ - 1 + (k + (ch_dim3 << 1)) * ch_dim2) * ch_dim1]
157 = wa1[i__ - 2] * (cc[m1 + (i__ - 1 + (k * 5 + 1) *
158 cc_dim2) * cc_dim1] + tr11 * (cc[m1 + (i__ - 1 + (k *
159 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
160 5 + 2) * cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__
161 - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic
162 - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]) - (ti11 * (
163 cc[m1 + (i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[
164 m1 + (ic + (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 *
165 (cc[m1 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] +
166 cc[m1 + (ic + (k * 5 + 4) * cc_dim2) * cc_dim1]))) -
167 wa1[i__ - 1] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2)
168 * cc_dim1] + tr11 * (cc[m1 + (i__ + (k * 5 + 3) *
169 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
170 cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__ + (k * 5
171 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
172 * cc_dim2) * cc_dim1]) + (ti11 * (cc[m1 + (i__ - 1 +
173 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
174 (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 * (cc[m1 + (
175 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
176 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1])));
177 ch[m2 + (i__ + (k + (ch_dim3 << 1)) * ch_dim2) * ch_dim1] =
178 wa1[i__ - 2] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2)
179 * cc_dim1] + tr11 * (cc[m1 + (i__ + (k * 5 + 3) *
180 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
181 cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__ + (k * 5
182 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
183 * cc_dim2) * cc_dim1]) + (ti11 * (cc[m1 + (i__ - 1 +
184 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
185 (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 * (cc[m1 + (
186 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
187 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]))) + wa1[
188 i__ - 1] * (cc[m1 + (i__ - 1 + (k * 5 + 1) * cc_dim2)
189 * cc_dim1] + tr11 * (cc[m1 + (i__ - 1 + (k * 5 + 3) *
190 cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k * 5 + 2) *
191 cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__ - 1 + (k *
192 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
193 5 + 4) * cc_dim2) * cc_dim1]) - (ti11 * (cc[m1 + (
194 i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic
195 + (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 * (cc[m1
196 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (
197 ic + (k * 5 + 4) * cc_dim2) * cc_dim1])));
198 ch[m2 + (i__ - 1 + (k + ch_dim3 * 3) * ch_dim2) * ch_dim1] =
199 wa2[i__ - 2] * (cc[m1 + (i__ - 1 + (k * 5 + 1) *
200 cc_dim2) * cc_dim1] + tr12 * (cc[m1 + (i__ - 1 + (k *
201 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
202 5 + 2) * cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__
203 - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic
204 - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]) - (ti12 * (
205 cc[m1 + (i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[
206 m1 + (ic + (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 *
207 (cc[m1 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] +
208 cc[m1 + (ic + (k * 5 + 4) * cc_dim2) * cc_dim1]))) -
209 wa2[i__ - 1] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2)
210 * cc_dim1] + tr12 * (cc[m1 + (i__ + (k * 5 + 3) *
211 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
212 cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__ + (k * 5
213 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
214 * cc_dim2) * cc_dim1]) + (ti12 * (cc[m1 + (i__ - 1 +
215 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
216 (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 * (cc[m1 + (
217 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
218 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1])));
219 ch[m2 + (i__ + (k + ch_dim3 * 3) * ch_dim2) * ch_dim1] = wa2[
220 i__ - 2] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2) *
221 cc_dim1] + tr12 * (cc[m1 + (i__ + (k * 5 + 3) *
222 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
223 cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__ + (k * 5
224 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
225 * cc_dim2) * cc_dim1]) + (ti12 * (cc[m1 + (i__ - 1 +
226 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
227 (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 * (cc[m1 + (
228 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
229 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]))) + wa2[
230 i__ - 1] * (cc[m1 + (i__ - 1 + (k * 5 + 1) * cc_dim2)
231 * cc_dim1] + tr12 * (cc[m1 + (i__ - 1 + (k * 5 + 3) *
232 cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k * 5 + 2) *
233 cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__ - 1 + (k *
234 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
235 5 + 4) * cc_dim2) * cc_dim1]) - (ti12 * (cc[m1 + (
236 i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic
237 + (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 * (cc[m1
238 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (
239 ic + (k * 5 + 4) * cc_dim2) * cc_dim1])));
240 ch[m2 + (i__ - 1 + (k + (ch_dim3 << 2)) * ch_dim2) * ch_dim1]
241 = wa3[i__ - 2] * (cc[m1 + (i__ - 1 + (k * 5 + 1) *
242 cc_dim2) * cc_dim1] + tr12 * (cc[m1 + (i__ - 1 + (k *
243 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
244 5 + 2) * cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__
245 - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic
246 - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]) + (ti12 * (
247 cc[m1 + (i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[
248 m1 + (ic + (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 *
249 (cc[m1 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] +
250 cc[m1 + (ic + (k * 5 + 4) * cc_dim2) * cc_dim1]))) -
251 wa3[i__ - 1] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2)
252 * cc_dim1] + tr12 * (cc[m1 + (i__ + (k * 5 + 3) *
253 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
254 cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__ + (k * 5
255 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
256 * cc_dim2) * cc_dim1]) - (ti12 * (cc[m1 + (i__ - 1 +
257 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
258 (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 * (cc[m1 + (
259 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
260 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1])));
261 ch[m2 + (i__ + (k + (ch_dim3 << 2)) * ch_dim2) * ch_dim1] =
262 wa3[i__ - 2] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2)
263 * cc_dim1] + tr12 * (cc[m1 + (i__ + (k * 5 + 3) *
264 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
265 cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__ + (k * 5
266 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
267 * cc_dim2) * cc_dim1]) - (ti12 * (cc[m1 + (i__ - 1 +
268 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
269 (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 * (cc[m1 + (
270 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
271 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]))) + wa3[
272 i__ - 1] * (cc[m1 + (i__ - 1 + (k * 5 + 1) * cc_dim2)
273 * cc_dim1] + tr12 * (cc[m1 + (i__ - 1 + (k * 5 + 3) *
274 cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k * 5 + 2) *
275 cc_dim2) * cc_dim1]) + tr11 * (cc[m1 + (i__ - 1 + (k *
276 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
277 5 + 4) * cc_dim2) * cc_dim1]) + (ti12 * (cc[m1 + (
278 i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic
279 + (k * 5 + 2) * cc_dim2) * cc_dim1]) - ti11 * (cc[m1
280 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (
281 ic + (k * 5 + 4) * cc_dim2) * cc_dim1])));
282 ch[m2 + (i__ - 1 + (k + ch_dim3 * 5) * ch_dim2) * ch_dim1] =
283 wa4[i__ - 2] * (cc[m1 + (i__ - 1 + (k * 5 + 1) *
284 cc_dim2) * cc_dim1] + tr11 * (cc[m1 + (i__ - 1 + (k *
285 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
286 5 + 2) * cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__
287 - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic
288 - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]) + (ti11 * (
289 cc[m1 + (i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[
290 m1 + (ic + (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 *
291 (cc[m1 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] +
292 cc[m1 + (ic + (k * 5 + 4) * cc_dim2) * cc_dim1]))) -
293 wa4[i__ - 1] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2)
294 * cc_dim1] + tr11 * (cc[m1 + (i__ + (k * 5 + 3) *
295 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
296 cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__ + (k * 5
297 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
298 * cc_dim2) * cc_dim1]) - (ti11 * (cc[m1 + (i__ - 1 +
299 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
300 (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 * (cc[m1 + (
301 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
302 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1])));
303 ch[m2 + (i__ + (k + ch_dim3 * 5) * ch_dim2) * ch_dim1] = wa4[
304 i__ - 2] * (cc[m1 + (i__ + (k * 5 + 1) * cc_dim2) *
305 cc_dim1] + tr11 * (cc[m1 + (i__ + (k * 5 + 3) *
306 cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 2) *
307 cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__ + (k * 5
308 + 5) * cc_dim2) * cc_dim1] - cc[m1 + (ic + (k * 5 + 4)
309 * cc_dim2) * cc_dim1]) - (ti11 * (cc[m1 + (i__ - 1 +
310 (k * 5 + 3) * cc_dim2) * cc_dim1] - cc[m1 + (ic - 1 +
311 (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 * (cc[m1 + (
312 i__ - 1 + (k * 5 + 5) * cc_dim2) * cc_dim1] - cc[m1 +
313 (ic - 1 + (k * 5 + 4) * cc_dim2) * cc_dim1]))) + wa4[
314 i__ - 1] * (cc[m1 + (i__ - 1 + (k * 5 + 1) * cc_dim2)
315 * cc_dim1] + tr11 * (cc[m1 + (i__ - 1 + (k * 5 + 3) *
316 cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k * 5 + 2) *
317 cc_dim2) * cc_dim1]) + tr12 * (cc[m1 + (i__ - 1 + (k *
318 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (ic - 1 + (k *
319 5 + 4) * cc_dim2) * cc_dim1]) + (ti11 * (cc[m1 + (
320 i__ + (k * 5 + 3) * cc_dim2) * cc_dim1] + cc[m1 + (ic
321 + (k * 5 + 2) * cc_dim2) * cc_dim1]) + ti12 * (cc[m1
322 + (i__ + (k * 5 + 5) * cc_dim2) * cc_dim1] + cc[m1 + (
323 ic + (k * 5 + 4) * cc_dim2) * cc_dim1])));
324/* L1002: */
325 }
326/* L102: */
327 }
328/* L103: */
329 }
330 return 0;
331} /* mradb5_ */
332