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