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