daxpy Subroutine

public subroutine daxpy(N, Da, Dx, Incx, Dy, Incy)

NAME

daxpy(3f) - [M_odepack::matrix] Compute a constant times a vector plus a vector.

SYNOPSIS

subroutine daxpy(N,Da,Dx,Incx,Dy,Incy)
integer,intent(in)          :: N
real(kind=dp),intent(in)    :: Da
real(kind=dp),intent(in)    :: Dx(*)
integer,intent(in)          :: Incx
real(kind=dp),intent(inout) :: Dy(*)
integer,intent(in)          :: Incy

DESCRIPTION

daxpy(3f) computes a constant times a vector plus a vector. It uses unrolled loops for increments equal to one.

Overwrite double precision DY with double precision DADX + DY. For I = 0 to N-1, replace DY(LY+IINCY) with DADX(LX+IINCX) + DY(LY+I*INCY),

where LX = 1 if INCX .GE. 0, else LX = 1+(1-N)*INCX, and LY is defined in a similar way using INCY.

INPUT OPTIONS

N

number of elements in input vector(s)

DA

double precision scalar multiplier

DX

double precision vector with N elements

INCX

storage spacing between elements of DX

DY

double precision vector with N elements

INCY

storage spacing between elements of DY

RETURNS

DY

double precision result (unchanged if N .LE. 0)

REFERENCES

B L A S Subprogram

C. L. Lawson, R. J. Hanson, D. R. Kincaid and F. T. Krogh, Basic linear algebra subprograms for Fortran usage, Algorithm No. 539, Transactions on Mathematical Software 5, 3 (September 1979), pp. 308-323.

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: N
real(kind=dp), intent(in) :: Da
real(kind=dp), intent(in) :: Dx(*)
integer, intent(in) :: Incx
real(kind=dp), intent(inout) :: Dy(*)
integer, intent(in) :: Incy

Called by

proc~~daxpy~2~~CalledByGraph proc~daxpy~2 M_odepack::daxpy proc~dainvg~2 M_odepack::dainvg proc~dainvg~2->proc~daxpy~2 proc~dgefa~2 M_odepack::dgefa proc~dainvg~2->proc~dgefa~2 proc~dgesl M_odepack::dgesl proc~dainvg~2->proc~dgesl proc~dgefa~2->proc~daxpy~2 proc~dgesl->proc~daxpy~2 proc~dprepji~2 M_odepack::dprepji proc~dprepji~2->proc~daxpy~2 proc~dprepji~2->proc~dgefa~2 proc~dprepj~2 M_odepack::dprepj proc~dprepj~2->proc~daxpy~2 proc~dprepj~2->proc~dgefa~2 proc~dprja~2 M_odepack::dprja proc~dprja~2->proc~daxpy~2 proc~dprja~2->proc~dgefa~2 proc~dsolsy~2 M_odepack::dsolsy proc~dsolsy~2->proc~daxpy~2 proc~dsolsy~2->proc~dgesl proc~dstodpk M_odepack::dstodpk proc~dstodpk->proc~daxpy~2 proc~dstoka~2 M_odepack::dstoka proc~dstoka~2->proc~daxpy~2 proc~daigbt~2 M_odepack::daigbt proc~daigbt~2->proc~dgefa~2 proc~daigbt~2->proc~dgesl proc~dlsodi~2 M_odepack::dlsodi proc~dlsodi~2->proc~dainvg~2 proc~dlsodkr~2 M_odepack::dlsodkr proc~dlsodkr~2->proc~dstoka~2 proc~dlsodpk~2 M_odepack::dlsodpk proc~dlsodpk~2->proc~dstodpk proc~dpjibt~2 M_odepack::dpjibt proc~dpjibt~2->proc~dgefa~2 proc~dpjibt~2->proc~dgesl proc~dslsbt~2 M_odepack::dslsbt proc~dslsbt~2->proc~dgesl proc~dlsoibt~2 M_odepack::dlsoibt proc~dlsoibt~2->proc~daigbt~2