Mercurial > hg > octave-lyh
changeset 10508:9a5e2d13fa5a
fix typos in blas-xtra
author | Jaroslav Hajek <highegg@gmail.com> |
---|---|
date | Sun, 11 Apr 2010 14:33:00 +0200 |
parents | 424795ef82b8 |
children | ddbd812d09aa |
files | libcruft/ChangeLog libcruft/blas-xtra/cmatm3.f libcruft/blas-xtra/dmatm3.f libcruft/blas-xtra/smatm3.f libcruft/blas-xtra/xsdot.f libcruft/blas-xtra/zmatm3.f |
diffstat | 6 files changed, 14 insertions(+), 4 deletions(-) [+] |
line wrap: on
line diff
--- a/libcruft/ChangeLog +++ b/libcruft/ChangeLog @@ -1,3 +1,9 @@ +2010-04-11 Jaroslav Hajek <highegg@gmail.com> + + * blas-xtra/cmatm3.f, blas-xtra/zmatm3.f, + blas-xtra/dmatm3.f, blas-xtra/xsdot.f: Fix typos. Add missing + decls. + 2010-03-21 John W. Eaton <jwe@octave.org> * Makefile.am (libcruft_la_LIBADD): Include ../libgnu/libgnu.la
--- a/libcruft/blas-xtra/cmatm3.f +++ b/libcruft/blas-xtra/cmatm3.f @@ -35,9 +35,10 @@ complex a(m*k,np),b(k*n,np) complex c(m*n,np) - complex cdot,one,zero + complex cdotu,one,zero parameter (one = 1e0, zero = 0e0) external cdotu,cgemv,cgemm + integer i c quick return if possible. if (np <= 0) return
--- a/libcruft/blas-xtra/dmatm3.f +++ b/libcruft/blas-xtra/dmatm3.f @@ -35,9 +35,10 @@ double precision a(m*k,np),b(k*n,np) double precision c(m*n,np) - double precision sdot,one,zero + double precision ddot,one,zero parameter (one = 1d0, zero = 0d0) external ddot,dgemv,dgemm + integer i c quick return if possible. if (np <= 0) return
--- a/libcruft/blas-xtra/smatm3.f +++ b/libcruft/blas-xtra/smatm3.f @@ -38,6 +38,7 @@ real sdot,one,zero parameter (one = 1e0, zero = 0e0) external sdot,sgemv,sgemm + integer i c quick return if possible. if (np <= 0) return
--- a/libcruft/blas-xtra/xsdot.f +++ b/libcruft/blas-xtra/xsdot.f @@ -1,5 +1,5 @@ subroutine xsdot (n, dx, incx, dy, incy, retval) - real ddot, dx(*), dy(*), retval + real ddot, dx(*), dy(*), retval, sdot integer n, incx, incy retval = sdot (n, dx, incx, dy, incy) return
--- a/libcruft/blas-xtra/zmatm3.f +++ b/libcruft/blas-xtra/zmatm3.f @@ -35,9 +35,10 @@ double complex a(m*k,np),b(k*n,np) double complex c(m*n,np) - double complex cdot,one,zero + double complex zdotu,one,zero parameter (one = 1d0, zero = 0d0) external zdotu,zgemv,zgemm + integer i c quick return if possible. if (np <= 0) return