Rename *ll* and *ul* to ll and ul in $defint
[maxima.git] / share / fftpack5 / lisp / mradbg.lisp
blobfb4268c16ba72164d3663fc897370a4d26708682
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 mradbg (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 ((idij 0) (is 0) (ar2h 0.0d0) (ai2 0.0d0) (ar2 0.0d0) (ds2 0.0d0)
31 (dc2 0.0d0) (ik 0) (ar1h 0.0d0) (lc 0) (l 0) (ai1 0.0d0) (ar1 0.0d0)
32 (ic 0) (j2 0) (jc 0) (j 0) (m1 0) (m2 0) (i 0) (k 0) (ipph 0)
33 (ipp2 0) (nbd 0) (idp2 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 idp2 nbd ipp2 ipph k i m2 m1 j
38 jc j2 ic l lc ik is idij))
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 idp2 (f2cl-lib:int-add ido 2))
48 (setf nbd (the f2cl-lib:integer4 (truncate (- ido 1) 2)))
49 (setf ipp2 (f2cl-lib:int-add ip 2))
50 (setf ipph (the f2cl-lib:integer4 (truncate (+ ip 1) 2)))
51 (if (< ido l1) (go label103))
52 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
53 ((> k l1) nil)
54 (tagbody
55 (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
56 ((> i ido) nil)
57 (tagbody
58 (setf m2 m2s)
59 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
60 ((> m1 m1d) nil)
61 (tagbody
62 (setf m2 (f2cl-lib:int-add m2 im2))
63 (setf (f2cl-lib:fref ch-%data%
64 (m2 i k 1)
65 ((1 in2) (1 ido) (1 l1) (1 ip))
66 ch-%offset%)
67 (f2cl-lib:fref cc-%data%
68 (m1 i 1 k)
69 ((1 in1) (1 ido) (1 ip) (1 l1))
70 cc-%offset%))
71 label1001))
72 label101))
73 label102))
74 (go label106)
75 label103
76 (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
77 ((> i ido) nil)
78 (tagbody
79 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
80 ((> k l1) nil)
81 (tagbody
82 (setf m2 m2s)
83 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
84 ((> m1 m1d) nil)
85 (tagbody
86 (setf m2 (f2cl-lib:int-add m2 im2))
87 (setf (f2cl-lib:fref ch-%data%
88 (m2 i k 1)
89 ((1 in2) (1 ido) (1 l1) (1 ip))
90 ch-%offset%)
91 (f2cl-lib:fref cc-%data%
92 (m1 i 1 k)
93 ((1 in1) (1 ido) (1 ip) (1 l1))
94 cc-%offset%))
95 label1004))
96 label104))
97 label105))
98 label106
99 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
100 ((> j ipph) nil)
101 (tagbody
102 (setf jc (f2cl-lib:int-sub ipp2 j))
103 (setf j2 (f2cl-lib:int-add j j))
104 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
105 ((> k l1) nil)
106 (tagbody
107 (setf m2 m2s)
108 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
109 ((> m1 m1d) nil)
110 (tagbody
111 (setf m2 (f2cl-lib:int-add m2 im2))
112 (setf (f2cl-lib:fref ch-%data%
113 (m2 1 k j)
114 ((1 in2) (1 ido) (1 l1) (1 ip))
115 ch-%offset%)
117 (f2cl-lib:fref cc-%data%
118 (m1 ido (f2cl-lib:int-sub j2 2) k)
119 ((1 in1) (1 ido) (1 ip) (1 l1))
120 cc-%offset%)
121 (f2cl-lib:fref cc-%data%
122 (m1 ido (f2cl-lib:int-sub j2 2) k)
123 ((1 in1) (1 ido) (1 ip) (1 l1))
124 cc-%offset%)))
125 (setf (f2cl-lib:fref ch-%data%
126 (m2 1 k jc)
127 ((1 in2) (1 ido) (1 l1) (1 ip))
128 ch-%offset%)
130 (f2cl-lib:fref cc-%data%
131 (m1 1 (f2cl-lib:int-sub j2 1) k)
132 ((1 in1) (1 ido) (1 ip) (1 l1))
133 cc-%offset%)
134 (f2cl-lib:fref cc-%data%
135 (m1 1 (f2cl-lib:int-sub j2 1) k)
136 ((1 in1) (1 ido) (1 ip) (1 l1))
137 cc-%offset%)))
138 label1007))
139 label107))
140 label108))
141 (if (= ido 1) (go label116))
142 (if (< nbd l1) (go label112))
143 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
144 ((> j ipph) nil)
145 (tagbody
146 (setf jc (f2cl-lib:int-sub ipp2 j))
147 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
148 ((> k l1) nil)
149 (tagbody
150 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
151 ((> i ido) nil)
152 (tagbody
153 (setf ic (f2cl-lib:int-sub idp2 i))
154 (setf m2 m2s)
155 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
156 ((> m1 m1d) nil)
157 (tagbody
158 (setf m2 (f2cl-lib:int-add m2 im2))
159 (setf (f2cl-lib:fref ch-%data%
160 (m2 (f2cl-lib:int-sub i 1) k j)
161 ((1 in2) (1 ido) (1 l1) (1 ip))
162 ch-%offset%)
164 (f2cl-lib:fref cc-%data%
165 (m1 (f2cl-lib:int-sub i 1)
166 (f2cl-lib:int-sub
167 (f2cl-lib:int-mul 2 j)
170 ((1 in1) (1 ido) (1 ip) (1 l1))
171 cc-%offset%)
172 (f2cl-lib:fref cc-%data%
173 (m1 (f2cl-lib:int-sub ic 1)
174 (f2cl-lib:int-sub
175 (f2cl-lib:int-mul 2 j)
178 ((1 in1) (1 ido) (1 ip) (1 l1))
179 cc-%offset%)))
180 (setf (f2cl-lib:fref ch-%data%
181 (m2 (f2cl-lib:int-sub i 1) k jc)
182 ((1 in2) (1 ido) (1 l1) (1 ip))
183 ch-%offset%)
185 (f2cl-lib:fref cc-%data%
186 (m1 (f2cl-lib:int-sub i 1)
187 (f2cl-lib:int-sub
188 (f2cl-lib:int-mul 2 j)
191 ((1 in1) (1 ido) (1 ip) (1 l1))
192 cc-%offset%)
193 (f2cl-lib:fref cc-%data%
194 (m1 (f2cl-lib:int-sub ic 1)
195 (f2cl-lib:int-sub
196 (f2cl-lib:int-mul 2 j)
199 ((1 in1) (1 ido) (1 ip) (1 l1))
200 cc-%offset%)))
201 (setf (f2cl-lib:fref ch-%data%
202 (m2 i k j)
203 ((1 in2) (1 ido) (1 l1) (1 ip))
204 ch-%offset%)
206 (f2cl-lib:fref cc-%data%
207 (m1 i
208 (f2cl-lib:int-sub
209 (f2cl-lib:int-mul 2 j)
212 ((1 in1) (1 ido) (1 ip) (1 l1))
213 cc-%offset%)
214 (f2cl-lib:fref cc-%data%
215 (m1 ic
216 (f2cl-lib:int-sub
217 (f2cl-lib:int-mul 2 j)
220 ((1 in1) (1 ido) (1 ip) (1 l1))
221 cc-%offset%)))
222 (setf (f2cl-lib:fref ch-%data%
223 (m2 i k jc)
224 ((1 in2) (1 ido) (1 l1) (1 ip))
225 ch-%offset%)
227 (f2cl-lib:fref cc-%data%
228 (m1 i
229 (f2cl-lib:int-sub
230 (f2cl-lib:int-mul 2 j)
233 ((1 in1) (1 ido) (1 ip) (1 l1))
234 cc-%offset%)
235 (f2cl-lib:fref cc-%data%
236 (m1 ic
237 (f2cl-lib:int-sub
238 (f2cl-lib:int-mul 2 j)
241 ((1 in1) (1 ido) (1 ip) (1 l1))
242 cc-%offset%)))
243 label1009))
244 label109))
245 label110))
246 label111))
247 (go label116)
248 label112
249 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
250 ((> j ipph) nil)
251 (tagbody
252 (setf jc (f2cl-lib:int-sub ipp2 j))
253 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
254 ((> i ido) nil)
255 (tagbody
256 (setf ic (f2cl-lib:int-sub idp2 i))
257 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
258 ((> k l1) nil)
259 (tagbody
260 (setf m2 m2s)
261 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
262 ((> m1 m1d) nil)
263 (tagbody
264 (setf m2 (f2cl-lib:int-add m2 im2))
265 (setf (f2cl-lib:fref ch-%data%
266 (m2 (f2cl-lib:int-sub i 1) k j)
267 ((1 in2) (1 ido) (1 l1) (1 ip))
268 ch-%offset%)
270 (f2cl-lib:fref cc-%data%
271 (m1 (f2cl-lib:int-sub i 1)
272 (f2cl-lib:int-sub
273 (f2cl-lib:int-mul 2 j)
276 ((1 in1) (1 ido) (1 ip) (1 l1))
277 cc-%offset%)
278 (f2cl-lib:fref cc-%data%
279 (m1 (f2cl-lib:int-sub ic 1)
280 (f2cl-lib:int-sub
281 (f2cl-lib:int-mul 2 j)
284 ((1 in1) (1 ido) (1 ip) (1 l1))
285 cc-%offset%)))
286 (setf (f2cl-lib:fref ch-%data%
287 (m2 (f2cl-lib:int-sub i 1) k jc)
288 ((1 in2) (1 ido) (1 l1) (1 ip))
289 ch-%offset%)
291 (f2cl-lib:fref cc-%data%
292 (m1 (f2cl-lib:int-sub i 1)
293 (f2cl-lib:int-sub
294 (f2cl-lib:int-mul 2 j)
297 ((1 in1) (1 ido) (1 ip) (1 l1))
298 cc-%offset%)
299 (f2cl-lib:fref cc-%data%
300 (m1 (f2cl-lib:int-sub ic 1)
301 (f2cl-lib:int-sub
302 (f2cl-lib:int-mul 2 j)
305 ((1 in1) (1 ido) (1 ip) (1 l1))
306 cc-%offset%)))
307 (setf (f2cl-lib:fref ch-%data%
308 (m2 i k j)
309 ((1 in2) (1 ido) (1 l1) (1 ip))
310 ch-%offset%)
312 (f2cl-lib:fref cc-%data%
313 (m1 i
314 (f2cl-lib:int-sub
315 (f2cl-lib:int-mul 2 j)
318 ((1 in1) (1 ido) (1 ip) (1 l1))
319 cc-%offset%)
320 (f2cl-lib:fref cc-%data%
321 (m1 ic
322 (f2cl-lib:int-sub
323 (f2cl-lib:int-mul 2 j)
326 ((1 in1) (1 ido) (1 ip) (1 l1))
327 cc-%offset%)))
328 (setf (f2cl-lib:fref ch-%data%
329 (m2 i k jc)
330 ((1 in2) (1 ido) (1 l1) (1 ip))
331 ch-%offset%)
333 (f2cl-lib:fref cc-%data%
334 (m1 i
335 (f2cl-lib:int-sub
336 (f2cl-lib:int-mul 2 j)
339 ((1 in1) (1 ido) (1 ip) (1 l1))
340 cc-%offset%)
341 (f2cl-lib:fref cc-%data%
342 (m1 ic
343 (f2cl-lib:int-sub
344 (f2cl-lib:int-mul 2 j)
347 ((1 in1) (1 ido) (1 ip) (1 l1))
348 cc-%offset%)))
349 label1013))
350 label113))
351 label114))
352 label115))
353 label116
354 (setf ar1 1.0d0)
355 (setf ai1 0.0d0)
356 (f2cl-lib:fdo (l 2 (f2cl-lib:int-add l 1))
357 ((> l ipph) nil)
358 (tagbody
359 (setf lc (f2cl-lib:int-sub ipp2 l))
360 (setf ar1h (- (* dcp ar1) (* dsp ai1)))
361 (setf ai1 (+ (* dcp ai1) (* dsp ar1)))
362 (setf ar1 ar1h)
363 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
364 ((> ik idl1) nil)
365 (tagbody
366 (setf m2 m2s)
367 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
368 ((> m1 m1d) nil)
369 (tagbody
370 (setf m2 (f2cl-lib:int-add m2 im2))
371 (setf (f2cl-lib:fref c2-%data%
372 (m1 ik l)
373 ((1 in1) (1 idl1) (1 ip))
374 c2-%offset%)
376 (f2cl-lib:fref ch2-%data%
377 (m2 ik 1)
378 ((1 in2) (1 idl1) (1 ip))
379 ch2-%offset%)
380 (* ar1
381 (f2cl-lib:fref ch2-%data%
382 (m2 ik 2)
383 ((1 in2) (1 idl1) (1 ip))
384 ch2-%offset%))))
385 (setf (f2cl-lib:fref c2-%data%
386 (m1 ik lc)
387 ((1 in1) (1 idl1) (1 ip))
388 c2-%offset%)
389 (* ai1
390 (f2cl-lib:fref ch2-%data%
391 (m2 ik ip)
392 ((1 in2) (1 idl1) (1 ip))
393 ch2-%offset%)))
394 label1017))
395 label117))
396 (setf dc2 ar1)
397 (setf ds2 ai1)
398 (setf ar2 ar1)
399 (setf ai2 ai1)
400 (f2cl-lib:fdo (j 3 (f2cl-lib:int-add j 1))
401 ((> j ipph) nil)
402 (tagbody
403 (setf jc (f2cl-lib:int-sub ipp2 j))
404 (setf ar2h (- (* dc2 ar2) (* ds2 ai2)))
405 (setf ai2 (+ (* dc2 ai2) (* ds2 ar2)))
406 (setf ar2 ar2h)
407 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
408 ((> ik idl1) nil)
409 (tagbody
410 (setf m2 m2s)
411 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
412 ((> m1 m1d) nil)
413 (tagbody
414 (setf m2 (f2cl-lib:int-add m2 im2))
415 (setf (f2cl-lib:fref c2-%data%
416 (m1 ik l)
417 ((1 in1) (1 idl1) (1 ip))
418 c2-%offset%)
420 (f2cl-lib:fref c2-%data%
421 (m1 ik l)
422 ((1 in1) (1 idl1) (1 ip))
423 c2-%offset%)
424 (* ar2
425 (f2cl-lib:fref ch2-%data%
426 (m2 ik j)
427 ((1 in2) (1 idl1) (1 ip))
428 ch2-%offset%))))
429 (setf (f2cl-lib:fref c2-%data%
430 (m1 ik lc)
431 ((1 in1) (1 idl1) (1 ip))
432 c2-%offset%)
434 (f2cl-lib:fref c2-%data%
435 (m1 ik lc)
436 ((1 in1) (1 idl1) (1 ip))
437 c2-%offset%)
438 (* ai2
439 (f2cl-lib:fref ch2-%data%
440 (m2 ik jc)
441 ((1 in2) (1 idl1) (1 ip))
442 ch2-%offset%))))
443 label1018))
444 label118))
445 label119))
446 label120))
447 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
448 ((> j ipph) nil)
449 (tagbody
450 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
451 ((> ik idl1) nil)
452 (tagbody
453 (setf m2 m2s)
454 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
455 ((> m1 m1d) nil)
456 (tagbody
457 (setf m2 (f2cl-lib:int-add m2 im2))
458 (setf (f2cl-lib:fref ch2-%data%
459 (m2 ik 1)
460 ((1 in2) (1 idl1) (1 ip))
461 ch2-%offset%)
463 (f2cl-lib:fref ch2-%data%
464 (m2 ik 1)
465 ((1 in2) (1 idl1) (1 ip))
466 ch2-%offset%)
467 (f2cl-lib:fref ch2-%data%
468 (m2 ik j)
469 ((1 in2) (1 idl1) (1 ip))
470 ch2-%offset%)))
471 label1021))
472 label121))
473 label122))
474 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
475 ((> j ipph) nil)
476 (tagbody
477 (setf jc (f2cl-lib:int-sub ipp2 j))
478 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
479 ((> k l1) nil)
480 (tagbody
481 (setf m2 m2s)
482 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
483 ((> m1 m1d) nil)
484 (tagbody
485 (setf m2 (f2cl-lib:int-add m2 im2))
486 (setf (f2cl-lib:fref ch-%data%
487 (m2 1 k j)
488 ((1 in2) (1 ido) (1 l1) (1 ip))
489 ch-%offset%)
491 (f2cl-lib:fref c1-%data%
492 (m1 1 k j)
493 ((1 in1) (1 ido) (1 l1) (1 ip))
494 c1-%offset%)
495 (f2cl-lib:fref c1-%data%
496 (m1 1 k jc)
497 ((1 in1) (1 ido) (1 l1) (1 ip))
498 c1-%offset%)))
499 (setf (f2cl-lib:fref ch-%data%
500 (m2 1 k jc)
501 ((1 in2) (1 ido) (1 l1) (1 ip))
502 ch-%offset%)
504 (f2cl-lib:fref c1-%data%
505 (m1 1 k j)
506 ((1 in1) (1 ido) (1 l1) (1 ip))
507 c1-%offset%)
508 (f2cl-lib:fref c1-%data%
509 (m1 1 k jc)
510 ((1 in1) (1 ido) (1 l1) (1 ip))
511 c1-%offset%)))
512 label1023))
513 label123))
514 label124))
515 (if (= ido 1) (go label132))
516 (if (< nbd l1) (go label128))
517 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
518 ((> j ipph) nil)
519 (tagbody
520 (setf jc (f2cl-lib:int-sub ipp2 j))
521 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
522 ((> k l1) nil)
523 (tagbody
524 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
525 ((> i ido) nil)
526 (tagbody
527 (setf m2 m2s)
528 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
529 ((> m1 m1d) nil)
530 (tagbody
531 (setf m2 (f2cl-lib:int-add m2 im2))
532 (setf (f2cl-lib:fref ch-%data%
533 (m2 (f2cl-lib:int-sub i 1) k j)
534 ((1 in2) (1 ido) (1 l1) (1 ip))
535 ch-%offset%)
537 (f2cl-lib:fref c1-%data%
538 (m1 (f2cl-lib:int-sub i 1) k j)
539 ((1 in1) (1 ido) (1 l1) (1 ip))
540 c1-%offset%)
541 (f2cl-lib:fref c1-%data%
542 (m1 i k jc)
543 ((1 in1) (1 ido) (1 l1) (1 ip))
544 c1-%offset%)))
545 (setf (f2cl-lib:fref ch-%data%
546 (m2 (f2cl-lib:int-sub i 1) k jc)
547 ((1 in2) (1 ido) (1 l1) (1 ip))
548 ch-%offset%)
550 (f2cl-lib:fref c1-%data%
551 (m1 (f2cl-lib:int-sub i 1) k j)
552 ((1 in1) (1 ido) (1 l1) (1 ip))
553 c1-%offset%)
554 (f2cl-lib:fref c1-%data%
555 (m1 i k jc)
556 ((1 in1) (1 ido) (1 l1) (1 ip))
557 c1-%offset%)))
558 (setf (f2cl-lib:fref ch-%data%
559 (m2 i k j)
560 ((1 in2) (1 ido) (1 l1) (1 ip))
561 ch-%offset%)
563 (f2cl-lib:fref c1-%data%
564 (m1 i k j)
565 ((1 in1) (1 ido) (1 l1) (1 ip))
566 c1-%offset%)
567 (f2cl-lib:fref c1-%data%
568 (m1 (f2cl-lib:int-sub i 1) k jc)
569 ((1 in1) (1 ido) (1 l1) (1 ip))
570 c1-%offset%)))
571 (setf (f2cl-lib:fref ch-%data%
572 (m2 i k jc)
573 ((1 in2) (1 ido) (1 l1) (1 ip))
574 ch-%offset%)
576 (f2cl-lib:fref c1-%data%
577 (m1 i k j)
578 ((1 in1) (1 ido) (1 l1) (1 ip))
579 c1-%offset%)
580 (f2cl-lib:fref c1-%data%
581 (m1 (f2cl-lib:int-sub i 1) k jc)
582 ((1 in1) (1 ido) (1 l1) (1 ip))
583 c1-%offset%)))
584 label1025))
585 label125))
586 label126))
587 label127))
588 (go label132)
589 label128
590 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
591 ((> j ipph) nil)
592 (tagbody
593 (setf jc (f2cl-lib:int-sub ipp2 j))
594 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
595 ((> i ido) nil)
596 (tagbody
597 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
598 ((> k l1) nil)
599 (tagbody
600 (setf m2 m2s)
601 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
602 ((> m1 m1d) nil)
603 (tagbody
604 (setf m2 (f2cl-lib:int-add m2 im2))
605 (setf (f2cl-lib:fref ch-%data%
606 (m2 (f2cl-lib:int-sub i 1) k j)
607 ((1 in2) (1 ido) (1 l1) (1 ip))
608 ch-%offset%)
610 (f2cl-lib:fref c1-%data%
611 (m1 (f2cl-lib:int-sub i 1) k j)
612 ((1 in1) (1 ido) (1 l1) (1 ip))
613 c1-%offset%)
614 (f2cl-lib:fref c1-%data%
615 (m1 i k jc)
616 ((1 in1) (1 ido) (1 l1) (1 ip))
617 c1-%offset%)))
618 (setf (f2cl-lib:fref ch-%data%
619 (m2 (f2cl-lib:int-sub i 1) k jc)
620 ((1 in2) (1 ido) (1 l1) (1 ip))
621 ch-%offset%)
623 (f2cl-lib:fref c1-%data%
624 (m1 (f2cl-lib:int-sub i 1) k j)
625 ((1 in1) (1 ido) (1 l1) (1 ip))
626 c1-%offset%)
627 (f2cl-lib:fref c1-%data%
628 (m1 i k jc)
629 ((1 in1) (1 ido) (1 l1) (1 ip))
630 c1-%offset%)))
631 (setf (f2cl-lib:fref ch-%data%
632 (m2 i k j)
633 ((1 in2) (1 ido) (1 l1) (1 ip))
634 ch-%offset%)
636 (f2cl-lib:fref c1-%data%
637 (m1 i k j)
638 ((1 in1) (1 ido) (1 l1) (1 ip))
639 c1-%offset%)
640 (f2cl-lib:fref c1-%data%
641 (m1 (f2cl-lib:int-sub i 1) k jc)
642 ((1 in1) (1 ido) (1 l1) (1 ip))
643 c1-%offset%)))
644 (setf (f2cl-lib:fref ch-%data%
645 (m2 i k jc)
646 ((1 in2) (1 ido) (1 l1) (1 ip))
647 ch-%offset%)
649 (f2cl-lib:fref c1-%data%
650 (m1 i k j)
651 ((1 in1) (1 ido) (1 l1) (1 ip))
652 c1-%offset%)
653 (f2cl-lib:fref c1-%data%
654 (m1 (f2cl-lib:int-sub i 1) k jc)
655 ((1 in1) (1 ido) (1 l1) (1 ip))
656 c1-%offset%)))
657 label1029))
658 label129))
659 label130))
660 label131))
661 label132
662 (if (= ido 1) (go end_label))
663 (f2cl-lib:fdo (ik 1 (f2cl-lib:int-add ik 1))
664 ((> ik idl1) nil)
665 (tagbody
666 (setf m2 m2s)
667 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
668 ((> m1 m1d) nil)
669 (tagbody
670 (setf m2 (f2cl-lib:int-add m2 im2))
671 (setf (f2cl-lib:fref c2-%data%
672 (m1 ik 1)
673 ((1 in1) (1 idl1) (1 ip))
674 c2-%offset%)
675 (f2cl-lib:fref ch2-%data%
676 (m2 ik 1)
677 ((1 in2) (1 idl1) (1 ip))
678 ch2-%offset%))
679 label1033))
680 label133))
681 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
682 ((> j ip) nil)
683 (tagbody
684 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
685 ((> k l1) nil)
686 (tagbody
687 (setf m2 m2s)
688 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
689 ((> m1 m1d) nil)
690 (tagbody
691 (setf m2 (f2cl-lib:int-add m2 im2))
692 (setf (f2cl-lib:fref c1-%data%
693 (m1 1 k j)
694 ((1 in1) (1 ido) (1 l1) (1 ip))
695 c1-%offset%)
696 (f2cl-lib:fref ch-%data%
697 (m2 1 k j)
698 ((1 in2) (1 ido) (1 l1) (1 ip))
699 ch-%offset%))
700 label1034))
701 label134))
702 label135))
703 (if (> nbd l1) (go label139))
704 (setf is (f2cl-lib:int-sub ido))
705 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
706 ((> j ip) nil)
707 (tagbody
708 (setf is (f2cl-lib:int-add is ido))
709 (setf idij is)
710 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
711 ((> i ido) nil)
712 (tagbody
713 (setf idij (f2cl-lib:int-add idij 2))
714 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
715 ((> k l1) nil)
716 (tagbody
717 (setf m2 m2s)
718 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
719 ((> m1 m1d) nil)
720 (tagbody
721 (setf m2 (f2cl-lib:int-add m2 im2))
722 (setf (f2cl-lib:fref c1-%data%
723 (m1 (f2cl-lib:int-sub i 1) k j)
724 ((1 in1) (1 ido) (1 l1) (1 ip))
725 c1-%offset%)
728 (f2cl-lib:fref wa-%data%
729 ((f2cl-lib:int-sub idij 1))
730 ((1 ido))
731 wa-%offset%)
732 (f2cl-lib:fref ch-%data%
733 (m2 (f2cl-lib:int-sub i 1) k j)
734 ((1 in2) (1 ido) (1 l1) (1 ip))
735 ch-%offset%))
737 (f2cl-lib:fref wa-%data%
738 (idij)
739 ((1 ido))
740 wa-%offset%)
741 (f2cl-lib:fref ch-%data%
742 (m2 i k j)
743 ((1 in2) (1 ido) (1 l1) (1 ip))
744 ch-%offset%))))
745 (setf (f2cl-lib:fref c1-%data%
746 (m1 i k j)
747 ((1 in1) (1 ido) (1 l1) (1 ip))
748 c1-%offset%)
751 (f2cl-lib:fref wa-%data%
752 ((f2cl-lib:int-sub idij 1))
753 ((1 ido))
754 wa-%offset%)
755 (f2cl-lib:fref ch-%data%
756 (m2 i k j)
757 ((1 in2) (1 ido) (1 l1) (1 ip))
758 ch-%offset%))
760 (f2cl-lib:fref wa-%data%
761 (idij)
762 ((1 ido))
763 wa-%offset%)
764 (f2cl-lib:fref ch-%data%
765 (m2 (f2cl-lib:int-sub i 1) k j)
766 ((1 in2) (1 ido) (1 l1) (1 ip))
767 ch-%offset%))))
768 label1036))
769 label136))
770 label137))
771 label138))
772 (go label143)
773 label139
774 (setf is (f2cl-lib:int-sub ido))
775 (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
776 ((> j ip) nil)
777 (tagbody
778 (setf is (f2cl-lib:int-add is ido))
779 (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
780 ((> k l1) nil)
781 (tagbody
782 (setf idij is)
783 (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
784 ((> i ido) nil)
785 (tagbody
786 (setf idij (f2cl-lib:int-add idij 2))
787 (setf m2 m2s)
788 (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
789 ((> m1 m1d) nil)
790 (tagbody
791 (setf m2 (f2cl-lib:int-add m2 im2))
792 (setf (f2cl-lib:fref c1-%data%
793 (m1 (f2cl-lib:int-sub i 1) k j)
794 ((1 in1) (1 ido) (1 l1) (1 ip))
795 c1-%offset%)
798 (f2cl-lib:fref wa-%data%
799 ((f2cl-lib:int-sub idij 1))
800 ((1 ido))
801 wa-%offset%)
802 (f2cl-lib:fref ch-%data%
803 (m2 (f2cl-lib:int-sub i 1) k j)
804 ((1 in2) (1 ido) (1 l1) (1 ip))
805 ch-%offset%))
807 (f2cl-lib:fref wa-%data%
808 (idij)
809 ((1 ido))
810 wa-%offset%)
811 (f2cl-lib:fref ch-%data%
812 (m2 i k j)
813 ((1 in2) (1 ido) (1 l1) (1 ip))
814 ch-%offset%))))
815 (setf (f2cl-lib:fref c1-%data%
816 (m1 i k j)
817 ((1 in1) (1 ido) (1 l1) (1 ip))
818 c1-%offset%)
821 (f2cl-lib:fref wa-%data%
822 ((f2cl-lib:int-sub idij 1))
823 ((1 ido))
824 wa-%offset%)
825 (f2cl-lib:fref ch-%data%
826 (m2 i k j)
827 ((1 in2) (1 ido) (1 l1) (1 ip))
828 ch-%offset%))
830 (f2cl-lib:fref wa-%data%
831 (idij)
832 ((1 ido))
833 wa-%offset%)
834 (f2cl-lib:fref ch-%data%
835 (m2 (f2cl-lib:int-sub i 1) k j)
836 ((1 in2) (1 ido) (1 l1) (1 ip))
837 ch-%offset%))))
838 label1040))
839 label140))
840 label141))
841 label142))
842 label143
843 (go end_label)
844 end_label
845 (return
846 (values nil nil nil nil nil nil nil nil nil nil nil nil nil nil nil)))))
848 (in-package #:cl-user)
849 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
850 (eval-when (:load-toplevel :compile-toplevel :execute)
851 (setf (gethash 'fortran-to-lisp::mradbg
852 fortran-to-lisp::*f2cl-function-info*)
853 (fortran-to-lisp::make-f2cl-finfo
854 :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
855 (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
856 (fortran-to-lisp::integer4) (array double-float (*))
857 (array double-float (*)) (array double-float (*))
858 (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
859 (array double-float (*)) (array double-float (*))
860 (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
861 (array double-float (*)))
862 :return-values '(nil nil nil nil nil nil nil nil nil nil nil nil nil
863 nil nil)
864 :calls 'nil)))