share/fftpack5/lisp/c1f3kf.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 c1f3kf (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              (sn 0.0d0))
  32         (declare (type (f2cl-lib:integer4) k i)
  33                  (type (double-float) sn tr2 cr2 ti2 ci2 cr3 ci3 dr2 dr3 di2
  34                                       di3))
  35         (if (> ido 1) (go label102))
  36         (setf sn (/ 1.0d0 (f2cl-lib:freal (f2cl-lib:int-mul 3 l1))))
  37         (if (= na 1) (go label106))
  38         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  39                       ((> k l1) nil)
  40           (tagbody
  41             (setf tr2
  42                     (+
  43                      (f2cl-lib:fref cc-%data%
  44                                     (1 k 1 2)
  45                                     ((1 in1) (1 l1) (1 ido) (1 3))
  46                                     cc-%offset%)
  47                      (f2cl-lib:fref cc-%data%
  48                                     (1 k 1 3)
  49                                     ((1 in1) (1 l1) (1 ido) (1 3))
  50                                     cc-%offset%)))
  51             (setf cr2
  52                     (+
  53                      (f2cl-lib:fref cc-%data%
  54                                     (1 k 1 1)
  55                                     ((1 in1) (1 l1) (1 ido) (1 3))
  56                                     cc-%offset%)
  57                      (* taur tr2)))
  58             (setf (f2cl-lib:fref cc-%data%
  59                                  (1 k 1 1)
  60                                  ((1 in1) (1 l1) (1 ido) (1 3))
  61                                  cc-%offset%)
  62                     (* sn
  63                        (+
  64                         (f2cl-lib:fref cc-%data%
  65                                        (1 k 1 1)
  66                                        ((1 in1) (1 l1) (1 ido) (1 3))
  67                                        cc-%offset%)
  68                         tr2)))
  69             (setf ti2
  70                     (+
  71                      (f2cl-lib:fref cc-%data%
  72                                     (2 k 1 2)
  73                                     ((1 in1) (1 l1) (1 ido) (1 3))
  74                                     cc-%offset%)
  75                      (f2cl-lib:fref cc-%data%
  76                                     (2 k 1 3)
  77                                     ((1 in1) (1 l1) (1 ido) (1 3))
  78                                     cc-%offset%)))
  79             (setf ci2
  80                     (+
  81                      (f2cl-lib:fref cc-%data%
  82                                     (2 k 1 1)
  83                                     ((1 in1) (1 l1) (1 ido) (1 3))
  84                                     cc-%offset%)
  85                      (* taur ti2)))
  86             (setf (f2cl-lib:fref cc-%data%
  87                                  (2 k 1 1)
  88                                  ((1 in1) (1 l1) (1 ido) (1 3))
  89                                  cc-%offset%)
  90                     (* sn
  91                        (+
  92                         (f2cl-lib:fref cc-%data%
  93                                        (2 k 1 1)
  94                                        ((1 in1) (1 l1) (1 ido) (1 3))
  95                                        cc-%offset%)
  96                         ti2)))
  97             (setf cr3
  98                     (* taui
  99                        (-
 100                         (f2cl-lib:fref cc-%data%
 101                                        (1 k 1 2)
 102                                        ((1 in1) (1 l1) (1 ido) (1 3))
 103                                        cc-%offset%)
 104                         (f2cl-lib:fref cc-%data%
 105                                        (1 k 1 3)
 106                                        ((1 in1) (1 l1) (1 ido) (1 3))
 107                                        cc-%offset%))))
 108             (setf ci3
 109                     (* taui
 110                        (-
 111                         (f2cl-lib:fref cc-%data%
 112                                        (2 k 1 2)
 113                                        ((1 in1) (1 l1) (1 ido) (1 3))
 114                                        cc-%offset%)
 115                         (f2cl-lib:fref cc-%data%
 116                                        (2 k 1 3)
 117                                        ((1 in1) (1 l1) (1 ido) (1 3))
 118                                        cc-%offset%))))
 119             (setf (f2cl-lib:fref cc-%data%
 120                                  (1 k 1 2)
 121                                  ((1 in1) (1 l1) (1 ido) (1 3))
 122                                  cc-%offset%)
 123                     (* sn (- cr2 ci3)))
 124             (setf (f2cl-lib:fref cc-%data%
 125                                  (1 k 1 3)
 126                                  ((1 in1) (1 l1) (1 ido) (1 3))
 127                                  cc-%offset%)
 128                     (* sn (+ cr2 ci3)))
 129             (setf (f2cl-lib:fref cc-%data%
 130                                  (2 k 1 2)
 131                                  ((1 in1) (1 l1) (1 ido) (1 3))
 132                                  cc-%offset%)
 133                     (* sn (+ ci2 cr3)))
 134             (setf (f2cl-lib:fref cc-%data%
 135                                  (2 k 1 3)
 136                                  ((1 in1) (1 l1) (1 ido) (1 3))
 137                                  cc-%offset%)
 138                     (* sn (- ci2 cr3)))
 139            label101))
 140         (go end_label)
 141        label106
 142         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 143                       ((> k l1) nil)
 144           (tagbody
 145             (setf tr2
 146                     (+
 147                      (f2cl-lib:fref cc-%data%
 148                                     (1 k 1 2)
 149                                     ((1 in1) (1 l1) (1 ido) (1 3))
 150                                     cc-%offset%)
 151                      (f2cl-lib:fref cc-%data%
 152                                     (1 k 1 3)
 153                                     ((1 in1) (1 l1) (1 ido) (1 3))
 154                                     cc-%offset%)))
 155             (setf cr2
 156                     (+
 157                      (f2cl-lib:fref cc-%data%
 158                                     (1 k 1 1)
 159                                     ((1 in1) (1 l1) (1 ido) (1 3))
 160                                     cc-%offset%)
 161                      (* taur tr2)))
 162             (setf (f2cl-lib:fref ch-%data%
 163                                  (1 k 1 1)
 164                                  ((1 in2) (1 l1) (1 3) (1 ido))
 165                                  ch-%offset%)
 166                     (* sn
 167                        (+
 168                         (f2cl-lib:fref cc-%data%
 169                                        (1 k 1 1)
 170                                        ((1 in1) (1 l1) (1 ido) (1 3))
 171                                        cc-%offset%)
 172                         tr2)))
 173             (setf ti2
 174                     (+
 175                      (f2cl-lib:fref cc-%data%
 176                                     (2 k 1 2)
 177                                     ((1 in1) (1 l1) (1 ido) (1 3))
 178                                     cc-%offset%)
 179                      (f2cl-lib:fref cc-%data%
 180                                     (2 k 1 3)
 181                                     ((1 in1) (1 l1) (1 ido) (1 3))
 182                                     cc-%offset%)))
 183             (setf ci2
 184                     (+
 185                      (f2cl-lib:fref cc-%data%
 186                                     (2 k 1 1)
 187                                     ((1 in1) (1 l1) (1 ido) (1 3))
 188                                     cc-%offset%)
 189                      (* taur ti2)))
 190             (setf (f2cl-lib:fref ch-%data%
 191                                  (2 k 1 1)
 192                                  ((1 in2) (1 l1) (1 3) (1 ido))
 193                                  ch-%offset%)
 194                     (* sn
 195                        (+
 196                         (f2cl-lib:fref cc-%data%
 197                                        (2 k 1 1)
 198                                        ((1 in1) (1 l1) (1 ido) (1 3))
 199                                        cc-%offset%)
 200                         ti2)))
 201             (setf cr3
 202                     (* taui
 203                        (-
 204                         (f2cl-lib:fref cc-%data%
 205                                        (1 k 1 2)
 206                                        ((1 in1) (1 l1) (1 ido) (1 3))
 207                                        cc-%offset%)
 208                         (f2cl-lib:fref cc-%data%
 209                                        (1 k 1 3)
 210                                        ((1 in1) (1 l1) (1 ido) (1 3))
 211                                        cc-%offset%))))
 212             (setf ci3
 213                     (* taui
 214                        (-
 215                         (f2cl-lib:fref cc-%data%
 216                                        (2 k 1 2)
 217                                        ((1 in1) (1 l1) (1 ido) (1 3))
 218                                        cc-%offset%)
 219                         (f2cl-lib:fref cc-%data%
 220                                        (2 k 1 3)
 221                                        ((1 in1) (1 l1) (1 ido) (1 3))
 222                                        cc-%offset%))))
 223             (setf (f2cl-lib:fref ch-%data%
 224                                  (1 k 2 1)
 225                                  ((1 in2) (1 l1) (1 3) (1 ido))
 226                                  ch-%offset%)
 227                     (* sn (- cr2 ci3)))
 228             (setf (f2cl-lib:fref ch-%data%
 229                                  (1 k 3 1)
 230                                  ((1 in2) (1 l1) (1 3) (1 ido))
 231                                  ch-%offset%)
 232                     (* sn (+ cr2 ci3)))
 233             (setf (f2cl-lib:fref ch-%data%
 234                                  (2 k 2 1)
 235                                  ((1 in2) (1 l1) (1 3) (1 ido))
 236                                  ch-%offset%)
 237                     (* sn (+ ci2 cr3)))
 238             (setf (f2cl-lib:fref ch-%data%
 239                                  (2 k 3 1)
 240                                  ((1 in2) (1 l1) (1 3) (1 ido))
 241                                  ch-%offset%)
 242                     (* sn (- ci2 cr3)))
 243            label107))
 244         (go end_label)
 245        label102
 246         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 247                       ((> k l1) nil)
 248           (tagbody
 249             (setf tr2
 250                     (+
 251                      (f2cl-lib:fref cc-%data%
 252                                     (1 k 1 2)
 253                                     ((1 in1) (1 l1) (1 ido) (1 3))
 254                                     cc-%offset%)
 255                      (f2cl-lib:fref cc-%data%
 256                                     (1 k 1 3)
 257                                     ((1 in1) (1 l1) (1 ido) (1 3))
 258                                     cc-%offset%)))
 259             (setf cr2
 260                     (+
 261                      (f2cl-lib:fref cc-%data%
 262                                     (1 k 1 1)
 263                                     ((1 in1) (1 l1) (1 ido) (1 3))
 264                                     cc-%offset%)
 265                      (* taur tr2)))
 266             (setf (f2cl-lib:fref ch-%data%
 267                                  (1 k 1 1)
 268                                  ((1 in2) (1 l1) (1 3) (1 ido))
 269                                  ch-%offset%)
 270                     (+
 271                      (f2cl-lib:fref cc-%data%
 272                                     (1 k 1 1)
 273                                     ((1 in1) (1 l1) (1 ido) (1 3))
 274                                     cc-%offset%)
 275                      tr2))
 276             (setf ti2
 277                     (+
 278                      (f2cl-lib:fref cc-%data%
 279                                     (2 k 1 2)
 280                                     ((1 in1) (1 l1) (1 ido) (1 3))
 281                                     cc-%offset%)
 282                      (f2cl-lib:fref cc-%data%
 283                                     (2 k 1 3)
 284                                     ((1 in1) (1 l1) (1 ido) (1 3))
 285                                     cc-%offset%)))
 286             (setf ci2
 287                     (+
 288                      (f2cl-lib:fref cc-%data%
 289                                     (2 k 1 1)
 290                                     ((1 in1) (1 l1) (1 ido) (1 3))
 291                                     cc-%offset%)
 292                      (* taur ti2)))
 293             (setf (f2cl-lib:fref ch-%data%
 294                                  (2 k 1 1)
 295                                  ((1 in2) (1 l1) (1 3) (1 ido))
 296                                  ch-%offset%)
 297                     (+
 298                      (f2cl-lib:fref cc-%data%
 299                                     (2 k 1 1)
 300                                     ((1 in1) (1 l1) (1 ido) (1 3))
 301                                     cc-%offset%)
 302                      ti2))
 303             (setf cr3
 304                     (* taui
 305                        (-
 306                         (f2cl-lib:fref cc-%data%
 307                                        (1 k 1 2)
 308                                        ((1 in1) (1 l1) (1 ido) (1 3))
 309                                        cc-%offset%)
 310                         (f2cl-lib:fref cc-%data%
 311                                        (1 k 1 3)
 312                                        ((1 in1) (1 l1) (1 ido) (1 3))
 313                                        cc-%offset%))))
 314             (setf ci3
 315                     (* taui
 316                        (-
 317                         (f2cl-lib:fref cc-%data%
 318                                        (2 k 1 2)
 319                                        ((1 in1) (1 l1) (1 ido) (1 3))
 320                                        cc-%offset%)
 321                         (f2cl-lib:fref cc-%data%
 322                                        (2 k 1 3)
 323                                        ((1 in1) (1 l1) (1 ido) (1 3))
 324                                        cc-%offset%))))
 325             (setf (f2cl-lib:fref ch-%data%
 326                                  (1 k 2 1)
 327                                  ((1 in2) (1 l1) (1 3) (1 ido))
 328                                  ch-%offset%)
 329                     (- cr2 ci3))
 330             (setf (f2cl-lib:fref ch-%data%
 331                                  (1 k 3 1)
 332                                  ((1 in2) (1 l1) (1 3) (1 ido))
 333                                  ch-%offset%)
 334                     (+ cr2 ci3))
 335             (setf (f2cl-lib:fref ch-%data%
 336                                  (2 k 2 1)
 337                                  ((1 in2) (1 l1) (1 3) (1 ido))
 338                                  ch-%offset%)
 339                     (+ ci2 cr3))
 340             (setf (f2cl-lib:fref ch-%data%
 341                                  (2 k 3 1)
 342                                  ((1 in2) (1 l1) (1 3) (1 ido))
 343                                  ch-%offset%)
 344                     (- ci2 cr3))
 345            label103))
 346         (f2cl-lib:fdo (i 2 (f2cl-lib:int-add i 1))
 347                       ((> i ido) nil)
 348           (tagbody
 349             (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 350                           ((> k l1) nil)
 351               (tagbody
 352                 (setf tr2
 353                         (+
 354                          (f2cl-lib:fref cc-%data%
 355                                         (1 k i 2)
 356                                         ((1 in1) (1 l1) (1 ido) (1 3))
 357                                         cc-%offset%)
 358                          (f2cl-lib:fref cc-%data%
 359                                         (1 k i 3)
 360                                         ((1 in1) (1 l1) (1 ido) (1 3))
 361                                         cc-%offset%)))
 362                 (setf cr2
 363                         (+
 364                          (f2cl-lib:fref cc-%data%
 365                                         (1 k i 1)
 366                                         ((1 in1) (1 l1) (1 ido) (1 3))
 367                                         cc-%offset%)
 368                          (* taur tr2)))
 369                 (setf (f2cl-lib:fref ch-%data%
 370                                      (1 k 1 i)
 371                                      ((1 in2) (1 l1) (1 3) (1 ido))
 372                                      ch-%offset%)
 373                         (+
 374                          (f2cl-lib:fref cc-%data%
 375                                         (1 k i 1)
 376                                         ((1 in1) (1 l1) (1 ido) (1 3))
 377                                         cc-%offset%)
 378                          tr2))
 379                 (setf ti2
 380                         (+
 381                          (f2cl-lib:fref cc-%data%
 382                                         (2 k i 2)
 383                                         ((1 in1) (1 l1) (1 ido) (1 3))
 384                                         cc-%offset%)
 385                          (f2cl-lib:fref cc-%data%
 386                                         (2 k i 3)
 387                                         ((1 in1) (1 l1) (1 ido) (1 3))
 388                                         cc-%offset%)))
 389                 (setf ci2
 390                         (+
 391                          (f2cl-lib:fref cc-%data%
 392                                         (2 k i 1)
 393                                         ((1 in1) (1 l1) (1 ido) (1 3))
 394                                         cc-%offset%)
 395                          (* taur ti2)))
 396                 (setf (f2cl-lib:fref ch-%data%
 397                                      (2 k 1 i)
 398                                      ((1 in2) (1 l1) (1 3) (1 ido))
 399                                      ch-%offset%)
 400                         (+
 401                          (f2cl-lib:fref cc-%data%
 402                                         (2 k i 1)
 403                                         ((1 in1) (1 l1) (1 ido) (1 3))
 404                                         cc-%offset%)
 405                          ti2))
 406                 (setf cr3
 407                         (* taui
 408                            (-
 409                             (f2cl-lib:fref cc-%data%
 410                                            (1 k i 2)
 411                                            ((1 in1) (1 l1) (1 ido) (1 3))
 412                                            cc-%offset%)
 413                             (f2cl-lib:fref cc-%data%
 414                                            (1 k i 3)
 415                                            ((1 in1) (1 l1) (1 ido) (1 3))
 416                                            cc-%offset%))))
 417                 (setf ci3
 418                         (* taui
 419                            (-
 420                             (f2cl-lib:fref cc-%data%
 421                                            (2 k i 2)
 422                                            ((1 in1) (1 l1) (1 ido) (1 3))
 423                                            cc-%offset%)
 424                             (f2cl-lib:fref cc-%data%
 425                                            (2 k i 3)
 426                                            ((1 in1) (1 l1) (1 ido) (1 3))
 427                                            cc-%offset%))))
 428                 (setf dr2 (- cr2 ci3))
 429                 (setf dr3 (+ cr2 ci3))
 430                 (setf di2 (+ ci2 cr3))
 431                 (setf di3 (- ci2 cr3))
 432                 (setf (f2cl-lib:fref ch-%data%
 433                                      (2 k 2 i)
 434                                      ((1 in2) (1 l1) (1 3) (1 ido))
 435                                      ch-%offset%)
 436                         (-
 437                          (*
 438                           (f2cl-lib:fref wa-%data%
 439                                          (i 1 1)
 440                                          ((1 ido) (1 2) (1 2))
 441                                          wa-%offset%)
 442                           di2)
 443                          (*
 444                           (f2cl-lib:fref wa-%data%
 445                                          (i 1 2)
 446                                          ((1 ido) (1 2) (1 2))
 447                                          wa-%offset%)
 448                           dr2)))
 449                 (setf (f2cl-lib:fref ch-%data%
 450                                      (1 k 2 i)
 451                                      ((1 in2) (1 l1) (1 3) (1 ido))
 452                                      ch-%offset%)
 453                         (+
 454                          (*
 455                           (f2cl-lib:fref wa-%data%
 456                                          (i 1 1)
 457                                          ((1 ido) (1 2) (1 2))
 458                                          wa-%offset%)
 459                           dr2)
 460                          (*
 461                           (f2cl-lib:fref wa-%data%
 462                                          (i 1 2)
 463                                          ((1 ido) (1 2) (1 2))
 464                                          wa-%offset%)
 465                           di2)))
 466                 (setf (f2cl-lib:fref ch-%data%
 467                                      (2 k 3 i)
 468                                      ((1 in2) (1 l1) (1 3) (1 ido))
 469                                      ch-%offset%)
 470                         (-
 471                          (*
 472                           (f2cl-lib:fref wa-%data%
 473                                          (i 2 1)
 474                                          ((1 ido) (1 2) (1 2))
 475                                          wa-%offset%)
 476                           di3)
 477                          (*
 478                           (f2cl-lib:fref wa-%data%
 479                                          (i 2 2)
 480                                          ((1 ido) (1 2) (1 2))
 481                                          wa-%offset%)
 482                           dr3)))
 483                 (setf (f2cl-lib:fref ch-%data%
 484                                      (1 k 3 i)
 485                                      ((1 in2) (1 l1) (1 3) (1 ido))
 486                                      ch-%offset%)
 487                         (+
 488                          (*
 489                           (f2cl-lib:fref wa-%data%
 490                                          (i 2 1)
 491                                          ((1 ido) (1 2) (1 2))
 492                                          wa-%offset%)
 493                           dr3)
 494                          (*
 495                           (f2cl-lib:fref wa-%data%
 496                                          (i 2 2)
 497                                          ((1 ido) (1 2) (1 2))
 498                                          wa-%offset%)
 499                           di3)))
 500                label104))
 501            label105))
 502         (go end_label)
 503        end_label
 504         (return (values nil nil nil nil nil nil nil nil))))))
 505
 506 (in-package #:cl-user)
 507 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
 508 (eval-when (:load-toplevel :compile-toplevel :execute)
 509   (setf (gethash 'fortran-to-lisp::c1f3kf
 510                  fortran-to-lisp::*f2cl-function-info*)
 511           (fortran-to-lisp::make-f2cl-finfo
 512            :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 513                         (fortran-to-lisp::integer4) (array double-float (*))
 514                         (fortran-to-lisp::integer4) (array double-float (*))
 515                         (fortran-to-lisp::integer4) (array double-float (*)))
 516            :return-values '(nil nil nil nil nil nil nil nil)
 517            :calls 'nil)))
 518