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* ((one (f2cl-lib:cmplx
1.0 0.0)) (zero (f2cl-lib:cmplx
0.0 0.0)))
21 (declare (type (f2cl-lib:complex16
) one
)
22 (type (f2cl-lib:complex16
) zero
)
24 (defun zgbmv (trans m n kl ku alpha a lda x incx beta y incy
)
25 (declare (type (array f2cl-lib
:complex16
(*)) y x a
)
26 (type (f2cl-lib:complex16
) beta alpha
)
27 (type (f2cl-lib:integer4
) incy incx lda ku kl n m
)
28 (type (simple-string *) trans
))
29 (f2cl-lib:with-multi-array-data
30 ((trans character trans-%data% trans-%offset%
)
31 (a f2cl-lib
:complex16 a-%data% a-%offset%
)
32 (x f2cl-lib
:complex16 x-%data% x-%offset%
)
33 (y f2cl-lib
:complex16 y-%data% y-%offset%
))
34 (prog ((noconj nil
) (i 0) (info 0) (ix 0) (iy 0) (j 0) (jx 0) (jy 0)
35 (k 0) (kup1 0) (kx 0) (ky 0) (lenx 0) (leny 0) (temp #C
(0.0
0.0)))
36 (declare (type f2cl-lib
:logical noconj
)
37 (type (f2cl-lib:integer4
) i info ix iy j jx jy k kup1 kx ky
39 (type (f2cl-lib:complex16
) temp
))
42 ((and (not (lsame trans
"N"))
43 (not (lsame trans
"T"))
44 (not (lsame trans
"C")))
54 ((< lda
(f2cl-lib:int-add kl ku
1))
62 (xerbla "ZGBMV " info
)
64 (if (or (= m
0) (= n
0) (and (= alpha zero
) (= beta one
)))
66 (setf noconj
(lsame trans
"T"))
81 (f2cl-lib:int-sub lenx
1)
90 (f2cl-lib:int-sub leny
1)
98 (f2cl-lib:fdo
(i 1 (f2cl-lib:int-add i
1))
101 (setf (f2cl-lib:fref y-%data%
(i) ((1 *)) y-%offset%
)
105 (f2cl-lib:fdo
(i 1 (f2cl-lib:int-add i
1))
108 (setf (f2cl-lib:fref y-%data%
(i) ((1 *)) y-%offset%
)
110 (f2cl-lib:fref y-%data%
119 (f2cl-lib:fdo
(i 1 (f2cl-lib:int-add i
1))
122 (setf (f2cl-lib:fref y-%data%
(iy) ((1 *)) y-%offset%
)
124 (setf iy
(f2cl-lib:int-add iy incy
))
127 (f2cl-lib:fdo
(i 1 (f2cl-lib:int-add i
1))
130 (setf (f2cl-lib:fref y-%data%
(iy) ((1 *)) y-%offset%
)
132 (f2cl-lib:fref y-%data%
136 (setf iy
(f2cl-lib:int-add iy incy
))
138 (if (= alpha zero
) (go end_label
))
139 (setf kup1
(f2cl-lib:int-add ku
1))
145 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
149 ((/= (f2cl-lib:fref x
(jx) ((1 *))) zero
)
152 (f2cl-lib:fref x-%data%
156 (setf k
(f2cl-lib:int-sub kup1 j
))
158 (max (the f2cl-lib
:integer4
1)
159 (the f2cl-lib
:integer4
163 (f2cl-lib:int-add i
1))
165 (min (the f2cl-lib
:integer4 m
)
166 (the f2cl-lib
:integer4
167 (f2cl-lib:int-add j kl
))))
170 (setf (f2cl-lib:fref y-%data%
(i) ((1 *)) y-%offset%
)
172 (f2cl-lib:fref y-%data%
177 (f2cl-lib:fref a-%data%
178 ((f2cl-lib:int-add k i
) j
)
182 (setf jx
(f2cl-lib:int-add jx incx
))
185 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
189 ((/= (f2cl-lib:fref x
(jx) ((1 *))) zero
)
192 (f2cl-lib:fref x-%data%
197 (setf k
(f2cl-lib:int-sub kup1 j
))
199 (max (the f2cl-lib
:integer4
1)
200 (the f2cl-lib
:integer4
204 (f2cl-lib:int-add i
1))
206 (min (the f2cl-lib
:integer4 m
)
207 (the f2cl-lib
:integer4
208 (f2cl-lib:int-add j kl
))))
211 (setf (f2cl-lib:fref y-%data%
(iy) ((1 *)) y-%offset%
)
213 (f2cl-lib:fref y-%data%
218 (f2cl-lib:fref a-%data%
219 ((f2cl-lib:int-add k i
) j
)
222 (setf iy
(f2cl-lib:int-add iy incy
))
224 (setf jx
(f2cl-lib:int-add jx incx
))
225 (if (> j ku
) (setf ky
(f2cl-lib:int-add ky incy
)))
231 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
235 (setf k
(f2cl-lib:int-sub kup1 j
))
239 (max (the f2cl-lib
:integer4
1)
240 (the f2cl-lib
:integer4
244 (f2cl-lib:int-add i
1))
246 (min (the f2cl-lib
:integer4 m
)
247 (the f2cl-lib
:integer4
248 (f2cl-lib:int-add j kl
))))
254 (f2cl-lib:fref a-%data%
255 ((f2cl-lib:int-add k i
) j
)
258 (f2cl-lib:fref x-%data%
265 (max (the f2cl-lib
:integer4
1)
266 (the f2cl-lib
:integer4
270 (f2cl-lib:int-add i
1))
272 (min (the f2cl-lib
:integer4 m
)
273 (the f2cl-lib
:integer4
274 (f2cl-lib:int-add j kl
))))
281 (f2cl-lib:fref a-%data%
282 ((f2cl-lib:int-add k i
) j
)
285 (f2cl-lib:fref x-%data%
290 (setf (f2cl-lib:fref y-%data%
(jy) ((1 *)) y-%offset%
)
291 (+ (f2cl-lib:fref y-%data%
(jy) ((1 *)) y-%offset%
)
293 (setf jy
(f2cl-lib:int-add jy incy
))
296 (f2cl-lib:fdo
(j 1 (f2cl-lib:int-add j
1))
301 (setf k
(f2cl-lib:int-sub kup1 j
))
305 (max (the f2cl-lib
:integer4
1)
306 (the f2cl-lib
:integer4
310 (f2cl-lib:int-add i
1))
312 (min (the f2cl-lib
:integer4 m
)
313 (the f2cl-lib
:integer4
314 (f2cl-lib:int-add j kl
))))
320 (f2cl-lib:fref a-%data%
321 ((f2cl-lib:int-add k i
) j
)
324 (f2cl-lib:fref x-%data%
328 (setf ix
(f2cl-lib:int-add ix incx
))
332 (max (the f2cl-lib
:integer4
1)
333 (the f2cl-lib
:integer4
337 (f2cl-lib:int-add i
1))
339 (min (the f2cl-lib
:integer4 m
)
340 (the f2cl-lib
:integer4
341 (f2cl-lib:int-add j kl
))))
348 (f2cl-lib:fref a-%data%
349 ((f2cl-lib:int-add k i
) j
)
352 (f2cl-lib:fref x-%data%
356 (setf ix
(f2cl-lib:int-add ix incx
))
358 (setf (f2cl-lib:fref y-%data%
(jy) ((1 *)) y-%offset%
)
359 (+ (f2cl-lib:fref y-%data%
(jy) ((1 *)) y-%offset%
)
361 (setf jy
(f2cl-lib:int-add jy incy
))
362 (if (> j ku
) (setf kx
(f2cl-lib:int-add kx incx
)))
367 (values nil nil nil nil nil nil nil nil nil nil nil nil nil
))))))
369 (in-package #-gcl
#:cl-user
#+gcl
"CL-USER")
370 #+#.
(cl:if
(cl:find-package
'#:f2cl
) '(and) '(or))
371 (eval-when (:load-toplevel
:compile-toplevel
:execute
)
372 (setf (gethash 'fortran-to-lisp
::zgbmv fortran-to-lisp
::*f2cl-function-info
*)
373 (fortran-to-lisp::make-f2cl-finfo
374 :arg-types
'((simple-string) (fortran-to-lisp::integer4
)
375 (fortran-to-lisp::integer4
) (fortran-to-lisp::integer4
)
376 (fortran-to-lisp::integer4
)
377 (fortran-to-lisp::complex16
)
378 (array fortran-to-lisp
::complex16
(*))
379 (fortran-to-lisp::integer4
)
380 (array fortran-to-lisp
::complex16
(*))
381 (fortran-to-lisp::integer4
)
382 (fortran-to-lisp::complex16
)
383 (array fortran-to-lisp
::complex16
(*))
384 (fortran-to-lisp::integer4
))
385 :return-values
'(nil nil nil nil nil nil nil nil nil nil nil nil
387 :calls
'(fortran-to-lisp::xerbla fortran-to-lisp
::lsame
))))