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 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
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
))
40 (f2cl-lib:int-add
(f2cl-lib:int-mul
(f2cl-lib:int-sub m
1) im1
)
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
)))
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))
56 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
59 (setf m2
(f2cl-lib:int-add m2 im2
))
60 (setf (f2cl-lib:fref ch2-%data%
62 ((1 in2
) (1 idl1
) (1 ip
))
64 (f2cl-lib:fref c2-%data%
66 ((1 in1
) (1 idl1
) (1 ip
))
70 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
73 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
77 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
80 (setf m2
(f2cl-lib:int-add m2 im2
))
81 (setf (f2cl-lib:fref ch-%data%
83 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
85 (f2cl-lib:fref c1-%data%
87 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
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))
97 (setf is
(f2cl-lib:int-add is ido
))
99 (f2cl-lib:fdo
(i 3 (f2cl-lib:int-add i
2))
102 (setf idij
(f2cl-lib:int-add idij
2))
103 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
107 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
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
))
117 (f2cl-lib:fref wa-%data%
118 ((f2cl-lib:int-sub idij
1))
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
))
126 (f2cl-lib:fref wa-%data%
130 (f2cl-lib:fref c1-%data%
132 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
134 (setf (f2cl-lib:fref ch-%data%
136 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
140 (f2cl-lib:fref wa-%data%
141 ((f2cl-lib:int-sub idij
1))
144 (f2cl-lib:fref c1-%data%
146 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
149 (f2cl-lib:fref wa-%data%
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
))
163 (setf is
(f2cl-lib:int-sub ido
))
164 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
167 (setf is
(f2cl-lib:int-add is ido
))
168 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
172 (f2cl-lib:fdo
(i 3 (f2cl-lib:int-add i
2))
175 (setf idij
(f2cl-lib:int-add idij
2))
177 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
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
))
187 (f2cl-lib:fref wa-%data%
188 ((f2cl-lib:int-sub idij
1))
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
))
196 (f2cl-lib:fref wa-%data%
200 (f2cl-lib:fref c1-%data%
202 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
204 (setf (f2cl-lib:fref ch-%data%
206 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
210 (f2cl-lib:fref wa-%data%
211 ((f2cl-lib:int-sub idij
1))
214 (f2cl-lib:fref c1-%data%
216 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
219 (f2cl-lib:fref wa-%data%
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
))
232 (if (< nbd l1
) (go label115
))
233 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
236 (setf jc
(f2cl-lib:int-sub ipp2 j
))
237 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
240 (f2cl-lib:fdo
(i 3 (f2cl-lib:int-add i
2))
244 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
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
))
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
))
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
))
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
))
266 (f2cl-lib:fref ch-%data%
268 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
270 (f2cl-lib:fref ch-%data%
272 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
274 (setf (f2cl-lib:fref c1-%data%
276 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
279 (f2cl-lib:fref ch-%data%
281 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
283 (f2cl-lib:fref ch-%data%
285 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
287 (setf (f2cl-lib:fref c1-%data%
289 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
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
))
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
))
306 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
309 (setf jc
(f2cl-lib:int-sub ipp2 j
))
310 (f2cl-lib:fdo
(i 3 (f2cl-lib:int-add i
2))
313 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
317 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
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
))
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
))
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
))
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
))
339 (f2cl-lib:fref ch-%data%
341 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
343 (f2cl-lib:fref ch-%data%
345 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
347 (setf (f2cl-lib:fref c1-%data%
349 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
352 (f2cl-lib:fref ch-%data%
354 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
356 (f2cl-lib:fref ch-%data%
358 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
360 (setf (f2cl-lib:fref c1-%data%
362 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
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
))
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
))
379 (f2cl-lib:fdo
(ik 1 (f2cl-lib:int-add ik
1))
383 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
386 (setf m2
(f2cl-lib:int-add m2 im2
))
387 (setf (f2cl-lib:fref c2-%data%
389 ((1 in1
) (1 idl1
) (1 ip
))
391 (f2cl-lib:fref ch2-%data%
393 ((1 in2
) (1 idl1
) (1 ip
))
398 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
401 (setf jc
(f2cl-lib:int-sub ipp2 j
))
402 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
406 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
409 (setf m2
(f2cl-lib:int-add m2 im2
))
410 (setf (f2cl-lib:fref c1-%data%
412 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
415 (f2cl-lib:fref ch-%data%
417 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
419 (f2cl-lib:fref ch-%data%
421 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
423 (setf (f2cl-lib:fref c1-%data%
425 ((1 in1
) (1 ido
) (1 l1
) (1 ip
))
428 (f2cl-lib:fref ch-%data%
430 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
432 (f2cl-lib:fref ch-%data%
434 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
441 (f2cl-lib:fdo
(l 2 (f2cl-lib:int-add l
1))
444 (setf lc
(f2cl-lib:int-sub ipp2 l
))
445 (setf ar1h
(- (* dcp ar1
) (* dsp ai1
)))
446 (setf ai1
(+ (* dcp ai1
) (* dsp ar1
)))
448 (f2cl-lib:fdo
(ik 1 (f2cl-lib:int-add ik
1))
452 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
455 (setf m2
(f2cl-lib:int-add m2 im2
))
456 (setf (f2cl-lib:fref ch2-%data%
458 ((1 in2
) (1 idl1
) (1 ip
))
461 (f2cl-lib:fref c2-%data%
463 ((1 in1
) (1 idl1
) (1 ip
))
466 (f2cl-lib:fref c2-%data%
468 ((1 in1
) (1 idl1
) (1 ip
))
470 (setf (f2cl-lib:fref ch2-%data%
472 ((1 in2
) (1 idl1
) (1 ip
))
475 (f2cl-lib:fref c2-%data%
477 ((1 in1
) (1 idl1
) (1 ip
))
485 (f2cl-lib:fdo
(j 3 (f2cl-lib:int-add j
1))
488 (setf jc
(f2cl-lib:int-sub ipp2 j
))
489 (setf ar2h
(- (* dc2 ar2
) (* ds2 ai2
)))
490 (setf ai2
(+ (* dc2 ai2
) (* ds2 ar2
)))
492 (f2cl-lib:fdo
(ik 1 (f2cl-lib:int-add ik
1))
496 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
499 (setf m2
(f2cl-lib:int-add m2 im2
))
500 (setf (f2cl-lib:fref ch2-%data%
502 ((1 in2
) (1 idl1
) (1 ip
))
505 (f2cl-lib:fref ch2-%data%
507 ((1 in2
) (1 idl1
) (1 ip
))
510 (f2cl-lib:fref c2-%data%
512 ((1 in1
) (1 idl1
) (1 ip
))
514 (setf (f2cl-lib:fref ch2-%data%
516 ((1 in2
) (1 idl1
) (1 ip
))
519 (f2cl-lib:fref ch2-%data%
521 ((1 in2
) (1 idl1
) (1 ip
))
524 (f2cl-lib:fref c2-%data%
526 ((1 in1
) (1 idl1
) (1 ip
))
532 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
535 (f2cl-lib:fdo
(ik 1 (f2cl-lib:int-add ik
1))
539 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
542 (setf m2
(f2cl-lib:int-add m2 im2
))
543 (setf (f2cl-lib:fref ch2-%data%
545 ((1 in2
) (1 idl1
) (1 ip
))
548 (f2cl-lib:fref ch2-%data%
550 ((1 in2
) (1 idl1
) (1 ip
))
552 (f2cl-lib:fref c2-%data%
554 ((1 in1
) (1 idl1
) (1 ip
))
559 (if (< ido l1
) (go label132
))
560 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
563 (f2cl-lib:fdo
(i 1 (f2cl-lib:int-add i
1))
567 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
570 (setf m2
(f2cl-lib:int-add m2 im2
))
571 (setf (f2cl-lib:fref cc-%data%
573 ((1 in1
) (1 ido
) (1 ip
) (1 l1
))
575 (f2cl-lib:fref ch-%data%
577 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
584 (f2cl-lib:fdo
(i 1 (f2cl-lib:int-add i
1))
587 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
591 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
594 (setf m2
(f2cl-lib:int-add m2 im2
))
595 (setf (f2cl-lib:fref cc-%data%
597 ((1 in1
) (1 ido
) (1 ip
) (1 l1
))
599 (f2cl-lib:fref ch-%data%
601 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
607 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
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))
616 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
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
))
624 (f2cl-lib:fref ch-%data%
626 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
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
))
632 (f2cl-lib:fref ch-%data%
634 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
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))
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))
649 (f2cl-lib:fdo
(i 3 (f2cl-lib:int-add i
2))
652 (setf ic
(f2cl-lib:int-sub idp2 i
))
654 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
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
))
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
))
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
))
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
))
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
))
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
))
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
))
691 (f2cl-lib:fref ch-%data%
693 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
695 (f2cl-lib:fref ch-%data%
697 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
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
))
704 (f2cl-lib:fref ch-%data%
706 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
708 (f2cl-lib:fref ch-%data%
710 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
718 (f2cl-lib:fdo
(j 2 (f2cl-lib:int-add j
1))
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))
726 (setf ic
(f2cl-lib:int-sub idp2 i
))
727 (f2cl-lib:fdo
(k 1 (f2cl-lib:int-add k
1))
731 (f2cl-lib:fdo
(m1 1 (f2cl-lib:int-add m1 im1
))
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
))
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
))
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
))
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
))
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
))
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
))
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
))
768 (f2cl-lib:fref ch-%data%
770 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
772 (f2cl-lib:fref ch-%data%
774 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
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
))
781 (f2cl-lib:fref ch-%data%
783 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
785 (f2cl-lib:fref ch-%data%
787 ((1 in2
) (1 ido
) (1 l1
) (1 ip
))
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