SUBROUTINE COPY ( N, DX, IX, DY, IY )

Argument Definitions (+ indicates altered content)
INTEGER            IX,         IY,         N
REAL*8             DX(*),     +DY(*)
Description
                                                                        * 
  Copyright (C) 1996 Leif Laaksonen, Dage Sundholm                      * 
  Copyright (C) 1996-2010 Jacek Kobus               * 
                                                                        * 
  This program is free software; you can redistribute it and/or modify  * 
  it under the terms of the GNU General Public License version 2 as     * 
  published by the Free Software Foundation.                            * 
                                                                        * 
 
### blas ### 
 
    This is simplified replacement for BLAS routines: dcopy, daxpy, dscal, 
    ddot, mxv. Note that ix and iy are ignored and are set to 1.
Source file:blas.f
Local Variables (+ indicates altered content)
INTEGER           +I

SUBROUTINE AXPY ( N, DA, DX, IX, DY, IY )

Argument Definitions (+ indicates altered content)
INTEGER            IX,         IY,         N
REAL*8             DA,         DX(*),     +DY(*)
Source file:blas.f
Local Variables (+ indicates altered content)
INTEGER           +I

SUBROUTINE SCAL ( N, DA, DX, IX )

Argument Definitions (+ indicates altered content)
INTEGER            IX,         N
REAL*8             DA,        +DX(*)
Source file:blas.f
Local Variables (+ indicates altered content)
INTEGER           +I

REAL FUNCTION DOT ( N, DX, IX, DY, IY )

Argument Definitions (+ indicates altered content)
INTEGER            IX,         IY,         N
REAL*8             DX(*),      DY(*)
Source file:blas.f
Local Variables (+ indicates altered content)
INTEGER           +I

SUBROUTINE GEMV ( NR1, NC, DX, DV, DVR )

Argument Definitions (+ indicates altered content)
INTEGER            NC,         NR1
REAL*8             DV(*),     +DVR(*),     DX(NR1,*)
Description
    Multiplies the matrix DX(nr,nc) times the vector DV(nc) and stores 
    the result in the vector DY(nr) by means of calling DGEMV 
 
    subroutine dgemv (trans,nr1,nc,alpha,dx,nr2,dv,incr1,beta, 
   &     dvr,incr2) 
    integer*4  incr1,incr2,inr,nc,nr1,nr2 
    real*8  alpha,beta,dx(nr1,nc),dv(nc),dvr(nr1),s 
    character*1 trans 
 
    do inr=1,nr1 
       s=0.d0 
       do inc=1,nc 
          s=s+dx(inr,inc)*dv(inc) 
       enddo 
       dvr(inr)=s 
    enddo 
 
    return 
    end 
 
Multiplies the matrix DX(nr,nc) times the vector DV(nc) and stores 
the result in the vector DY(nr) (simplified version of DGEMV)
Source file:blas.f
Local Variables (+ indicates altered content)
INTEGER           +INC,       +INR
REAL*8            +S