1DROTM(1) BLAS routine DROTM(1)
2
3
4
6 DROTM - APPLY THE MODIFIED GIVENS TRANSFORMATION, H, TO THE 2 BY N
7 MATRIX (DX**T) , WHERE **T INDICATES TRANSPOSE
8
10 SUBROUTINE DROTM(N,DX,INCX,DY,INCY,DPARAM)
11
12 INTEGER INCX,INCY,N
13
14 DOUBLE PRECISION
15 DPARAM(5),DX(1),DY(1)
16
18 APPLY THE MODIFIED GIVENS TRANSFORMATION, H, TO THE 2 BY N MATRIX
19 (DY**T)
20
21 DX(LX+I*INCX), I = 0 TO N-1, WHERE LX = 1 IF INCX .GE. 0, ELSE
22 LX = (-INCX)*N, AND SIMILARLY FOR SY USING LY AND INCY.
23 WITH DPARAM(1)=DFLAG, H HAS ONE OF THE FOLLOWING FORMS..
24
25 DFLAG=-1.D0 DFLAG=0.D0 DFLAG=1.D0 DFLAG=-2.D0
26
27 (DH11 DH12) (1.D0 DH12) (DH11 1.D0) (1.D0 0.D0)
28 H=( ) ( ) ( ) ( )
29 (DH21 DH22), (DH21 1.D0), (-1.D0 DH22), (0.D0 1.D0).
30 SEE DROTMG FOR A DESCRIPTION OF DATA STORAGE IN DPARAM.
31
32
34 N (input) INTEGER
35 number of elements in input vector(s)
36
37 DX (input/output) DOUBLE PRECISION array, dimension N
38 double precision vector with 5 elements
39
40 INCX (input) INTEGER
41 storage spacing between elements of DX
42
43 DY (input/output) DOUBLE PRECISION array, dimension N
44 double precision vector with N elements
45
46 INCY (input) INTEGER
47 storage spacing between elements of DY
48
49 DPARAM (input/output) DOUBLE PRECISION array, dimension 5
50 DPARAM(1)=DFLAG DPARAM(2)=DH11 DPARAM(3)=DH21 DPARAM(4)=DH12
51 DPARAM(5)=DH22
52
53
54
55BLAS routine November 2008 DROTM(1)