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