← Back to branch summary

~ubuntu-branches/debian/squeeze/maxima/squeeze

« back to all changes in this revision

Viewing changes to src/numerical/slatec/zwrsk.lisp

  • Committer: Bazaar Package Importer
  • Author(s): Camm Maguire
  • Date: 2006-10-18 14:52:42 UTC
  • mto: (1.1.5 upstream)
  • mto: This revision was merged to the branch mainline in revision 4.
  • Revision ID: james.westby@ubuntu.com-20061018145242-vzyrm5hmxr8kiosf
Import upstream version 5.10.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/numerical/slatec/zwrsk.lisp
1
 
;;; Compiled by f2cl version 2.0 beta 2002-05-06
 
1
;;; Compiled by f2cl version 2.0 beta Date: 2006/01/31 15:11:05 
 
2
;;; Using Lisp CMU Common Lisp Snapshot 2006-01 (19C)
2
3
;;; 
3
4
;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
4
5
;;;           (:coerce-assigns :as-needed) (:array-type ':simple-array)
5
6
;;;           (:array-slicing nil) (:declare-common nil)
6
7
;;;           (:float-format double-float))
7
8
 
8
 
(in-package "SLATEC")
 
9
(in-package :slatec)
9
10
 
10
11
 
11
12
(defun zwrsk (zrr zri fnu kode n yr yi nz cwr cwi tol elim alim)
12
 
  (declare (type (simple-array double-float (*)) yr yi)
13
 
           (type f2cl-lib:integer4 kode n nz)
14
 
           (type (simple-array double-float (*)) cwr cwi)
15
 
           (type double-float zrr zri fnu tol elim alim))
 
13
  (declare (type (simple-array double-float (*)) cwi cwr)
 
14
           (type (simple-array double-float (*)) yi yr)
 
15
           (type f2cl-lib:integer4 nz n kode)
 
16
           (type double-float alim elim tol fnu zri zrr))
16
17
  (prog ((i 0) (nw 0) (act 0.0) (acw 0.0) (ascle 0.0) (cinui 0.0) (cinur 0.0)
17
18
         (csclr 0.0) (cti 0.0) (ctr 0.0) (c1i 0.0) (c1r 0.0) (c2i 0.0)
18
19
         (c2r 0.0) (pti 0.0) (ptr 0.0) (ract 0.0) (sti 0.0) (str 0.0))
19
 
    (declare
20
 
     (type double-float str sti ract ptr pti c2r c2i c1r c1i ctr cti csclr
21
 
      cinur cinui ascle acw act)
22
 
     (type f2cl-lib:integer4 nw i))
 
20
    (declare (type double-float str sti ract ptr pti c2r c2i c1r c1i ctr cti
 
21
                                csclr cinur cinui ascle acw act)
 
22
             (type f2cl-lib:integer4 nw i))
23
23
    (setf nz 0)
24
24
    (multiple-value-bind
25
 
        (var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-7 var-8 var-9 var-10)
 
25
          (var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-7 var-8 var-9 var-10)
26
26
        (zbknu zrr zri fnu kode 2 cwr cwi nw tol elim alim)
27
 
      (declare
28
 
       (ignore var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-8 var-9 var-10))
 
27
      (declare (ignore var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-8 var-9
 
28
                       var-10))
29
29
      (setf nw var-7))
30
30
    (if (/= nw 0) (go label50))
31
31
    (zrati zrr zri fnu n yr yi tol)