reference/drotmf.f · 三字经/OpenBLAS

Fetch the repository succeeded.

richgiteeai

Watch

Unwatch Watching Releases Only Ignoring

1 Star 0 Fork 0

三字经/OpenBLAS

Code Issues 0 Pull Requests 0 Wiki Insights Pipelines

Service

Tencent Cloud Serverless

Create your Gitee Account

Explore and code with more than 14 million developers,Free private repositories !:)

Already have an account? Sign in

文件

develop

Branches (57)

Tags (69)

Manage

develop

gh-pages

release-0.3.0

issue5267

revert-5565-jn/makefile-rule-dynamic

revert-5643-neov2_param

revert-5479-forklock

revert-4741-Pthread_Scalability_Improvement

gemm_batch

revert-5180-openmp_use_cmake

revert-5251-issue5250

timeout-riscv-ci

revert-4920-issue4917

revert-4628-CI-c910v-mips64-loongarch64

risc-v

revert-4153-dgemv

jenkins

revert-4103-windows_perf

utest-riscv

revert-3518-elbrus

v0.3.33

v0.3.32

v0.3.31

v0.3.30

v0.3.29

v0.3.28

v0.3.27

v0.3.26

v0.3.25

v0.3.24

v0.3.23

v0.3.22

v0.3.21

v0.3.20

v0.3.19

v0.3.18

v0.3.17

v0.3.16

v0.3.15

v0.3.14

Clone or Download

Clone/Download

HTTPS SSH SVN SVN+SSH Download ZIP

Prompt

To download the code, please copy the following command and execute it in the terminal

To ensure that your submitted code identity is correctly recognized by Gitee, please execute the following command.

git config --global user.name userName 
git config --global user.email userEmail

When using the SSH protocol for the first time to clone or push code, follow the prompts below to complete the SSH configuration.

1 Generate RSA keys.

2 Obtain the content of the RSA public key and configure it in SSH Public Keys

To use SVN on Gitee, please visit the usage guide

When using the HTTPS protocol, the command line will prompt for account and password verification as follows. For security reasons, Gitee recommends configure and use personal access tokens instead of login passwords for cloning, pushing, and other operations.

Username for 'https://gitee.com': userName

Password for 'https://userName@gitee.com': # Private Token

Branches 57

Tags 69

OpenBLAS

reference

drotmf.f

drotmf.f 3.18 KB

SUBROUTINE DROTMF (N,DX,INCX,DY,INCY,DPARAM)
C
C APPLY THE MODIFIED GIVENS TRANSFORMATION, H, TO THE 2 BY N MATRIX
C
C (DX**T) , WHERE **T INDICATES TRANSPOSE. THE ELEMENTS OF DX ARE IN
C (DY**T)
C
C DX(LX+I*INCX), I = 0 TO N-1, WHERE LX = 1 IF INCX .GE. 0, ELSE
C LX = (-INCX)*N, AND SIMILARLY FOR SY USING LY AND INCY.
C WITH DPARAM(1)=DFLAG, H HAS ONE OF THE FOLLOWING FORMS..
C
C DFLAG=-1.D0 DFLAG=0.D0 DFLAG=1.D0 DFLAG=-2.D0
C
C (DH11 DH12) (1.D0 DH12) (DH11 1.D0) (1.D0 0.D0)
C H=( ) ( ) ( ) ( )
C (DH21 DH22), (DH21 1.D0), (-1.D0 DH22), (0.D0 1.D0).
C SEE DROTMG FOR A DESCRIPTION OF DATA STORAGE IN DPARAM.
C
 DOUBLE PRECISION DFLAG,DH12,DH22,DX,TWO,Z,DH11,DH21,
 1 DPARAM,DY,W,ZERO
 DIMENSION DX(1),DY(1),DPARAM(5)
 DATA ZERO,TWO/0.D0,2.D0/
C
 DFLAG=DPARAM(1)
 IF(N .LE. 0 .OR.(DFLAG+TWO.EQ.ZERO)) GO TO 140
 IF(.NOT.(INCX.EQ.INCY.AND. INCX .GT.0)) GO TO 70
C
 NSTEPS=N*INCX
 IF(DFLAG) 50,10,30
 10 CONTINUE
 DH12=DPARAM(4)
 DH21=DPARAM(3)
 DO 20 I=1,NSTEPS,INCX
 W=DX(I)
 Z=DY(I)
 DX(I)=W+Z*DH12
 DY(I)=W*DH21+Z
 20 CONTINUE
 GO TO 140
 30 CONTINUE
 DH11=DPARAM(2)
 DH22=DPARAM(5)
 DO 40 I=1,NSTEPS,INCX
 W=DX(I)
 Z=DY(I)
 DX(I)=W*DH11+Z
 DY(I)=-W+DH22*Z
 40 CONTINUE
 GO TO 140
 50 CONTINUE
 DH11=DPARAM(2)
 DH12=DPARAM(4)
 DH21=DPARAM(3)
 DH22=DPARAM(5)
 DO 60 I=1,NSTEPS,INCX
 W=DX(I)
 Z=DY(I)
 DX(I)=W*DH11+Z*DH12
 DY(I)=W*DH21+Z*DH22
 60 CONTINUE
 GO TO 140
 70 CONTINUE
 KX=1
 KY=1
 IF(INCX .LT. 0) KX=1+(1-N)*INCX
 IF(INCY .LT. 0) KY=1+(1-N)*INCY
C
 IF(DFLAG)120,80,100
 80 CONTINUE
 DH12=DPARAM(4)
 DH21=DPARAM(3)
 DO 90 I=1,N
 W=DX(KX)
 Z=DY(KY)
 DX(KX)=W+Z*DH12
 DY(KY)=W*DH21+Z
 KX=KX+INCX
 KY=KY+INCY
 90 CONTINUE
 GO TO 140
 100 CONTINUE
 DH11=DPARAM(2)
 DH22=DPARAM(5)
 DO 110 I=1,N
 W=DX(KX)
 Z=DY(KY)
 DX(KX)=W*DH11+Z
 DY(KY)=-W+DH22*Z
 KX=KX+INCX
 KY=KY+INCY
 110 CONTINUE
 GO TO 140
 120 CONTINUE
 DH11=DPARAM(2)
 DH12=DPARAM(4)
 DH21=DPARAM(3)
 DH22=DPARAM(5)
 DO 130 I=1,N
 W=DX(KX)
 Z=DY(KY)
 DX(KX)=W*DH11+Z*DH12
 DY(KY)=W*DH21+Z*DH22
 KX=KX+INCX
 KY=KY+INCY
 130 CONTINUE
 140 CONTINUE
 RETURN
 END

Copy Edit Raw Blame History

Zhang Xianyi authored 2011年01月24日 22:54 +08:00 . Import GotoBLAS2 1.13 BSD version codes.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108

SUBROUTINEDROTMF(N,DX,INCX,DY,INCY,DPARAM)
C
CAPPLYTHEMODIFIEDGIVENSTRANSFORMATION,H,TOTHE2BYNMATRIX
C
C(DX**T),WHERE**TINDICATESTRANSPOSE.THEELEMENTSOFDXAREIN
C(DY**T)
C
CDX(LX+I*INCX),I=0TON-1,WHERELX=1IFINCX.GE.0,ELSE
CLX=(-INCX)*N,ANDSIMILARLYFORSYUSINGLYANDINCY.
CWITHDPARAM(1)=DFLAG,HHASONEOFTHEFOLLOWINGFORMS..
C
CDFLAG=-1.D0DFLAG=0.D0DFLAG=1.D0DFLAG=-2.D0
C
C(DH11DH12)(1.D0DH12)(DH111.D0)(1.D00.D0)
CH=()()()()
C(DH21DH22),(DH211.D0),(-1.D0DH22),(0.D01.D0).
CSEEDROTMGFORADESCRIPTIONOFDATASTORAGEINDPARAM.
C
DOUBLE PRECISIONDFLAG,DH12,DH22,DX,TWO,Z,DH11,DH21,
1DPARAM,DY,W,ZERO
DIMENSIONDX(1),DY(1),DPARAM(5)
DATAZERO,TWO/0.D0,2.D0/
C
DFLAG=DPARAM(1)
IF(N.LE.0.OR.(DFLAG+TWO.EQ.ZERO))GO TO140
IF(.NOT.(INCX.EQ.INCY.AND.INCX.GT.0))GO TO70
C
NSTEPS=N*INCX
IF(DFLAG)50,10,30
10CONTINUE
DH12=DPARAM(4)
DH21=DPARAM(3)
DO20I=1,NSTEPS,INCX
W=DX(I)
Z=DY(I)
DX(I)=W+Z*DH12
DY(I)=W*DH21+Z
20CONTINUE
GO TO140
30CONTINUE
DH11=DPARAM(2)
DH22=DPARAM(5)
DO40I=1,NSTEPS,INCX
W=DX(I)
Z=DY(I)
DX(I)=W*DH11+Z
DY(I)=-W+DH22*Z
40CONTINUE
GO TO140
50CONTINUE
DH11=DPARAM(2)
DH12=DPARAM(4)
DH21=DPARAM(3)
DH22=DPARAM(5)
DO60I=1,NSTEPS,INCX
W=DX(I)
Z=DY(I)
DX(I)=W*DH11+Z*DH12
DY(I)=W*DH21+Z*DH22
60CONTINUE
GO TO140
70CONTINUE
KX=1
KY=1
IF(INCX.LT.0)KX=1+(1-N)*INCX
IF(INCY.LT.0)KY=1+(1-N)*INCY
C
IF(DFLAG)120,80,100
80CONTINUE
DH12=DPARAM(4)
DH21=DPARAM(3)
DO90I=1,N
W=DX(KX)
Z=DY(KY)
DX(KX)=W+Z*DH12
DY(KY)=W*DH21+Z
KX=KX+INCX
KY=KY+INCY
90CONTINUE
GO TO140
100CONTINUE
DH11=DPARAM(2)
DH22=DPARAM(5)
DO110I=1,N
W=DX(KX)
Z=DY(KY)
DX(KX)=W*DH11+Z
DY(KY)=-W+DH22*Z
KX=KX+INCX
KY=KY+INCY
110CONTINUE
GO TO140
120CONTINUE
DH11=DPARAM(2)
DH12=DPARAM(4)
DH21=DPARAM(3)
DH22=DPARAM(5)
DO130I=1,N
W=DX(KX)
Z=DY(KY)
DX(KX)=W*DH11+Z*DH12
DY(KY)=W*DH21+Z*DH22
KX=KX+INCX
KY=KY+INCY
130CONTINUE
140CONTINUE
RETURN
END