Rename *ll* and *ul* to ll and ul in $defint
[maxima.git] / share / fftpack5 / lisp / mradfg.lisp
blob7496d4c1359835dfbd0d64ed2d6037a7af9b7318
1 ;;; Compiled by f2cl version:
2 ;;; ("f2cl1.l,v 95098eb54f13 2013/04/01 00:45:16 toy $"
3 ;;; "f2cl2.l,v 95098eb54f13 2013/04/01 00:45:16 toy $"
4 ;;; "f2cl3.l,v 96616d88fb7e 2008/02/22 22:19:34 rtoy $"
5 ;;; "f2cl4.l,v 96616d88fb7e 2008/02/22 22:19:34 rtoy $"
6 ;;; "f2cl5.l,v 95098eb54f13 2013/04/01 00:45:16 toy $"
7 ;;; "f2cl6.l,v 1d5cbacbb977 2008/08/24 00:56:27 rtoy $"
8 ;;; "macros.l,v 1409c1352feb 2013/03/24 20:44:50 toy $")
10 ;;; Using Lisp CMU Common Lisp snapshot-2020-04 (21D Unicode)
11 ;;;
12 ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
13 ;;; (:coerce-assigns :as-needed) (:array-type ':array)
14 ;;; (:array-slicing t) (:declare-common nil)
15 ;;; (:float-format single-float))
17 (in-package "FFTPACK5")
20 (defun mradfg (m ido ip l1 idl1 cc c1 c2 im1 in1 ch ch2 im2 in2 wa)
21 (declare (type (array double-float (*)) wa ch2 ch c2 c1 cc)
22 (type (f2cl-lib:integer4) in2 im2 in1 im1 idl1 l1 ip ido m))
23 (f2cl-lib:with-multi-array-data
24 ((cc double-float cc-%data% cc-%offset%)
25 (c1 double-float c1-%data% c1-%offset%)
26 (c2 double-float c2-%data% c2-%offset%)
27 (ch double-float ch-%data% ch-%offset%)
28 (ch2 double-float ch2-%data% ch2-%offset%)
29 (wa double-float wa-%data% wa-%offset%))
30 (prog ((ic 0) (j2 0) (ar2h 0.0d0) (ai2 0.0d0) (ar2 0.0d0) (ds2 0.0d0)
31 (dc2 0.0d0) (ar1h 0.0d0) (lc 0) (l 0) (ai1 0.0d0) (ar1 0.0d0) (jc 0)
32 (i 0) (idij 0) (is 0) (k 0) (j 0) (m1 0) (m2 0) (ik 0) (nbd 0)
33 (idp2 0) (ipp2 0) (ipph 0) (dsp 0.0d0) (dcp 0.0d0) (arg 0.0d0)
34 (tpi 0.0d0) (m2s 0) (m1d 0))
35 (declare (type (double-float) tpi arg dcp dsp ar1 ai1 ar1h dc2 ds2 ar2
36 ai2 ar2h)
37 (type (f2cl-lib:integer4) m1d m2s ipph ipp2 idp2 nbd ik m2 m1 j
38 k is idij i jc l lc j2 ic))
39 (setf m1d
40 (f2cl-lib:int-add (f2cl-lib:int-mul (f2cl-lib:int-sub m 1) im1)
41 1))
42 (setf m2s (f2cl-lib:int-sub 1 im2))
43 (setf tpi (* 2.0d0 4.0d0 (atan 1.0d0)))
44 (setf arg (/ tpi (f2cl-lib:ffloat ip)))
45 (setf dcp (cos arg))
46 (setf dsp (sin arg))
47 (setf ipph (the f2cl-lib:integer4 (truncate (+ ip 1) 2)))
48 (setf ipp2 (f2cl-lib:int-add ip 2))
49 (setf idp2 (f2cl-lib:int-add ido 2))
50 (setf nbd (the f2cl-lib:integer4 (truncate (- ido 1) 2)))
51 (if (= ido 1) (go label119))
52 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
53 ((> ik idl1) nil)
54 (tagbody
55 (setf m2 m2s)
56 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
57 ((> m1 m1d) nil)
58 (tagbody
59 (setf m2 (f2cl-lib:int-add m2 im2))
60 (setf (f2cl-lib:fref ch2-%data%
61 (m2 ik 1)
62 ((1 in2) (1 idl1) (1 ip))
63 ch2-%offset%)
64 (f2cl-lib:fref c2-%data%
65 (m1 ik 1)
66 ((1 in1) (1 idl1) (1 ip))
67 c2-%offset%))
68 label1001))
69 label101))
70 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
71 ((> j ip) nil)
72 (tagbody
73 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
74 ((> k l1) nil)
75 (tagbody
76 (setf m2 m2s)
77 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
78 ((> m1 m1d) nil)
79 (tagbody
80 (setf m2 (f2cl-lib:int-add m2 im2))
81 (setf (f2cl-lib:fref ch-%data%
82 (m2 1 k j)
83 ((1 in2) (1 ido) (1 l1) (1 ip))
84 ch-%offset%)
85 (f2cl-lib:fref c1-%data%
86 (m1 1 k j)
87 ((1 in1) (1 ido) (1 l1) (1 ip))
88 c1-%offset%))
89 label1002))
90 label102))
91 label103))
92 (if (> nbd l1) (go label107))
93 (setf is (f2cl-lib:int-sub ido))
94 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
95 ((> j ip) nil)
96 (tagbody
97 (setf is (f2cl-lib:int-add is ido))
98 (setf idij is)
99 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
100 ((> i ido) nil)
101 (tagbody
102 (setf idij (f2cl-lib:int-add idij 2))
103 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
104 ((> k l1) nil)
105 (tagbody
106 (setf m2 m2s)
107 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
108 ((> m1 m1d) nil)
109 (tagbody
110 (setf m2 (f2cl-lib:int-add m2 im2))
111 (setf (f2cl-lib:fref ch-%data%
112 (m2 (f2cl-lib:int-sub i 1) k j)
113 ((1 in2) (1 ido) (1 l1) (1 ip))
114 ch-%offset%)
117 (f2cl-lib:fref wa-%data%
118 ((f2cl-lib:int-sub idij 1))
119 ((1 ido))
120 wa-%offset%)
121 (f2cl-lib:fref c1-%data%
122 (m1 (f2cl-lib:int-sub i 1) k j)
123 ((1 in1) (1 ido) (1 l1) (1 ip))
124 c1-%offset%))
126 (f2cl-lib:fref wa-%data%
127 (idij)
128 ((1 ido))
129 wa-%offset%)
130 (f2cl-lib:fref c1-%data%
131 (m1 i k j)
132 ((1 in1) (1 ido) (1 l1) (1 ip))
133 c1-%offset%))))
134 (setf (f2cl-lib:fref ch-%data%
135 (m2 i k j)
136 ((1 in2) (1 ido) (1 l1) (1 ip))
137 ch-%offset%)
140 (f2cl-lib:fref wa-%data%
141 ((f2cl-lib:int-sub idij 1))
142 ((1 ido))
143 wa-%offset%)
144 (f2cl-lib:fref c1-%data%
145 (m1 i k j)
146 ((1 in1) (1 ido) (1 l1) (1 ip))
147 c1-%offset%))
149 (f2cl-lib:fref wa-%data%
150 (idij)
151 ((1 ido))
152 wa-%offset%)
153 (f2cl-lib:fref c1-%data%
154 (m1 (f2cl-lib:int-sub i 1) k j)
155 ((1 in1) (1 ido) (1 l1) (1 ip))
156 c1-%offset%))))
157 label1004))
158 label104))
159 label105))
160 label106))
161 (go label111)
162 label107
163 (setf is (f2cl-lib:int-sub ido))
164 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
165 ((> j ip) nil)
166 (tagbody
167 (setf is (f2cl-lib:int-add is ido))
168 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
169 ((> k l1) nil)
170 (tagbody
171 (setf idij is)
172 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
173 ((> i ido) nil)
174 (tagbody
175 (setf idij (f2cl-lib:int-add idij 2))
176 (setf m2 m2s)
177 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
178 ((> m1 m1d) nil)
179 (tagbody
180 (setf m2 (f2cl-lib:int-add m2 im2))
181 (setf (f2cl-lib:fref ch-%data%
182 (m2 (f2cl-lib:int-sub i 1) k j)
183 ((1 in2) (1 ido) (1 l1) (1 ip))
184 ch-%offset%)
187 (f2cl-lib:fref wa-%data%
188 ((f2cl-lib:int-sub idij 1))
189 ((1 ido))
190 wa-%offset%)
191 (f2cl-lib:fref c1-%data%
192 (m1 (f2cl-lib:int-sub i 1) k j)
193 ((1 in1) (1 ido) (1 l1) (1 ip))
194 c1-%offset%))
196 (f2cl-lib:fref wa-%data%
197 (idij)
198 ((1 ido))
199 wa-%offset%)
200 (f2cl-lib:fref c1-%data%
201 (m1 i k j)
202 ((1 in1) (1 ido) (1 l1) (1 ip))
203 c1-%offset%))))
204 (setf (f2cl-lib:fref ch-%data%
205 (m2 i k j)
206 ((1 in2) (1 ido) (1 l1) (1 ip))
207 ch-%offset%)
210 (f2cl-lib:fref wa-%data%
211 ((f2cl-lib:int-sub idij 1))
212 ((1 ido))
213 wa-%offset%)
214 (f2cl-lib:fref c1-%data%
215 (m1 i k j)
216 ((1 in1) (1 ido) (1 l1) (1 ip))
217 c1-%offset%))
219 (f2cl-lib:fref wa-%data%
220 (idij)
221 ((1 ido))
222 wa-%offset%)
223 (f2cl-lib:fref c1-%data%
224 (m1 (f2cl-lib:int-sub i 1) k j)
225 ((1 in1) (1 ido) (1 l1) (1 ip))
226 c1-%offset%))))
227 label1008))
228 label108))
229 label109))
230 label110))
231 label111
232 (if (< nbd l1) (go label115))
233 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
234 ((> j ipph) nil)
235 (tagbody
236 (setf jc (f2cl-lib:int-sub ipp2 j))
237 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
238 ((> k l1) nil)
239 (tagbody
240 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
241 ((> i ido) nil)
242 (tagbody
243 (setf m2 m2s)
244 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
245 ((> m1 m1d) nil)
246 (tagbody
247 (setf m2 (f2cl-lib:int-add m2 im2))
248 (setf (f2cl-lib:fref c1-%data%
249 (m1 (f2cl-lib:int-sub i 1) k j)
250 ((1 in1) (1 ido) (1 l1) (1 ip))
251 c1-%offset%)
253 (f2cl-lib:fref ch-%data%
254 (m2 (f2cl-lib:int-sub i 1) k j)
255 ((1 in2) (1 ido) (1 l1) (1 ip))
256 ch-%offset%)
257 (f2cl-lib:fref ch-%data%
258 (m2 (f2cl-lib:int-sub i 1) k jc)
259 ((1 in2) (1 ido) (1 l1) (1 ip))
260 ch-%offset%)))
261 (setf (f2cl-lib:fref c1-%data%
262 (m1 (f2cl-lib:int-sub i 1) k jc)
263 ((1 in1) (1 ido) (1 l1) (1 ip))
264 c1-%offset%)
266 (f2cl-lib:fref ch-%data%
267 (m2 i k j)
268 ((1 in2) (1 ido) (1 l1) (1 ip))
269 ch-%offset%)
270 (f2cl-lib:fref ch-%data%
271 (m2 i k jc)
272 ((1 in2) (1 ido) (1 l1) (1 ip))
273 ch-%offset%)))
274 (setf (f2cl-lib:fref c1-%data%
275 (m1 i k j)
276 ((1 in1) (1 ido) (1 l1) (1 ip))
277 c1-%offset%)
279 (f2cl-lib:fref ch-%data%
280 (m2 i k j)
281 ((1 in2) (1 ido) (1 l1) (1 ip))
282 ch-%offset%)
283 (f2cl-lib:fref ch-%data%
284 (m2 i k jc)
285 ((1 in2) (1 ido) (1 l1) (1 ip))
286 ch-%offset%)))
287 (setf (f2cl-lib:fref c1-%data%
288 (m1 i k jc)
289 ((1 in1) (1 ido) (1 l1) (1 ip))
290 c1-%offset%)
292 (f2cl-lib:fref ch-%data%
293 (m2 (f2cl-lib:int-sub i 1) k jc)
294 ((1 in2) (1 ido) (1 l1) (1 ip))
295 ch-%offset%)
296 (f2cl-lib:fref ch-%data%
297 (m2 (f2cl-lib:int-sub i 1) k j)
298 ((1 in2) (1 ido) (1 l1) (1 ip))
299 ch-%offset%)))
300 label1012))
301 label112))
302 label113))
303 label114))
304 (go label121)
305 label115
306 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
307 ((> j ipph) nil)
308 (tagbody
309 (setf jc (f2cl-lib:int-sub ipp2 j))
310 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
311 ((> i ido) nil)
312 (tagbody
313 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
314 ((> k l1) nil)
315 (tagbody
316 (setf m2 m2s)
317 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
318 ((> m1 m1d) nil)
319 (tagbody
320 (setf m2 (f2cl-lib:int-add m2 im2))
321 (setf (f2cl-lib:fref c1-%data%
322 (m1 (f2cl-lib:int-sub i 1) k j)
323 ((1 in1) (1 ido) (1 l1) (1 ip))
324 c1-%offset%)
326 (f2cl-lib:fref ch-%data%
327 (m2 (f2cl-lib:int-sub i 1) k j)
328 ((1 in2) (1 ido) (1 l1) (1 ip))
329 ch-%offset%)
330 (f2cl-lib:fref ch-%data%
331 (m2 (f2cl-lib:int-sub i 1) k jc)
332 ((1 in2) (1 ido) (1 l1) (1 ip))
333 ch-%offset%)))
334 (setf (f2cl-lib:fref c1-%data%
335 (m1 (f2cl-lib:int-sub i 1) k jc)
336 ((1 in1) (1 ido) (1 l1) (1 ip))
337 c1-%offset%)
339 (f2cl-lib:fref ch-%data%
340 (m2 i k j)
341 ((1 in2) (1 ido) (1 l1) (1 ip))
342 ch-%offset%)
343 (f2cl-lib:fref ch-%data%
344 (m2 i k jc)
345 ((1 in2) (1 ido) (1 l1) (1 ip))
346 ch-%offset%)))
347 (setf (f2cl-lib:fref c1-%data%
348 (m1 i k j)
349 ((1 in1) (1 ido) (1 l1) (1 ip))
350 c1-%offset%)
352 (f2cl-lib:fref ch-%data%
353 (m2 i k j)
354 ((1 in2) (1 ido) (1 l1) (1 ip))
355 ch-%offset%)
356 (f2cl-lib:fref ch-%data%
357 (m2 i k jc)
358 ((1 in2) (1 ido) (1 l1) (1 ip))
359 ch-%offset%)))
360 (setf (f2cl-lib:fref c1-%data%
361 (m1 i k jc)
362 ((1 in1) (1 ido) (1 l1) (1 ip))
363 c1-%offset%)
365 (f2cl-lib:fref ch-%data%
366 (m2 (f2cl-lib:int-sub i 1) k jc)
367 ((1 in2) (1 ido) (1 l1) (1 ip))
368 ch-%offset%)
369 (f2cl-lib:fref ch-%data%
370 (m2 (f2cl-lib:int-sub i 1) k j)
371 ((1 in2) (1 ido) (1 l1) (1 ip))
372 ch-%offset%)))
373 label1016))
374 label116))
375 label117))
376 label118))
377 (go label121)
378 label119
379 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
380 ((> ik idl1) nil)
381 (tagbody
382 (setf m2 m2s)
383 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
384 ((> m1 m1d) nil)
385 (tagbody
386 (setf m2 (f2cl-lib:int-add m2 im2))
387 (setf (f2cl-lib:fref c2-%data%
388 (m1 ik 1)
389 ((1 in1) (1 idl1) (1 ip))
390 c2-%offset%)
391 (f2cl-lib:fref ch2-%data%
392 (m2 ik 1)
393 ((1 in2) (1 idl1) (1 ip))
394 ch2-%offset%))
395 label1020))
396 label120))
397 label121
398 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
399 ((> j ipph) nil)
400 (tagbody
401 (setf jc (f2cl-lib:int-sub ipp2 j))
402 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
403 ((> k l1) nil)
404 (tagbody
405 (setf m2 m2s)
406 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
407 ((> m1 m1d) nil)
408 (tagbody
409 (setf m2 (f2cl-lib:int-add m2 im2))
410 (setf (f2cl-lib:fref c1-%data%
411 (m1 1 k j)
412 ((1 in1) (1 ido) (1 l1) (1 ip))
413 c1-%offset%)
415 (f2cl-lib:fref ch-%data%
416 (m2 1 k j)
417 ((1 in2) (1 ido) (1 l1) (1 ip))
418 ch-%offset%)
419 (f2cl-lib:fref ch-%data%
420 (m2 1 k jc)
421 ((1 in2) (1 ido) (1 l1) (1 ip))
422 ch-%offset%)))
423 (setf (f2cl-lib:fref c1-%data%
424 (m1 1 k jc)
425 ((1 in1) (1 ido) (1 l1) (1 ip))
426 c1-%offset%)
428 (f2cl-lib:fref ch-%data%
429 (m2 1 k jc)
430 ((1 in2) (1 ido) (1 l1) (1 ip))
431 ch-%offset%)
432 (f2cl-lib:fref ch-%data%
433 (m2 1 k j)
434 ((1 in2) (1 ido) (1 l1) (1 ip))
435 ch-%offset%)))
436 label1022))
437 label122))
438 label123))
439 (setf ar1 1.0d0)
440 (setf ai1 0.0d0)
441 (f2cl-lib:fdo (l 2 (f2cl-lib:int-add l 1))
442 ((> l ipph) nil)
443 (tagbody
444 (setf lc (f2cl-lib:int-sub ipp2 l))
445 (setf ar1h (- (* dcp ar1) (* dsp ai1)))
446 (setf ai1 (+ (* dcp ai1) (* dsp ar1)))
447 (setf ar1 ar1h)
448 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
449 ((> ik idl1) nil)
450 (tagbody
451 (setf m2 m2s)
452 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
453 ((> m1 m1d) nil)
454 (tagbody
455 (setf m2 (f2cl-lib:int-add m2 im2))
456 (setf (f2cl-lib:fref ch2-%data%
457 (m2 ik l)
458 ((1 in2) (1 idl1) (1 ip))
459 ch2-%offset%)
461 (f2cl-lib:fref c2-%data%
462 (m1 ik 1)
463 ((1 in1) (1 idl1) (1 ip))
464 c2-%offset%)
465 (* ar1
466 (f2cl-lib:fref c2-%data%
467 (m1 ik 2)
468 ((1 in1) (1 idl1) (1 ip))
469 c2-%offset%))))
470 (setf (f2cl-lib:fref ch2-%data%
471 (m2 ik lc)
472 ((1 in2) (1 idl1) (1 ip))
473 ch2-%offset%)
474 (* ai1
475 (f2cl-lib:fref c2-%data%
476 (m1 ik ip)
477 ((1 in1) (1 idl1) (1 ip))
478 c2-%offset%)))
479 label1024))
480 label124))
481 (setf dc2 ar1)
482 (setf ds2 ai1)
483 (setf ar2 ar1)
484 (setf ai2 ai1)
485 (f2cl-lib:fdo (j 3 (f2cl-lib:int-add j 1))
486 ((> j ipph) nil)
487 (tagbody
488 (setf jc (f2cl-lib:int-sub ipp2 j))
489 (setf ar2h (- (* dc2 ar2) (* ds2 ai2)))
490 (setf ai2 (+ (* dc2 ai2) (* ds2 ar2)))
491 (setf ar2 ar2h)
492 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
493 ((> ik idl1) nil)
494 (tagbody
495 (setf m2 m2s)
496 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
497 ((> m1 m1d) nil)
498 (tagbody
499 (setf m2 (f2cl-lib:int-add m2 im2))
500 (setf (f2cl-lib:fref ch2-%data%
501 (m2 ik l)
502 ((1 in2) (1 idl1) (1 ip))
503 ch2-%offset%)
505 (f2cl-lib:fref ch2-%data%
506 (m2 ik l)
507 ((1 in2) (1 idl1) (1 ip))
508 ch2-%offset%)
509 (* ar2
510 (f2cl-lib:fref c2-%data%
511 (m1 ik j)
512 ((1 in1) (1 idl1) (1 ip))
513 c2-%offset%))))
514 (setf (f2cl-lib:fref ch2-%data%
515 (m2 ik lc)
516 ((1 in2) (1 idl1) (1 ip))
517 ch2-%offset%)
519 (f2cl-lib:fref ch2-%data%
520 (m2 ik lc)
521 ((1 in2) (1 idl1) (1 ip))
522 ch2-%offset%)
523 (* ai2
524 (f2cl-lib:fref c2-%data%
525 (m1 ik jc)
526 ((1 in1) (1 idl1) (1 ip))
527 c2-%offset%))))
528 label1025))
529 label125))
530 label126))
531 label127))
532 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
533 ((> j ipph) nil)
534 (tagbody
535 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
536 ((> ik idl1) nil)
537 (tagbody
538 (setf m2 m2s)
539 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
540 ((> m1 m1d) nil)
541 (tagbody
542 (setf m2 (f2cl-lib:int-add m2 im2))
543 (setf (f2cl-lib:fref ch2-%data%
544 (m2 ik 1)
545 ((1 in2) (1 idl1) (1 ip))
546 ch2-%offset%)
548 (f2cl-lib:fref ch2-%data%
549 (m2 ik 1)
550 ((1 in2) (1 idl1) (1 ip))
551 ch2-%offset%)
552 (f2cl-lib:fref c2-%data%
553 (m1 ik j)
554 ((1 in1) (1 idl1) (1 ip))
555 c2-%offset%)))
556 label1028))
557 label128))
558 label129))
559 (if (< ido l1) (go label132))
560 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
561 ((> k l1) nil)
562 (tagbody
563 (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
564 ((> i ido) nil)
565 (tagbody
566 (setf m2 m2s)
567 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
568 ((> m1 m1d) nil)
569 (tagbody
570 (setf m2 (f2cl-lib:int-add m2 im2))
571 (setf (f2cl-lib:fref cc-%data%
572 (m1 i 1 k)
573 ((1 in1) (1 ido) (1 ip) (1 l1))
574 cc-%offset%)
575 (f2cl-lib:fref ch-%data%
576 (m2 i k 1)
577 ((1 in2) (1 ido) (1 l1) (1 ip))
578 ch-%offset%))
579 label1030))
580 label130))
581 label131))
582 (go label135)
583 label132
584 (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
585 ((> i ido) nil)
586 (tagbody
587 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
588 ((> k l1) nil)
589 (tagbody
590 (setf m2 m2s)
591 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
592 ((> m1 m1d) nil)
593 (tagbody
594 (setf m2 (f2cl-lib:int-add m2 im2))
595 (setf (f2cl-lib:fref cc-%data%
596 (m1 i 1 k)
597 ((1 in1) (1 ido) (1 ip) (1 l1))
598 cc-%offset%)
599 (f2cl-lib:fref ch-%data%
600 (m2 i k 1)
601 ((1 in2) (1 ido) (1 l1) (1 ip))
602 ch-%offset%))
603 label1033))
604 label133))
605 label134))
606 label135
607 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
608 ((> j ipph) nil)
609 (tagbody
610 (setf jc (f2cl-lib:int-sub ipp2 j))
611 (setf j2 (f2cl-lib:int-add j j))
612 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
613 ((> k l1) nil)
614 (tagbody
615 (setf m2 m2s)
616 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
617 ((> m1 m1d) nil)
618 (tagbody
619 (setf m2 (f2cl-lib:int-add m2 im2))
620 (setf (f2cl-lib:fref cc-%data%
621 (m1 ido (f2cl-lib:int-sub j2 2) k)
622 ((1 in1) (1 ido) (1 ip) (1 l1))
623 cc-%offset%)
624 (f2cl-lib:fref ch-%data%
625 (m2 1 k j)
626 ((1 in2) (1 ido) (1 l1) (1 ip))
627 ch-%offset%))
628 (setf (f2cl-lib:fref cc-%data%
629 (m1 1 (f2cl-lib:int-sub j2 1) k)
630 ((1 in1) (1 ido) (1 ip) (1 l1))
631 cc-%offset%)
632 (f2cl-lib:fref ch-%data%
633 (m2 1 k jc)
634 ((1 in2) (1 ido) (1 l1) (1 ip))
635 ch-%offset%))
636 label1036))
637 label136))
638 label137))
639 (if (= ido 1) (go end_label))
640 (if (< nbd l1) (go label141))
641 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
642 ((> j ipph) nil)
643 (tagbody
644 (setf jc (f2cl-lib:int-sub ipp2 j))
645 (setf j2 (f2cl-lib:int-add j j))
646 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
647 ((> k l1) nil)
648 (tagbody
649 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
650 ((> i ido) nil)
651 (tagbody
652 (setf ic (f2cl-lib:int-sub idp2 i))
653 (setf m2 m2s)
654 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
655 ((> m1 m1d) nil)
656 (tagbody
657 (setf m2 (f2cl-lib:int-add m2 im2))
658 (setf (f2cl-lib:fref cc-%data%
659 (m1 (f2cl-lib:int-sub i 1)
660 (f2cl-lib:int-sub j2 1) k)
661 ((1 in1) (1 ido) (1 ip) (1 l1))
662 cc-%offset%)
664 (f2cl-lib:fref ch-%data%
665 (m2 (f2cl-lib:int-sub i 1) k j)
666 ((1 in2) (1 ido) (1 l1) (1 ip))
667 ch-%offset%)
668 (f2cl-lib:fref ch-%data%
669 (m2 (f2cl-lib:int-sub i 1) k jc)
670 ((1 in2) (1 ido) (1 l1) (1 ip))
671 ch-%offset%)))
672 (setf (f2cl-lib:fref cc-%data%
673 (m1 (f2cl-lib:int-sub ic 1)
674 (f2cl-lib:int-sub j2 2) k)
675 ((1 in1) (1 ido) (1 ip) (1 l1))
676 cc-%offset%)
678 (f2cl-lib:fref ch-%data%
679 (m2 (f2cl-lib:int-sub i 1) k j)
680 ((1 in2) (1 ido) (1 l1) (1 ip))
681 ch-%offset%)
682 (f2cl-lib:fref ch-%data%
683 (m2 (f2cl-lib:int-sub i 1) k jc)
684 ((1 in2) (1 ido) (1 l1) (1 ip))
685 ch-%offset%)))
686 (setf (f2cl-lib:fref cc-%data%
687 (m1 i (f2cl-lib:int-sub j2 1) k)
688 ((1 in1) (1 ido) (1 ip) (1 l1))
689 cc-%offset%)
691 (f2cl-lib:fref ch-%data%
692 (m2 i k j)
693 ((1 in2) (1 ido) (1 l1) (1 ip))
694 ch-%offset%)
695 (f2cl-lib:fref ch-%data%
696 (m2 i k jc)
697 ((1 in2) (1 ido) (1 l1) (1 ip))
698 ch-%offset%)))
699 (setf (f2cl-lib:fref cc-%data%
700 (m1 ic (f2cl-lib:int-sub j2 2) k)
701 ((1 in1) (1 ido) (1 ip) (1 l1))
702 cc-%offset%)
704 (f2cl-lib:fref ch-%data%
705 (m2 i k jc)
706 ((1 in2) (1 ido) (1 l1) (1 ip))
707 ch-%offset%)
708 (f2cl-lib:fref ch-%data%
709 (m2 i k j)
710 ((1 in2) (1 ido) (1 l1) (1 ip))
711 ch-%offset%)))
712 label1038))
713 label138))
714 label139))
715 label140))
716 (go end_label)
717 label141
718 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
719 ((> j ipph) nil)
720 (tagbody
721 (setf jc (f2cl-lib:int-sub ipp2 j))
722 (setf j2 (f2cl-lib:int-add j j))
723 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
724 ((> i ido) nil)
725 (tagbody
726 (setf ic (f2cl-lib:int-sub idp2 i))
727 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
728 ((> k l1) nil)
729 (tagbody
730 (setf m2 m2s)
731 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
732 ((> m1 m1d) nil)
733 (tagbody
734 (setf m2 (f2cl-lib:int-add m2 im2))
735 (setf (f2cl-lib:fref cc-%data%
736 (m1 (f2cl-lib:int-sub i 1)
737 (f2cl-lib:int-sub j2 1) k)
738 ((1 in1) (1 ido) (1 ip) (1 l1))
739 cc-%offset%)
741 (f2cl-lib:fref ch-%data%
742 (m2 (f2cl-lib:int-sub i 1) k j)
743 ((1 in2) (1 ido) (1 l1) (1 ip))
744 ch-%offset%)
745 (f2cl-lib:fref ch-%data%
746 (m2 (f2cl-lib:int-sub i 1) k jc)
747 ((1 in2) (1 ido) (1 l1) (1 ip))
748 ch-%offset%)))
749 (setf (f2cl-lib:fref cc-%data%
750 (m1 (f2cl-lib:int-sub ic 1)
751 (f2cl-lib:int-sub j2 2) k)
752 ((1 in1) (1 ido) (1 ip) (1 l1))
753 cc-%offset%)
755 (f2cl-lib:fref ch-%data%
756 (m2 (f2cl-lib:int-sub i 1) k j)
757 ((1 in2) (1 ido) (1 l1) (1 ip))
758 ch-%offset%)
759 (f2cl-lib:fref ch-%data%
760 (m2 (f2cl-lib:int-sub i 1) k jc)
761 ((1 in2) (1 ido) (1 l1) (1 ip))
762 ch-%offset%)))
763 (setf (f2cl-lib:fref cc-%data%
764 (m1 i (f2cl-lib:int-sub j2 1) k)
765 ((1 in1) (1 ido) (1 ip) (1 l1))
766 cc-%offset%)
768 (f2cl-lib:fref ch-%data%
769 (m2 i k j)
770 ((1 in2) (1 ido) (1 l1) (1 ip))
771 ch-%offset%)
772 (f2cl-lib:fref ch-%data%
773 (m2 i k jc)
774 ((1 in2) (1 ido) (1 l1) (1 ip))
775 ch-%offset%)))
776 (setf (f2cl-lib:fref cc-%data%
777 (m1 ic (f2cl-lib:int-sub j2 2) k)
778 ((1 in1) (1 ido) (1 ip) (1 l1))
779 cc-%offset%)
781 (f2cl-lib:fref ch-%data%
782 (m2 i k jc)
783 ((1 in2) (1 ido) (1 l1) (1 ip))
784 ch-%offset%)
785 (f2cl-lib:fref ch-%data%
786 (m2 i k j)
787 ((1 in2) (1 ido) (1 l1) (1 ip))
788 ch-%offset%)))
789 label1042))
790 label142))
791 label143))
792 label144))
793 (go end_label)
794 end_label
795 (return
796 (values nil nil nil nil nil nil nil nil nil nil nil nil nil nil nil)))))
798 (in-package #:cl-user)
799 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
800 (eval-when (:load-toplevel :compile-toplevel :execute)
801 (setf (gethash 'fortran-to-lisp::mradfg
802 fortran-to-lisp::*f2cl-function-info*)
803 (fortran-to-lisp::make-f2cl-finfo
804 :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
805 (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
806 (fortran-to-lisp::integer4) (array double-float (*))
807 (array double-float (*)) (array double-float (*))
808 (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
809 (array double-float (*)) (array double-float (*))
810 (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
811 (array double-float (*)))
812 :return-values '(nil nil nil nil nil nil nil nil nil nil nil nil nil
813 nil nil)
814 :calls 'nil)))