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