NLSL: zaypx.f90 Source File

NLSL

- Main Page
- Modules
- Data Types List
- Files

- File List
- File Members

zaypx.f90
Go to the documentation of this file.
1 c Version 1.4 10/10/94
2 c**********************************************************************
3 c
4 c complex double precision vector scale and add
5 c ---------------------------------------------
6 c
7 c This subroutine will update a complex double precision
8 c vector Y by the formula,
9 c
10 c Y=x+a*Y ,
11 c
12 c where a is a complex double precision constant and X is a
13 c complex double presion vector.
14 c
15 c
16 c**********************************************************************
17 c
18 subroutine zaypx(x,y,scale,ndim)
19 c
20 use nlsdim
21 c*djs use rnddbl
22 c
23 integer ndim
24 double complex scale
25 c
26 double complex x,y
27 dimension x(mxdim),y(mxdim)
28 c
29 integer iel
30 c
31 c*djs double complex CZERO
32 c*djs parameter (CZERO=(0.0D0,0.0D0))
33 c
34 c######################################################################
35 c
36 do iel=1,ndim
37 y(iel)=x(iel)+scale*y(iel)
38 c*djs if (abs(y(iel)).lt.RNDOFF) y(iel)=CZERO
39 end do
40 c
41 return
42 end
zaypx
subroutine zaypx(x, y, scale, ndim)
Definition: zaypx.f90:19
nlsdim
Definition: nlsdim.f90:36
nlsdim::mxdim
integer, parameter mxdim
Definition: nlsdim.f90:39