share/fftpack5/lisp/mradb5.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 mradb5 (m ido l1 cc im1 in1 ch im2 in2 wa1 wa2 wa3 wa4)
  21   (declare (type (array double-float (*)) wa4 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        (wa4 double-float wa4-%data% wa4-%offset%))
  30     (prog ((ic 0) (i 0) (idp2 0) (m1 0) (m2 0) (k 0) (ti12 0.0d0) (tr12 0.0d0)
  31            (ti11 0.0d0) (tr11 0.0d0) (arg 0.0d0) (m2s 0) (m1d 0))
  32       (declare (type (double-float) arg tr11 ti11 tr12 ti12)
  33                (type (f2cl-lib:integer4) m1d m2s k m2 m1 idp2 i ic))
  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       (setf arg (/ (* 2.0d0 4.0d0 (atan 1.0d0)) 5.0d0))
  39       (setf tr11 (cos arg))
  40       (setf ti11 (sin arg))
  41       (setf tr12 (cos (* 2.0d0 arg)))
  42       (setf ti12 (sin (* 2.0d0 arg)))
  43       (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  44                     ((> k l1) nil)
  45         (tagbody
  46           (setf m2 m2s)
  47           (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
  48                         ((> m1 m1d) nil)
  49             (tagbody
  50               (setf m2 (f2cl-lib:int-add m2 im2))
  51               (setf (f2cl-lib:fref ch-%data%
  52                                    (m2 1 k 1)
  53                                    ((1 in2) (1 ido) (1 l1) (1 5))
  54                                    ch-%offset%)
  55                       (+
  56                        (f2cl-lib:fref cc-%data%
  57                                       (m1 1 1 k)
  58                                       ((1 in1) (1 ido) (1 5) (1 l1))
  59                                       cc-%offset%)
  60                        (* 2.0d0
  61                           (f2cl-lib:fref cc-%data%
  62                                          (m1 ido 2 k)
  63                                          ((1 in1) (1 ido) (1 5) (1 l1))
  64                                          cc-%offset%))
  65                        (* 2.0d0
  66                           (f2cl-lib:fref cc-%data%
  67                                          (m1 ido 4 k)
  68                                          ((1 in1) (1 ido) (1 5) (1 l1))
  69                                          cc-%offset%))))
  70               (setf (f2cl-lib:fref ch-%data%
  71                                    (m2 1 k 2)
  72                                    ((1 in2) (1 ido) (1 l1) (1 5))
  73                                    ch-%offset%)
  74                       (-
  75                        (+
  76                         (f2cl-lib:fref cc-%data%
  77                                        (m1 1 1 k)
  78                                        ((1 in1) (1 ido) (1 5) (1 l1))
  79                                        cc-%offset%)
  80                         (* tr11
  81                            2.0d0
  82                            (f2cl-lib:fref cc-%data%
  83                                           (m1 ido 2 k)
  84                                           ((1 in1) (1 ido) (1 5) (1 l1))
  85                                           cc-%offset%))
  86                         (* tr12
  87                            2.0d0
  88                            (f2cl-lib:fref cc-%data%
  89                                           (m1 ido 4 k)
  90                                           ((1 in1) (1 ido) (1 5) (1 l1))
  91                                           cc-%offset%)))
  92                        (+
  93                         (* ti11
  94                            2.0d0
  95                            (f2cl-lib:fref cc-%data%
  96                                           (m1 1 3 k)
  97                                           ((1 in1) (1 ido) (1 5) (1 l1))
  98                                           cc-%offset%))
  99                         (* ti12
 100                            2.0d0
 101                            (f2cl-lib:fref cc-%data%
 102                                           (m1 1 5 k)
 103                                           ((1 in1) (1 ido) (1 5) (1 l1))
 104                                           cc-%offset%)))))
 105               (setf (f2cl-lib:fref ch-%data%
 106                                    (m2 1 k 3)
 107                                    ((1 in2) (1 ido) (1 l1) (1 5))
 108                                    ch-%offset%)
 109                       (-
 110                        (+
 111                         (f2cl-lib:fref cc-%data%
 112                                        (m1 1 1 k)
 113                                        ((1 in1) (1 ido) (1 5) (1 l1))
 114                                        cc-%offset%)
 115                         (* tr12
 116                            2.0d0
 117                            (f2cl-lib:fref cc-%data%
 118                                           (m1 ido 2 k)
 119                                           ((1 in1) (1 ido) (1 5) (1 l1))
 120                                           cc-%offset%))
 121                         (* tr11
 122                            2.0d0
 123                            (f2cl-lib:fref cc-%data%
 124                                           (m1 ido 4 k)
 125                                           ((1 in1) (1 ido) (1 5) (1 l1))
 126                                           cc-%offset%)))
 127                        (+
 128                         (* ti12
 129                            2.0d0
 130                            (f2cl-lib:fref cc-%data%
 131                                           (m1 1 3 k)
 132                                           ((1 in1) (1 ido) (1 5) (1 l1))
 133                                           cc-%offset%))
 134                         (* (- ti11)
 135                            2.0d0
 136                            (f2cl-lib:fref cc-%data%
 137                                           (m1 1 5 k)
 138                                           ((1 in1) (1 ido) (1 5) (1 l1))
 139                                           cc-%offset%)))))
 140               (setf (f2cl-lib:fref ch-%data%
 141                                    (m2 1 k 4)
 142                                    ((1 in2) (1 ido) (1 l1) (1 5))
 143                                    ch-%offset%)
 144                       (+
 145                        (f2cl-lib:fref cc-%data%
 146                                       (m1 1 1 k)
 147                                       ((1 in1) (1 ido) (1 5) (1 l1))
 148                                       cc-%offset%)
 149                        (* tr12
 150                           2.0d0
 151                           (f2cl-lib:fref cc-%data%
 152                                          (m1 ido 2 k)
 153                                          ((1 in1) (1 ido) (1 5) (1 l1))
 154                                          cc-%offset%))
 155                        (* tr11
 156                           2.0d0
 157                           (f2cl-lib:fref cc-%data%
 158                                          (m1 ido 4 k)
 159                                          ((1 in1) (1 ido) (1 5) (1 l1))
 160                                          cc-%offset%))
 161                        (+
 162                         (* ti12
 163                            2.0d0
 164                            (f2cl-lib:fref cc-%data%
 165                                           (m1 1 3 k)
 166                                           ((1 in1) (1 ido) (1 5) (1 l1))
 167                                           cc-%offset%))
 168                         (* (- ti11)
 169                            2.0d0
 170                            (f2cl-lib:fref cc-%data%
 171                                           (m1 1 5 k)
 172                                           ((1 in1) (1 ido) (1 5) (1 l1))
 173                                           cc-%offset%)))))
 174               (setf (f2cl-lib:fref ch-%data%
 175                                    (m2 1 k 5)
 176                                    ((1 in2) (1 ido) (1 l1) (1 5))
 177                                    ch-%offset%)
 178                       (+
 179                        (f2cl-lib:fref cc-%data%
 180                                       (m1 1 1 k)
 181                                       ((1 in1) (1 ido) (1 5) (1 l1))
 182                                       cc-%offset%)
 183                        (* tr11
 184                           2.0d0
 185                           (f2cl-lib:fref cc-%data%
 186                                          (m1 ido 2 k)
 187                                          ((1 in1) (1 ido) (1 5) (1 l1))
 188                                          cc-%offset%))
 189                        (* tr12
 190                           2.0d0
 191                           (f2cl-lib:fref cc-%data%
 192                                          (m1 ido 4 k)
 193                                          ((1 in1) (1 ido) (1 5) (1 l1))
 194                                          cc-%offset%))
 195                        (+
 196                         (* ti11
 197                            2.0d0
 198                            (f2cl-lib:fref cc-%data%
 199                                           (m1 1 3 k)
 200                                           ((1 in1) (1 ido) (1 5) (1 l1))
 201                                           cc-%offset%))
 202                         (* ti12
 203                            2.0d0
 204                            (f2cl-lib:fref cc-%data%
 205                                           (m1 1 5 k)
 206                                           ((1 in1) (1 ido) (1 5) (1 l1))
 207                                           cc-%offset%)))))
 208              label1001))
 209          label101))
 210       (if (= ido 1) (go end_label))
 211       (setf idp2 (f2cl-lib:int-add ido 2))
 212       (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
 213                     ((> k l1) nil)
 214         (tagbody
 215           (f2cl-lib:fdo (i 3 (f2cl-lib:int-add i 2))
 216                         ((> i ido) nil)
 217             (tagbody
 218               (setf ic (f2cl-lib:int-sub idp2 i))
 219               (setf m2 m2s)
 220               (f2cl-lib:fdo (m1 1 (f2cl-lib:int-add m1 im1))
 221                             ((> m1 m1d) nil)
 222                 (tagbody
 223                   (setf m2 (f2cl-lib:int-add m2 im2))
 224                   (setf (f2cl-lib:fref ch-%data%
 225                                        (m2 (f2cl-lib:int-sub i 1) k 1)
 226                                        ((1 in2) (1 ido) (1 l1) (1 5))
 227                                        ch-%offset%)
 228                           (+
 229                            (f2cl-lib:fref cc-%data%
 230                                           (m1 (f2cl-lib:int-sub i 1) 1 k)
 231                                           ((1 in1) (1 ido) (1 5) (1 l1))
 232                                           cc-%offset%)
 233                            (+
 234                             (f2cl-lib:fref cc-%data%
 235                                            (m1 (f2cl-lib:int-sub i 1) 3 k)
 236                                            ((1 in1) (1 ido) (1 5) (1 l1))
 237                                            cc-%offset%)
 238                             (f2cl-lib:fref cc-%data%
 239                                            (m1 (f2cl-lib:int-sub ic 1) 2 k)
 240                                            ((1 in1) (1 ido) (1 5) (1 l1))
 241                                            cc-%offset%))
 242                            (+
 243                             (f2cl-lib:fref cc-%data%
 244                                            (m1 (f2cl-lib:int-sub i 1) 5 k)
 245                                            ((1 in1) (1 ido) (1 5) (1 l1))
 246                                            cc-%offset%)
 247                             (f2cl-lib:fref cc-%data%
 248                                            (m1 (f2cl-lib:int-sub ic 1) 4 k)
 249                                            ((1 in1) (1 ido) (1 5) (1 l1))
 250                                            cc-%offset%))))
 251                   (setf (f2cl-lib:fref ch-%data%
 252                                        (m2 i k 1)
 253                                        ((1 in2) (1 ido) (1 l1) (1 5))
 254                                        ch-%offset%)
 255                           (+
 256                            (f2cl-lib:fref cc-%data%
 257                                           (m1 i 1 k)
 258                                           ((1 in1) (1 ido) (1 5) (1 l1))
 259                                           cc-%offset%)
 260                            (-
 261                             (f2cl-lib:fref cc-%data%
 262                                            (m1 i 3 k)
 263                                            ((1 in1) (1 ido) (1 5) (1 l1))
 264                                            cc-%offset%)
 265                             (f2cl-lib:fref cc-%data%
 266                                            (m1 ic 2 k)
 267                                            ((1 in1) (1 ido) (1 5) (1 l1))
 268                                            cc-%offset%))
 269                            (-
 270                             (f2cl-lib:fref cc-%data%
 271                                            (m1 i 5 k)
 272                                            ((1 in1) (1 ido) (1 5) (1 l1))
 273                                            cc-%offset%)
 274                             (f2cl-lib:fref cc-%data%
 275                                            (m1 ic 4 k)
 276                                            ((1 in1) (1 ido) (1 5) (1 l1))
 277                                            cc-%offset%))))
 278                   (setf (f2cl-lib:fref ch-%data%
 279                                        (m2 (f2cl-lib:int-sub i 1) k 2)
 280                                        ((1 in2) (1 ido) (1 l1) (1 5))
 281                                        ch-%offset%)
 282                           (-
 283                            (*
 284                             (f2cl-lib:fref wa1-%data%
 285                                            ((f2cl-lib:int-sub i 2))
 286                                            ((1 ido))
 287                                            wa1-%offset%)
 288                             (-
 289                              (+
 290                               (f2cl-lib:fref cc-%data%
 291                                              (m1 (f2cl-lib:int-sub i 1) 1 k)
 292                                              ((1 in1) (1 ido) (1 5) (1 l1))
 293                                              cc-%offset%)
 294                               (* tr11
 295                                  (+
 296                                   (f2cl-lib:fref cc-%data%
 297                                                  (m1 (f2cl-lib:int-sub i 1) 3
 298                                                   k)
 299                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 300                                                  cc-%offset%)
 301                                   (f2cl-lib:fref cc-%data%
 302                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 303                                                   k)
 304                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 305                                                  cc-%offset%)))
 306                               (* tr12
 307                                  (+
 308                                   (f2cl-lib:fref cc-%data%
 309                                                  (m1 (f2cl-lib:int-sub i 1) 5
 310                                                   k)
 311                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 312                                                  cc-%offset%)
 313                                   (f2cl-lib:fref cc-%data%
 314                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 315                                                   k)
 316                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 317                                                  cc-%offset%))))
 318                              (+
 319                               (* ti11
 320                                  (+
 321                                   (f2cl-lib:fref cc-%data%
 322                                                  (m1 i 3 k)
 323                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 324                                                  cc-%offset%)
 325                                   (f2cl-lib:fref cc-%data%
 326                                                  (m1 ic 2 k)
 327                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 328                                                  cc-%offset%)))
 329                               (* ti12
 330                                  (+
 331                                   (f2cl-lib:fref cc-%data%
 332                                                  (m1 i 5 k)
 333                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 334                                                  cc-%offset%)
 335                                   (f2cl-lib:fref cc-%data%
 336                                                  (m1 ic 4 k)
 337                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 338                                                  cc-%offset%))))))
 339                            (*
 340                             (f2cl-lib:fref wa1-%data%
 341                                            ((f2cl-lib:int-sub i 1))
 342                                            ((1 ido))
 343                                            wa1-%offset%)
 344                             (+
 345                              (f2cl-lib:fref cc-%data%
 346                                             (m1 i 1 k)
 347                                             ((1 in1) (1 ido) (1 5) (1 l1))
 348                                             cc-%offset%)
 349                              (* tr11
 350                                 (-
 351                                  (f2cl-lib:fref cc-%data%
 352                                                 (m1 i 3 k)
 353                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 354                                                 cc-%offset%)
 355                                  (f2cl-lib:fref cc-%data%
 356                                                 (m1 ic 2 k)
 357                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 358                                                 cc-%offset%)))
 359                              (* tr12
 360                                 (-
 361                                  (f2cl-lib:fref cc-%data%
 362                                                 (m1 i 5 k)
 363                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 364                                                 cc-%offset%)
 365                                  (f2cl-lib:fref cc-%data%
 366                                                 (m1 ic 4 k)
 367                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 368                                                 cc-%offset%)))
 369                              (+
 370                               (* ti11
 371                                  (-
 372                                   (f2cl-lib:fref cc-%data%
 373                                                  (m1 (f2cl-lib:int-sub i 1) 3
 374                                                   k)
 375                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 376                                                  cc-%offset%)
 377                                   (f2cl-lib:fref cc-%data%
 378                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 379                                                   k)
 380                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 381                                                  cc-%offset%)))
 382                               (* ti12
 383                                  (-
 384                                   (f2cl-lib:fref cc-%data%
 385                                                  (m1 (f2cl-lib:int-sub i 1) 5
 386                                                   k)
 387                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 388                                                  cc-%offset%)
 389                                   (f2cl-lib:fref cc-%data%
 390                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 391                                                   k)
 392                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 393                                                  cc-%offset%))))))))
 394                   (setf (f2cl-lib:fref ch-%data%
 395                                        (m2 i k 2)
 396                                        ((1 in2) (1 ido) (1 l1) (1 5))
 397                                        ch-%offset%)
 398                           (+
 399                            (*
 400                             (f2cl-lib:fref wa1-%data%
 401                                            ((f2cl-lib:int-sub i 2))
 402                                            ((1 ido))
 403                                            wa1-%offset%)
 404                             (+
 405                              (f2cl-lib:fref cc-%data%
 406                                             (m1 i 1 k)
 407                                             ((1 in1) (1 ido) (1 5) (1 l1))
 408                                             cc-%offset%)
 409                              (* tr11
 410                                 (-
 411                                  (f2cl-lib:fref cc-%data%
 412                                                 (m1 i 3 k)
 413                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 414                                                 cc-%offset%)
 415                                  (f2cl-lib:fref cc-%data%
 416                                                 (m1 ic 2 k)
 417                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 418                                                 cc-%offset%)))
 419                              (* tr12
 420                                 (-
 421                                  (f2cl-lib:fref cc-%data%
 422                                                 (m1 i 5 k)
 423                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 424                                                 cc-%offset%)
 425                                  (f2cl-lib:fref cc-%data%
 426                                                 (m1 ic 4 k)
 427                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 428                                                 cc-%offset%)))
 429                              (+
 430                               (* ti11
 431                                  (-
 432                                   (f2cl-lib:fref cc-%data%
 433                                                  (m1 (f2cl-lib:int-sub i 1) 3
 434                                                   k)
 435                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 436                                                  cc-%offset%)
 437                                   (f2cl-lib:fref cc-%data%
 438                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 439                                                   k)
 440                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 441                                                  cc-%offset%)))
 442                               (* ti12
 443                                  (-
 444                                   (f2cl-lib:fref cc-%data%
 445                                                  (m1 (f2cl-lib:int-sub i 1) 5
 446                                                   k)
 447                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 448                                                  cc-%offset%)
 449                                   (f2cl-lib:fref cc-%data%
 450                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 451                                                   k)
 452                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 453                                                  cc-%offset%))))))
 454                            (*
 455                             (f2cl-lib:fref wa1-%data%
 456                                            ((f2cl-lib:int-sub i 1))
 457                                            ((1 ido))
 458                                            wa1-%offset%)
 459                             (-
 460                              (+
 461                               (f2cl-lib:fref cc-%data%
 462                                              (m1 (f2cl-lib:int-sub i 1) 1 k)
 463                                              ((1 in1) (1 ido) (1 5) (1 l1))
 464                                              cc-%offset%)
 465                               (* tr11
 466                                  (+
 467                                   (f2cl-lib:fref cc-%data%
 468                                                  (m1 (f2cl-lib:int-sub i 1) 3
 469                                                   k)
 470                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 471                                                  cc-%offset%)
 472                                   (f2cl-lib:fref cc-%data%
 473                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 474                                                   k)
 475                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 476                                                  cc-%offset%)))
 477                               (* tr12
 478                                  (+
 479                                   (f2cl-lib:fref cc-%data%
 480                                                  (m1 (f2cl-lib:int-sub i 1) 5
 481                                                   k)
 482                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 483                                                  cc-%offset%)
 484                                   (f2cl-lib:fref cc-%data%
 485                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 486                                                   k)
 487                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 488                                                  cc-%offset%))))
 489                              (+
 490                               (* ti11
 491                                  (+
 492                                   (f2cl-lib:fref cc-%data%
 493                                                  (m1 i 3 k)
 494                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 495                                                  cc-%offset%)
 496                                   (f2cl-lib:fref cc-%data%
 497                                                  (m1 ic 2 k)
 498                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 499                                                  cc-%offset%)))
 500                               (* ti12
 501                                  (+
 502                                   (f2cl-lib:fref cc-%data%
 503                                                  (m1 i 5 k)
 504                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 505                                                  cc-%offset%)
 506                                   (f2cl-lib:fref cc-%data%
 507                                                  (m1 ic 4 k)
 508                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 509                                                  cc-%offset%))))))))
 510                   (setf (f2cl-lib:fref ch-%data%
 511                                        (m2 (f2cl-lib:int-sub i 1) k 3)
 512                                        ((1 in2) (1 ido) (1 l1) (1 5))
 513                                        ch-%offset%)
 514                           (-
 515                            (*
 516                             (f2cl-lib:fref wa2-%data%
 517                                            ((f2cl-lib:int-sub i 2))
 518                                            ((1 ido))
 519                                            wa2-%offset%)
 520                             (-
 521                              (+
 522                               (f2cl-lib:fref cc-%data%
 523                                              (m1 (f2cl-lib:int-sub i 1) 1 k)
 524                                              ((1 in1) (1 ido) (1 5) (1 l1))
 525                                              cc-%offset%)
 526                               (* tr12
 527                                  (+
 528                                   (f2cl-lib:fref cc-%data%
 529                                                  (m1 (f2cl-lib:int-sub i 1) 3
 530                                                   k)
 531                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 532                                                  cc-%offset%)
 533                                   (f2cl-lib:fref cc-%data%
 534                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 535                                                   k)
 536                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 537                                                  cc-%offset%)))
 538                               (* tr11
 539                                  (+
 540                                   (f2cl-lib:fref cc-%data%
 541                                                  (m1 (f2cl-lib:int-sub i 1) 5
 542                                                   k)
 543                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 544                                                  cc-%offset%)
 545                                   (f2cl-lib:fref cc-%data%
 546                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 547                                                   k)
 548                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 549                                                  cc-%offset%))))
 550                              (-
 551                               (* ti12
 552                                  (+
 553                                   (f2cl-lib:fref cc-%data%
 554                                                  (m1 i 3 k)
 555                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 556                                                  cc-%offset%)
 557                                   (f2cl-lib:fref cc-%data%
 558                                                  (m1 ic 2 k)
 559                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 560                                                  cc-%offset%)))
 561                               (* ti11
 562                                  (+
 563                                   (f2cl-lib:fref cc-%data%
 564                                                  (m1 i 5 k)
 565                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 566                                                  cc-%offset%)
 567                                   (f2cl-lib:fref cc-%data%
 568                                                  (m1 ic 4 k)
 569                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 570                                                  cc-%offset%))))))
 571                            (*
 572                             (f2cl-lib:fref wa2-%data%
 573                                            ((f2cl-lib:int-sub i 1))
 574                                            ((1 ido))
 575                                            wa2-%offset%)
 576                             (+
 577                              (f2cl-lib:fref cc-%data%
 578                                             (m1 i 1 k)
 579                                             ((1 in1) (1 ido) (1 5) (1 l1))
 580                                             cc-%offset%)
 581                              (* tr12
 582                                 (-
 583                                  (f2cl-lib:fref cc-%data%
 584                                                 (m1 i 3 k)
 585                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 586                                                 cc-%offset%)
 587                                  (f2cl-lib:fref cc-%data%
 588                                                 (m1 ic 2 k)
 589                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 590                                                 cc-%offset%)))
 591                              (* tr11
 592                                 (-
 593                                  (f2cl-lib:fref cc-%data%
 594                                                 (m1 i 5 k)
 595                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 596                                                 cc-%offset%)
 597                                  (f2cl-lib:fref cc-%data%
 598                                                 (m1 ic 4 k)
 599                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 600                                                 cc-%offset%)))
 601                              (-
 602                               (* ti12
 603                                  (-
 604                                   (f2cl-lib:fref cc-%data%
 605                                                  (m1 (f2cl-lib:int-sub i 1) 3
 606                                                   k)
 607                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 608                                                  cc-%offset%)
 609                                   (f2cl-lib:fref cc-%data%
 610                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 611                                                   k)
 612                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 613                                                  cc-%offset%)))
 614                               (* ti11
 615                                  (-
 616                                   (f2cl-lib:fref cc-%data%
 617                                                  (m1 (f2cl-lib:int-sub i 1) 5
 618                                                   k)
 619                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 620                                                  cc-%offset%)
 621                                   (f2cl-lib:fref cc-%data%
 622                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 623                                                   k)
 624                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 625                                                  cc-%offset%))))))))
 626                   (setf (f2cl-lib:fref ch-%data%
 627                                        (m2 i k 3)
 628                                        ((1 in2) (1 ido) (1 l1) (1 5))
 629                                        ch-%offset%)
 630                           (+
 631                            (*
 632                             (f2cl-lib:fref wa2-%data%
 633                                            ((f2cl-lib:int-sub i 2))
 634                                            ((1 ido))
 635                                            wa2-%offset%)
 636                             (+
 637                              (f2cl-lib:fref cc-%data%
 638                                             (m1 i 1 k)
 639                                             ((1 in1) (1 ido) (1 5) (1 l1))
 640                                             cc-%offset%)
 641                              (* tr12
 642                                 (-
 643                                  (f2cl-lib:fref cc-%data%
 644                                                 (m1 i 3 k)
 645                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 646                                                 cc-%offset%)
 647                                  (f2cl-lib:fref cc-%data%
 648                                                 (m1 ic 2 k)
 649                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 650                                                 cc-%offset%)))
 651                              (* tr11
 652                                 (-
 653                                  (f2cl-lib:fref cc-%data%
 654                                                 (m1 i 5 k)
 655                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 656                                                 cc-%offset%)
 657                                  (f2cl-lib:fref cc-%data%
 658                                                 (m1 ic 4 k)
 659                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 660                                                 cc-%offset%)))
 661                              (-
 662                               (* ti12
 663                                  (-
 664                                   (f2cl-lib:fref cc-%data%
 665                                                  (m1 (f2cl-lib:int-sub i 1) 3
 666                                                   k)
 667                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 668                                                  cc-%offset%)
 669                                   (f2cl-lib:fref cc-%data%
 670                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 671                                                   k)
 672                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 673                                                  cc-%offset%)))
 674                               (* ti11
 675                                  (-
 676                                   (f2cl-lib:fref cc-%data%
 677                                                  (m1 (f2cl-lib:int-sub i 1) 5
 678                                                   k)
 679                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 680                                                  cc-%offset%)
 681                                   (f2cl-lib:fref cc-%data%
 682                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 683                                                   k)
 684                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 685                                                  cc-%offset%))))))
 686                            (*
 687                             (f2cl-lib:fref wa2-%data%
 688                                            ((f2cl-lib:int-sub i 1))
 689                                            ((1 ido))
 690                                            wa2-%offset%)
 691                             (-
 692                              (+
 693                               (f2cl-lib:fref cc-%data%
 694                                              (m1 (f2cl-lib:int-sub i 1) 1 k)
 695                                              ((1 in1) (1 ido) (1 5) (1 l1))
 696                                              cc-%offset%)
 697                               (* tr12
 698                                  (+
 699                                   (f2cl-lib:fref cc-%data%
 700                                                  (m1 (f2cl-lib:int-sub i 1) 3
 701                                                   k)
 702                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 703                                                  cc-%offset%)
 704                                   (f2cl-lib:fref cc-%data%
 705                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 706                                                   k)
 707                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 708                                                  cc-%offset%)))
 709                               (* tr11
 710                                  (+
 711                                   (f2cl-lib:fref cc-%data%
 712                                                  (m1 (f2cl-lib:int-sub i 1) 5
 713                                                   k)
 714                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 715                                                  cc-%offset%)
 716                                   (f2cl-lib:fref cc-%data%
 717                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 718                                                   k)
 719                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 720                                                  cc-%offset%))))
 721                              (-
 722                               (* ti12
 723                                  (+
 724                                   (f2cl-lib:fref cc-%data%
 725                                                  (m1 i 3 k)
 726                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 727                                                  cc-%offset%)
 728                                   (f2cl-lib:fref cc-%data%
 729                                                  (m1 ic 2 k)
 730                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 731                                                  cc-%offset%)))
 732                               (* ti11
 733                                  (+
 734                                   (f2cl-lib:fref cc-%data%
 735                                                  (m1 i 5 k)
 736                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 737                                                  cc-%offset%)
 738                                   (f2cl-lib:fref cc-%data%
 739                                                  (m1 ic 4 k)
 740                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 741                                                  cc-%offset%))))))))
 742                   (setf (f2cl-lib:fref ch-%data%
 743                                        (m2 (f2cl-lib:int-sub i 1) k 4)
 744                                        ((1 in2) (1 ido) (1 l1) (1 5))
 745                                        ch-%offset%)
 746                           (-
 747                            (*
 748                             (f2cl-lib:fref wa3-%data%
 749                                            ((f2cl-lib:int-sub i 2))
 750                                            ((1 ido))
 751                                            wa3-%offset%)
 752                             (+
 753                              (f2cl-lib:fref cc-%data%
 754                                             (m1 (f2cl-lib:int-sub i 1) 1 k)
 755                                             ((1 in1) (1 ido) (1 5) (1 l1))
 756                                             cc-%offset%)
 757                              (* tr12
 758                                 (+
 759                                  (f2cl-lib:fref cc-%data%
 760                                                 (m1 (f2cl-lib:int-sub i 1) 3 k)
 761                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 762                                                 cc-%offset%)
 763                                  (f2cl-lib:fref cc-%data%
 764                                                 (m1 (f2cl-lib:int-sub ic 1) 2
 765                                                  k)
 766                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 767                                                 cc-%offset%)))
 768                              (* tr11
 769                                 (+
 770                                  (f2cl-lib:fref cc-%data%
 771                                                 (m1 (f2cl-lib:int-sub i 1) 5 k)
 772                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 773                                                 cc-%offset%)
 774                                  (f2cl-lib:fref cc-%data%
 775                                                 (m1 (f2cl-lib:int-sub ic 1) 4
 776                                                  k)
 777                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 778                                                 cc-%offset%)))
 779                              (-
 780                               (* ti12
 781                                  (+
 782                                   (f2cl-lib:fref cc-%data%
 783                                                  (m1 i 3 k)
 784                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 785                                                  cc-%offset%)
 786                                   (f2cl-lib:fref cc-%data%
 787                                                  (m1 ic 2 k)
 788                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 789                                                  cc-%offset%)))
 790                               (* ti11
 791                                  (+
 792                                   (f2cl-lib:fref cc-%data%
 793                                                  (m1 i 5 k)
 794                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 795                                                  cc-%offset%)
 796                                   (f2cl-lib:fref cc-%data%
 797                                                  (m1 ic 4 k)
 798                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 799                                                  cc-%offset%))))))
 800                            (*
 801                             (f2cl-lib:fref wa3-%data%
 802                                            ((f2cl-lib:int-sub i 1))
 803                                            ((1 ido))
 804                                            wa3-%offset%)
 805                             (-
 806                              (+
 807                               (f2cl-lib:fref cc-%data%
 808                                              (m1 i 1 k)
 809                                              ((1 in1) (1 ido) (1 5) (1 l1))
 810                                              cc-%offset%)
 811                               (* tr12
 812                                  (-
 813                                   (f2cl-lib:fref cc-%data%
 814                                                  (m1 i 3 k)
 815                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 816                                                  cc-%offset%)
 817                                   (f2cl-lib:fref cc-%data%
 818                                                  (m1 ic 2 k)
 819                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 820                                                  cc-%offset%)))
 821                               (* tr11
 822                                  (-
 823                                   (f2cl-lib:fref cc-%data%
 824                                                  (m1 i 5 k)
 825                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 826                                                  cc-%offset%)
 827                                   (f2cl-lib:fref cc-%data%
 828                                                  (m1 ic 4 k)
 829                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 830                                                  cc-%offset%))))
 831                              (-
 832                               (* ti12
 833                                  (-
 834                                   (f2cl-lib:fref cc-%data%
 835                                                  (m1 (f2cl-lib:int-sub i 1) 3
 836                                                   k)
 837                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 838                                                  cc-%offset%)
 839                                   (f2cl-lib:fref cc-%data%
 840                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 841                                                   k)
 842                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 843                                                  cc-%offset%)))
 844                               (* ti11
 845                                  (-
 846                                   (f2cl-lib:fref cc-%data%
 847                                                  (m1 (f2cl-lib:int-sub i 1) 5
 848                                                   k)
 849                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 850                                                  cc-%offset%)
 851                                   (f2cl-lib:fref cc-%data%
 852                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 853                                                   k)
 854                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 855                                                  cc-%offset%))))))))
 856                   (setf (f2cl-lib:fref ch-%data%
 857                                        (m2 i k 4)
 858                                        ((1 in2) (1 ido) (1 l1) (1 5))
 859                                        ch-%offset%)
 860                           (+
 861                            (*
 862                             (f2cl-lib:fref wa3-%data%
 863                                            ((f2cl-lib:int-sub i 2))
 864                                            ((1 ido))
 865                                            wa3-%offset%)
 866                             (-
 867                              (+
 868                               (f2cl-lib:fref cc-%data%
 869                                              (m1 i 1 k)
 870                                              ((1 in1) (1 ido) (1 5) (1 l1))
 871                                              cc-%offset%)
 872                               (* tr12
 873                                  (-
 874                                   (f2cl-lib:fref cc-%data%
 875                                                  (m1 i 3 k)
 876                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 877                                                  cc-%offset%)
 878                                   (f2cl-lib:fref cc-%data%
 879                                                  (m1 ic 2 k)
 880                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 881                                                  cc-%offset%)))
 882                               (* tr11
 883                                  (-
 884                                   (f2cl-lib:fref cc-%data%
 885                                                  (m1 i 5 k)
 886                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 887                                                  cc-%offset%)
 888                                   (f2cl-lib:fref cc-%data%
 889                                                  (m1 ic 4 k)
 890                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 891                                                  cc-%offset%))))
 892                              (-
 893                               (* ti12
 894                                  (-
 895                                   (f2cl-lib:fref cc-%data%
 896                                                  (m1 (f2cl-lib:int-sub i 1) 3
 897                                                   k)
 898                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 899                                                  cc-%offset%)
 900                                   (f2cl-lib:fref cc-%data%
 901                                                  (m1 (f2cl-lib:int-sub ic 1) 2
 902                                                   k)
 903                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 904                                                  cc-%offset%)))
 905                               (* ti11
 906                                  (-
 907                                   (f2cl-lib:fref cc-%data%
 908                                                  (m1 (f2cl-lib:int-sub i 1) 5
 909                                                   k)
 910                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 911                                                  cc-%offset%)
 912                                   (f2cl-lib:fref cc-%data%
 913                                                  (m1 (f2cl-lib:int-sub ic 1) 4
 914                                                   k)
 915                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 916                                                  cc-%offset%))))))
 917                            (*
 918                             (f2cl-lib:fref wa3-%data%
 919                                            ((f2cl-lib:int-sub i 1))
 920                                            ((1 ido))
 921                                            wa3-%offset%)
 922                             (+
 923                              (f2cl-lib:fref cc-%data%
 924                                             (m1 (f2cl-lib:int-sub i 1) 1 k)
 925                                             ((1 in1) (1 ido) (1 5) (1 l1))
 926                                             cc-%offset%)
 927                              (* tr12
 928                                 (+
 929                                  (f2cl-lib:fref cc-%data%
 930                                                 (m1 (f2cl-lib:int-sub i 1) 3 k)
 931                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 932                                                 cc-%offset%)
 933                                  (f2cl-lib:fref cc-%data%
 934                                                 (m1 (f2cl-lib:int-sub ic 1) 2
 935                                                  k)
 936                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 937                                                 cc-%offset%)))
 938                              (* tr11
 939                                 (+
 940                                  (f2cl-lib:fref cc-%data%
 941                                                 (m1 (f2cl-lib:int-sub i 1) 5 k)
 942                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 943                                                 cc-%offset%)
 944                                  (f2cl-lib:fref cc-%data%
 945                                                 (m1 (f2cl-lib:int-sub ic 1) 4
 946                                                  k)
 947                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 948                                                 cc-%offset%)))
 949                              (-
 950                               (* ti12
 951                                  (+
 952                                   (f2cl-lib:fref cc-%data%
 953                                                  (m1 i 3 k)
 954                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 955                                                  cc-%offset%)
 956                                   (f2cl-lib:fref cc-%data%
 957                                                  (m1 ic 2 k)
 958                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 959                                                  cc-%offset%)))
 960                               (* ti11
 961                                  (+
 962                                   (f2cl-lib:fref cc-%data%
 963                                                  (m1 i 5 k)
 964                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 965                                                  cc-%offset%)
 966                                   (f2cl-lib:fref cc-%data%
 967                                                  (m1 ic 4 k)
 968                                                  ((1 in1) (1 ido) (1 5) (1 l1))
 969                                                  cc-%offset%))))))))
 970                   (setf (f2cl-lib:fref ch-%data%
 971                                        (m2 (f2cl-lib:int-sub i 1) k 5)
 972                                        ((1 in2) (1 ido) (1 l1) (1 5))
 973                                        ch-%offset%)
 974                           (-
 975                            (*
 976                             (f2cl-lib:fref wa4-%data%
 977                                            ((f2cl-lib:int-sub i 2))
 978                                            ((1 ido))
 979                                            wa4-%offset%)
 980                             (+
 981                              (f2cl-lib:fref cc-%data%
 982                                             (m1 (f2cl-lib:int-sub i 1) 1 k)
 983                                             ((1 in1) (1 ido) (1 5) (1 l1))
 984                                             cc-%offset%)
 985                              (* tr11
 986                                 (+
 987                                  (f2cl-lib:fref cc-%data%
 988                                                 (m1 (f2cl-lib:int-sub i 1) 3 k)
 989                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 990                                                 cc-%offset%)
 991                                  (f2cl-lib:fref cc-%data%
 992                                                 (m1 (f2cl-lib:int-sub ic 1) 2
 993                                                  k)
 994                                                 ((1 in1) (1 ido) (1 5) (1 l1))
 995                                                 cc-%offset%)))
 996                              (* tr12
 997                                 (+
 998                                  (f2cl-lib:fref cc-%data%
 999                                                 (m1 (f2cl-lib:int-sub i 1) 5 k)
1000                                                 ((1 in1) (1 ido) (1 5) (1 l1))
1001                                                 cc-%offset%)
1002                                  (f2cl-lib:fref cc-%data%
1003                                                 (m1 (f2cl-lib:int-sub ic 1) 4
1004                                                  k)
1005                                                 ((1 in1) (1 ido) (1 5) (1 l1))
1006                                                 cc-%offset%)))
1007                              (+
1008                               (* ti11
1009                                  (+
1010                                   (f2cl-lib:fref cc-%data%
1011                                                  (m1 i 3 k)
1012                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1013                                                  cc-%offset%)
1014                                   (f2cl-lib:fref cc-%data%
1015                                                  (m1 ic 2 k)
1016                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1017                                                  cc-%offset%)))
1018                               (* ti12
1019                                  (+
1020                                   (f2cl-lib:fref cc-%data%
1021                                                  (m1 i 5 k)
1022                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1023                                                  cc-%offset%)
1024                                   (f2cl-lib:fref cc-%data%
1025                                                  (m1 ic 4 k)
1026                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1027                                                  cc-%offset%))))))
1028                            (*
1029                             (f2cl-lib:fref wa4-%data%
1030                                            ((f2cl-lib:int-sub i 1))
1031                                            ((1 ido))
1032                                            wa4-%offset%)
1033                             (-
1034                              (+
1035                               (f2cl-lib:fref cc-%data%
1036                                              (m1 i 1 k)
1037                                              ((1 in1) (1 ido) (1 5) (1 l1))
1038                                              cc-%offset%)
1039                               (* tr11
1040                                  (-
1041                                   (f2cl-lib:fref cc-%data%
1042                                                  (m1 i 3 k)
1043                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1044                                                  cc-%offset%)
1045                                   (f2cl-lib:fref cc-%data%
1046                                                  (m1 ic 2 k)
1047                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1048                                                  cc-%offset%)))
1049                               (* tr12
1050                                  (-
1051                                   (f2cl-lib:fref cc-%data%
1052                                                  (m1 i 5 k)
1053                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1054                                                  cc-%offset%)
1055                                   (f2cl-lib:fref cc-%data%
1056                                                  (m1 ic 4 k)
1057                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1058                                                  cc-%offset%))))
1059                              (+
1060                               (* ti11
1061                                  (-
1062                                   (f2cl-lib:fref cc-%data%
1063                                                  (m1 (f2cl-lib:int-sub i 1) 3
1064                                                   k)
1065                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1066                                                  cc-%offset%)
1067                                   (f2cl-lib:fref cc-%data%
1068                                                  (m1 (f2cl-lib:int-sub ic 1) 2
1069                                                   k)
1070                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1071                                                  cc-%offset%)))
1072                               (* ti12
1073                                  (-
1074                                   (f2cl-lib:fref cc-%data%
1075                                                  (m1 (f2cl-lib:int-sub i 1) 5
1076                                                   k)
1077                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1078                                                  cc-%offset%)
1079                                   (f2cl-lib:fref cc-%data%
1080                                                  (m1 (f2cl-lib:int-sub ic 1) 4
1081                                                   k)
1082                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1083                                                  cc-%offset%))))))))
1084                   (setf (f2cl-lib:fref ch-%data%
1085                                        (m2 i k 5)
1086                                        ((1 in2) (1 ido) (1 l1) (1 5))
1087                                        ch-%offset%)
1088                           (+
1089                            (*
1090                             (f2cl-lib:fref wa4-%data%
1091                                            ((f2cl-lib:int-sub i 2))
1092                                            ((1 ido))
1093                                            wa4-%offset%)
1094                             (-
1095                              (+
1096                               (f2cl-lib:fref cc-%data%
1097                                              (m1 i 1 k)
1098                                              ((1 in1) (1 ido) (1 5) (1 l1))
1099                                              cc-%offset%)
1100                               (* tr11
1101                                  (-
1102                                   (f2cl-lib:fref cc-%data%
1103                                                  (m1 i 3 k)
1104                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1105                                                  cc-%offset%)
1106                                   (f2cl-lib:fref cc-%data%
1107                                                  (m1 ic 2 k)
1108                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1109                                                  cc-%offset%)))
1110                               (* tr12
1111                                  (-
1112                                   (f2cl-lib:fref cc-%data%
1113                                                  (m1 i 5 k)
1114                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1115                                                  cc-%offset%)
1116                                   (f2cl-lib:fref cc-%data%
1117                                                  (m1 ic 4 k)
1118                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1119                                                  cc-%offset%))))
1120                              (+
1121                               (* ti11
1122                                  (-
1123                                   (f2cl-lib:fref cc-%data%
1124                                                  (m1 (f2cl-lib:int-sub i 1) 3
1125                                                   k)
1126                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1127                                                  cc-%offset%)
1128                                   (f2cl-lib:fref cc-%data%
1129                                                  (m1 (f2cl-lib:int-sub ic 1) 2
1130                                                   k)
1131                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1132                                                  cc-%offset%)))
1133                               (* ti12
1134                                  (-
1135                                   (f2cl-lib:fref cc-%data%
1136                                                  (m1 (f2cl-lib:int-sub i 1) 5
1137                                                   k)
1138                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1139                                                  cc-%offset%)
1140                                   (f2cl-lib:fref cc-%data%
1141                                                  (m1 (f2cl-lib:int-sub ic 1) 4
1142                                                   k)
1143                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1144                                                  cc-%offset%))))))
1145                            (*
1146                             (f2cl-lib:fref wa4-%data%
1147                                            ((f2cl-lib:int-sub i 1))
1148                                            ((1 ido))
1149                                            wa4-%offset%)
1150                             (+
1151                              (f2cl-lib:fref cc-%data%
1152                                             (m1 (f2cl-lib:int-sub i 1) 1 k)
1153                                             ((1 in1) (1 ido) (1 5) (1 l1))
1154                                             cc-%offset%)
1155                              (* tr11
1156                                 (+
1157                                  (f2cl-lib:fref cc-%data%
1158                                                 (m1 (f2cl-lib:int-sub i 1) 3 k)
1159                                                 ((1 in1) (1 ido) (1 5) (1 l1))
1160                                                 cc-%offset%)
1161                                  (f2cl-lib:fref cc-%data%
1162                                                 (m1 (f2cl-lib:int-sub ic 1) 2
1163                                                  k)
1164                                                 ((1 in1) (1 ido) (1 5) (1 l1))
1165                                                 cc-%offset%)))
1166                              (* tr12
1167                                 (+
1168                                  (f2cl-lib:fref cc-%data%
1169                                                 (m1 (f2cl-lib:int-sub i 1) 5 k)
1170                                                 ((1 in1) (1 ido) (1 5) (1 l1))
1171                                                 cc-%offset%)
1172                                  (f2cl-lib:fref cc-%data%
1173                                                 (m1 (f2cl-lib:int-sub ic 1) 4
1174                                                  k)
1175                                                 ((1 in1) (1 ido) (1 5) (1 l1))
1176                                                 cc-%offset%)))
1177                              (+
1178                               (* ti11
1179                                  (+
1180                                   (f2cl-lib:fref cc-%data%
1181                                                  (m1 i 3 k)
1182                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1183                                                  cc-%offset%)
1184                                   (f2cl-lib:fref cc-%data%
1185                                                  (m1 ic 2 k)
1186                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1187                                                  cc-%offset%)))
1188                               (* ti12
1189                                  (+
1190                                   (f2cl-lib:fref cc-%data%
1191                                                  (m1 i 5 k)
1192                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1193                                                  cc-%offset%)
1194                                   (f2cl-lib:fref cc-%data%
1195                                                  (m1 ic 4 k)
1196                                                  ((1 in1) (1 ido) (1 5) (1 l1))
1197                                                  cc-%offset%))))))))
1198                  label1002))
1199              label102))
1200          label103))
1201       (go end_label)
1202      end_label
1203       (return (values nil nil nil nil nil nil nil nil nil nil nil nil nil)))))
1204
1205 (in-package #:cl-user)
1206 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
1207 (eval-when (:load-toplevel :compile-toplevel :execute)
1208   (setf (gethash 'fortran-to-lisp::mradb5
1209                  fortran-to-lisp::*f2cl-function-info*)
1210           (fortran-to-lisp::make-f2cl-finfo
1211            :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
1212                         (fortran-to-lisp::integer4) (array double-float (*))
1213                         (fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
1214                         (array double-float (*)) (fortran-to-lisp::integer4)
1215                         (fortran-to-lisp::integer4) (array double-float (*))
1216                         (array double-float (*)) (array double-float (*))
1217                         (array double-float (*)))
1218            :return-values '(nil nil nil nil nil nil nil nil nil nil nil nil
1219                             nil)
1220            :calls 'nil)))
1221