src/numerical/slatec/dqwgts.lisp

   1 ;;; Compiled by f2cl version:
   2 ;;; ("f2cl1.l,v 46c1f6a93b0d 2012/05/03 04:40:28 toy $"
   3 ;;;  "f2cl2.l,v 96616d88fb7e 2008/02/22 22:19:34 rtoy $"
   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 46c1f6a93b0d 2012/05/03 04:40:28 toy $"
   7 ;;;  "f2cl6.l,v 1d5cbacbb977 2008/08/24 00:56:27 rtoy $"
   8 ;;;  "macros.l,v fceac530ef0c 2011/11/26 04:02:26 toy $")
   9
  10 ;;; Using Lisp CMU Common Lisp snapshot-2012-04 (20C 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 double-float))
  16
  17 (in-package :slatec)
  18
  19
  20 (defun dqwgts (x a b alfa beta integr)
  21   (declare (type (f2cl-lib:integer4) integr)
  22            (type (double-float) beta alfa b a x))
  23   (prog ((bmx 0.0) (xma 0.0) (dqwgts 0.0))
  24     (declare (type (double-float) dqwgts xma bmx))
  25     (setf xma (- x a))
  26     (setf bmx (- b x))
  27     (setf dqwgts (* (expt xma alfa) (expt bmx beta)))
  28     (f2cl-lib:computed-goto (label40 label10 label20 label30) integr)
  29    label10
  30     (setf dqwgts (* dqwgts (f2cl-lib:flog xma)))
  31     (go label40)
  32    label20
  33     (setf dqwgts (* dqwgts (f2cl-lib:flog bmx)))
  34     (go label40)
  35    label30
  36     (setf dqwgts (* dqwgts (f2cl-lib:flog xma) (f2cl-lib:flog bmx)))
  37    label40
  38     (go end_label)
  39    end_label
  40     (return (values dqwgts nil nil nil nil nil nil))))
  41
  42 (in-package #:cl-user)
  43 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
  44 (eval-when (:load-toplevel :compile-toplevel :execute)
  45   (setf (gethash 'fortran-to-lisp::dqwgts
  46                  fortran-to-lisp::*f2cl-function-info*)
  47           (fortran-to-lisp::make-f2cl-finfo
  48            :arg-types '((double-float) (double-float) (double-float)
  49                         (double-float) (double-float)
  50                         (fortran-to-lisp::integer4))
  51            :return-values '(nil nil nil nil nil nil)
  52            :calls 'nil)))
  53