;;; Compiled by f2cl version:
;;; ("f2cl1.l,v 95098eb54f13 2013/04/01 00:45:16 toy $"
;;;  "f2cl2.l,v 95098eb54f13 2013/04/01 00:45:16 toy $"
;;;  "f2cl3.l,v 96616d88fb7e 2008/02/22 22:19:34 rtoy $"
;;;  "f2cl4.l,v 96616d88fb7e 2008/02/22 22:19:34 rtoy $"
;;;  "f2cl5.l,v 95098eb54f13 2013/04/01 00:45:16 toy $"
;;;  "f2cl6.l,v 1d5cbacbb977 2008/08/24 00:56:27 rtoy $"
;;;  "macros.l,v 1409c1352feb 2013/03/24 20:44:50 toy $")

;;; Using Lisp CMU Common Lisp snapshot-2020-04 (21D Unicode)
;;; 
;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
;;;           (:coerce-assigns :as-needed) (:array-type ':array)
;;;           (:array-slicing t) (:declare-common nil)
;;;           (:float-format single-float))

(in-package "FFTPACK5")


(defun mrftb1 (m im n in c ch wa fac)
  (declare (type (array double-float (*)) fac)
           (type (array double-float (*)) wa ch c)
           (type (f2cl-lib:integer4) in n im m))
  (f2cl-lib:with-multi-array-data
      ((c double-float c-%data% c-%offset%)
       (ch double-float ch-%data% ch-%offset%)
       (wa double-float wa-%data% wa-%offset%)
       (fac double-float fac-%data% fac-%offset%))
    (prog ((ix4 0) (ix3 0) (ix2 0) (idl1 0) (ido 0) (l2 0) (iw 0) (l1 0) (j 0)
           (i 0) (m2 0) (nl 0) (modn 0) (halfm 0.0d0) (half 0.0d0) (ip 0)
           (k1 0) (na 0) (nf 0))
      (declare (type (double-float) half halfm)
               (type (f2cl-lib:integer4) nf na k1 ip modn nl m2 i j l1 iw l2
                                         ido idl1 ix2 ix3 ix4))
      (setf nf
              (f2cl-lib:int
               (f2cl-lib:fref fac-%data% (2) ((1 15)) fac-%offset%)))
      (setf na 0)
      (f2cl-lib:fdo (k1 1 (f2cl-lib:int-add k1 1))
                    ((> k1 nf) nil)
        (tagbody
          (setf ip
                  (f2cl-lib:int
                   (f2cl-lib:fref fac-%data%
                                  ((f2cl-lib:int-add k1 2))
                                  ((1 15))
                                  fac-%offset%)))
          (setf na (f2cl-lib:int-sub 1 na))
          (if (<= ip 5) (go label10))
          (if (= k1 nf) (go label10))
          (setf na (f2cl-lib:int-sub 1 na))
         label10))
      (setf half 0.5d0)
      (setf halfm -0.5d0)
      (setf modn (mod n 2))
      (setf nl (f2cl-lib:int-sub n 2))
      (if (/= modn 0) (setf nl (f2cl-lib:int-sub n 1)))
      (if (= na 0) (go label120))
      (setf m2 (f2cl-lib:int-sub 1 im))
      (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
                    ((> i m) nil)
        (tagbody
          (setf m2 (f2cl-lib:int-add m2 im))
          (setf (f2cl-lib:fref ch-%data% (i 1) ((1 m) (1 *)) ch-%offset%)
                  (f2cl-lib:fref c-%data% (m2 1) ((1 in) (1 *)) c-%offset%))
          (setf (f2cl-lib:fref ch-%data% (i n) ((1 m) (1 *)) ch-%offset%)
                  (f2cl-lib:fref c-%data% (m2 n) ((1 in) (1 *)) c-%offset%))
         label117))
      (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 2))
                    ((> j nl) nil)
        (tagbody
          (setf m2 (f2cl-lib:int-sub 1 im))
          (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
                        ((> i m) nil)
            (tagbody
              (setf m2 (f2cl-lib:int-add m2 im))
              (setf (f2cl-lib:fref ch-%data% (i j) ((1 m) (1 *)) ch-%offset%)
                      (* half
                         (f2cl-lib:fref c-%data%
                                        (m2 j)
                                        ((1 in) (1 *))
                                        c-%offset%)))
              (setf (f2cl-lib:fref ch-%data%
                                   (i (f2cl-lib:int-add j 1))
                                   ((1 m) (1 *))
                                   ch-%offset%)
                      (* halfm
                         (f2cl-lib:fref c-%data%
                                        (m2 (f2cl-lib:int-add j 1))
                                        ((1 in) (1 *))
                                        c-%offset%)))
             label118))))
     label118
      (go label124)
     label120
      (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 2))
                    ((> j nl) nil)
        (tagbody
          (setf m2 (f2cl-lib:int-sub 1 im))
          (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
                        ((> i m) nil)
            (tagbody
              (setf m2 (f2cl-lib:int-add m2 im))
              (setf (f2cl-lib:fref c-%data% (m2 j) ((1 in) (1 *)) c-%offset%)
                      (* half
                         (f2cl-lib:fref c-%data%
                                        (m2 j)
                                        ((1 in) (1 *))
                                        c-%offset%)))
              (setf (f2cl-lib:fref c-%data%
                                   (m2 (f2cl-lib:int-add j 1))
                                   ((1 in) (1 *))
                                   c-%offset%)
                      (* halfm
                         (f2cl-lib:fref c-%data%
                                        (m2 (f2cl-lib:int-add j 1))
                                        ((1 in) (1 *))
                                        c-%offset%)))
             label122))))
     label122
     label124
      (setf l1 1)
      (setf iw 1)
      (f2cl-lib:fdo (k1 1 (f2cl-lib:int-add k1 1))
                    ((> k1 nf) nil)
        (tagbody
          (setf ip
                  (f2cl-lib:int
                   (f2cl-lib:fref fac-%data%
                                  ((f2cl-lib:int-add k1 2))
                                  ((1 15))
                                  fac-%offset%)))
          (setf l2 (f2cl-lib:int-mul ip l1))
          (setf ido (the f2cl-lib:integer4 (truncate n l2)))
          (setf idl1 (f2cl-lib:int-mul ido l1))
          (if (/= ip 4) (go label103))
          (setf ix2 (f2cl-lib:int-add iw ido))
          (setf ix3 (f2cl-lib:int-add ix2 ido))
          (if (/= na 0) (go label101))
          (mradb4 m ido l1 c im in ch 1 m
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix2)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix3)
                                 ((1 n))
                                 wa-%offset%))
          (go label102)
         label101
          (mradb4 m ido l1 ch 1 m c im in
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix2)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix3)
                                 ((1 n))
                                 wa-%offset%))
         label102
          (setf na (f2cl-lib:int-sub 1 na))
          (go label115)
         label103
          (if (/= ip 2) (go label106))
          (if (/= na 0) (go label104))
          (mradb2 m ido l1 c im in ch 1 m
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%))
          (go label105)
         label104
          (mradb2 m ido l1 ch 1 m c im in
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%))
         label105
          (setf na (f2cl-lib:int-sub 1 na))
          (go label115)
         label106
          (if (/= ip 3) (go label109))
          (setf ix2 (f2cl-lib:int-add iw ido))
          (if (/= na 0) (go label107))
          (mradb3 m ido l1 c im in ch 1 m
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix2)
                                 ((1 n))
                                 wa-%offset%))
          (go label108)
         label107
          (mradb3 m ido l1 ch 1 m c im in
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix2)
                                 ((1 n))
                                 wa-%offset%))
         label108
          (setf na (f2cl-lib:int-sub 1 na))
          (go label115)
         label109
          (if (/= ip 5) (go label112))
          (setf ix2 (f2cl-lib:int-add iw ido))
          (setf ix3 (f2cl-lib:int-add ix2 ido))
          (setf ix4 (f2cl-lib:int-add ix3 ido))
          (if (/= na 0) (go label110))
          (mradb5 m ido l1 c im in ch 1 m
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix2)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix3)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix4)
                                 ((1 n))
                                 wa-%offset%))
          (go label111)
         label110
          (mradb5 m ido l1 ch 1 m c im in
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix2)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix3)
                                 ((1 n))
                                 wa-%offset%)
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (ix4)
                                 ((1 n))
                                 wa-%offset%))
         label111
          (setf na (f2cl-lib:int-sub 1 na))
          (go label115)
         label112
          (if (/= na 0) (go label113))
          (mradbg m ido ip l1 idl1 c c c im in ch ch 1 m
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%))
          (go label114)
         label113
          (mradbg m ido ip l1 idl1 ch ch ch 1 m c c im in
           (f2cl-lib:array-slice wa-%data%
                                 double-float
                                 (iw)
                                 ((1 n))
                                 wa-%offset%))
         label114
          (if (= ido 1) (setf na (f2cl-lib:int-sub 1 na)))
         label115
          (setf l1 l2)
          (setf iw
                  (f2cl-lib:int-add iw
                                    (f2cl-lib:int-mul (f2cl-lib:int-sub ip 1)
                                                      ido)))
         label116))
      (go end_label)
     end_label
      (return (values nil nil nil nil nil nil nil nil)))))

(in-package #-gcl #:cl-user #+gcl "CL-USER")
#+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
(eval-when (:load-toplevel :compile-toplevel :execute)
  (setf (gethash 'fortran-to-lisp::mrftb1
                 fortran-to-lisp::*f2cl-function-info*)
          (fortran-to-lisp::make-f2cl-finfo
           :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
                        (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
                        (array double-float (*)) (array double-float (*))
                        (array double-float (*)) (array double-float (*)))
           :return-values '(nil nil nil nil nil nil nil nil)
           :calls '(fortran-to-lisp::mradbg fortran-to-lisp::mradb5
                    fortran-to-lisp::mradb3 fortran-to-lisp::mradb2
                    fortran-to-lisp::mradb4))))

