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