share/fftpack5/lisp/cmfgkb.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 cmfgkb (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) (wai 0.0d0) (war 0.0d0)
  30            (idlj 0) (lc 0) (l 0) (jc 0) (j 0) (m1 0) (m2 0) (ki 0) (ipph 0)
  31            (ipp2 0) (m2s 0) (m1d 0))
  32       (declare (type (double-float) war wai 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                            (f2cl-lib:fref wa-%data%
 215                                           (1 (f2cl-lib:int-sub l 1) 2)
 216                                           ((1 ido)
 217                                            (1
 218                                             (f2cl-lib:int-add ip
 219                                                               (f2cl-lib:int-sub
 220                                                                1)))
 221                                            (1 2))
 222                                           wa-%offset%)
 223                            (f2cl-lib:fref ch1-%data%
 224                                           (1 m2 ki ip)
 225                                           ((1 2) (1 in2) (1 lid) (1 ip))
 226                                           ch1-%offset%)))
 227                   (setf (f2cl-lib:fref cc1-%data%
 228                                        (2 m1 ki l)
 229                                        ((1 2) (1 in1) (1 lid) (1 ip))
 230                                        cc1-%offset%)
 231                           (+
 232                            (f2cl-lib:fref ch1-%data%
 233                                           (2 m2 ki 1)
 234                                           ((1 2) (1 in2) (1 lid) (1 ip))
 235                                           ch1-%offset%)
 236                            (*
 237                             (f2cl-lib:fref wa-%data%
 238                                            (1 (f2cl-lib:int-sub l 1) 1)
 239                                            ((1 ido)
 240                                             (1
 241                                              (f2cl-lib:int-add ip
 242                                                                (f2cl-lib:int-sub
 243                                                                 1)))
 244                                             (1 2))
 245                                            wa-%offset%)
 246                             (f2cl-lib:fref ch1-%data%
 247                                            (2 m2 ki 2)
 248                                            ((1 2) (1 in2) (1 lid) (1 ip))
 249                                            ch1-%offset%))))
 250                   (setf (f2cl-lib:fref cc1-%data%
 251                                        (2 m1 ki lc)
 252                                        ((1 2) (1 in1) (1 lid) (1 ip))
 253                                        cc1-%offset%)
 254                           (*
 255                            (f2cl-lib:fref wa-%data%
 256                                           (1 (f2cl-lib:int-sub l 1) 2)
 257                                           ((1 ido)
 258                                            (1
 259                                             (f2cl-lib:int-add ip
 260                                                               (f2cl-lib:int-sub
 261                                                                1)))
 262                                            (1 2))
 263                                           wa-%offset%)
 264                            (f2cl-lib:fref ch1-%data%
 265                                           (2 m2 ki ip)
 266                                           ((1 2) (1 in2) (1 lid) (1 ip))
 267                                           ch1-%offset%)))
 268                  label113))))
 269          label113
 270           (f2cl-lib:fdo (j 3 (f2cl-lib:int-add j 1))
 271                         ((> j ipph) nil)
 272             (tagbody
 273               (setf jc (f2cl-lib:int-sub ipp2 j))
 274               (setf idlj
 275                       (mod
 276                        (f2cl-lib:int-mul (f2cl-lib:int-sub l 1)
 277                                          (f2cl-lib:int-sub j 1))
 278                        ip))
 279               (setf war
 280                       (f2cl-lib:fref wa-%data%
 281                                      (1 idlj 1)
 282                                      ((1 ido)
 283                                       (1
 284                                        (f2cl-lib:int-add ip
 285                                                          (f2cl-lib:int-sub 1)))
 286                                       (1 2))
 287                                      wa-%offset%))
 288               (setf wai
 289                       (f2cl-lib:fref wa-%data%
 290                                      (1 idlj 2)
 291                                      ((1 ido)
 292                                       (1
 293                                        (f2cl-lib:int-add ip
 294                                                          (f2cl-lib:int-sub 1)))
 295                                       (1 2))
 296                                      wa-%offset%))
 297               (f2cl-lib:fdo (ki 1 (f2cl-lib:int-add ki 1))
 298                             ((> ki lid) nil)
 299                 (tagbody
 300                   (setf m2 m2s)
 301                   (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 302                                 ((> m1 m1d) nil)
 303                     (tagbody
 304                       (setf m2 (f2cl-lib:int-add m2 im2))
 305                       (setf (f2cl-lib:fref cc1-%data%
 306                                            (1 m1 ki l)
 307                                            ((1 2) (1 in1) (1 lid) (1 ip))
 308                                            cc1-%offset%)
 309                               (+
 310                                (f2cl-lib:fref cc1-%data%
 311                                               (1 m1 ki l)
 312                                               ((1 2) (1 in1) (1 lid) (1 ip))
 313                                               cc1-%offset%)
 314                                (* war
 315                                   (f2cl-lib:fref ch1-%data%
 316                                                  (1 m2 ki j)
 317                                                  ((1 2) (1 in2) (1 lid) (1 ip))
 318                                                  ch1-%offset%))))
 319                       (setf (f2cl-lib:fref cc1-%data%
 320                                            (1 m1 ki lc)
 321                                            ((1 2) (1 in1) (1 lid) (1 ip))
 322                                            cc1-%offset%)
 323                               (+
 324                                (f2cl-lib:fref cc1-%data%
 325                                               (1 m1 ki lc)
 326                                               ((1 2) (1 in1) (1 lid) (1 ip))
 327                                               cc1-%offset%)
 328                                (* wai
 329                                   (f2cl-lib:fref ch1-%data%
 330                                                  (1 m2 ki jc)
 331                                                  ((1 2) (1 in2) (1 lid) (1 ip))
 332                                                  ch1-%offset%))))
 333                       (setf (f2cl-lib:fref cc1-%data%
 334                                            (2 m1 ki l)
 335                                            ((1 2) (1 in1) (1 lid) (1 ip))
 336                                            cc1-%offset%)
 337                               (+
 338                                (f2cl-lib:fref cc1-%data%
 339                                               (2 m1 ki l)
 340                                               ((1 2) (1 in1) (1 lid) (1 ip))
 341                                               cc1-%offset%)
 342                                (* war
 343                                   (f2cl-lib:fref ch1-%data%
 344                                                  (2 m2 ki j)
 345                                                  ((1 2) (1 in2) (1 lid) (1 ip))
 346                                                  ch1-%offset%))))
 347                       (setf (f2cl-lib:fref cc1-%data%
 348                                            (2 m1 ki lc)
 349                                            ((1 2) (1 in1) (1 lid) (1 ip))
 350                                            cc1-%offset%)
 351                               (+
 352                                (f2cl-lib:fref cc1-%data%
 353                                               (2 m1 ki lc)
 354                                               ((1 2) (1 in1) (1 lid) (1 ip))
 355                                               cc1-%offset%)
 356                                (* wai
 357                                   (f2cl-lib:fref ch1-%data%
 358                                                  (2 m2 ki jc)
 359                                                  ((1 2) (1 in2) (1 lid) (1 ip))
 360                                                  ch1-%offset%))))
 361                      label114))))
 362              label114
 363              label115))
 364          label116))
 365       (if (or (> ido 1) (= na 1)) (go label136))
 366       (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
 367                     ((> j ipph) nil)
 368         (tagbody
 369           (setf jc (f2cl-lib:int-sub ipp2 j))
 370           (f2cl-lib:fdo (ki 1 (f2cl-lib:int-add ki 1))
 371                         ((> ki lid) nil)
 372             (tagbody
 373               (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 374                             ((> m1 m1d) nil)
 375                 (tagbody
 376                   (setf chold1
 377                           (-
 378                            (f2cl-lib:fref cc1-%data%
 379                                           (1 m1 ki j)
 380                                           ((1 2) (1 in1) (1 lid) (1 ip))
 381                                           cc1-%offset%)
 382                            (f2cl-lib:fref cc1-%data%
 383                                           (2 m1 ki jc)
 384                                           ((1 2) (1 in1) (1 lid) (1 ip))
 385                                           cc1-%offset%)))
 386                   (setf chold2
 387                           (+
 388                            (f2cl-lib:fref cc1-%data%
 389                                           (1 m1 ki j)
 390                                           ((1 2) (1 in1) (1 lid) (1 ip))
 391                                           cc1-%offset%)
 392                            (f2cl-lib:fref cc1-%data%
 393                                           (2 m1 ki jc)
 394                                           ((1 2) (1 in1) (1 lid) (1 ip))
 395                                           cc1-%offset%)))
 396                   (setf (f2cl-lib:fref cc1-%data%
 397                                        (1 m1 ki j)
 398                                        ((1 2) (1 in1) (1 lid) (1 ip))
 399                                        cc1-%offset%)
 400                           chold1)
 401                   (setf (f2cl-lib:fref cc1-%data%
 402                                        (2 m1 ki jc)
 403                                        ((1 2) (1 in1) (1 lid) (1 ip))
 404                                        cc1-%offset%)
 405                           (-
 406                            (f2cl-lib:fref cc1-%data%
 407                                           (2 m1 ki j)
 408                                           ((1 2) (1 in1) (1 lid) (1 ip))
 409                                           cc1-%offset%)
 410                            (f2cl-lib:fref cc1-%data%
 411                                           (1 m1 ki jc)
 412                                           ((1 2) (1 in1) (1 lid) (1 ip))
 413                                           cc1-%offset%)))
 414                   (setf (f2cl-lib:fref cc1-%data%
 415                                        (2 m1 ki j)
 416                                        ((1 2) (1 in1) (1 lid) (1 ip))
 417                                        cc1-%offset%)
 418                           (+
 419                            (f2cl-lib:fref cc1-%data%
 420                                           (2 m1 ki j)
 421                                           ((1 2) (1 in1) (1 lid) (1 ip))
 422                                           cc1-%offset%)
 423                            (f2cl-lib:fref cc1-%data%
 424                                           (1 m1 ki jc)
 425                                           ((1 2) (1 in1) (1 lid) (1 ip))
 426                                           cc1-%offset%)))
 427                   (setf (f2cl-lib:fref cc1-%data%
 428                                        (1 m1 ki jc)
 429                                        ((1 2) (1 in1) (1 lid) (1 ip))
 430                                        cc1-%offset%)
 431                           chold2)
 432                  label119))))
 433          label119
 434          label120))
 435       (go end_label)
 436      label136
 437       (f2cl-lib:fdo (ki 1 (f2cl-lib:int-add ki 1))
 438                     ((> ki lid) nil)
 439         (tagbody
 440           (setf m2 m2s)
 441           (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 442                         ((> m1 m1d) nil)
 443             (tagbody
 444               (setf m2 (f2cl-lib:int-add m2 im2))
 445               (setf (f2cl-lib:fref ch1-%data%
 446                                    (1 m2 ki 1)
 447                                    ((1 2) (1 in2) (1 lid) (1 ip))
 448                                    ch1-%offset%)
 449                       (f2cl-lib:fref cc1-%data%
 450                                      (1 m1 ki 1)
 451                                      ((1 2) (1 in1) (1 lid) (1 ip))
 452                                      cc1-%offset%))
 453               (setf (f2cl-lib:fref ch1-%data%
 454                                    (2 m2 ki 1)
 455                                    ((1 2) (1 in2) (1 lid) (1 ip))
 456                                    ch1-%offset%)
 457                       (f2cl-lib:fref cc1-%data%
 458                                      (2 m1 ki 1)
 459                                      ((1 2) (1 in1) (1 lid) (1 ip))
 460                                      cc1-%offset%))
 461              label137))))
 462      label137
 463       (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
 464                     ((> j ipph) nil)
 465         (tagbody
 466           (setf jc (f2cl-lib:int-sub ipp2 j))
 467           (f2cl-lib:fdo (ki 1 (f2cl-lib:int-add ki 1))
 468                         ((> ki lid) nil)
 469             (tagbody
 470               (setf m2 m2s)
 471               (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 472                             ((> m1 m1d) nil)
 473                 (tagbody
 474                   (setf m2 (f2cl-lib:int-add m2 im2))
 475                   (setf (f2cl-lib:fref ch1-%data%
 476                                        (1 m2 ki j)
 477                                        ((1 2) (1 in2) (1 lid) (1 ip))
 478                                        ch1-%offset%)
 479                           (-
 480                            (f2cl-lib:fref cc1-%data%
 481                                           (1 m1 ki j)
 482                                           ((1 2) (1 in1) (1 lid) (1 ip))
 483                                           cc1-%offset%)
 484                            (f2cl-lib:fref cc1-%data%
 485                                           (2 m1 ki jc)
 486                                           ((1 2) (1 in1) (1 lid) (1 ip))
 487                                           cc1-%offset%)))
 488                   (setf (f2cl-lib:fref ch1-%data%
 489                                        (1 m2 ki jc)
 490                                        ((1 2) (1 in2) (1 lid) (1 ip))
 491                                        ch1-%offset%)
 492                           (+
 493                            (f2cl-lib:fref cc1-%data%
 494                                           (1 m1 ki j)
 495                                           ((1 2) (1 in1) (1 lid) (1 ip))
 496                                           cc1-%offset%)
 497                            (f2cl-lib:fref cc1-%data%
 498                                           (2 m1 ki jc)
 499                                           ((1 2) (1 in1) (1 lid) (1 ip))
 500                                           cc1-%offset%)))
 501                   (setf (f2cl-lib:fref ch1-%data%
 502                                        (2 m2 ki jc)
 503                                        ((1 2) (1 in2) (1 lid) (1 ip))
 504                                        ch1-%offset%)
 505                           (-
 506                            (f2cl-lib:fref cc1-%data%
 507                                           (2 m1 ki j)
 508                                           ((1 2) (1 in1) (1 lid) (1 ip))
 509                                           cc1-%offset%)
 510                            (f2cl-lib:fref cc1-%data%
 511                                           (1 m1 ki jc)
 512                                           ((1 2) (1 in1) (1 lid) (1 ip))
 513                                           cc1-%offset%)))
 514                   (setf (f2cl-lib:fref ch1-%data%
 515                                        (2 m2 ki j)
 516                                        ((1 2) (1 in2) (1 lid) (1 ip))
 517                                        ch1-%offset%)
 518                           (+
 519                            (f2cl-lib:fref cc1-%data%
 520                                           (2 m1 ki j)
 521                                           ((1 2) (1 in1) (1 lid) (1 ip))
 522                                           cc1-%offset%)
 523                            (f2cl-lib:fref cc1-%data%
 524                                           (1 m1 ki jc)
 525                                           ((1 2) (1 in1) (1 lid) (1 ip))
 526                                           cc1-%offset%)))
 527                  label134))))
 528          label134
 529          label135))
 530       (if (= ido 1) (go end_label))
 531       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
 532                     ((> i ido) nil)
 533         (tagbody
 534           (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 535                         ((> k l1) nil)
 536             (tagbody
 537               (setf m2 m2s)
 538               (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 539                             ((> m1 m1d) nil)
 540                 (tagbody
 541                   (setf m2 (f2cl-lib:int-add m2 im2))
 542                   (setf (f2cl-lib:fref cc-%data%
 543                                        (1 m1 k 1 i)
 544                                        ((1 2) (1 in1) (1 l1) (1 ip) (1 ido))
 545                                        cc-%offset%)
 546                           (f2cl-lib:fref ch-%data%
 547                                          (1 m2 k i 1)
 548                                          ((1 2) (1 in2) (1 l1) (1 ido) (1 ip))
 549                                          ch-%offset%))
 550                   (setf (f2cl-lib:fref cc-%data%
 551                                        (2 m1 k 1 i)
 552                                        ((1 2) (1 in1) (1 l1) (1 ip) (1 ido))
 553                                        cc-%offset%)
 554                           (f2cl-lib:fref ch-%data%
 555                                          (2 m2 k i 1)
 556                                          ((1 2) (1 in2) (1 l1) (1 ido) (1 ip))
 557                                          ch-%offset%))
 558                  label130))))
 559          label130
 560          label131))
 561       (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
 562                     ((> j ip) nil)
 563         (tagbody
 564           (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 565                         ((> k l1) nil)
 566             (tagbody
 567               (setf m2 m2s)
 568               (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 569                             ((> m1 m1d) nil)
 570                 (tagbody
 571                   (setf m2 (f2cl-lib:int-add m2 im2))
 572                   (setf (f2cl-lib:fref cc-%data%
 573                                        (1 m1 k j 1)
 574                                        ((1 2) (1 in1) (1 l1) (1 ip) (1 ido))
 575                                        cc-%offset%)
 576                           (f2cl-lib:fref ch-%data%
 577                                          (1 m2 k 1 j)
 578                                          ((1 2) (1 in2) (1 l1) (1 ido) (1 ip))
 579                                          ch-%offset%))
 580                   (setf (f2cl-lib:fref cc-%data%
 581                                        (2 m1 k j 1)
 582                                        ((1 2) (1 in1) (1 l1) (1 ip) (1 ido))
 583                                        cc-%offset%)
 584                           (f2cl-lib:fref ch-%data%
 585                                          (2 m2 k 1 j)
 586                                          ((1 2) (1 in2) (1 l1) (1 ido) (1 ip))
 587                                          ch-%offset%))
 588                  label122))))
 589          label122
 590          label123))
 591       (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
 592                     ((> j ip) nil)
 593         (tagbody
 594           (f2cl-lib:fdo (i 2 (f2cl-lib:int-add i 1))
 595                         ((> i ido) nil)
 596             (tagbody
 597               (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 598                             ((> k l1) nil)
 599                 (tagbody
 600                   (setf m2 m2s)
 601                   (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 602                                 ((> m1 m1d) nil)
 603                     (tagbody
 604                       (setf m2 (f2cl-lib:int-add m2 im2))
 605                       (setf (f2cl-lib:fref cc-%data%
 606                                            (1 m1 k j i)
 607                                            ((1 2) (1 in1) (1 l1) (1 ip)
 608                                             (1 ido))
 609                                            cc-%offset%)
 610                               (-
 611                                (*
 612                                 (f2cl-lib:fref wa-%data%
 613                                                (i (f2cl-lib:int-sub j 1) 1)
 614                                                ((1 ido)
 615                                                 (1
 616                                                  (f2cl-lib:int-add ip
 617                                                                    (f2cl-lib:int-sub
 618                                                                     1)))
 619                                                 (1 2))
 620                                                wa-%offset%)
 621                                 (f2cl-lib:fref ch-%data%
 622                                                (1 m2 k i j)
 623                                                ((1 2) (1 in2) (1 l1) (1 ido)
 624                                                 (1 ip))
 625                                                ch-%offset%))
 626                                (*
 627                                 (f2cl-lib:fref wa-%data%
 628                                                (i (f2cl-lib:int-sub j 1) 2)
 629                                                ((1 ido)
 630                                                 (1
 631                                                  (f2cl-lib:int-add ip
 632                                                                    (f2cl-lib:int-sub
 633                                                                     1)))
 634                                                 (1 2))
 635                                                wa-%offset%)
 636                                 (f2cl-lib:fref ch-%data%
 637                                                (2 m2 k i j)
 638                                                ((1 2) (1 in2) (1 l1) (1 ido)
 639                                                 (1 ip))
 640                                                ch-%offset%))))
 641                       (setf (f2cl-lib:fref cc-%data%
 642                                            (2 m1 k j i)
 643                                            ((1 2) (1 in1) (1 l1) (1 ip)
 644                                             (1 ido))
 645                                            cc-%offset%)
 646                               (+
 647                                (*
 648                                 (f2cl-lib:fref wa-%data%
 649                                                (i (f2cl-lib:int-sub j 1) 1)
 650                                                ((1 ido)
 651                                                 (1
 652                                                  (f2cl-lib:int-add ip
 653                                                                    (f2cl-lib:int-sub
 654                                                                     1)))
 655                                                 (1 2))
 656                                                wa-%offset%)
 657                                 (f2cl-lib:fref ch-%data%
 658                                                (2 m2 k i j)
 659                                                ((1 2) (1 in2) (1 l1) (1 ido)
 660                                                 (1 ip))
 661                                                ch-%offset%))
 662                                (*
 663                                 (f2cl-lib:fref wa-%data%
 664                                                (i (f2cl-lib:int-sub j 1) 2)
 665                                                ((1 ido)
 666                                                 (1
 667                                                  (f2cl-lib:int-add ip
 668                                                                    (f2cl-lib:int-sub
 669                                                                     1)))
 670                                                 (1 2))
 671                                                wa-%offset%)
 672                                 (f2cl-lib:fref ch-%data%
 673                                                (1 m2 k i j)
 674                                                ((1 2) (1 in2) (1 l1) (1 ido)
 675                                                 (1 ip))
 676                                                ch-%offset%))))
 677                      label124))))
 678              label124
 679              label125))
 680          label126))
 681       (go end_label)
 682      end_label
 683       (return
 684        (values nil nil nil nil nil nil nil nil nil nil nil nil nil nil nil)))))
 685
 686 (in-package #:cl-user)
 687 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
 688 (eval-when (:load-toplevel :compile-toplevel :execute)
 689   (setf (gethash 'fortran-to-lisp::cmfgkb
 690                  fortran-to-lisp::*f2cl-function-info*)
 691           (fortran-to-lisp::make-f2cl-finfo
 692            :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 693                         (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 694                         (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 695                         (array double-float (*)) (array double-float (*))
 696                         (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 697                         (array double-float (*)) (array double-float (*))
 698                         (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
 699                         (array double-float (*)))
 700            :return-values '(nil nil nil nil nil nil nil nil nil nil nil nil nil
 701                             nil nil)
 702            :calls 'nil)))
 703