share/fftpack5/lisp/c1f3kb.lisp

   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 $")
   9
  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))
  16
  17 (in-package "FFTPACK5")
  18
  19
  20 (let ((taur -0.5d0) (taui 0.866025403784439d0))
  21   (declare (type (double-float) taur taui))
  22   (defun c1f3kb (ido l1 na cc in1 ch in2 wa)
  23     (declare (type (array double-float (*)) wa ch cc)
  24              (type (f2cl-lib:integer4) in2 in1 na l1 ido))
  25     (f2cl-lib:with-multi-array-data
  26         ((cc double-float cc-%data% cc-%offset%)
  27          (ch double-float ch-%data% ch-%offset%)
  28          (wa double-float wa-%data% wa-%offset%))
  29       (prog ((di3 0.0d0) (di2 0.0d0) (dr3 0.0d0) (dr2 0.0d0) (i 0) (ci3 0.0d0)
  30              (cr3 0.0d0) (ci2 0.0d0) (ti2 0.0d0) (cr2 0.0d0) (tr2 0.0d0) (k 0))
  31         (declare (type (f2cl-lib:integer4) k i)
  32                  (type (double-float) tr2 cr2 ti2 ci2 cr3 ci3 dr2 dr3 di2 di3))
  33         (if (or (> ido 1) (= na 1)) (go label102))
  34         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  35                       ((> k l1) nil)
  36           (tagbody
  37             (setf tr2
  38                     (+
  39                      (f2cl-lib:fref cc-%data%
  40                                     (1 k 1 2)
  41                                     ((1 in1) (1 l1) (1 ido) (1 3))
  42                                     cc-%offset%)
  43                      (f2cl-lib:fref cc-%data%
  44                                     (1 k 1 3)
  45                                     ((1 in1) (1 l1) (1 ido) (1 3))
  46                                     cc-%offset%)))
  47             (setf cr2
  48                     (+
  49                      (f2cl-lib:fref cc-%data%
  50                                     (1 k 1 1)
  51                                     ((1 in1) (1 l1) (1 ido) (1 3))
  52                                     cc-%offset%)
  53                      (* taur tr2)))
  54             (setf (f2cl-lib:fref cc-%data%
  55                                  (1 k 1 1)
  56                                  ((1 in1) (1 l1) (1 ido) (1 3))
  57                                  cc-%offset%)
  58                     (+
  59                      (f2cl-lib:fref cc-%data%
  60                                     (1 k 1 1)
  61                                     ((1 in1) (1 l1) (1 ido) (1 3))
  62                                     cc-%offset%)
  63                      tr2))
  64             (setf ti2
  65                     (+
  66                      (f2cl-lib:fref cc-%data%
  67                                     (2 k 1 2)
  68                                     ((1 in1) (1 l1) (1 ido) (1 3))
  69                                     cc-%offset%)
  70                      (f2cl-lib:fref cc-%data%
  71                                     (2 k 1 3)
  72                                     ((1 in1) (1 l1) (1 ido) (1 3))
  73                                     cc-%offset%)))
  74             (setf ci2
  75                     (+
  76                      (f2cl-lib:fref cc-%data%
  77                                     (2 k 1 1)
  78                                     ((1 in1) (1 l1) (1 ido) (1 3))
  79                                     cc-%offset%)
  80                      (* taur ti2)))
  81             (setf (f2cl-lib:fref cc-%data%
  82                                  (2 k 1 1)
  83                                  ((1 in1) (1 l1) (1 ido) (1 3))
  84                                  cc-%offset%)
  85                     (+
  86                      (f2cl-lib:fref cc-%data%
  87                                     (2 k 1 1)
  88                                     ((1 in1) (1 l1) (1 ido) (1 3))
  89                                     cc-%offset%)
  90                      ti2))
  91             (setf cr3
  92                     (* taui
  93                        (-
  94                         (f2cl-lib:fref cc-%data%
  95                                        (1 k 1 2)
  96                                        ((1 in1) (1 l1) (1 ido) (1 3))
  97                                        cc-%offset%)
  98                         (f2cl-lib:fref cc-%data%
  99                                        (1 k 1 3)
 100                                        ((1 in1) (1 l1) (1 ido) (1 3))
 101                                        cc-%offset%))))
 102             (setf ci3
 103                     (* taui
 104                        (-
 105                         (f2cl-lib:fref cc-%data%
 106                                        (2 k 1 2)
 107                                        ((1 in1) (1 l1) (1 ido) (1 3))
 108                                        cc-%offset%)
 109                         (f2cl-lib:fref cc-%data%
 110                                        (2 k 1 3)
 111                                        ((1 in1) (1 l1) (1 ido) (1 3))
 112                                        cc-%offset%))))
 113             (setf (f2cl-lib:fref cc-%data%
 114                                  (1 k 1 2)
 115                                  ((1 in1) (1 l1) (1 ido) (1 3))
 116                                  cc-%offset%)
 117                     (- cr2 ci3))
 118             (setf (f2cl-lib:fref cc-%data%
 119                                  (1 k 1 3)
 120                                  ((1 in1) (1 l1) (1 ido) (1 3))
 121                                  cc-%offset%)
 122                     (+ cr2 ci3))
 123             (setf (f2cl-lib:fref cc-%data%
 124                                  (2 k 1 2)
 125                                  ((1 in1) (1 l1) (1 ido) (1 3))
 126                                  cc-%offset%)
 127                     (+ ci2 cr3))
 128             (setf (f2cl-lib:fref cc-%data%
 129                                  (2 k 1 3)
 130                                  ((1 in1) (1 l1) (1 ido) (1 3))
 131                                  cc-%offset%)
 132                     (- ci2 cr3))
 133            label101))
 134         (go end_label)
 135        label102
 136         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 137                       ((> k l1) nil)
 138           (tagbody
 139             (setf tr2
 140                     (+
 141                      (f2cl-lib:fref cc-%data%
 142                                     (1 k 1 2)
 143                                     ((1 in1) (1 l1) (1 ido) (1 3))
 144                                     cc-%offset%)
 145                      (f2cl-lib:fref cc-%data%
 146                                     (1 k 1 3)
 147                                     ((1 in1) (1 l1) (1 ido) (1 3))
 148                                     cc-%offset%)))
 149             (setf cr2
 150                     (+
 151                      (f2cl-lib:fref cc-%data%
 152                                     (1 k 1 1)
 153                                     ((1 in1) (1 l1) (1 ido) (1 3))
 154                                     cc-%offset%)
 155                      (* taur tr2)))
 156             (setf (f2cl-lib:fref ch-%data%
 157                                  (1 k 1 1)
 158                                  ((1 in2) (1 l1) (1 3) (1 ido))
 159                                  ch-%offset%)
 160                     (+
 161                      (f2cl-lib:fref cc-%data%
 162                                     (1 k 1 1)
 163                                     ((1 in1) (1 l1) (1 ido) (1 3))
 164                                     cc-%offset%)
 165                      tr2))
 166             (setf ti2
 167                     (+
 168                      (f2cl-lib:fref cc-%data%
 169                                     (2 k 1 2)
 170                                     ((1 in1) (1 l1) (1 ido) (1 3))
 171                                     cc-%offset%)
 172                      (f2cl-lib:fref cc-%data%
 173                                     (2 k 1 3)
 174                                     ((1 in1) (1 l1) (1 ido) (1 3))
 175                                     cc-%offset%)))
 176             (setf ci2
 177                     (+
 178                      (f2cl-lib:fref cc-%data%
 179                                     (2 k 1 1)
 180                                     ((1 in1) (1 l1) (1 ido) (1 3))
 181                                     cc-%offset%)
 182                      (* taur ti2)))
 183             (setf (f2cl-lib:fref ch-%data%
 184                                  (2 k 1 1)
 185                                  ((1 in2) (1 l1) (1 3) (1 ido))
 186                                  ch-%offset%)
 187                     (+
 188                      (f2cl-lib:fref cc-%data%
 189                                     (2 k 1 1)
 190                                     ((1 in1) (1 l1) (1 ido) (1 3))
 191                                     cc-%offset%)
 192                      ti2))
 193             (setf cr3
 194                     (* taui
 195                        (-
 196                         (f2cl-lib:fref cc-%data%
 197                                        (1 k 1 2)
 198                                        ((1 in1) (1 l1) (1 ido) (1 3))
 199                                        cc-%offset%)
 200                         (f2cl-lib:fref cc-%data%
 201                                        (1 k 1 3)
 202                                        ((1 in1) (1 l1) (1 ido) (1 3))
 203                                        cc-%offset%))))
 204             (setf ci3
 205                     (* taui
 206                        (-
 207                         (f2cl-lib:fref cc-%data%
 208                                        (2 k 1 2)
 209                                        ((1 in1) (1 l1) (1 ido) (1 3))
 210                                        cc-%offset%)
 211                         (f2cl-lib:fref cc-%data%
 212                                        (2 k 1 3)
 213                                        ((1 in1) (1 l1) (1 ido) (1 3))
 214                                        cc-%offset%))))
 215             (setf (f2cl-lib:fref ch-%data%
 216                                  (1 k 2 1)
 217                                  ((1 in2) (1 l1) (1 3) (1 ido))
 218                                  ch-%offset%)
 219                     (- cr2 ci3))
 220             (setf (f2cl-lib:fref ch-%data%
 221                                  (1 k 3 1)
 222                                  ((1 in2) (1 l1) (1 3) (1 ido))
 223                                  ch-%offset%)
 224                     (+ cr2 ci3))
 225             (setf (f2cl-lib:fref ch-%data%
 226                                  (2 k 2 1)
 227                                  ((1 in2) (1 l1) (1 3) (1 ido))
 228                                  ch-%offset%)
 229                     (+ ci2 cr3))
 230             (setf (f2cl-lib:fref ch-%data%
 231                                  (2 k 3 1)
 232                                  ((1 in2) (1 l1) (1 3) (1 ido))
 233                                  ch-%offset%)
 234                     (- ci2 cr3))
 235            label103))
 236         (if (= ido 1) (go end_label))
 237         (f2cl-lib:fdo (i 2 (f2cl-lib:int-add i 1))
 238                       ((> i ido) nil)
 239           (tagbody
 240             (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 241                           ((> k l1) nil)
 242               (tagbody
 243                 (setf tr2
 244                         (+
 245                          (f2cl-lib:fref cc-%data%
 246                                         (1 k i 2)
 247                                         ((1 in1) (1 l1) (1 ido) (1 3))
 248                                         cc-%offset%)
 249                          (f2cl-lib:fref cc-%data%
 250                                         (1 k i 3)
 251                                         ((1 in1) (1 l1) (1 ido) (1 3))
 252                                         cc-%offset%)))
 253                 (setf cr2
 254                         (+
 255                          (f2cl-lib:fref cc-%data%
 256                                         (1 k i 1)
 257                                         ((1 in1) (1 l1) (1 ido) (1 3))
 258                                         cc-%offset%)
 259                          (* taur tr2)))
 260                 (setf (f2cl-lib:fref ch-%data%
 261                                      (1 k 1 i)
 262                                      ((1 in2) (1 l1) (1 3) (1 ido))
 263                                      ch-%offset%)
 264                         (+
 265                          (f2cl-lib:fref cc-%data%
 266                                         (1 k i 1)
 267                                         ((1 in1) (1 l1) (1 ido) (1 3))
 268                                         cc-%offset%)
 269                          tr2))
 270                 (setf ti2
 271                         (+
 272                          (f2cl-lib:fref cc-%data%
 273                                         (2 k i 2)
 274                                         ((1 in1) (1 l1) (1 ido) (1 3))
 275                                         cc-%offset%)
 276                          (f2cl-lib:fref cc-%data%
 277                                         (2 k i 3)
 278                                         ((1 in1) (1 l1) (1 ido) (1 3))
 279                                         cc-%offset%)))
 280                 (setf ci2
 281                         (+
 282                          (f2cl-lib:fref cc-%data%
 283                                         (2 k i 1)
 284                                         ((1 in1) (1 l1) (1 ido) (1 3))
 285                                         cc-%offset%)
 286                          (* taur ti2)))
 287                 (setf (f2cl-lib:fref ch-%data%
 288                                      (2 k 1 i)
 289                                      ((1 in2) (1 l1) (1 3) (1 ido))
 290                                      ch-%offset%)
 291                         (+
 292                          (f2cl-lib:fref cc-%data%
 293                                         (2 k i 1)
 294                                         ((1 in1) (1 l1) (1 ido) (1 3))
 295                                         cc-%offset%)
 296                          ti2))
 297                 (setf cr3
 298                         (* taui
 299                            (-
 300                             (f2cl-lib:fref cc-%data%
 301                                            (1 k i 2)
 302                                            ((1 in1) (1 l1) (1 ido) (1 3))
 303                                            cc-%offset%)
 304                             (f2cl-lib:fref cc-%data%
 305                                            (1 k i 3)
 306                                            ((1 in1) (1 l1) (1 ido) (1 3))
 307                                            cc-%offset%))))
 308                 (setf ci3
 309                         (* taui
 310                            (-
 311                             (f2cl-lib:fref cc-%data%
 312                                            (2 k i 2)
 313                                            ((1 in1) (1 l1) (1 ido) (1 3))
 314                                            cc-%offset%)
 315                             (f2cl-lib:fref cc-%data%
 316                                            (2 k i 3)
 317                                            ((1 in1) (1 l1) (1 ido) (1 3))
 318                                            cc-%offset%))))
 319                 (setf dr2 (- cr2 ci3))
 320                 (setf dr3 (+ cr2 ci3))
 321                 (setf di2 (+ ci2 cr3))
 322                 (setf di3 (- ci2 cr3))
 323                 (setf (f2cl-lib:fref ch-%data%
 324                                      (2 k 2 i)
 325                                      ((1 in2) (1 l1) (1 3) (1 ido))
 326                                      ch-%offset%)
 327                         (+
 328                          (*
 329                           (f2cl-lib:fref wa-%data%
 330                                          (i 1 1)
 331                                          ((1 ido) (1 2) (1 2))
 332                                          wa-%offset%)
 333                           di2)
 334                          (*
 335                           (f2cl-lib:fref wa-%data%
 336                                          (i 1 2)
 337                                          ((1 ido) (1 2) (1 2))
 338                                          wa-%offset%)
 339                           dr2)))
 340                 (setf (f2cl-lib:fref ch-%data%
 341                                      (1 k 2 i)
 342                                      ((1 in2) (1 l1) (1 3) (1 ido))
 343                                      ch-%offset%)
 344                         (-
 345                          (*
 346                           (f2cl-lib:fref wa-%data%
 347                                          (i 1 1)
 348                                          ((1 ido) (1 2) (1 2))
 349                                          wa-%offset%)
 350                           dr2)
 351                          (*
 352                           (f2cl-lib:fref wa-%data%
 353                                          (i 1 2)
 354                                          ((1 ido) (1 2) (1 2))
 355                                          wa-%offset%)
 356                           di2)))
 357                 (setf (f2cl-lib:fref ch-%data%
 358                                      (2 k 3 i)
 359                                      ((1 in2) (1 l1) (1 3) (1 ido))
 360                                      ch-%offset%)
 361                         (+
 362                          (*
 363                           (f2cl-lib:fref wa-%data%
 364                                          (i 2 1)
 365                                          ((1 ido) (1 2) (1 2))
 366                                          wa-%offset%)
 367                           di3)
 368                          (*
 369                           (f2cl-lib:fref wa-%data%
 370                                          (i 2 2)
 371                                          ((1 ido) (1 2) (1 2))
 372                                          wa-%offset%)
 373                           dr3)))
 374                 (setf (f2cl-lib:fref ch-%data%
 375                                      (1 k 3 i)
 376                                      ((1 in2) (1 l1) (1 3) (1 ido))
 377                                      ch-%offset%)
 378                         (-
 379                          (*
 380                           (f2cl-lib:fref wa-%data%
 381                                          (i 2 1)
 382                                          ((1 ido) (1 2) (1 2))
 383                                          wa-%offset%)
 384                           dr3)
 385                          (*
 386                           (f2cl-lib:fref wa-%data%
 387                                          (i 2 2)
 388                                          ((1 ido) (1 2) (1 2))
 389                                          wa-%offset%)
 390                           di3)))
 391                label104))
 392            label105))
 393         (go end_label)
 394        end_label
 395         (return (values nil nil nil nil nil nil nil nil))))))
 396
 397 (in-package #-gcl #:cl-user #+gcl "CL-USER")
 398 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
 399 (eval-when (:load-toplevel :compile-toplevel :execute)
 400   (setf (gethash 'fortran-to-lisp::c1f3kb
 401                  fortran-to-lisp::*f2cl-function-info*)
 402           (fortran-to-lisp::make-f2cl-finfo
 403            :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 404                         (fortran-to-lisp::integer4) (array double-float (*))
 405                         (fortran-to-lisp::integer4) (array double-float (*))
 406                         (fortran-to-lisp::integer4) (array double-float (*)))
 407            :return-values '(nil nil nil nil nil nil nil nil)
 408            :calls 'nil)))
 409