16 SUBROUTINE daxpby (N,DZ,DA,DX,DB,DY)
17 USE stel_kinds,
ONLY: rprec, dp
55 REAL(rprec),
PARAMETER :: ZERO = 0, one = 1
56 REAL(rprec) DA, DB, DX(N), DY(N), DZ(N)
70 ELSE IF (db.EQ.one)
THEN
75 ELSE IF (db.EQ.-one)
THEN
86 ELSE IF (da.EQ.one)
THEN
92 ELSE IF (db.EQ.one)
THEN
97 ELSE IF (db.EQ.-one)
THEN
100 dz(i) = dx(i) - dy(i)
105 dz(i) = dx(i) + db * dy(i)
108 ELSE IF (da.EQ.-one)
THEN
114 ELSE IF (db.EQ.one)
THEN
117 dz(i) = -dx(i) + dy(i)
119 ELSE IF (db.EQ.-one)
THEN
122 dz(i) = -dx(i) - dy(i)
127 dz(i) = -dx(i) + db * dy(i)
136 ELSE IF (db.EQ.one)
THEN
139 dz(i) = da * dx(i) + dy(i)
141 ELSE IF (db.EQ.-one)
THEN
144 dz(i) = da * dx(i) - dy(i)
149 dz(i) = da * dx(i) + db * dy(i)