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