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 (defun r1f3kf (ido l1 cc in1 ch in2 wa1 wa2
)
21 (declare (type (array double-float
(*)) wa2 wa1 ch cc
)
22 (type (f2cl-lib:integer4
) in2 in1 l1 ido
))
23 (f2cl-lib:with-multi-array-data
24 ((cc double-float cc-%data% cc-%offset%
)
25 (ch double-float ch-%data% ch-%offset%
)
26 (wa1 double-float wa1-%data% wa1-%offset%
)
27 (wa2 double-float wa2-%data% wa2-%offset%
))
28 (prog ((ic 0) (i 0) (idp2 0) (k 0) (taui 0.0d0
) (taur 0.0d0
) (arg 0.0d0
))
29 (declare (type (double-float) arg taur taui
)
30 (type (f2cl-lib:integer4
) k idp2 i ic
))
31 (setf arg
(/ (* 2.0d0
4.0d0
(atan 1.0d0
)) 3.0d0
))
34 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
37 (setf (f2cl-lib:fref ch-%data%
39 ((1 in2
) (1 ido
) (1 3) (1 l1
))
42 (f2cl-lib:fref cc-%data%
44 ((1 in1
) (1 ido
) (1 l1
) (1 3))
47 (f2cl-lib:fref cc-%data%
49 ((1 in1
) (1 ido
) (1 l1
) (1 3))
51 (f2cl-lib:fref cc-%data%
53 ((1 in1
) (1 ido
) (1 l1
) (1 3))
55 (setf (f2cl-lib:fref ch-%data%
57 ((1 in2
) (1 ido
) (1 3) (1 l1
))
61 (f2cl-lib:fref cc-%data%
63 ((1 in1
) (1 ido
) (1 l1
) (1 3))
65 (f2cl-lib:fref cc-%data%
67 ((1 in1
) (1 ido
) (1 l1
) (1 3))
69 (setf (f2cl-lib:fref ch-%data%
71 ((1 in2
) (1 ido
) (1 3) (1 l1
))
74 (f2cl-lib:fref cc-%data%
76 ((1 in1
) (1 ido
) (1 l1
) (1 3))
80 (f2cl-lib:fref cc-%data%
82 ((1 in1
) (1 ido
) (1 l1
) (1 3))
84 (f2cl-lib:fref cc-%data%
86 ((1 in1
) (1 ido
) (1 l1
) (1 3))
89 (if (= ido
1) (go end_label
))
90 (setf idp2
(f2cl-lib:int-add ido
2))
91 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
94 (f2cl-lib:fdo
(i 3 (f2cl-lib:int-add i
2))
97 (setf ic
(f2cl-lib:int-sub idp2 i
))
98 (setf (f2cl-lib:fref ch-%data%
99 (1 (f2cl-lib:int-sub i
1) 1 k
)
100 ((1 in2
) (1 ido
) (1 3) (1 l1
))
103 (f2cl-lib:fref cc-%data%
104 (1 (f2cl-lib:int-sub i
1) k
1)
105 ((1 in1
) (1 ido
) (1 l1
) (1 3))
109 (f2cl-lib:fref wa1-%data%
110 ((f2cl-lib:int-sub i
2))
113 (f2cl-lib:fref cc-%data%
114 (1 (f2cl-lib:int-sub i
1) k
2)
115 ((1 in1
) (1 ido
) (1 l1
) (1 3))
118 (f2cl-lib:fref wa1-%data%
119 ((f2cl-lib:int-sub i
1))
122 (f2cl-lib:fref cc-%data%
124 ((1 in1
) (1 ido
) (1 l1
) (1 3))
128 (f2cl-lib:fref wa2-%data%
129 ((f2cl-lib:int-sub i
2))
132 (f2cl-lib:fref cc-%data%
133 (1 (f2cl-lib:int-sub i
1) k
3)
134 ((1 in1
) (1 ido
) (1 l1
) (1 3))
137 (f2cl-lib:fref wa2-%data%
138 ((f2cl-lib:int-sub i
1))
141 (f2cl-lib:fref cc-%data%
143 ((1 in1
) (1 ido
) (1 l1
) (1 3))
145 (setf (f2cl-lib:fref ch-%data%
147 ((1 in2
) (1 ido
) (1 3) (1 l1
))
150 (f2cl-lib:fref cc-%data%
152 ((1 in1
) (1 ido
) (1 l1
) (1 3))
157 (f2cl-lib:fref wa1-%data%
158 ((f2cl-lib:int-sub i
2))
161 (f2cl-lib:fref cc-%data%
163 ((1 in1
) (1 ido
) (1 l1
) (1 3))
166 (f2cl-lib:fref wa1-%data%
167 ((f2cl-lib:int-sub i
1))
170 (f2cl-lib:fref cc-%data%
171 (1 (f2cl-lib:int-sub i
1) k
2)
172 ((1 in1
) (1 ido
) (1 l1
) (1 3))
176 (f2cl-lib:fref wa2-%data%
177 ((f2cl-lib:int-sub i
2))
180 (f2cl-lib:fref cc-%data%
182 ((1 in1
) (1 ido
) (1 l1
) (1 3))
185 (f2cl-lib:fref wa2-%data%
186 ((f2cl-lib:int-sub i
1))
189 (f2cl-lib:fref cc-%data%
190 (1 (f2cl-lib:int-sub i
1) k
3)
191 ((1 in1
) (1 ido
) (1 l1
) (1 3))
193 (setf (f2cl-lib:fref ch-%data%
194 (1 (f2cl-lib:int-sub i
1) 3 k
)
195 ((1 in2
) (1 ido
) (1 3) (1 l1
))
198 (f2cl-lib:fref cc-%data%
199 (1 (f2cl-lib:int-sub i
1) k
1)
200 ((1 in1
) (1 ido
) (1 l1
) (1 3))
205 (f2cl-lib:fref wa1-%data%
206 ((f2cl-lib:int-sub i
2))
209 (f2cl-lib:fref cc-%data%
210 (1 (f2cl-lib:int-sub i
1) k
2)
211 ((1 in1
) (1 ido
) (1 l1
) (1 3))
214 (f2cl-lib:fref wa1-%data%
215 ((f2cl-lib:int-sub i
1))
218 (f2cl-lib:fref cc-%data%
220 ((1 in1
) (1 ido
) (1 l1
) (1 3))
224 (f2cl-lib:fref wa2-%data%
225 ((f2cl-lib:int-sub i
2))
228 (f2cl-lib:fref cc-%data%
229 (1 (f2cl-lib:int-sub i
1) k
3)
230 ((1 in1
) (1 ido
) (1 l1
) (1 3))
233 (f2cl-lib:fref wa2-%data%
234 ((f2cl-lib:int-sub i
1))
237 (f2cl-lib:fref cc-%data%
239 ((1 in1
) (1 ido
) (1 l1
) (1 3))
244 (f2cl-lib:fref wa1-%data%
245 ((f2cl-lib:int-sub i
2))
248 (f2cl-lib:fref cc-%data%
250 ((1 in1
) (1 ido
) (1 l1
) (1 3))
253 (f2cl-lib:fref wa1-%data%
254 ((f2cl-lib:int-sub i
1))
257 (f2cl-lib:fref cc-%data%
258 (1 (f2cl-lib:int-sub i
1) k
2)
259 ((1 in1
) (1 ido
) (1 l1
) (1 3))
263 (f2cl-lib:fref wa2-%data%
264 ((f2cl-lib:int-sub i
2))
267 (f2cl-lib:fref cc-%data%
269 ((1 in1
) (1 ido
) (1 l1
) (1 3))
272 (f2cl-lib:fref wa2-%data%
273 ((f2cl-lib:int-sub i
1))
276 (f2cl-lib:fref cc-%data%
277 (1 (f2cl-lib:int-sub i
1) k
3)
278 ((1 in1
) (1 ido
) (1 l1
) (1 3))
280 (setf (f2cl-lib:fref ch-%data%
281 (1 (f2cl-lib:int-sub ic
1) 2 k
)
282 ((1 in2
) (1 ido
) (1 3) (1 l1
))
286 (f2cl-lib:fref cc-%data%
287 (1 (f2cl-lib:int-sub i
1) k
1)
288 ((1 in1
) (1 ido
) (1 l1
) (1 3))
293 (f2cl-lib:fref wa1-%data%
294 ((f2cl-lib:int-sub i
2))
297 (f2cl-lib:fref cc-%data%
298 (1 (f2cl-lib:int-sub i
1) k
2)
299 ((1 in1
) (1 ido
) (1 l1
) (1 3))
302 (f2cl-lib:fref wa1-%data%
303 ((f2cl-lib:int-sub i
1))
306 (f2cl-lib:fref cc-%data%
308 ((1 in1
) (1 ido
) (1 l1
) (1 3))
312 (f2cl-lib:fref wa2-%data%
313 ((f2cl-lib:int-sub i
2))
316 (f2cl-lib:fref cc-%data%
317 (1 (f2cl-lib:int-sub i
1) k
3)
318 ((1 in1
) (1 ido
) (1 l1
) (1 3))
321 (f2cl-lib:fref wa2-%data%
322 ((f2cl-lib:int-sub i
1))
325 (f2cl-lib:fref cc-%data%
327 ((1 in1
) (1 ido
) (1 l1
) (1 3))
332 (f2cl-lib:fref wa1-%data%
333 ((f2cl-lib:int-sub i
2))
336 (f2cl-lib:fref cc-%data%
338 ((1 in1
) (1 ido
) (1 l1
) (1 3))
341 (f2cl-lib:fref wa1-%data%
342 ((f2cl-lib:int-sub i
1))
345 (f2cl-lib:fref cc-%data%
346 (1 (f2cl-lib:int-sub i
1) k
2)
347 ((1 in1
) (1 ido
) (1 l1
) (1 3))
351 (f2cl-lib:fref wa2-%data%
352 ((f2cl-lib:int-sub i
2))
355 (f2cl-lib:fref cc-%data%
357 ((1 in1
) (1 ido
) (1 l1
) (1 3))
360 (f2cl-lib:fref wa2-%data%
361 ((f2cl-lib:int-sub i
1))
364 (f2cl-lib:fref cc-%data%
365 (1 (f2cl-lib:int-sub i
1) k
3)
366 ((1 in1
) (1 ido
) (1 l1
) (1 3))
368 (setf (f2cl-lib:fref ch-%data%
370 ((1 in2
) (1 ido
) (1 3) (1 l1
))
373 (f2cl-lib:fref cc-%data%
375 ((1 in1
) (1 ido
) (1 l1
) (1 3))
381 (f2cl-lib:fref wa1-%data%
382 ((f2cl-lib:int-sub i
2))
385 (f2cl-lib:fref cc-%data%
387 ((1 in1
) (1 ido
) (1 l1
) (1 3))
390 (f2cl-lib:fref wa1-%data%
391 ((f2cl-lib:int-sub i
1))
394 (f2cl-lib:fref cc-%data%
395 (1 (f2cl-lib:int-sub i
1) k
2)
396 ((1 in1
) (1 ido
) (1 l1
) (1 3))
400 (f2cl-lib:fref wa2-%data%
401 ((f2cl-lib:int-sub i
2))
404 (f2cl-lib:fref cc-%data%
406 ((1 in1
) (1 ido
) (1 l1
) (1 3))
409 (f2cl-lib:fref wa2-%data%
410 ((f2cl-lib:int-sub i
1))
413 (f2cl-lib:fref cc-%data%
414 (1 (f2cl-lib:int-sub i
1) k
3)
415 ((1 in1
) (1 ido
) (1 l1
) (1 3))
421 (f2cl-lib:fref wa2-%data%
422 ((f2cl-lib:int-sub i
2))
425 (f2cl-lib:fref cc-%data%
426 (1 (f2cl-lib:int-sub i
1) k
3)
427 ((1 in1
) (1 ido
) (1 l1
) (1 3))
430 (f2cl-lib:fref wa2-%data%
431 ((f2cl-lib:int-sub i
1))
434 (f2cl-lib:fref cc-%data%
436 ((1 in1
) (1 ido
) (1 l1
) (1 3))
440 (f2cl-lib:fref wa1-%data%
441 ((f2cl-lib:int-sub i
2))
444 (f2cl-lib:fref cc-%data%
445 (1 (f2cl-lib:int-sub i
1) k
2)
446 ((1 in1
) (1 ido
) (1 l1
) (1 3))
449 (f2cl-lib:fref wa1-%data%
450 ((f2cl-lib:int-sub i
1))
453 (f2cl-lib:fref cc-%data%
455 ((1 in1
) (1 ido
) (1 l1
) (1 3))
457 (setf (f2cl-lib:fref ch-%data%
459 ((1 in2
) (1 ido
) (1 3) (1 l1
))
466 (f2cl-lib:fref wa2-%data%
467 ((f2cl-lib:int-sub i
2))
470 (f2cl-lib:fref cc-%data%
471 (1 (f2cl-lib:int-sub i
1) k
3)
472 ((1 in1
) (1 ido
) (1 l1
) (1 3))
475 (f2cl-lib:fref wa2-%data%
476 ((f2cl-lib:int-sub i
1))
479 (f2cl-lib:fref cc-%data%
481 ((1 in1
) (1 ido
) (1 l1
) (1 3))
485 (f2cl-lib:fref wa1-%data%
486 ((f2cl-lib:int-sub i
2))
489 (f2cl-lib:fref cc-%data%
490 (1 (f2cl-lib:int-sub i
1) k
2)
491 ((1 in1
) (1 ido
) (1 l1
) (1 3))
494 (f2cl-lib:fref wa1-%data%
495 ((f2cl-lib:int-sub i
1))
498 (f2cl-lib:fref cc-%data%
500 ((1 in1
) (1 ido
) (1 l1
) (1 3))
503 (f2cl-lib:fref cc-%data%
505 ((1 in1
) (1 ido
) (1 l1
) (1 3))
511 (f2cl-lib:fref wa1-%data%
512 ((f2cl-lib:int-sub i
2))
515 (f2cl-lib:fref cc-%data%
517 ((1 in1
) (1 ido
) (1 l1
) (1 3))
520 (f2cl-lib:fref wa1-%data%
521 ((f2cl-lib:int-sub i
1))
524 (f2cl-lib:fref cc-%data%
525 (1 (f2cl-lib:int-sub i
1) k
2)
526 ((1 in1
) (1 ido
) (1 l1
) (1 3))
530 (f2cl-lib:fref wa2-%data%
531 ((f2cl-lib:int-sub i
2))
534 (f2cl-lib:fref cc-%data%
536 ((1 in1
) (1 ido
) (1 l1
) (1 3))
539 (f2cl-lib:fref wa2-%data%
540 ((f2cl-lib:int-sub i
1))
543 (f2cl-lib:fref cc-%data%
544 (1 (f2cl-lib:int-sub i
1) k
3)
545 ((1 in1
) (1 ido
) (1 l1
) (1 3))
551 (return (values nil nil nil nil nil nil nil nil
)))))
553 (in-package #:cl-user
)
554 #+#.
(cl:if
(cl:find-package
'#:f2cl
) '(and) '(or))
555 (eval-when (:load-toplevel
:compile-toplevel
:execute
)
556 (setf (gethash 'fortran-to-lisp
::r1f3kf
557 fortran-to-lisp
::*f2cl-function-info
*)
558 (fortran-to-lisp::make-f2cl-finfo
559 :arg-types
'((fortran-to-lisp::integer4
) (fortran-to-lisp::integer4
)
560 (array double-float
(*)) (fortran-to-lisp::integer4
)
561 (array double-float
(*)) (fortran-to-lisp::integer4
)
562 (array double-float
(*)) (array double-float
(*)))
563 :return-values
'(nil nil nil nil nil nil nil nil
)