share/fftpack5/lisp/cmf4kf.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 (defun cmf4kf (lot ido l1 na cc im1 in1 ch im2 in2 wa)
  21   (declare (type (array double-float (*)) wa ch cc)
  22            (type (f2cl-lib:integer4) in2 im2 in1 im1 na l1 ido lot))
  23   (f2cl-lib:with-multi-array-data
  24       ((cc double-float cc-%data% cc-%offset%)
  25        (ch double-float ch-%data% ch-%offset%)
  26        (wa double-float wa-%data% wa-%offset%))
  27     (prog ((ci4 0.0d0) (ci2 0.0d0) (cr4 0.0d0) (cr2 0.0d0) (ci3 0.0d0)
  28            (cr3 0.0d0) (i 0) (m2 0) (tr3 0.0d0) (ti4 0.0d0) (tr2 0.0d0)
  29            (tr1 0.0d0) (ti3 0.0d0) (tr4 0.0d0) (ti2 0.0d0) (ti1 0.0d0) (m1 0)
  30            (k 0) (sn 0.0d0) (m2s 0) (m1d 0))
  31       (declare (type (f2cl-lib:integer4) m1d m2s k m1 m2 i)
  32                (type (double-float) sn ti1 ti2 tr4 ti3 tr1 tr2 ti4 tr3 cr3 ci3
  33                                     cr2 cr4 ci2 ci4))
  34       (setf m1d
  35               (f2cl-lib:int-add (f2cl-lib:int-mul (f2cl-lib:int-sub lot 1) im1)
  36                                 1))
  37       (setf m2s (f2cl-lib:int-sub 1 im2))
  38       (if (> ido 1) (go label102))
  39       (setf sn (/ 1.0d0 (f2cl-lib:freal (f2cl-lib:int-mul 4 l1))))
  40       (if (= na 1) (go label106))
  41       (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  42                     ((> k l1) nil)
  43         (tagbody
  44           (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
  45                         ((> m1 m1d) nil)
  46             (tagbody
  47               (setf ti1
  48                       (-
  49                        (f2cl-lib:fref cc-%data%
  50                                       (2 m1 k 1 1)
  51                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  52                                       cc-%offset%)
  53                        (f2cl-lib:fref cc-%data%
  54                                       (2 m1 k 1 3)
  55                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  56                                       cc-%offset%)))
  57               (setf ti2
  58                       (+
  59                        (f2cl-lib:fref cc-%data%
  60                                       (2 m1 k 1 1)
  61                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  62                                       cc-%offset%)
  63                        (f2cl-lib:fref cc-%data%
  64                                       (2 m1 k 1 3)
  65                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  66                                       cc-%offset%)))
  67               (setf tr4
  68                       (-
  69                        (f2cl-lib:fref cc-%data%
  70                                       (2 m1 k 1 2)
  71                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  72                                       cc-%offset%)
  73                        (f2cl-lib:fref cc-%data%
  74                                       (2 m1 k 1 4)
  75                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  76                                       cc-%offset%)))
  77               (setf ti3
  78                       (+
  79                        (f2cl-lib:fref cc-%data%
  80                                       (2 m1 k 1 2)
  81                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  82                                       cc-%offset%)
  83                        (f2cl-lib:fref cc-%data%
  84                                       (2 m1 k 1 4)
  85                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  86                                       cc-%offset%)))
  87               (setf tr1
  88                       (-
  89                        (f2cl-lib:fref cc-%data%
  90                                       (1 m1 k 1 1)
  91                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  92                                       cc-%offset%)
  93                        (f2cl-lib:fref cc-%data%
  94                                       (1 m1 k 1 3)
  95                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
  96                                       cc-%offset%)))
  97               (setf tr2
  98                       (+
  99                        (f2cl-lib:fref cc-%data%
 100                                       (1 m1 k 1 1)
 101                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 102                                       cc-%offset%)
 103                        (f2cl-lib:fref cc-%data%
 104                                       (1 m1 k 1 3)
 105                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 106                                       cc-%offset%)))
 107               (setf ti4
 108                       (-
 109                        (f2cl-lib:fref cc-%data%
 110                                       (1 m1 k 1 4)
 111                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 112                                       cc-%offset%)
 113                        (f2cl-lib:fref cc-%data%
 114                                       (1 m1 k 1 2)
 115                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 116                                       cc-%offset%)))
 117               (setf tr3
 118                       (+
 119                        (f2cl-lib:fref cc-%data%
 120                                       (1 m1 k 1 2)
 121                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 122                                       cc-%offset%)
 123                        (f2cl-lib:fref cc-%data%
 124                                       (1 m1 k 1 4)
 125                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 126                                       cc-%offset%)))
 127               (setf (f2cl-lib:fref cc-%data%
 128                                    (1 m1 k 1 1)
 129                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 130                                    cc-%offset%)
 131                       (* sn (+ tr2 tr3)))
 132               (setf (f2cl-lib:fref cc-%data%
 133                                    (1 m1 k 1 3)
 134                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 135                                    cc-%offset%)
 136                       (* sn (- tr2 tr3)))
 137               (setf (f2cl-lib:fref cc-%data%
 138                                    (2 m1 k 1 1)
 139                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 140                                    cc-%offset%)
 141                       (* sn (+ ti2 ti3)))
 142               (setf (f2cl-lib:fref cc-%data%
 143                                    (2 m1 k 1 3)
 144                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 145                                    cc-%offset%)
 146                       (* sn (- ti2 ti3)))
 147               (setf (f2cl-lib:fref cc-%data%
 148                                    (1 m1 k 1 2)
 149                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 150                                    cc-%offset%)
 151                       (* sn (+ tr1 tr4)))
 152               (setf (f2cl-lib:fref cc-%data%
 153                                    (1 m1 k 1 4)
 154                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 155                                    cc-%offset%)
 156                       (* sn (- tr1 tr4)))
 157               (setf (f2cl-lib:fref cc-%data%
 158                                    (2 m1 k 1 2)
 159                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 160                                    cc-%offset%)
 161                       (* sn (+ ti1 ti4)))
 162               (setf (f2cl-lib:fref cc-%data%
 163                                    (2 m1 k 1 4)
 164                                    ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 165                                    cc-%offset%)
 166                       (* sn (- ti1 ti4)))
 167              label101))))
 168      label101
 169       (go end_label)
 170      label106
 171       (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 172                     ((> k l1) nil)
 173         (tagbody
 174           (setf m2 m2s)
 175           (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 176                         ((> m1 m1d) nil)
 177             (tagbody
 178               (setf m2 (f2cl-lib:int-add m2 im2))
 179               (setf ti1
 180                       (-
 181                        (f2cl-lib:fref cc-%data%
 182                                       (2 m1 k 1 1)
 183                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 184                                       cc-%offset%)
 185                        (f2cl-lib:fref cc-%data%
 186                                       (2 m1 k 1 3)
 187                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 188                                       cc-%offset%)))
 189               (setf ti2
 190                       (+
 191                        (f2cl-lib:fref cc-%data%
 192                                       (2 m1 k 1 1)
 193                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 194                                       cc-%offset%)
 195                        (f2cl-lib:fref cc-%data%
 196                                       (2 m1 k 1 3)
 197                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 198                                       cc-%offset%)))
 199               (setf tr4
 200                       (-
 201                        (f2cl-lib:fref cc-%data%
 202                                       (2 m1 k 1 2)
 203                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 204                                       cc-%offset%)
 205                        (f2cl-lib:fref cc-%data%
 206                                       (2 m1 k 1 4)
 207                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 208                                       cc-%offset%)))
 209               (setf ti3
 210                       (+
 211                        (f2cl-lib:fref cc-%data%
 212                                       (2 m1 k 1 2)
 213                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 214                                       cc-%offset%)
 215                        (f2cl-lib:fref cc-%data%
 216                                       (2 m1 k 1 4)
 217                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 218                                       cc-%offset%)))
 219               (setf tr1
 220                       (-
 221                        (f2cl-lib:fref cc-%data%
 222                                       (1 m1 k 1 1)
 223                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 224                                       cc-%offset%)
 225                        (f2cl-lib:fref cc-%data%
 226                                       (1 m1 k 1 3)
 227                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 228                                       cc-%offset%)))
 229               (setf tr2
 230                       (+
 231                        (f2cl-lib:fref cc-%data%
 232                                       (1 m1 k 1 1)
 233                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 234                                       cc-%offset%)
 235                        (f2cl-lib:fref cc-%data%
 236                                       (1 m1 k 1 3)
 237                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 238                                       cc-%offset%)))
 239               (setf ti4
 240                       (-
 241                        (f2cl-lib:fref cc-%data%
 242                                       (1 m1 k 1 4)
 243                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 244                                       cc-%offset%)
 245                        (f2cl-lib:fref cc-%data%
 246                                       (1 m1 k 1 2)
 247                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 248                                       cc-%offset%)))
 249               (setf tr3
 250                       (+
 251                        (f2cl-lib:fref cc-%data%
 252                                       (1 m1 k 1 2)
 253                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 254                                       cc-%offset%)
 255                        (f2cl-lib:fref cc-%data%
 256                                       (1 m1 k 1 4)
 257                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 258                                       cc-%offset%)))
 259               (setf (f2cl-lib:fref ch-%data%
 260                                    (1 m2 k 1 1)
 261                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 262                                    ch-%offset%)
 263                       (* sn (+ tr2 tr3)))
 264               (setf (f2cl-lib:fref ch-%data%
 265                                    (1 m2 k 3 1)
 266                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 267                                    ch-%offset%)
 268                       (* sn (- tr2 tr3)))
 269               (setf (f2cl-lib:fref ch-%data%
 270                                    (2 m2 k 1 1)
 271                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 272                                    ch-%offset%)
 273                       (* sn (+ ti2 ti3)))
 274               (setf (f2cl-lib:fref ch-%data%
 275                                    (2 m2 k 3 1)
 276                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 277                                    ch-%offset%)
 278                       (* sn (- ti2 ti3)))
 279               (setf (f2cl-lib:fref ch-%data%
 280                                    (1 m2 k 2 1)
 281                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 282                                    ch-%offset%)
 283                       (* sn (+ tr1 tr4)))
 284               (setf (f2cl-lib:fref ch-%data%
 285                                    (1 m2 k 4 1)
 286                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 287                                    ch-%offset%)
 288                       (* sn (- tr1 tr4)))
 289               (setf (f2cl-lib:fref ch-%data%
 290                                    (2 m2 k 2 1)
 291                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 292                                    ch-%offset%)
 293                       (* sn (+ ti1 ti4)))
 294               (setf (f2cl-lib:fref ch-%data%
 295                                    (2 m2 k 4 1)
 296                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 297                                    ch-%offset%)
 298                       (* sn (- ti1 ti4)))
 299              label107))))
 300      label107
 301       (go end_label)
 302      label102
 303       (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 304                     ((> k l1) nil)
 305         (tagbody
 306           (setf m2 m2s)
 307           (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 308                         ((> m1 m1d) nil)
 309             (tagbody
 310               (setf m2 (f2cl-lib:int-add m2 im2))
 311               (setf ti1
 312                       (-
 313                        (f2cl-lib:fref cc-%data%
 314                                       (2 m1 k 1 1)
 315                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 316                                       cc-%offset%)
 317                        (f2cl-lib:fref cc-%data%
 318                                       (2 m1 k 1 3)
 319                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 320                                       cc-%offset%)))
 321               (setf ti2
 322                       (+
 323                        (f2cl-lib:fref cc-%data%
 324                                       (2 m1 k 1 1)
 325                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 326                                       cc-%offset%)
 327                        (f2cl-lib:fref cc-%data%
 328                                       (2 m1 k 1 3)
 329                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 330                                       cc-%offset%)))
 331               (setf tr4
 332                       (-
 333                        (f2cl-lib:fref cc-%data%
 334                                       (2 m1 k 1 2)
 335                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 336                                       cc-%offset%)
 337                        (f2cl-lib:fref cc-%data%
 338                                       (2 m1 k 1 4)
 339                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 340                                       cc-%offset%)))
 341               (setf ti3
 342                       (+
 343                        (f2cl-lib:fref cc-%data%
 344                                       (2 m1 k 1 2)
 345                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 346                                       cc-%offset%)
 347                        (f2cl-lib:fref cc-%data%
 348                                       (2 m1 k 1 4)
 349                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 350                                       cc-%offset%)))
 351               (setf tr1
 352                       (-
 353                        (f2cl-lib:fref cc-%data%
 354                                       (1 m1 k 1 1)
 355                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 356                                       cc-%offset%)
 357                        (f2cl-lib:fref cc-%data%
 358                                       (1 m1 k 1 3)
 359                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 360                                       cc-%offset%)))
 361               (setf tr2
 362                       (+
 363                        (f2cl-lib:fref cc-%data%
 364                                       (1 m1 k 1 1)
 365                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 366                                       cc-%offset%)
 367                        (f2cl-lib:fref cc-%data%
 368                                       (1 m1 k 1 3)
 369                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 370                                       cc-%offset%)))
 371               (setf ti4
 372                       (-
 373                        (f2cl-lib:fref cc-%data%
 374                                       (1 m1 k 1 4)
 375                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 376                                       cc-%offset%)
 377                        (f2cl-lib:fref cc-%data%
 378                                       (1 m1 k 1 2)
 379                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 380                                       cc-%offset%)))
 381               (setf tr3
 382                       (+
 383                        (f2cl-lib:fref cc-%data%
 384                                       (1 m1 k 1 2)
 385                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 386                                       cc-%offset%)
 387                        (f2cl-lib:fref cc-%data%
 388                                       (1 m1 k 1 4)
 389                                       ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 390                                       cc-%offset%)))
 391               (setf (f2cl-lib:fref ch-%data%
 392                                    (1 m2 k 1 1)
 393                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 394                                    ch-%offset%)
 395                       (+ tr2 tr3))
 396               (setf (f2cl-lib:fref ch-%data%
 397                                    (1 m2 k 3 1)
 398                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 399                                    ch-%offset%)
 400                       (- tr2 tr3))
 401               (setf (f2cl-lib:fref ch-%data%
 402                                    (2 m2 k 1 1)
 403                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 404                                    ch-%offset%)
 405                       (+ ti2 ti3))
 406               (setf (f2cl-lib:fref ch-%data%
 407                                    (2 m2 k 3 1)
 408                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 409                                    ch-%offset%)
 410                       (- ti2 ti3))
 411               (setf (f2cl-lib:fref ch-%data%
 412                                    (1 m2 k 2 1)
 413                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 414                                    ch-%offset%)
 415                       (+ tr1 tr4))
 416               (setf (f2cl-lib:fref ch-%data%
 417                                    (1 m2 k 4 1)
 418                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 419                                    ch-%offset%)
 420                       (- tr1 tr4))
 421               (setf (f2cl-lib:fref ch-%data%
 422                                    (2 m2 k 2 1)
 423                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 424                                    ch-%offset%)
 425                       (+ ti1 ti4))
 426               (setf (f2cl-lib:fref ch-%data%
 427                                    (2 m2 k 4 1)
 428                                    ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 429                                    ch-%offset%)
 430                       (- ti1 ti4))
 431              label103))))
 432      label103
 433       (f2cl-lib:fdo (i 2 (f2cl-lib:int-add i 1))
 434                     ((> i ido) nil)
 435         (tagbody
 436           (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 437                         ((> k l1) nil)
 438             (tagbody
 439               (setf m2 m2s)
 440               (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 441                             ((> m1 m1d) nil)
 442                 (tagbody
 443                   (setf m2 (f2cl-lib:int-add m2 im2))
 444                   (setf ti1
 445                           (-
 446                            (f2cl-lib:fref cc-%data%
 447                                           (2 m1 k i 1)
 448                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 449                                           cc-%offset%)
 450                            (f2cl-lib:fref cc-%data%
 451                                           (2 m1 k i 3)
 452                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 453                                           cc-%offset%)))
 454                   (setf ti2
 455                           (+
 456                            (f2cl-lib:fref cc-%data%
 457                                           (2 m1 k i 1)
 458                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 459                                           cc-%offset%)
 460                            (f2cl-lib:fref cc-%data%
 461                                           (2 m1 k i 3)
 462                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 463                                           cc-%offset%)))
 464                   (setf ti3
 465                           (+
 466                            (f2cl-lib:fref cc-%data%
 467                                           (2 m1 k i 2)
 468                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 469                                           cc-%offset%)
 470                            (f2cl-lib:fref cc-%data%
 471                                           (2 m1 k i 4)
 472                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 473                                           cc-%offset%)))
 474                   (setf tr4
 475                           (-
 476                            (f2cl-lib:fref cc-%data%
 477                                           (2 m1 k i 2)
 478                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 479                                           cc-%offset%)
 480                            (f2cl-lib:fref cc-%data%
 481                                           (2 m1 k i 4)
 482                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 483                                           cc-%offset%)))
 484                   (setf tr1
 485                           (-
 486                            (f2cl-lib:fref cc-%data%
 487                                           (1 m1 k i 1)
 488                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 489                                           cc-%offset%)
 490                            (f2cl-lib:fref cc-%data%
 491                                           (1 m1 k i 3)
 492                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 493                                           cc-%offset%)))
 494                   (setf tr2
 495                           (+
 496                            (f2cl-lib:fref cc-%data%
 497                                           (1 m1 k i 1)
 498                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 499                                           cc-%offset%)
 500                            (f2cl-lib:fref cc-%data%
 501                                           (1 m1 k i 3)
 502                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 503                                           cc-%offset%)))
 504                   (setf ti4
 505                           (-
 506                            (f2cl-lib:fref cc-%data%
 507                                           (1 m1 k i 4)
 508                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 509                                           cc-%offset%)
 510                            (f2cl-lib:fref cc-%data%
 511                                           (1 m1 k i 2)
 512                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 513                                           cc-%offset%)))
 514                   (setf tr3
 515                           (+
 516                            (f2cl-lib:fref cc-%data%
 517                                           (1 m1 k i 2)
 518                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 519                                           cc-%offset%)
 520                            (f2cl-lib:fref cc-%data%
 521                                           (1 m1 k i 4)
 522                                           ((1 2) (1 in1) (1 l1) (1 ido) (1 4))
 523                                           cc-%offset%)))
 524                   (setf (f2cl-lib:fref ch-%data%
 525                                        (1 m2 k 1 i)
 526                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 527                                        ch-%offset%)
 528                           (+ tr2 tr3))
 529                   (setf cr3 (- tr2 tr3))
 530                   (setf (f2cl-lib:fref ch-%data%
 531                                        (2 m2 k 1 i)
 532                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 533                                        ch-%offset%)
 534                           (+ ti2 ti3))
 535                   (setf ci3 (- ti2 ti3))
 536                   (setf cr2 (+ tr1 tr4))
 537                   (setf cr4 (- tr1 tr4))
 538                   (setf ci2 (+ ti1 ti4))
 539                   (setf ci4 (- ti1 ti4))
 540                   (setf (f2cl-lib:fref ch-%data%
 541                                        (1 m2 k 2 i)
 542                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 543                                        ch-%offset%)
 544                           (+
 545                            (*
 546                             (f2cl-lib:fref wa-%data%
 547                                            (i 1 1)
 548                                            ((1 ido) (1 3) (1 2))
 549                                            wa-%offset%)
 550                             cr2)
 551                            (*
 552                             (f2cl-lib:fref wa-%data%
 553                                            (i 1 2)
 554                                            ((1 ido) (1 3) (1 2))
 555                                            wa-%offset%)
 556                             ci2)))
 557                   (setf (f2cl-lib:fref ch-%data%
 558                                        (2 m2 k 2 i)
 559                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 560                                        ch-%offset%)
 561                           (-
 562                            (*
 563                             (f2cl-lib:fref wa-%data%
 564                                            (i 1 1)
 565                                            ((1 ido) (1 3) (1 2))
 566                                            wa-%offset%)
 567                             ci2)
 568                            (*
 569                             (f2cl-lib:fref wa-%data%
 570                                            (i 1 2)
 571                                            ((1 ido) (1 3) (1 2))
 572                                            wa-%offset%)
 573                             cr2)))
 574                   (setf (f2cl-lib:fref ch-%data%
 575                                        (1 m2 k 3 i)
 576                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 577                                        ch-%offset%)
 578                           (+
 579                            (*
 580                             (f2cl-lib:fref wa-%data%
 581                                            (i 2 1)
 582                                            ((1 ido) (1 3) (1 2))
 583                                            wa-%offset%)
 584                             cr3)
 585                            (*
 586                             (f2cl-lib:fref wa-%data%
 587                                            (i 2 2)
 588                                            ((1 ido) (1 3) (1 2))
 589                                            wa-%offset%)
 590                             ci3)))
 591                   (setf (f2cl-lib:fref ch-%data%
 592                                        (2 m2 k 3 i)
 593                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 594                                        ch-%offset%)
 595                           (-
 596                            (*
 597                             (f2cl-lib:fref wa-%data%
 598                                            (i 2 1)
 599                                            ((1 ido) (1 3) (1 2))
 600                                            wa-%offset%)
 601                             ci3)
 602                            (*
 603                             (f2cl-lib:fref wa-%data%
 604                                            (i 2 2)
 605                                            ((1 ido) (1 3) (1 2))
 606                                            wa-%offset%)
 607                             cr3)))
 608                   (setf (f2cl-lib:fref ch-%data%
 609                                        (1 m2 k 4 i)
 610                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 611                                        ch-%offset%)
 612                           (+
 613                            (*
 614                             (f2cl-lib:fref wa-%data%
 615                                            (i 3 1)
 616                                            ((1 ido) (1 3) (1 2))
 617                                            wa-%offset%)
 618                             cr4)
 619                            (*
 620                             (f2cl-lib:fref wa-%data%
 621                                            (i 3 2)
 622                                            ((1 ido) (1 3) (1 2))
 623                                            wa-%offset%)
 624                             ci4)))
 625                   (setf (f2cl-lib:fref ch-%data%
 626                                        (2 m2 k 4 i)
 627                                        ((1 2) (1 in2) (1 l1) (1 4) (1 ido))
 628                                        ch-%offset%)
 629                           (-
 630                            (*
 631                             (f2cl-lib:fref wa-%data%
 632                                            (i 3 1)
 633                                            ((1 ido) (1 3) (1 2))
 634                                            wa-%offset%)
 635                             ci4)
 636                            (*
 637                             (f2cl-lib:fref wa-%data%
 638                                            (i 3 2)
 639                                            ((1 ido) (1 3) (1 2))
 640                                            wa-%offset%)
 641                             cr4)))
 642                  label104))))
 643          label104
 644          label105))
 645       (go end_label)
 646      end_label
 647       (return (values nil nil nil nil nil nil nil nil nil nil nil)))))
 648
 649 (in-package #:cl-user)
 650 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
 651 (eval-when (:load-toplevel :compile-toplevel :execute)
 652   (setf (gethash 'fortran-to-lisp::cmf4kf
 653                  fortran-to-lisp::*f2cl-function-info*)
 654           (fortran-to-lisp::make-f2cl-finfo
 655            :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 656                         (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 657                         (array double-float (*)) (fortran-to-lisp::integer4)
 658                         (fortran-to-lisp::integer4) (array double-float (*))
 659                         (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 660                         (array double-float (*)))
 661            :return-values '(nil nil nil nil nil nil nil nil nil nil nil)
 662            :calls 'nil)))
 663