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