share/fftpack5/lisp/c1f5kf.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 ((tr11 0.3090169943749474d0)
  21       (ti11 -0.9510565162951536d0)
  22       (tr12 -0.8090169943749475d0)
  23       (ti12 -0.5877852522924731d0))
  24   (declare (type (double-float) tr11 ti11 tr12 ti12))
  25   (defun c1f5kf (ido l1 na cc in1 ch in2 wa)
  26     (declare (type (array double-float (*)) wa ch cc)
  27              (type (f2cl-lib:integer4) in2 in1 na l1 ido))
  28     (f2cl-lib:with-multi-array-data
  29         ((cc double-float cc-%data% cc-%offset%)
  30          (ch double-float ch-%data% ch-%offset%)
  31          (wa double-float wa-%data% wa-%offset%))
  32       (prog ((di2 0.0d0) (di5 0.0d0) (dr2 0.0d0) (dr5 0.0d0) (di4 0.0d0)
  33              (di3 0.0d0) (dr4 0.0d0) (dr3 0.0d0) (i 0) (ci4 0.0d0) (cr4 0.0d0)
  34              (ci5 0.0d0) (cr5 0.0d0) (ci3 0.0d0) (cr3 0.0d0) (ci2 0.0d0)
  35              (cr2 0.0d0) (chold2 0.0d0) (chold1 0.0d0) (tr3 0.0d0) (tr4 0.0d0)
  36              (tr2 0.0d0) (tr5 0.0d0) (ti3 0.0d0) (ti4 0.0d0) (ti2 0.0d0)
  37              (ti5 0.0d0) (k 0) (sn 0.0d0))
  38         (declare (type (f2cl-lib:integer4) k i)
  39                  (type (double-float) sn ti5 ti2 ti4 ti3 tr5 tr2 tr4 tr3 chold1
  40                                       chold2 cr2 ci2 cr3 ci3 cr5 ci5 cr4 ci4
  41                                       dr3 dr4 di3 di4 dr5 dr2 di5 di2))
  42         (if (> ido 1) (go label102))
  43         (setf sn (/ 1.0d0 (f2cl-lib:freal (f2cl-lib:int-mul 5 l1))))
  44         (if (= na 1) (go label106))
  45         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  46                       ((> k l1) nil)
  47           (tagbody
  48             (setf ti5
  49                     (-
  50                      (f2cl-lib:fref cc-%data%
  51                                     (2 k 1 2)
  52                                     ((1 in1) (1 l1) (1 ido) (1 5))
  53                                     cc-%offset%)
  54                      (f2cl-lib:fref cc-%data%
  55                                     (2 k 1 5)
  56                                     ((1 in1) (1 l1) (1 ido) (1 5))
  57                                     cc-%offset%)))
  58             (setf ti2
  59                     (+
  60                      (f2cl-lib:fref cc-%data%
  61                                     (2 k 1 2)
  62                                     ((1 in1) (1 l1) (1 ido) (1 5))
  63                                     cc-%offset%)
  64                      (f2cl-lib:fref cc-%data%
  65                                     (2 k 1 5)
  66                                     ((1 in1) (1 l1) (1 ido) (1 5))
  67                                     cc-%offset%)))
  68             (setf ti4
  69                     (-
  70                      (f2cl-lib:fref cc-%data%
  71                                     (2 k 1 3)
  72                                     ((1 in1) (1 l1) (1 ido) (1 5))
  73                                     cc-%offset%)
  74                      (f2cl-lib:fref cc-%data%
  75                                     (2 k 1 4)
  76                                     ((1 in1) (1 l1) (1 ido) (1 5))
  77                                     cc-%offset%)))
  78             (setf ti3
  79                     (+
  80                      (f2cl-lib:fref cc-%data%
  81                                     (2 k 1 3)
  82                                     ((1 in1) (1 l1) (1 ido) (1 5))
  83                                     cc-%offset%)
  84                      (f2cl-lib:fref cc-%data%
  85                                     (2 k 1 4)
  86                                     ((1 in1) (1 l1) (1 ido) (1 5))
  87                                     cc-%offset%)))
  88             (setf tr5
  89                     (-
  90                      (f2cl-lib:fref cc-%data%
  91                                     (1 k 1 2)
  92                                     ((1 in1) (1 l1) (1 ido) (1 5))
  93                                     cc-%offset%)
  94                      (f2cl-lib:fref cc-%data%
  95                                     (1 k 1 5)
  96                                     ((1 in1) (1 l1) (1 ido) (1 5))
  97                                     cc-%offset%)))
  98             (setf tr2
  99                     (+
 100                      (f2cl-lib:fref cc-%data%
 101                                     (1 k 1 2)
 102                                     ((1 in1) (1 l1) (1 ido) (1 5))
 103                                     cc-%offset%)
 104                      (f2cl-lib:fref cc-%data%
 105                                     (1 k 1 5)
 106                                     ((1 in1) (1 l1) (1 ido) (1 5))
 107                                     cc-%offset%)))
 108             (setf tr4
 109                     (-
 110                      (f2cl-lib:fref cc-%data%
 111                                     (1 k 1 3)
 112                                     ((1 in1) (1 l1) (1 ido) (1 5))
 113                                     cc-%offset%)
 114                      (f2cl-lib:fref cc-%data%
 115                                     (1 k 1 4)
 116                                     ((1 in1) (1 l1) (1 ido) (1 5))
 117                                     cc-%offset%)))
 118             (setf tr3
 119                     (+
 120                      (f2cl-lib:fref cc-%data%
 121                                     (1 k 1 3)
 122                                     ((1 in1) (1 l1) (1 ido) (1 5))
 123                                     cc-%offset%)
 124                      (f2cl-lib:fref cc-%data%
 125                                     (1 k 1 4)
 126                                     ((1 in1) (1 l1) (1 ido) (1 5))
 127                                     cc-%offset%)))
 128             (setf chold1
 129                     (* sn
 130                        (+
 131                         (f2cl-lib:fref cc-%data%
 132                                        (1 k 1 1)
 133                                        ((1 in1) (1 l1) (1 ido) (1 5))
 134                                        cc-%offset%)
 135                         tr2
 136                         tr3)))
 137             (setf chold2
 138                     (* sn
 139                        (+
 140                         (f2cl-lib:fref cc-%data%
 141                                        (2 k 1 1)
 142                                        ((1 in1) (1 l1) (1 ido) (1 5))
 143                                        cc-%offset%)
 144                         ti2
 145                         ti3)))
 146             (setf cr2
 147                     (+
 148                      (f2cl-lib:fref cc-%data%
 149                                     (1 k 1 1)
 150                                     ((1 in1) (1 l1) (1 ido) (1 5))
 151                                     cc-%offset%)
 152                      (* tr11 tr2)
 153                      (* tr12 tr3)))
 154             (setf ci2
 155                     (+
 156                      (f2cl-lib:fref cc-%data%
 157                                     (2 k 1 1)
 158                                     ((1 in1) (1 l1) (1 ido) (1 5))
 159                                     cc-%offset%)
 160                      (* tr11 ti2)
 161                      (* tr12 ti3)))
 162             (setf cr3
 163                     (+
 164                      (f2cl-lib:fref cc-%data%
 165                                     (1 k 1 1)
 166                                     ((1 in1) (1 l1) (1 ido) (1 5))
 167                                     cc-%offset%)
 168                      (* tr12 tr2)
 169                      (* tr11 tr3)))
 170             (setf ci3
 171                     (+
 172                      (f2cl-lib:fref cc-%data%
 173                                     (2 k 1 1)
 174                                     ((1 in1) (1 l1) (1 ido) (1 5))
 175                                     cc-%offset%)
 176                      (* tr12 ti2)
 177                      (* tr11 ti3)))
 178             (setf (f2cl-lib:fref cc-%data%
 179                                  (1 k 1 1)
 180                                  ((1 in1) (1 l1) (1 ido) (1 5))
 181                                  cc-%offset%)
 182                     chold1)
 183             (setf (f2cl-lib:fref cc-%data%
 184                                  (2 k 1 1)
 185                                  ((1 in1) (1 l1) (1 ido) (1 5))
 186                                  cc-%offset%)
 187                     chold2)
 188             (setf cr5 (+ (* ti11 tr5) (* ti12 tr4)))
 189             (setf ci5 (+ (* ti11 ti5) (* ti12 ti4)))
 190             (setf cr4 (- (* ti12 tr5) (* ti11 tr4)))
 191             (setf ci4 (- (* ti12 ti5) (* ti11 ti4)))
 192             (setf (f2cl-lib:fref cc-%data%
 193                                  (1 k 1 2)
 194                                  ((1 in1) (1 l1) (1 ido) (1 5))
 195                                  cc-%offset%)
 196                     (* sn (- cr2 ci5)))
 197             (setf (f2cl-lib:fref cc-%data%
 198                                  (1 k 1 5)
 199                                  ((1 in1) (1 l1) (1 ido) (1 5))
 200                                  cc-%offset%)
 201                     (* sn (+ cr2 ci5)))
 202             (setf (f2cl-lib:fref cc-%data%
 203                                  (2 k 1 2)
 204                                  ((1 in1) (1 l1) (1 ido) (1 5))
 205                                  cc-%offset%)
 206                     (* sn (+ ci2 cr5)))
 207             (setf (f2cl-lib:fref cc-%data%
 208                                  (2 k 1 3)
 209                                  ((1 in1) (1 l1) (1 ido) (1 5))
 210                                  cc-%offset%)
 211                     (* sn (+ ci3 cr4)))
 212             (setf (f2cl-lib:fref cc-%data%
 213                                  (1 k 1 3)
 214                                  ((1 in1) (1 l1) (1 ido) (1 5))
 215                                  cc-%offset%)
 216                     (* sn (- cr3 ci4)))
 217             (setf (f2cl-lib:fref cc-%data%
 218                                  (1 k 1 4)
 219                                  ((1 in1) (1 l1) (1 ido) (1 5))
 220                                  cc-%offset%)
 221                     (* sn (+ cr3 ci4)))
 222             (setf (f2cl-lib:fref cc-%data%
 223                                  (2 k 1 4)
 224                                  ((1 in1) (1 l1) (1 ido) (1 5))
 225                                  cc-%offset%)
 226                     (* sn (- ci3 cr4)))
 227             (setf (f2cl-lib:fref cc-%data%
 228                                  (2 k 1 5)
 229                                  ((1 in1) (1 l1) (1 ido) (1 5))
 230                                  cc-%offset%)
 231                     (* sn (- ci2 cr5)))
 232            label101))
 233         (go end_label)
 234        label106
 235         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 236                       ((> k l1) nil)
 237           (tagbody
 238             (setf ti5
 239                     (-
 240                      (f2cl-lib:fref cc-%data%
 241                                     (2 k 1 2)
 242                                     ((1 in1) (1 l1) (1 ido) (1 5))
 243                                     cc-%offset%)
 244                      (f2cl-lib:fref cc-%data%
 245                                     (2 k 1 5)
 246                                     ((1 in1) (1 l1) (1 ido) (1 5))
 247                                     cc-%offset%)))
 248             (setf ti2
 249                     (+
 250                      (f2cl-lib:fref cc-%data%
 251                                     (2 k 1 2)
 252                                     ((1 in1) (1 l1) (1 ido) (1 5))
 253                                     cc-%offset%)
 254                      (f2cl-lib:fref cc-%data%
 255                                     (2 k 1 5)
 256                                     ((1 in1) (1 l1) (1 ido) (1 5))
 257                                     cc-%offset%)))
 258             (setf ti4
 259                     (-
 260                      (f2cl-lib:fref cc-%data%
 261                                     (2 k 1 3)
 262                                     ((1 in1) (1 l1) (1 ido) (1 5))
 263                                     cc-%offset%)
 264                      (f2cl-lib:fref cc-%data%
 265                                     (2 k 1 4)
 266                                     ((1 in1) (1 l1) (1 ido) (1 5))
 267                                     cc-%offset%)))
 268             (setf ti3
 269                     (+
 270                      (f2cl-lib:fref cc-%data%
 271                                     (2 k 1 3)
 272                                     ((1 in1) (1 l1) (1 ido) (1 5))
 273                                     cc-%offset%)
 274                      (f2cl-lib:fref cc-%data%
 275                                     (2 k 1 4)
 276                                     ((1 in1) (1 l1) (1 ido) (1 5))
 277                                     cc-%offset%)))
 278             (setf tr5
 279                     (-
 280                      (f2cl-lib:fref cc-%data%
 281                                     (1 k 1 2)
 282                                     ((1 in1) (1 l1) (1 ido) (1 5))
 283                                     cc-%offset%)
 284                      (f2cl-lib:fref cc-%data%
 285                                     (1 k 1 5)
 286                                     ((1 in1) (1 l1) (1 ido) (1 5))
 287                                     cc-%offset%)))
 288             (setf tr2
 289                     (+
 290                      (f2cl-lib:fref cc-%data%
 291                                     (1 k 1 2)
 292                                     ((1 in1) (1 l1) (1 ido) (1 5))
 293                                     cc-%offset%)
 294                      (f2cl-lib:fref cc-%data%
 295                                     (1 k 1 5)
 296                                     ((1 in1) (1 l1) (1 ido) (1 5))
 297                                     cc-%offset%)))
 298             (setf tr4
 299                     (-
 300                      (f2cl-lib:fref cc-%data%
 301                                     (1 k 1 3)
 302                                     ((1 in1) (1 l1) (1 ido) (1 5))
 303                                     cc-%offset%)
 304                      (f2cl-lib:fref cc-%data%
 305                                     (1 k 1 4)
 306                                     ((1 in1) (1 l1) (1 ido) (1 5))
 307                                     cc-%offset%)))
 308             (setf tr3
 309                     (+
 310                      (f2cl-lib:fref cc-%data%
 311                                     (1 k 1 3)
 312                                     ((1 in1) (1 l1) (1 ido) (1 5))
 313                                     cc-%offset%)
 314                      (f2cl-lib:fref cc-%data%
 315                                     (1 k 1 4)
 316                                     ((1 in1) (1 l1) (1 ido) (1 5))
 317                                     cc-%offset%)))
 318             (setf (f2cl-lib:fref ch-%data%
 319                                  (1 k 1 1)
 320                                  ((1 in2) (1 l1) (1 5) (1 ido))
 321                                  ch-%offset%)
 322                     (* sn
 323                        (+
 324                         (f2cl-lib:fref cc-%data%
 325                                        (1 k 1 1)
 326                                        ((1 in1) (1 l1) (1 ido) (1 5))
 327                                        cc-%offset%)
 328                         tr2
 329                         tr3)))
 330             (setf (f2cl-lib:fref ch-%data%
 331                                  (2 k 1 1)
 332                                  ((1 in2) (1 l1) (1 5) (1 ido))
 333                                  ch-%offset%)
 334                     (* sn
 335                        (+
 336                         (f2cl-lib:fref cc-%data%
 337                                        (2 k 1 1)
 338                                        ((1 in1) (1 l1) (1 ido) (1 5))
 339                                        cc-%offset%)
 340                         ti2
 341                         ti3)))
 342             (setf cr2
 343                     (+
 344                      (f2cl-lib:fref cc-%data%
 345                                     (1 k 1 1)
 346                                     ((1 in1) (1 l1) (1 ido) (1 5))
 347                                     cc-%offset%)
 348                      (* tr11 tr2)
 349                      (* tr12 tr3)))
 350             (setf ci2
 351                     (+
 352                      (f2cl-lib:fref cc-%data%
 353                                     (2 k 1 1)
 354                                     ((1 in1) (1 l1) (1 ido) (1 5))
 355                                     cc-%offset%)
 356                      (* tr11 ti2)
 357                      (* tr12 ti3)))
 358             (setf cr3
 359                     (+
 360                      (f2cl-lib:fref cc-%data%
 361                                     (1 k 1 1)
 362                                     ((1 in1) (1 l1) (1 ido) (1 5))
 363                                     cc-%offset%)
 364                      (* tr12 tr2)
 365                      (* tr11 tr3)))
 366             (setf ci3
 367                     (+
 368                      (f2cl-lib:fref cc-%data%
 369                                     (2 k 1 1)
 370                                     ((1 in1) (1 l1) (1 ido) (1 5))
 371                                     cc-%offset%)
 372                      (* tr12 ti2)
 373                      (* tr11 ti3)))
 374             (setf cr5 (+ (* ti11 tr5) (* ti12 tr4)))
 375             (setf ci5 (+ (* ti11 ti5) (* ti12 ti4)))
 376             (setf cr4 (- (* ti12 tr5) (* ti11 tr4)))
 377             (setf ci4 (- (* ti12 ti5) (* ti11 ti4)))
 378             (setf (f2cl-lib:fref ch-%data%
 379                                  (1 k 2 1)
 380                                  ((1 in2) (1 l1) (1 5) (1 ido))
 381                                  ch-%offset%)
 382                     (* sn (- cr2 ci5)))
 383             (setf (f2cl-lib:fref ch-%data%
 384                                  (1 k 5 1)
 385                                  ((1 in2) (1 l1) (1 5) (1 ido))
 386                                  ch-%offset%)
 387                     (* sn (+ cr2 ci5)))
 388             (setf (f2cl-lib:fref ch-%data%
 389                                  (2 k 2 1)
 390                                  ((1 in2) (1 l1) (1 5) (1 ido))
 391                                  ch-%offset%)
 392                     (* sn (+ ci2 cr5)))
 393             (setf (f2cl-lib:fref ch-%data%
 394                                  (2 k 3 1)
 395                                  ((1 in2) (1 l1) (1 5) (1 ido))
 396                                  ch-%offset%)
 397                     (* sn (+ ci3 cr4)))
 398             (setf (f2cl-lib:fref ch-%data%
 399                                  (1 k 3 1)
 400                                  ((1 in2) (1 l1) (1 5) (1 ido))
 401                                  ch-%offset%)
 402                     (* sn (- cr3 ci4)))
 403             (setf (f2cl-lib:fref ch-%data%
 404                                  (1 k 4 1)
 405                                  ((1 in2) (1 l1) (1 5) (1 ido))
 406                                  ch-%offset%)
 407                     (* sn (+ cr3 ci4)))
 408             (setf (f2cl-lib:fref ch-%data%
 409                                  (2 k 4 1)
 410                                  ((1 in2) (1 l1) (1 5) (1 ido))
 411                                  ch-%offset%)
 412                     (* sn (- ci3 cr4)))
 413             (setf (f2cl-lib:fref ch-%data%
 414                                  (2 k 5 1)
 415                                  ((1 in2) (1 l1) (1 5) (1 ido))
 416                                  ch-%offset%)
 417                     (* sn (- ci2 cr5)))
 418            label107))
 419         (go end_label)
 420        label102
 421         (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 422                       ((> k l1) nil)
 423           (tagbody
 424             (setf ti5
 425                     (-
 426                      (f2cl-lib:fref cc-%data%
 427                                     (2 k 1 2)
 428                                     ((1 in1) (1 l1) (1 ido) (1 5))
 429                                     cc-%offset%)
 430                      (f2cl-lib:fref cc-%data%
 431                                     (2 k 1 5)
 432                                     ((1 in1) (1 l1) (1 ido) (1 5))
 433                                     cc-%offset%)))
 434             (setf ti2
 435                     (+
 436                      (f2cl-lib:fref cc-%data%
 437                                     (2 k 1 2)
 438                                     ((1 in1) (1 l1) (1 ido) (1 5))
 439                                     cc-%offset%)
 440                      (f2cl-lib:fref cc-%data%
 441                                     (2 k 1 5)
 442                                     ((1 in1) (1 l1) (1 ido) (1 5))
 443                                     cc-%offset%)))
 444             (setf ti4
 445                     (-
 446                      (f2cl-lib:fref cc-%data%
 447                                     (2 k 1 3)
 448                                     ((1 in1) (1 l1) (1 ido) (1 5))
 449                                     cc-%offset%)
 450                      (f2cl-lib:fref cc-%data%
 451                                     (2 k 1 4)
 452                                     ((1 in1) (1 l1) (1 ido) (1 5))
 453                                     cc-%offset%)))
 454             (setf ti3
 455                     (+
 456                      (f2cl-lib:fref cc-%data%
 457                                     (2 k 1 3)
 458                                     ((1 in1) (1 l1) (1 ido) (1 5))
 459                                     cc-%offset%)
 460                      (f2cl-lib:fref cc-%data%
 461                                     (2 k 1 4)
 462                                     ((1 in1) (1 l1) (1 ido) (1 5))
 463                                     cc-%offset%)))
 464             (setf tr5
 465                     (-
 466                      (f2cl-lib:fref cc-%data%
 467                                     (1 k 1 2)
 468                                     ((1 in1) (1 l1) (1 ido) (1 5))
 469                                     cc-%offset%)
 470                      (f2cl-lib:fref cc-%data%
 471                                     (1 k 1 5)
 472                                     ((1 in1) (1 l1) (1 ido) (1 5))
 473                                     cc-%offset%)))
 474             (setf tr2
 475                     (+
 476                      (f2cl-lib:fref cc-%data%
 477                                     (1 k 1 2)
 478                                     ((1 in1) (1 l1) (1 ido) (1 5))
 479                                     cc-%offset%)
 480                      (f2cl-lib:fref cc-%data%
 481                                     (1 k 1 5)
 482                                     ((1 in1) (1 l1) (1 ido) (1 5))
 483                                     cc-%offset%)))
 484             (setf tr4
 485                     (-
 486                      (f2cl-lib:fref cc-%data%
 487                                     (1 k 1 3)
 488                                     ((1 in1) (1 l1) (1 ido) (1 5))
 489                                     cc-%offset%)
 490                      (f2cl-lib:fref cc-%data%
 491                                     (1 k 1 4)
 492                                     ((1 in1) (1 l1) (1 ido) (1 5))
 493                                     cc-%offset%)))
 494             (setf tr3
 495                     (+
 496                      (f2cl-lib:fref cc-%data%
 497                                     (1 k 1 3)
 498                                     ((1 in1) (1 l1) (1 ido) (1 5))
 499                                     cc-%offset%)
 500                      (f2cl-lib:fref cc-%data%
 501                                     (1 k 1 4)
 502                                     ((1 in1) (1 l1) (1 ido) (1 5))
 503                                     cc-%offset%)))
 504             (setf (f2cl-lib:fref ch-%data%
 505                                  (1 k 1 1)
 506                                  ((1 in2) (1 l1) (1 5) (1 ido))
 507                                  ch-%offset%)
 508                     (+
 509                      (f2cl-lib:fref cc-%data%
 510                                     (1 k 1 1)
 511                                     ((1 in1) (1 l1) (1 ido) (1 5))
 512                                     cc-%offset%)
 513                      tr2
 514                      tr3))
 515             (setf (f2cl-lib:fref ch-%data%
 516                                  (2 k 1 1)
 517                                  ((1 in2) (1 l1) (1 5) (1 ido))
 518                                  ch-%offset%)
 519                     (+
 520                      (f2cl-lib:fref cc-%data%
 521                                     (2 k 1 1)
 522                                     ((1 in1) (1 l1) (1 ido) (1 5))
 523                                     cc-%offset%)
 524                      ti2
 525                      ti3))
 526             (setf cr2
 527                     (+
 528                      (f2cl-lib:fref cc-%data%
 529                                     (1 k 1 1)
 530                                     ((1 in1) (1 l1) (1 ido) (1 5))
 531                                     cc-%offset%)
 532                      (* tr11 tr2)
 533                      (* tr12 tr3)))
 534             (setf ci2
 535                     (+
 536                      (f2cl-lib:fref cc-%data%
 537                                     (2 k 1 1)
 538                                     ((1 in1) (1 l1) (1 ido) (1 5))
 539                                     cc-%offset%)
 540                      (* tr11 ti2)
 541                      (* tr12 ti3)))
 542             (setf cr3
 543                     (+
 544                      (f2cl-lib:fref cc-%data%
 545                                     (1 k 1 1)
 546                                     ((1 in1) (1 l1) (1 ido) (1 5))
 547                                     cc-%offset%)
 548                      (* tr12 tr2)
 549                      (* tr11 tr3)))
 550             (setf ci3
 551                     (+
 552                      (f2cl-lib:fref cc-%data%
 553                                     (2 k 1 1)
 554                                     ((1 in1) (1 l1) (1 ido) (1 5))
 555                                     cc-%offset%)
 556                      (* tr12 ti2)
 557                      (* tr11 ti3)))
 558             (setf cr5 (+ (* ti11 tr5) (* ti12 tr4)))
 559             (setf ci5 (+ (* ti11 ti5) (* ti12 ti4)))
 560             (setf cr4 (- (* ti12 tr5) (* ti11 tr4)))
 561             (setf ci4 (- (* ti12 ti5) (* ti11 ti4)))
 562             (setf (f2cl-lib:fref ch-%data%
 563                                  (1 k 2 1)
 564                                  ((1 in2) (1 l1) (1 5) (1 ido))
 565                                  ch-%offset%)
 566                     (- cr2 ci5))
 567             (setf (f2cl-lib:fref ch-%data%
 568                                  (1 k 5 1)
 569                                  ((1 in2) (1 l1) (1 5) (1 ido))
 570                                  ch-%offset%)
 571                     (+ cr2 ci5))
 572             (setf (f2cl-lib:fref ch-%data%
 573                                  (2 k 2 1)
 574                                  ((1 in2) (1 l1) (1 5) (1 ido))
 575                                  ch-%offset%)
 576                     (+ ci2 cr5))
 577             (setf (f2cl-lib:fref ch-%data%
 578                                  (2 k 3 1)
 579                                  ((1 in2) (1 l1) (1 5) (1 ido))
 580                                  ch-%offset%)
 581                     (+ ci3 cr4))
 582             (setf (f2cl-lib:fref ch-%data%
 583                                  (1 k 3 1)
 584                                  ((1 in2) (1 l1) (1 5) (1 ido))
 585                                  ch-%offset%)
 586                     (- cr3 ci4))
 587             (setf (f2cl-lib:fref ch-%data%
 588                                  (1 k 4 1)
 589                                  ((1 in2) (1 l1) (1 5) (1 ido))
 590                                  ch-%offset%)
 591                     (+ cr3 ci4))
 592             (setf (f2cl-lib:fref ch-%data%
 593                                  (2 k 4 1)
 594                                  ((1 in2) (1 l1) (1 5) (1 ido))
 595                                  ch-%offset%)
 596                     (- ci3 cr4))
 597             (setf (f2cl-lib:fref ch-%data%
 598                                  (2 k 5 1)
 599                                  ((1 in2) (1 l1) (1 5) (1 ido))
 600                                  ch-%offset%)
 601                     (- ci2 cr5))
 602            label103))
 603         (f2cl-lib:fdo (i 2 (f2cl-lib:int-add i 1))
 604                       ((> i ido) nil)
 605           (tagbody
 606             (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 607                           ((> k l1) nil)
 608               (tagbody
 609                 (setf ti5
 610                         (-
 611                          (f2cl-lib:fref cc-%data%
 612                                         (2 k i 2)
 613                                         ((1 in1) (1 l1) (1 ido) (1 5))
 614                                         cc-%offset%)
 615                          (f2cl-lib:fref cc-%data%
 616                                         (2 k i 5)
 617                                         ((1 in1) (1 l1) (1 ido) (1 5))
 618                                         cc-%offset%)))
 619                 (setf ti2
 620                         (+
 621                          (f2cl-lib:fref cc-%data%
 622                                         (2 k i 2)
 623                                         ((1 in1) (1 l1) (1 ido) (1 5))
 624                                         cc-%offset%)
 625                          (f2cl-lib:fref cc-%data%
 626                                         (2 k i 5)
 627                                         ((1 in1) (1 l1) (1 ido) (1 5))
 628                                         cc-%offset%)))
 629                 (setf ti4
 630                         (-
 631                          (f2cl-lib:fref cc-%data%
 632                                         (2 k i 3)
 633                                         ((1 in1) (1 l1) (1 ido) (1 5))
 634                                         cc-%offset%)
 635                          (f2cl-lib:fref cc-%data%
 636                                         (2 k i 4)
 637                                         ((1 in1) (1 l1) (1 ido) (1 5))
 638                                         cc-%offset%)))
 639                 (setf ti3
 640                         (+
 641                          (f2cl-lib:fref cc-%data%
 642                                         (2 k i 3)
 643                                         ((1 in1) (1 l1) (1 ido) (1 5))
 644                                         cc-%offset%)
 645                          (f2cl-lib:fref cc-%data%
 646                                         (2 k i 4)
 647                                         ((1 in1) (1 l1) (1 ido) (1 5))
 648                                         cc-%offset%)))
 649                 (setf tr5
 650                         (-
 651                          (f2cl-lib:fref cc-%data%
 652                                         (1 k i 2)
 653                                         ((1 in1) (1 l1) (1 ido) (1 5))
 654                                         cc-%offset%)
 655                          (f2cl-lib:fref cc-%data%
 656                                         (1 k i 5)
 657                                         ((1 in1) (1 l1) (1 ido) (1 5))
 658                                         cc-%offset%)))
 659                 (setf tr2
 660                         (+
 661                          (f2cl-lib:fref cc-%data%
 662                                         (1 k i 2)
 663                                         ((1 in1) (1 l1) (1 ido) (1 5))
 664                                         cc-%offset%)
 665                          (f2cl-lib:fref cc-%data%
 666                                         (1 k i 5)
 667                                         ((1 in1) (1 l1) (1 ido) (1 5))
 668                                         cc-%offset%)))
 669                 (setf tr4
 670                         (-
 671                          (f2cl-lib:fref cc-%data%
 672                                         (1 k i 3)
 673                                         ((1 in1) (1 l1) (1 ido) (1 5))
 674                                         cc-%offset%)
 675                          (f2cl-lib:fref cc-%data%
 676                                         (1 k i 4)
 677                                         ((1 in1) (1 l1) (1 ido) (1 5))
 678                                         cc-%offset%)))
 679                 (setf tr3
 680                         (+
 681                          (f2cl-lib:fref cc-%data%
 682                                         (1 k i 3)
 683                                         ((1 in1) (1 l1) (1 ido) (1 5))
 684                                         cc-%offset%)
 685                          (f2cl-lib:fref cc-%data%
 686                                         (1 k i 4)
 687                                         ((1 in1) (1 l1) (1 ido) (1 5))
 688                                         cc-%offset%)))
 689                 (setf (f2cl-lib:fref ch-%data%
 690                                      (1 k 1 i)
 691                                      ((1 in2) (1 l1) (1 5) (1 ido))
 692                                      ch-%offset%)
 693                         (+
 694                          (f2cl-lib:fref cc-%data%
 695                                         (1 k i 1)
 696                                         ((1 in1) (1 l1) (1 ido) (1 5))
 697                                         cc-%offset%)
 698                          tr2
 699                          tr3))
 700                 (setf (f2cl-lib:fref ch-%data%
 701                                      (2 k 1 i)
 702                                      ((1 in2) (1 l1) (1 5) (1 ido))
 703                                      ch-%offset%)
 704                         (+
 705                          (f2cl-lib:fref cc-%data%
 706                                         (2 k i 1)
 707                                         ((1 in1) (1 l1) (1 ido) (1 5))
 708                                         cc-%offset%)
 709                          ti2
 710                          ti3))
 711                 (setf cr2
 712                         (+
 713                          (f2cl-lib:fref cc-%data%
 714                                         (1 k i 1)
 715                                         ((1 in1) (1 l1) (1 ido) (1 5))
 716                                         cc-%offset%)
 717                          (* tr11 tr2)
 718                          (* tr12 tr3)))
 719                 (setf ci2
 720                         (+
 721                          (f2cl-lib:fref cc-%data%
 722                                         (2 k i 1)
 723                                         ((1 in1) (1 l1) (1 ido) (1 5))
 724                                         cc-%offset%)
 725                          (* tr11 ti2)
 726                          (* tr12 ti3)))
 727                 (setf cr3
 728                         (+
 729                          (f2cl-lib:fref cc-%data%
 730                                         (1 k i 1)
 731                                         ((1 in1) (1 l1) (1 ido) (1 5))
 732                                         cc-%offset%)
 733                          (* tr12 tr2)
 734                          (* tr11 tr3)))
 735                 (setf ci3
 736                         (+
 737                          (f2cl-lib:fref cc-%data%
 738                                         (2 k i 1)
 739                                         ((1 in1) (1 l1) (1 ido) (1 5))
 740                                         cc-%offset%)
 741                          (* tr12 ti2)
 742                          (* tr11 ti3)))
 743                 (setf cr5 (+ (* ti11 tr5) (* ti12 tr4)))
 744                 (setf ci5 (+ (* ti11 ti5) (* ti12 ti4)))
 745                 (setf cr4 (- (* ti12 tr5) (* ti11 tr4)))
 746                 (setf ci4 (- (* ti12 ti5) (* ti11 ti4)))
 747                 (setf dr3 (- cr3 ci4))
 748                 (setf dr4 (+ cr3 ci4))
 749                 (setf di3 (+ ci3 cr4))
 750                 (setf di4 (- ci3 cr4))
 751                 (setf dr5 (+ cr2 ci5))
 752                 (setf dr2 (- cr2 ci5))
 753                 (setf di5 (- ci2 cr5))
 754                 (setf di2 (+ ci2 cr5))
 755                 (setf (f2cl-lib:fref ch-%data%
 756                                      (1 k 2 i)
 757                                      ((1 in2) (1 l1) (1 5) (1 ido))
 758                                      ch-%offset%)
 759                         (+
 760                          (*
 761                           (f2cl-lib:fref wa-%data%
 762                                          (i 1 1)
 763                                          ((1 ido) (1 4) (1 2))
 764                                          wa-%offset%)
 765                           dr2)
 766                          (*
 767                           (f2cl-lib:fref wa-%data%
 768                                          (i 1 2)
 769                                          ((1 ido) (1 4) (1 2))
 770                                          wa-%offset%)
 771                           di2)))
 772                 (setf (f2cl-lib:fref ch-%data%
 773                                      (2 k 2 i)
 774                                      ((1 in2) (1 l1) (1 5) (1 ido))
 775                                      ch-%offset%)
 776                         (-
 777                          (*
 778                           (f2cl-lib:fref wa-%data%
 779                                          (i 1 1)
 780                                          ((1 ido) (1 4) (1 2))
 781                                          wa-%offset%)
 782                           di2)
 783                          (*
 784                           (f2cl-lib:fref wa-%data%
 785                                          (i 1 2)
 786                                          ((1 ido) (1 4) (1 2))
 787                                          wa-%offset%)
 788                           dr2)))
 789                 (setf (f2cl-lib:fref ch-%data%
 790                                      (1 k 3 i)
 791                                      ((1 in2) (1 l1) (1 5) (1 ido))
 792                                      ch-%offset%)
 793                         (+
 794                          (*
 795                           (f2cl-lib:fref wa-%data%
 796                                          (i 2 1)
 797                                          ((1 ido) (1 4) (1 2))
 798                                          wa-%offset%)
 799                           dr3)
 800                          (*
 801                           (f2cl-lib:fref wa-%data%
 802                                          (i 2 2)
 803                                          ((1 ido) (1 4) (1 2))
 804                                          wa-%offset%)
 805                           di3)))
 806                 (setf (f2cl-lib:fref ch-%data%
 807                                      (2 k 3 i)
 808                                      ((1 in2) (1 l1) (1 5) (1 ido))
 809                                      ch-%offset%)
 810                         (-
 811                          (*
 812                           (f2cl-lib:fref wa-%data%
 813                                          (i 2 1)
 814                                          ((1 ido) (1 4) (1 2))
 815                                          wa-%offset%)
 816                           di3)
 817                          (*
 818                           (f2cl-lib:fref wa-%data%
 819                                          (i 2 2)
 820                                          ((1 ido) (1 4) (1 2))
 821                                          wa-%offset%)
 822                           dr3)))
 823                 (setf (f2cl-lib:fref ch-%data%
 824                                      (1 k 4 i)
 825                                      ((1 in2) (1 l1) (1 5) (1 ido))
 826                                      ch-%offset%)
 827                         (+
 828                          (*
 829                           (f2cl-lib:fref wa-%data%
 830                                          (i 3 1)
 831                                          ((1 ido) (1 4) (1 2))
 832                                          wa-%offset%)
 833                           dr4)
 834                          (*
 835                           (f2cl-lib:fref wa-%data%
 836                                          (i 3 2)
 837                                          ((1 ido) (1 4) (1 2))
 838                                          wa-%offset%)
 839                           di4)))
 840                 (setf (f2cl-lib:fref ch-%data%
 841                                      (2 k 4 i)
 842                                      ((1 in2) (1 l1) (1 5) (1 ido))
 843                                      ch-%offset%)
 844                         (-
 845                          (*
 846                           (f2cl-lib:fref wa-%data%
 847                                          (i 3 1)
 848                                          ((1 ido) (1 4) (1 2))
 849                                          wa-%offset%)
 850                           di4)
 851                          (*
 852                           (f2cl-lib:fref wa-%data%
 853                                          (i 3 2)
 854                                          ((1 ido) (1 4) (1 2))
 855                                          wa-%offset%)
 856                           dr4)))
 857                 (setf (f2cl-lib:fref ch-%data%
 858                                      (1 k 5 i)
 859                                      ((1 in2) (1 l1) (1 5) (1 ido))
 860                                      ch-%offset%)
 861                         (+
 862                          (*
 863                           (f2cl-lib:fref wa-%data%
 864                                          (i 4 1)
 865                                          ((1 ido) (1 4) (1 2))
 866                                          wa-%offset%)
 867                           dr5)
 868                          (*
 869                           (f2cl-lib:fref wa-%data%
 870                                          (i 4 2)
 871                                          ((1 ido) (1 4) (1 2))
 872                                          wa-%offset%)
 873                           di5)))
 874                 (setf (f2cl-lib:fref ch-%data%
 875                                      (2 k 5 i)
 876                                      ((1 in2) (1 l1) (1 5) (1 ido))
 877                                      ch-%offset%)
 878                         (-
 879                          (*
 880                           (f2cl-lib:fref wa-%data%
 881                                          (i 4 1)
 882                                          ((1 ido) (1 4) (1 2))
 883                                          wa-%offset%)
 884                           di5)
 885                          (*
 886                           (f2cl-lib:fref wa-%data%
 887                                          (i 4 2)
 888                                          ((1 ido) (1 4) (1 2))
 889                                          wa-%offset%)
 890                           dr5)))
 891                label104))
 892            label105))
 893         (go end_label)
 894        end_label
 895         (return (values nil nil nil nil nil nil nil nil))))))
 896
 897 (in-package #:cl-user)
 898 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
 899 (eval-when (:load-toplevel :compile-toplevel :execute)
 900   (setf (gethash 'fortran-to-lisp::c1f5kf
 901                  fortran-to-lisp::*f2cl-function-info*)
 902           (fortran-to-lisp::make-f2cl-finfo
 903            :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 904                         (fortran-to-lisp::integer4) (array double-float (*))
 905                         (fortran-to-lisp::integer4) (array double-float (*))
 906                         (fortran-to-lisp::integer4) (array double-float (*)))
 907            :return-values '(nil nil nil nil nil nil nil nil)
 908            :calls 'nil)))
 909