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