1 ;;; Compiled by f2cl version:
2 ;;; ("f2cl1.l,v 2edcbd958861 2012/05/30 03:34:52 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 3fe93de3be82 2012/05/06 02:17:14 toy $"
7 ;;; "f2cl6.l,v 1d5cbacbb977 2008/08/24 00:56:27 rtoy $"
8 ;;; "macros.l,v 3fe93de3be82 2012/05/06 02:17:14 toy $")
10 ;;; Using Lisp CMU Common Lisp 20d (20D Unicode)
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))
20 (let* ((zero (f2cl-lib:cmplx
0.0 0.0)))
21 (declare (type (f2cl-lib:complex16
) zero
) (ignorable zero
))
22 (defun zhpr (uplo n alpha x incx ap
)
23 (declare (type (array f2cl-lib
:complex16
(*)) ap x
)
24 (type (double-float) alpha
)
25 (type (f2cl-lib:integer4
) incx n
)
26 (type (simple-string *) uplo
))
27 (f2cl-lib:with-multi-array-data
28 ((uplo character uplo-%data% uplo-%offset%
)
29 (x f2cl-lib
:complex16 x-%data% x-%offset%
)
30 (ap f2cl-lib
:complex16 ap-%data% ap-%offset%
))
31 (prog ((i 0) (info 0) (ix 0) (j 0) (jx 0) (k 0) (kk 0) (kx 0)
33 (declare (type (f2cl-lib:integer4
) i info ix j jx k kk kx
)
34 (type (f2cl-lib:complex16
) temp
))
37 ((and (not (lsame uplo
"U")) (not (lsame uplo
"L")))
47 (if (or (= n
0) (= alpha
(f2cl-lib:dble zero
))) (go end_label
))
52 (f2cl-lib:int-mul
(f2cl-lib:int-sub n
1)
61 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
65 ((/= (f2cl-lib:fref x
(j) ((1 *))) zero
)
70 (f2cl-lib:fref x-%data%
76 (f2cl-lib:fdo
(i 1 (f2cl-lib:int-add i
1))
79 (f2cl-lib:int-sub
1)))
82 (setf (f2cl-lib:fref ap-%data%
87 (f2cl-lib:fref ap-%data%
92 (f2cl-lib:fref x-%data%
97 (setf k
(f2cl-lib:int-add k
1))
99 (setf (f2cl-lib:fref ap-%data%
101 (f2cl-lib:int-add kk j
)
108 (f2cl-lib:fref ap-%data%
110 (f2cl-lib:int-add kk j
)
116 (f2cl-lib:fref x-%data%
121 'f2cl-lib
:complex16
)))
123 (setf (f2cl-lib:fref ap-%data%
125 (f2cl-lib:int-add kk j
)
131 (f2cl-lib:fref ap-%data%
133 (f2cl-lib:int-add kk j
)
137 'f2cl-lib
:complex16
))))
138 (setf kk
(f2cl-lib:int-add kk j
))
142 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
146 ((/= (f2cl-lib:fref x
(jx) ((1 *))) zero
)
151 (f2cl-lib:fref x-%data%
155 'f2cl-lib
:complex16
))
157 (f2cl-lib:fdo
(k kk
(f2cl-lib:int-add k
1))
161 (f2cl-lib:int-sub
2)))
164 (setf (f2cl-lib:fref ap-%data%
169 (f2cl-lib:fref ap-%data%
174 (f2cl-lib:fref x-%data%
179 (setf ix
(f2cl-lib:int-add ix incx
))
181 (setf (f2cl-lib:fref ap-%data%
183 (f2cl-lib:int-add kk j
)
190 (f2cl-lib:fref ap-%data%
192 (f2cl-lib:int-add kk j
)
198 (f2cl-lib:fref x-%data%
203 'f2cl-lib
:complex16
)))
205 (setf (f2cl-lib:fref ap-%data%
207 (f2cl-lib:int-add kk j
)
213 (f2cl-lib:fref ap-%data%
215 (f2cl-lib:int-add kk j
)
219 'f2cl-lib
:complex16
))))
220 (setf jx
(f2cl-lib:int-add jx incx
))
221 (setf kk
(f2cl-lib:int-add kk j
))
226 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
230 ((/= (f2cl-lib:fref x
(j) ((1 *))) zero
)
235 (f2cl-lib:fref x-%data%
239 'f2cl-lib
:complex16
))
240 (setf (f2cl-lib:fref ap-%data%
(kk) ((1 *)) ap-%offset%
)
244 (f2cl-lib:fref ap-%data%
250 (f2cl-lib:fref x-%data%
254 'f2cl-lib
:complex16
))
255 (setf k
(f2cl-lib:int-add kk
1))
256 (f2cl-lib:fdo
(i (f2cl-lib:int-add j
1)
257 (f2cl-lib:int-add i
1))
260 (setf (f2cl-lib:fref ap-%data%
265 (f2cl-lib:fref ap-%data%
270 (f2cl-lib:fref x-%data%
275 (setf k
(f2cl-lib:int-add k
1))
278 (setf (f2cl-lib:fref ap-%data%
(kk) ((1 *)) ap-%offset%
)
281 (f2cl-lib:fref ap-%data%
285 'f2cl-lib
:complex16
))))
288 (f2cl-lib:int-sub
(f2cl-lib:int-add kk n
) j
)
293 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
297 ((/= (f2cl-lib:fref x
(jx) ((1 *))) zero
)
302 (f2cl-lib:fref x-%data%
306 'f2cl-lib
:complex16
))
307 (setf (f2cl-lib:fref ap-%data%
(kk) ((1 *)) ap-%offset%
)
311 (f2cl-lib:fref ap-%data%
317 (f2cl-lib:fref x-%data%
321 'f2cl-lib
:complex16
))
323 (f2cl-lib:fdo
(k (f2cl-lib:int-add kk
1)
324 (f2cl-lib:int-add k
1))
328 (f2cl-lib:int-sub j
)))
331 (setf ix
(f2cl-lib:int-add ix incx
))
332 (setf (f2cl-lib:fref ap-%data%
337 (f2cl-lib:fref ap-%data%
342 (f2cl-lib:fref x-%data%
349 (setf (f2cl-lib:fref ap-%data%
(kk) ((1 *)) ap-%offset%
)
352 (f2cl-lib:fref ap-%data%
356 'f2cl-lib
:complex16
))))
357 (setf jx
(f2cl-lib:int-add jx incx
))
360 (f2cl-lib:int-sub
(f2cl-lib:int-add kk n
) j
)
365 (return (values nil nil nil nil nil nil
))))))
367 (in-package #-gcl
#:cl-user
#+gcl
"CL-USER")
368 #+#.
(cl:if
(cl:find-package
'#:f2cl
) '(and) '(or))
369 (eval-when (:load-toplevel
:compile-toplevel
:execute
)
370 (setf (gethash 'fortran-to-lisp
::zhpr fortran-to-lisp
::*f2cl-function-info
*)
371 (fortran-to-lisp::make-f2cl-finfo
372 :arg-types
'((simple-string) (fortran-to-lisp::integer4
)
373 (double-float) (array fortran-to-lisp
::complex16
(*))
374 (fortran-to-lisp::integer4
)
375 (array fortran-to-lisp
::complex16
(*)))
376 :return-values
'(nil nil nil nil nil nil
)
377 :calls
'(fortran-to-lisp::xerbla fortran-to-lisp
::lsame
))))