share/colnew/lisp/dmzsol.lisp

   1 ;;; Compiled by f2cl version:
   2 ;;; ("f2cl1.l,v 1.221 2010/05/26 19:25:52 rtoy Exp $"
   3 ;;;  "f2cl2.l,v 1.37 2008/02/22 22:19:33 rtoy Exp $"
   4 ;;;  "f2cl3.l,v 1.6 2008/02/22 22:19:33 rtoy Exp $"
   5 ;;;  "f2cl4.l,v 1.7 2008/02/22 22:19:34 rtoy Exp $"
   6 ;;;  "f2cl5.l,v 1.204 2010/02/23 05:21:30 rtoy Exp $"
   7 ;;;  "f2cl6.l,v 1.48 2008/08/24 00:56:27 rtoy Exp $"
   8 ;;;  "macros.l,v 1.114 2010/05/17 01:42:14 rtoy Exp $")
   9
  10 ;;; Using Lisp CMU Common Lisp CVS Head 2010-05-25 18:21:07 (20A 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 :colnew)
  18
  19
  20 (defun dmzsol (kd mstar n v z dmz)
  21   (declare (type (array double-float (*)) z)
  22            (type (array double-float (*)) dmz v)
  23            (type (f2cl-lib:integer4) n mstar kd))
  24   (f2cl-lib:with-multi-array-data
  25       ((v double-float v-%data% v-%offset%)
  26        (dmz double-float dmz-%data% dmz-%offset%)
  27        (z double-float z-%data% z-%offset%))
  28     (prog ((l 0) (fact 0.0) (j 0) (i 0) (jz 0))
  29       (declare (type double-float fact) (type (f2cl-lib:integer4) jz i j l))
  30       (setf jz 1)
  31       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  32                     ((> i n) nil)
  33         (tagbody
  34           (f2cl-lib:fdo (j 1 (f2cl-lib:int-add j 1))
  35                         ((> j mstar) nil)
  36             (tagbody
  37               (setf fact (f2cl-lib:fref z-%data% (jz) ((1 1)) z-%offset%))
  38               (f2cl-lib:fdo (l 1 (f2cl-lib:int-add l 1))
  39                             ((> l kd) nil)
  40                 (tagbody
  41                   (setf (f2cl-lib:fref dmz-%data%
  42                                        (l i)
  43                                        ((1 kd) (1 1))
  44                                        dmz-%offset%)
  45                           (+
  46                            (f2cl-lib:fref dmz-%data%
  47                                           (l i)
  48                                           ((1 kd) (1 1))
  49                                           dmz-%offset%)
  50                            (* fact
  51                               (f2cl-lib:fref v-%data%
  52                                              (l jz)
  53                                              ((1 kd) (1 1))
  54                                              v-%offset%))))
  55                  label10))
  56               (setf jz (f2cl-lib:int-add jz 1))
  57              label20))
  58          label30))
  59       (go end_label)
  60      end_label
  61       (return (values nil nil nil nil nil nil)))))
  62
  63 (in-package #:cl-user)
  64 #+#.(cl:if (cl:find-package '#:f2cl) '(and) '(or))
  65 (eval-when (:load-toplevel :compile-toplevel :execute)
  66   (setf (gethash 'fortran-to-lisp::dmzsol
  67                  fortran-to-lisp::*f2cl-function-info*)
  68           (fortran-to-lisp::make-f2cl-finfo
  69            :arg-types '((fortran-to-lisp::integer4) (fortran-to-lisp::integer4)
  70                         (fortran-to-lisp::integer4) (array double-float (*))
  71                         (array double-float (1)) (array double-float (*)))
  72            :return-values '(nil nil nil nil nil nil)
  73            :calls 'nil)))
  74