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)
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 (let ((tr11 0.3090169943749474d0
)
21 (ti11 0.9510565162951536d0
)
22 (tr12 -
0.8090169943749475d0
)
23 (ti12 0.5877852522924731d0
))
24 (declare (type (double-float) tr11 ti11 tr12 ti12
))
25 (defun c1f5kb (ido l1 na cc in1 ch in2 wa
)
26 (declare (type (array double-float
(*)) wa ch cc
)
27 (type (f2cl-lib:integer4
) in2 in1 na l1 ido
))
28 (f2cl-lib:with-multi-array-data
29 ((cc double-float cc-%data% cc-%offset%
)
30 (ch double-float ch-%data% ch-%offset%
)
31 (wa double-float wa-%data% wa-%offset%
))
32 (prog ((di2 0.0d0
) (di5 0.0d0
) (dr2 0.0d0
) (dr5 0.0d0
) (di4 0.0d0
)
33 (di3 0.0d0
) (dr4 0.0d0
) (dr3 0.0d0
) (i 0) (ci4 0.0d0
) (cr4 0.0d0
)
34 (ci5 0.0d0
) (cr5 0.0d0
) (ci3 0.0d0
) (cr3 0.0d0
) (ci2 0.0d0
)
35 (cr2 0.0d0
) (chold2 0.0d0
) (chold1 0.0d0
) (tr3 0.0d0
) (tr4 0.0d0
)
36 (tr2 0.0d0
) (tr5 0.0d0
) (ti3 0.0d0
) (ti4 0.0d0
) (ti2 0.0d0
)
38 (declare (type (f2cl-lib:integer4
) k i
)
39 (type (double-float) ti5 ti2 ti4 ti3 tr5 tr2 tr4 tr3 chold1
40 chold2 cr2 ci2 cr3 ci3 cr5 ci5 cr4 ci4
41 dr3 dr4 di3 di4 dr5 dr2 di5 di2
))
42 (if (or (> ido
1) (= na
1)) (go label102
))
43 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
48 (f2cl-lib:fref cc-%data%
50 ((1 in1
) (1 l1
) (1 ido
) (1 5))
52 (f2cl-lib:fref cc-%data%
54 ((1 in1
) (1 l1
) (1 ido
) (1 5))
58 (f2cl-lib:fref cc-%data%
60 ((1 in1
) (1 l1
) (1 ido
) (1 5))
62 (f2cl-lib:fref cc-%data%
64 ((1 in1
) (1 l1
) (1 ido
) (1 5))
68 (f2cl-lib:fref cc-%data%
70 ((1 in1
) (1 l1
) (1 ido
) (1 5))
72 (f2cl-lib:fref cc-%data%
74 ((1 in1
) (1 l1
) (1 ido
) (1 5))
78 (f2cl-lib:fref cc-%data%
80 ((1 in1
) (1 l1
) (1 ido
) (1 5))
82 (f2cl-lib:fref cc-%data%
84 ((1 in1
) (1 l1
) (1 ido
) (1 5))
88 (f2cl-lib:fref cc-%data%
90 ((1 in1
) (1 l1
) (1 ido
) (1 5))
92 (f2cl-lib:fref cc-%data%
94 ((1 in1
) (1 l1
) (1 ido
) (1 5))
98 (f2cl-lib:fref cc-%data%
100 ((1 in1
) (1 l1
) (1 ido
) (1 5))
102 (f2cl-lib:fref cc-%data%
104 ((1 in1
) (1 l1
) (1 ido
) (1 5))
108 (f2cl-lib:fref cc-%data%
110 ((1 in1
) (1 l1
) (1 ido
) (1 5))
112 (f2cl-lib:fref cc-%data%
114 ((1 in1
) (1 l1
) (1 ido
) (1 5))
118 (f2cl-lib:fref cc-%data%
120 ((1 in1
) (1 l1
) (1 ido
) (1 5))
122 (f2cl-lib:fref cc-%data%
124 ((1 in1
) (1 l1
) (1 ido
) (1 5))
128 (f2cl-lib:fref cc-%data%
130 ((1 in1
) (1 l1
) (1 ido
) (1 5))
136 (f2cl-lib:fref cc-%data%
138 ((1 in1
) (1 l1
) (1 ido
) (1 5))
144 (f2cl-lib:fref cc-%data%
146 ((1 in1
) (1 l1
) (1 ido
) (1 5))
152 (f2cl-lib:fref cc-%data%
154 ((1 in1
) (1 l1
) (1 ido
) (1 5))
160 (f2cl-lib:fref cc-%data%
162 ((1 in1
) (1 l1
) (1 ido
) (1 5))
168 (f2cl-lib:fref cc-%data%
170 ((1 in1
) (1 l1
) (1 ido
) (1 5))
174 (setf (f2cl-lib:fref cc-%data%
176 ((1 in1
) (1 l1
) (1 ido
) (1 5))
179 (setf (f2cl-lib:fref cc-%data%
181 ((1 in1
) (1 l1
) (1 ido
) (1 5))
184 (setf cr5
(+ (* ti11 tr5
) (* ti12 tr4
)))
185 (setf ci5
(+ (* ti11 ti5
) (* ti12 ti4
)))
186 (setf cr4
(- (* ti12 tr5
) (* ti11 tr4
)))
187 (setf ci4
(- (* ti12 ti5
) (* ti11 ti4
)))
188 (setf (f2cl-lib:fref cc-%data%
190 ((1 in1
) (1 l1
) (1 ido
) (1 5))
193 (setf (f2cl-lib:fref cc-%data%
195 ((1 in1
) (1 l1
) (1 ido
) (1 5))
198 (setf (f2cl-lib:fref cc-%data%
200 ((1 in1
) (1 l1
) (1 ido
) (1 5))
203 (setf (f2cl-lib:fref cc-%data%
205 ((1 in1
) (1 l1
) (1 ido
) (1 5))
208 (setf (f2cl-lib:fref cc-%data%
210 ((1 in1
) (1 l1
) (1 ido
) (1 5))
213 (setf (f2cl-lib:fref cc-%data%
215 ((1 in1
) (1 l1
) (1 ido
) (1 5))
218 (setf (f2cl-lib:fref cc-%data%
220 ((1 in1
) (1 l1
) (1 ido
) (1 5))
223 (setf (f2cl-lib:fref cc-%data%
225 ((1 in1
) (1 l1
) (1 ido
) (1 5))
231 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
236 (f2cl-lib:fref cc-%data%
238 ((1 in1
) (1 l1
) (1 ido
) (1 5))
240 (f2cl-lib:fref cc-%data%
242 ((1 in1
) (1 l1
) (1 ido
) (1 5))
246 (f2cl-lib:fref cc-%data%
248 ((1 in1
) (1 l1
) (1 ido
) (1 5))
250 (f2cl-lib:fref cc-%data%
252 ((1 in1
) (1 l1
) (1 ido
) (1 5))
256 (f2cl-lib:fref cc-%data%
258 ((1 in1
) (1 l1
) (1 ido
) (1 5))
260 (f2cl-lib:fref cc-%data%
262 ((1 in1
) (1 l1
) (1 ido
) (1 5))
266 (f2cl-lib:fref cc-%data%
268 ((1 in1
) (1 l1
) (1 ido
) (1 5))
270 (f2cl-lib:fref cc-%data%
272 ((1 in1
) (1 l1
) (1 ido
) (1 5))
276 (f2cl-lib:fref cc-%data%
278 ((1 in1
) (1 l1
) (1 ido
) (1 5))
280 (f2cl-lib:fref cc-%data%
282 ((1 in1
) (1 l1
) (1 ido
) (1 5))
286 (f2cl-lib:fref cc-%data%
288 ((1 in1
) (1 l1
) (1 ido
) (1 5))
290 (f2cl-lib:fref cc-%data%
292 ((1 in1
) (1 l1
) (1 ido
) (1 5))
296 (f2cl-lib:fref cc-%data%
298 ((1 in1
) (1 l1
) (1 ido
) (1 5))
300 (f2cl-lib:fref cc-%data%
302 ((1 in1
) (1 l1
) (1 ido
) (1 5))
306 (f2cl-lib:fref cc-%data%
308 ((1 in1
) (1 l1
) (1 ido
) (1 5))
310 (f2cl-lib:fref cc-%data%
312 ((1 in1
) (1 l1
) (1 ido
) (1 5))
314 (setf (f2cl-lib:fref ch-%data%
316 ((1 in2
) (1 l1
) (1 5) (1 ido
))
319 (f2cl-lib:fref cc-%data%
321 ((1 in1
) (1 l1
) (1 ido
) (1 5))
325 (setf (f2cl-lib:fref ch-%data%
327 ((1 in2
) (1 l1
) (1 5) (1 ido
))
330 (f2cl-lib:fref cc-%data%
332 ((1 in1
) (1 l1
) (1 ido
) (1 5))
338 (f2cl-lib:fref cc-%data%
340 ((1 in1
) (1 l1
) (1 ido
) (1 5))
346 (f2cl-lib:fref cc-%data%
348 ((1 in1
) (1 l1
) (1 ido
) (1 5))
354 (f2cl-lib:fref cc-%data%
356 ((1 in1
) (1 l1
) (1 ido
) (1 5))
362 (f2cl-lib:fref cc-%data%
364 ((1 in1
) (1 l1
) (1 ido
) (1 5))
368 (setf cr5
(+ (* ti11 tr5
) (* ti12 tr4
)))
369 (setf ci5
(+ (* ti11 ti5
) (* ti12 ti4
)))
370 (setf cr4
(- (* ti12 tr5
) (* ti11 tr4
)))
371 (setf ci4
(- (* ti12 ti5
) (* ti11 ti4
)))
372 (setf (f2cl-lib:fref ch-%data%
374 ((1 in2
) (1 l1
) (1 5) (1 ido
))
377 (setf (f2cl-lib:fref ch-%data%
379 ((1 in2
) (1 l1
) (1 5) (1 ido
))
382 (setf (f2cl-lib:fref ch-%data%
384 ((1 in2
) (1 l1
) (1 5) (1 ido
))
387 (setf (f2cl-lib:fref ch-%data%
389 ((1 in2
) (1 l1
) (1 5) (1 ido
))
392 (setf (f2cl-lib:fref ch-%data%
394 ((1 in2
) (1 l1
) (1 5) (1 ido
))
397 (setf (f2cl-lib:fref ch-%data%
399 ((1 in2
) (1 l1
) (1 5) (1 ido
))
402 (setf (f2cl-lib:fref ch-%data%
404 ((1 in2
) (1 l1
) (1 5) (1 ido
))
407 (setf (f2cl-lib:fref ch-%data%
409 ((1 in2
) (1 l1
) (1 5) (1 ido
))
413 (if (= ido
1) (go end_label
))
414 (f2cl-lib:fdo
(i 2 (f2cl-lib:int-add i
1))
417 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
422 (f2cl-lib:fref cc-%data%
424 ((1 in1
) (1 l1
) (1 ido
) (1 5))
426 (f2cl-lib:fref cc-%data%
428 ((1 in1
) (1 l1
) (1 ido
) (1 5))
432 (f2cl-lib:fref cc-%data%
434 ((1 in1
) (1 l1
) (1 ido
) (1 5))
436 (f2cl-lib:fref cc-%data%
438 ((1 in1
) (1 l1
) (1 ido
) (1 5))
442 (f2cl-lib:fref cc-%data%
444 ((1 in1
) (1 l1
) (1 ido
) (1 5))
446 (f2cl-lib:fref cc-%data%
448 ((1 in1
) (1 l1
) (1 ido
) (1 5))
452 (f2cl-lib:fref cc-%data%
454 ((1 in1
) (1 l1
) (1 ido
) (1 5))
456 (f2cl-lib:fref cc-%data%
458 ((1 in1
) (1 l1
) (1 ido
) (1 5))
462 (f2cl-lib:fref cc-%data%
464 ((1 in1
) (1 l1
) (1 ido
) (1 5))
466 (f2cl-lib:fref cc-%data%
468 ((1 in1
) (1 l1
) (1 ido
) (1 5))
472 (f2cl-lib:fref cc-%data%
474 ((1 in1
) (1 l1
) (1 ido
) (1 5))
476 (f2cl-lib:fref cc-%data%
478 ((1 in1
) (1 l1
) (1 ido
) (1 5))
482 (f2cl-lib:fref cc-%data%
484 ((1 in1
) (1 l1
) (1 ido
) (1 5))
486 (f2cl-lib:fref cc-%data%
488 ((1 in1
) (1 l1
) (1 ido
) (1 5))
492 (f2cl-lib:fref cc-%data%
494 ((1 in1
) (1 l1
) (1 ido
) (1 5))
496 (f2cl-lib:fref cc-%data%
498 ((1 in1
) (1 l1
) (1 ido
) (1 5))
500 (setf (f2cl-lib:fref ch-%data%
502 ((1 in2
) (1 l1
) (1 5) (1 ido
))
505 (f2cl-lib:fref cc-%data%
507 ((1 in1
) (1 l1
) (1 ido
) (1 5))
511 (setf (f2cl-lib:fref ch-%data%
513 ((1 in2
) (1 l1
) (1 5) (1 ido
))
516 (f2cl-lib:fref cc-%data%
518 ((1 in1
) (1 l1
) (1 ido
) (1 5))
524 (f2cl-lib:fref cc-%data%
526 ((1 in1
) (1 l1
) (1 ido
) (1 5))
532 (f2cl-lib:fref cc-%data%
534 ((1 in1
) (1 l1
) (1 ido
) (1 5))
540 (f2cl-lib:fref cc-%data%
542 ((1 in1
) (1 l1
) (1 ido
) (1 5))
548 (f2cl-lib:fref cc-%data%
550 ((1 in1
) (1 l1
) (1 ido
) (1 5))
554 (setf cr5
(+ (* ti11 tr5
) (* ti12 tr4
)))
555 (setf ci5
(+ (* ti11 ti5
) (* ti12 ti4
)))
556 (setf cr4
(- (* ti12 tr5
) (* ti11 tr4
)))
557 (setf ci4
(- (* ti12 ti5
) (* ti11 ti4
)))
558 (setf dr3
(- cr3 ci4
))
559 (setf dr4
(+ cr3 ci4
))
560 (setf di3
(+ ci3 cr4
))
561 (setf di4
(- ci3 cr4
))
562 (setf dr5
(+ cr2 ci5
))
563 (setf dr2
(- cr2 ci5
))
564 (setf di5
(- ci2 cr5
))
565 (setf di2
(+ ci2 cr5
))
566 (setf (f2cl-lib:fref ch-%data%
568 ((1 in2
) (1 l1
) (1 5) (1 ido
))
572 (f2cl-lib:fref wa-%data%
574 ((1 ido
) (1 4) (1 2))
578 (f2cl-lib:fref wa-%data%
580 ((1 ido
) (1 4) (1 2))
583 (setf (f2cl-lib:fref ch-%data%
585 ((1 in2
) (1 l1
) (1 5) (1 ido
))
589 (f2cl-lib:fref wa-%data%
591 ((1 ido
) (1 4) (1 2))
595 (f2cl-lib:fref wa-%data%
597 ((1 ido
) (1 4) (1 2))
600 (setf (f2cl-lib:fref ch-%data%
602 ((1 in2
) (1 l1
) (1 5) (1 ido
))
606 (f2cl-lib:fref wa-%data%
608 ((1 ido
) (1 4) (1 2))
612 (f2cl-lib:fref wa-%data%
614 ((1 ido
) (1 4) (1 2))
617 (setf (f2cl-lib:fref ch-%data%
619 ((1 in2
) (1 l1
) (1 5) (1 ido
))
623 (f2cl-lib:fref wa-%data%
625 ((1 ido
) (1 4) (1 2))
629 (f2cl-lib:fref wa-%data%
631 ((1 ido
) (1 4) (1 2))
634 (setf (f2cl-lib:fref ch-%data%
636 ((1 in2
) (1 l1
) (1 5) (1 ido
))
640 (f2cl-lib:fref wa-%data%
642 ((1 ido
) (1 4) (1 2))
646 (f2cl-lib:fref wa-%data%
648 ((1 ido
) (1 4) (1 2))
651 (setf (f2cl-lib:fref ch-%data%
653 ((1 in2
) (1 l1
) (1 5) (1 ido
))
657 (f2cl-lib:fref wa-%data%
659 ((1 ido
) (1 4) (1 2))
663 (f2cl-lib:fref wa-%data%
665 ((1 ido
) (1 4) (1 2))
668 (setf (f2cl-lib:fref ch-%data%
670 ((1 in2
) (1 l1
) (1 5) (1 ido
))
674 (f2cl-lib:fref wa-%data%
676 ((1 ido
) (1 4) (1 2))
680 (f2cl-lib:fref wa-%data%
682 ((1 ido
) (1 4) (1 2))
685 (setf (f2cl-lib:fref ch-%data%
687 ((1 in2
) (1 l1
) (1 5) (1 ido
))
691 (f2cl-lib:fref wa-%data%
693 ((1 ido
) (1 4) (1 2))
697 (f2cl-lib:fref wa-%data%
699 ((1 ido
) (1 4) (1 2))
706 (return (values nil nil nil nil nil nil nil nil
))))))
708 (in-package #:cl-user
)
709 #+#.
(cl:if
(cl:find-package
'#:f2cl
) '(and) '(or))
710 (eval-when (:load-toplevel
:compile-toplevel
:execute
)
711 (setf (gethash 'fortran-to-lisp
::c1f5kb
712 fortran-to-lisp
::*f2cl-function-info
*)
713 (fortran-to-lisp::make-f2cl-finfo
714 :arg-types
'((fortran-to-lisp::integer4
) (fortran-to-lisp::integer4
)
715 (fortran-to-lisp::integer4
) (array double-float
(*))
716 (fortran-to-lisp::integer4
) (array double-float
(*))
717 (fortran-to-lisp::integer4
) (array double-float
(*)))
718 :return-values
'(nil nil nil nil nil nil nil nil
)