;;; Compiled by f2cl version 2.0 beta 2002-05-06
;;; 
;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
;;;           (:coerce-assigns :as-needed) (:array-type ':array)
;;;           (:array-slicing t) (:declare-common nil)
;;;           (:float-format single-float))

(in-package "QUADPACK")


(defun dqpsrt (limit last$ maxerr ermax elist iord nrmax)
  (declare (type (array f2cl-lib:integer4 (*)) iord)
   (type (array double-float (*)) elist) (type double-float ermax)
   (type f2cl-lib:integer4 nrmax maxerr last$ limit))
  (f2cl-lib:with-multi-array-data
      ((elist double-float elist-%data% elist-%offset%)
       (iord f2cl-lib:integer4 iord-%data% iord-%offset%))
    (prog ((i 0) (ibeg 0) (ido 0) (isucc 0) (j 0) (jbnd 0) (jupbn 0) (k 0)
           (errmax 0.0d0) (errmin 0.0d0))
      (declare (type double-float errmin errmax)
       (type f2cl-lib:integer4 k jupbn jbnd j isucc ido ibeg i))
      (if (> last$ 2) (go label10))
      (f2cl-lib:fset (f2cl-lib:fref iord-%data% (1) ((1 last$)) iord-%offset%)
                     1)
      (f2cl-lib:fset (f2cl-lib:fref iord-%data% (2) ((1 last$)) iord-%offset%)
                     2)
      (go label90)
     label10
      (setf errmax
              (f2cl-lib:fref elist-%data% (maxerr) ((1 last$)) elist-%offset%))
      (if (= nrmax 1) (go label30))
      (setf ido (f2cl-lib:int-sub nrmax 1))
      (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
                    ((> i ido) nil)
        (tagbody
          (setf isucc
                  (f2cl-lib:fref iord-%data%
                                 ((f2cl-lib:int-sub nrmax 1))
                                 ((1 last$))
                                 iord-%offset%))
          (if
           (<= errmax
               (f2cl-lib:fref elist-%data% (isucc) ((1 last$)) elist-%offset%))
           (go label30))
          (f2cl-lib:fset
           (f2cl-lib:fref iord-%data% (nrmax) ((1 last$)) iord-%offset%)
           isucc)
          (setf nrmax (f2cl-lib:int-sub nrmax 1))
         label20))
     label30
      (setf jupbn last$)
      (if (> last$ (+ (the f2cl-lib:integer4 (truncate limit 2)) 2))
          (setf jupbn (f2cl-lib:int-sub (f2cl-lib:int-add limit 3) last$)))
      (setf errmin
              (f2cl-lib:fref elist-%data% (last$) ((1 last$)) elist-%offset%))
      (setf jbnd (f2cl-lib:int-sub jupbn 1))
      (setf ibeg (f2cl-lib:int-add nrmax 1))
      (if (> ibeg jbnd) (go label50))
      (f2cl-lib:fdo (i ibeg (f2cl-lib:int-add i 1))
                    ((> i jbnd) nil)
        (tagbody
          (setf isucc
                  (f2cl-lib:fref iord-%data% (i) ((1 last$)) iord-%offset%))
          (if
           (>= errmax
               (f2cl-lib:fref elist-%data% (isucc) ((1 last$)) elist-%offset%))
           (go label60))
          (f2cl-lib:fset
           (f2cl-lib:fref iord-%data%
                          ((f2cl-lib:int-sub i 1))
                          ((1 last$))
                          iord-%offset%)
           isucc)
         label40))
     label50
      (f2cl-lib:fset
       (f2cl-lib:fref iord-%data% (jbnd) ((1 last$)) iord-%offset%)
       maxerr)
      (f2cl-lib:fset
       (f2cl-lib:fref iord-%data% (jupbn) ((1 last$)) iord-%offset%)
       last$)
      (go label90)
     label60
      (f2cl-lib:fset
       (f2cl-lib:fref iord-%data%
                      ((f2cl-lib:int-sub i 1))
                      ((1 last$))
                      iord-%offset%)
       maxerr)
      (setf k jbnd)
      (f2cl-lib:fdo (j i (f2cl-lib:int-add j 1))
                    ((> j jbnd) nil)
        (tagbody
          (setf isucc
                  (f2cl-lib:fref iord-%data% (k) ((1 last$)) iord-%offset%))
          (if
           (< errmin
              (f2cl-lib:fref elist-%data% (isucc) ((1 last$)) elist-%offset%))
           (go label80))
          (f2cl-lib:fset
           (f2cl-lib:fref iord-%data%
                          ((f2cl-lib:int-add k 1))
                          ((1 last$))
                          iord-%offset%)
           isucc)
          (setf k (f2cl-lib:int-sub k 1))
         label70))
      (f2cl-lib:fset (f2cl-lib:fref iord-%data% (i) ((1 last$)) iord-%offset%)
                     last$)
      (go label90)
     label80
      (f2cl-lib:fset
       (f2cl-lib:fref iord-%data%
                      ((f2cl-lib:int-add k 1))
                      ((1 last$))
                      iord-%offset%)
       last$)
     label90
      (setf maxerr
              (f2cl-lib:fref iord-%data% (nrmax) ((1 last$)) iord-%offset%))
      (setf ermax
              (f2cl-lib:fref elist-%data% (maxerr) ((1 last$)) elist-%offset%))
      (go end_label)
     end_label
      (return (values nil nil maxerr ermax nil nil nrmax)))))