reference/srotmgf.f · 三字经/OpenBLAS

代码拉取完成,页面将自动刷新

richgiteeai

Watch

不关注关注所有动态仅关注版本发行动态关注但不提醒动态

1 Star 0 Fork 0

三字经/OpenBLAS

代码 Issues 0 Pull Requests 0 Wiki 统计流水线

服务

加入 Gitee

与超过 1400万开发者一起发现、参与优秀开源项目,私有仓库也完全免费 :)

免费加入

已有帐号? 立即登录

文件

develop

分支 (57)

标签 (69)

管理

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

克隆/下载

HTTPS SSH SVN SVN+SSH 下载ZIP

提示

下载代码请复制以下命令到终端执行

为确保你提交的代码身份被 Gitee 正确识别,请执行以下命令完成配置

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

初次使用 SSH 协议进行代码克隆、推送等操作时,需按下述提示完成 SSH 配置

1 生成 RSA 密钥

2 获取 RSA 公钥内容,并配置到 SSH公钥中

在 Gitee 上使用 SVN,请访问使用指南

使用 HTTPS 协议时,命令行会出现如下账号密码验证步骤。基于安全考虑,Gitee 建议配置并使用私人令牌替代登录密码进行克隆、推送等操作

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

Password for 'https://userName@gitee.com': # 私人令牌

分支 57

标签 69

OpenBLAS

reference

srotmgf.f

srotmgf.f 4.64 KB

SUBROUTINE SROTMGF (SD1,SD2,SX1,SY1,SPARAM)
C
C CONSTRUCT THE MODIFIED GIVENS TRANSFORMATION MATRIX H WHICH ZEROS
C THE SECOND COMPONENT OF THE 2-VECTOR (SQRT(SD1)*SX1,SQRT(SD2)*
C SY2)**T.
C WITH SPARAM(1)=SFLAG, H HAS ONE OF THE FOLLOWING FORMS..
C
C SFLAG=-1.E0 SFLAG=0.E0 SFLAG=1.E0 SFLAG=-2.E0
C
C (SH11 SH12) (1.E0 SH12) (SH11 1.E0) (1.E0 0.E0)
C H=( ) ( ) ( ) ( )
C (SH21 SH22), (SH21 1.E0), (-1.E0 SH22), (0.E0 1.E0).
C LOCATIONS 2-4 OF SPARAM CONTAIN SH11,SH21,SH12, AND SH22
C RESPECTIVELY. (VALUES OF 1.E0, -1.E0, OR 0.E0 IMPLIED BY THE
C VALUE OF SPARAM(1) ARE NOT STORED IN SPARAM.)
C
C THE VALUES OF GAMSQ AND RGAMSQ SET IN THE DATA STATEMENT MAY BE
C INEXACT. THIS IS OK AS THEY ARE ONLY USED FOR TESTING THE SIZE
C OF SD1 AND SD2. ALL ACTUAL SCALING OF DATA IS DONE USING GAM.
C
 DIMENSION SPARAM(5)
C
 DATA ZERO,ONE,TWO /0.E0,1.E0,2.E0/
 DATA GAM,GAMSQ,RGAMSQ/4096.E0,1.67772E7,5.96046E-8/
 IF(.NOT. SD1 .LT. ZERO) GO TO 10
C GO ZERO-H-D-AND-SX1..
 GO TO 60
 10 CONTINUE
C CASE-SD1-NONNEGATIVE
 SP2=SD2*SY1
 IF(.NOT. SP2 .EQ. ZERO) GO TO 20
 SFLAG=-TWO
 GO TO 260
C REGULAR-CASE..
 20 CONTINUE
 SP1=SD1*SX1
 SQ2=SP2*SY1
 SQ1=SP1*SX1
C
 IF(.NOT. ABS(SQ1) .GT. ABS(SQ2)) GO TO 40
 SH21=-SY1/SX1
 SH12=SP2/SP1
C
 SU=ONE-SH12*SH21
C
 IF(.NOT. SU .LE. ZERO) GO TO 30
C GO ZERO-H-D-AND-SX1..
 GO TO 60
 30 CONTINUE
 SFLAG=ZERO
 SD1=SD1/SU
 SD2=SD2/SU
 SX1=SX1*SU
C GO SCALE-CHECK..
 GO TO 100
 40 CONTINUE
 IF(.NOT. SQ2 .LT. ZERO) GO TO 50
C GO ZERO-H-D-AND-SX1..
 GO TO 60
 50 CONTINUE
 SFLAG=ONE
 SH11=SP1/SP2
 SH22=SX1/SY1
 SU=ONE+SH11*SH22
 STEMP=SD2/SU
 SD2=SD1/SU
 SD1=STEMP
 SX1=SY1*SU
C GO SCALE-CHECK
 GO TO 100
C PROCEDURE..ZERO-H-D-AND-SX1..
 60 CONTINUE
 SFLAG=-ONE
 SH11=ZERO
 SH12=ZERO
 SH21=ZERO
 SH22=ZERO
C
 SD1=ZERO
 SD2=ZERO
 SX1=ZERO
C RETURN..
 GO TO 220
C PROCEDURE..FIX-H..
 70 CONTINUE
 IF(.NOT. SFLAG .GE. ZERO) GO TO 90
C
 IF(.NOT. SFLAG .EQ. ZERO) GO TO 80
 SH11=ONE
 SH22=ONE
 SFLAG=-ONE
 GO TO 90
 80 CONTINUE
 SH21=-ONE
 SH12=ONE
 SFLAG=-ONE
 90 CONTINUE
 GO TO IGO,(120,150,180,210)
C PROCEDURE..SCALE-CHECK
 100 CONTINUE
 110 CONTINUE
 IF(.NOT. SD1 .LE. RGAMSQ) GO TO 130
 IF(SD1 .EQ. ZERO) GO TO 160
 ASSIGN 120 TO IGO
C FIX-H..
 GO TO 70
 120 CONTINUE
 SD1=SD1*GAM**2
 SX1=SX1/GAM
 SH11=SH11/GAM
 SH12=SH12/GAM
 GO TO 110
 130 CONTINUE
 140 CONTINUE
 IF(.NOT. SD1 .GE. GAMSQ) GO TO 160
 ASSIGN 150 TO IGO
C FIX-H..
 GO TO 70
 150 CONTINUE
 SD1=SD1/GAM**2
 SX1=SX1*GAM
 SH11=SH11*GAM
 SH12=SH12*GAM
 GO TO 140
 160 CONTINUE
 170 CONTINUE
 IF(.NOT. ABS(SD2) .LE. RGAMSQ) GO TO 190
 IF(SD2 .EQ. ZERO) GO TO 220
 ASSIGN 180 TO IGO
C FIX-H..
 GO TO 70
 180 CONTINUE
 SD2=SD2*GAM**2
 SH21=SH21/GAM
 SH22=SH22/GAM
 GO TO 170
 190 CONTINUE
 200 CONTINUE
 IF(.NOT. ABS(SD2) .GE. GAMSQ) GO TO 220
 ASSIGN 210 TO IGO
C FIX-H..
 GO TO 70
 210 CONTINUE
 SD2=SD2/GAM**2
 SH21=SH21*GAM
 SH22=SH22*GAM
 GO TO 200
 220 CONTINUE
 IF(SFLAG)250,230,240
 230 CONTINUE
 SPARAM(3)=SH21
 SPARAM(4)=SH12
 GO TO 260
 240 CONTINUE
 SPARAM(2)=SH11
 SPARAM(5)=SH22
 GO TO 260
 250 CONTINUE
 SPARAM(2)=SH11
 SPARAM(3)=SH21
 SPARAM(4)=SH12
 SPARAM(5)=SH22
 260 CONTINUE
 SPARAM(1)=SFLAG
 RETURN
 END

一键复制编辑原始数据按行查看历史

Zhang Xianyi 提交于 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 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166

SUBROUTINESROTMGF(SD1,SD2,SX1,SY1,SPARAM)
C
CCONSTRUCTTHEMODIFIEDGIVENSTRANSFORMATIONMATRIXHWHICHZEROS
CTHESECONDCOMPONENTOFTHE2-VECTOR(SQRT(SD1)*SX1,SQRT(SD2)*
CSY2)**T.
CWITHSPARAM(1)=SFLAG,HHASONEOFTHEFOLLOWINGFORMS..
C
CSFLAG=-1.E0SFLAG=0.E0SFLAG=1.E0SFLAG=-2.E0
C
C(SH11SH12)(1.E0SH12)(SH111.E0)(1.E00.E0)
CH=()()()()
C(SH21SH22),(SH211.E0),(-1.E0SH22),(0.E01.E0).
CLOCATIONS2-4OFSPARAMCONTAINSH11,SH21,SH12,ANDSH22
CRESPECTIVELY.(VALUESOF1.E0,-1.E0,OR0.E0IMPLIEDBYTHE
CVALUEOFSPARAM(1)ARENOTSTOREDINSPARAM.)
C
CTHEVALUESOFGAMSQANDRGAMSQSETINTHEDATASTATEMENTMAYBE
CINEXACT.THISISOKASTHEYAREONLYUSEDFORTESTINGTHESIZE
COFSD1ANDSD2.ALLACTUALSCALINGOFDATAISDONEUSINGGAM.
C
DIMENSIONSPARAM(5)
C
DATAZERO,ONE,TWO/0.E0,1.E0,2.E0/
DATAGAM,GAMSQ,RGAMSQ/4096.E0,1.67772E7,5.96046E-8/
IF(.NOT.SD1.LT.ZERO)GO TO10
CGOZERO-H-D-AND-SX1..
GO TO60
10CONTINUE
CCASE-SD1-NONNEGATIVE
SP2=SD2*SY1
IF(.NOT.SP2.EQ.ZERO)GO TO20
SFLAG=-TWO
GO TO260
CREGULAR-CASE..
20CONTINUE
SP1=SD1*SX1
SQ2=SP2*SY1
SQ1=SP1*SX1
C
IF(.NOT.ABS(SQ1).GT.ABS(SQ2))GO TO40
SH21=-SY1/SX1
SH12=SP2/SP1
C
SU=ONE-SH12*SH21
C
IF(.NOT.SU.LE.ZERO)GO TO30
CGOZERO-H-D-AND-SX1..
GO TO60
30CONTINUE
SFLAG=ZERO
SD1=SD1/SU
SD2=SD2/SU
SX1=SX1*SU
CGOSCALE-CHECK..
GO TO100
40CONTINUE
IF(.NOT.SQ2.LT.ZERO)GO TO50
CGOZERO-H-D-AND-SX1..
GO TO60
50CONTINUE
SFLAG=ONE
SH11=SP1/SP2
SH22=SX1/SY1
SU=ONE+SH11*SH22
STEMP=SD2/SU
SD2=SD1/SU
SD1=STEMP
SX1=SY1*SU
CGOSCALE-CHECK
GO TO100
CPROCEDURE..ZERO-H-D-AND-SX1..
60CONTINUE
SFLAG=-ONE
SH11=ZERO
SH12=ZERO
SH21=ZERO
SH22=ZERO
C
SD1=ZERO
SD2=ZERO
SX1=ZERO
CRETURN..
GO TO220
CPROCEDURE..FIX-H..
70CONTINUE
IF(.NOT.SFLAG.GE.ZERO)GO TO90
C
IF(.NOT.SFLAG.EQ.ZERO)GO TO80
SH11=ONE
SH22=ONE
SFLAG=-ONE
GO TO90
80CONTINUE
SH21=-ONE
SH12=ONE
SFLAG=-ONE
90CONTINUE
GO TOIGO,(120,150,180,210)
CPROCEDURE..SCALE-CHECK
100CONTINUE
110CONTINUE
IF(.NOT.SD1.LE.RGAMSQ)GO TO130
IF(SD1.EQ.ZERO)GO TO160
ASSIGN120TOIGO
CFIX-H..
GO TO70
120CONTINUE
SD1=SD1*GAM**2
SX1=SX1/GAM
SH11=SH11/GAM
SH12=SH12/GAM
GO TO110
130CONTINUE
140CONTINUE
IF(.NOT.SD1.GE.GAMSQ)GO TO160
ASSIGN150TOIGO
CFIX-H..
GO TO70
150CONTINUE
SD1=SD1/GAM**2
SX1=SX1*GAM
SH11=SH11*GAM
SH12=SH12*GAM
GO TO140
160CONTINUE
170CONTINUE
IF(.NOT.ABS(SD2).LE.RGAMSQ)GO TO190
IF(SD2.EQ.ZERO)GO TO220
ASSIGN180TOIGO
CFIX-H..
GO TO70
180CONTINUE
SD2=SD2*GAM**2
SH21=SH21/GAM
SH22=SH22/GAM
GO TO170
190CONTINUE
200CONTINUE
IF(.NOT.ABS(SD2).GE.GAMSQ)GO TO220
ASSIGN210TOIGO
CFIX-H..
GO TO70
210CONTINUE
SD2=SD2/GAM**2
SH21=SH21*GAM
SH22=SH22*GAM
GO TO200
220CONTINUE
IF(SFLAG)250,230,240
230CONTINUE
SPARAM(3)=SH21
SPARAM(4)=SH12
GO TO260
240CONTINUE
SPARAM(2)=SH11
SPARAM(5)=SH22
GO TO260
250CONTINUE
SPARAM(2)=SH11
SPARAM(3)=SH21
SPARAM(4)=SH12
SPARAM(5)=SH22
260CONTINUE
SPARAM(1)=SFLAG
RETURN
END