reference/dswapf.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

dswapf.f

dswapf.f 1.27 KB

subroutine dswapf (n,dx,incx,dy,incy)
c
c interchanges two vectors.
c uses unrolled loops for increments equal one.
c jack dongarra, linpack, 3/11/78.
c modified 12/3/93, array(1) declarations changed to array(*)
c
 double precision dx(*),dy(*),dtemp
 integer i,incx,incy,ix,iy,m,mp1,n
c
 if(n.le.0)return
 if(incx.eq.1.and.incy.eq.1)go to 20
c
c code for unequal increments or equal increments not equal
c to 1
c
 ix = 1
 iy = 1
 if(incx.lt.0)ix = (-n+1)*incx + 1
 if(incy.lt.0)iy = (-n+1)*incy + 1
 do 10 i = 1,n
 dtemp = dx(ix)
 dx(ix) = dy(iy)
 dy(iy) = dtemp
 ix = ix + incx
 iy = iy + incy
 10 continue
 return
c
c code for both increments equal to 1
c
c
c clean-up loop
c
 20 m = mod(n,3)
 if( m .eq. 0 ) go to 40
 do 30 i = 1,m
 dtemp = dx(i)
 dx(i) = dy(i)
 dy(i) = dtemp
 30 continue
 if( n .lt. 3 ) return
 40 mp1 = m + 1
 do 50 i = mp1,n,3
 dtemp = dx(i)
 dx(i) = dy(i)
 dy(i) = dtemp
 dtemp = dx(i + 1)
 dx(i + 1) = dy(i + 1)
 dy(i + 1) = dtemp
 dtemp = dx(i + 2)
 dx(i + 2) = dy(i + 2)
 dy(i + 2) = dtemp
 50 continue
 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

subroutinedswapf(n,dx,incx,dy,incy)
c
cinterchangestwovectors.
cusesunrolledloopsforincrementsequalone.
cjackdongarra,linpack,3/11/78.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
double precisiondx(*),dy(*),dtemp
integeri,incx,incy,ix,iy,m,mp1,n
c
if(n.le.0)return
if(incx.eq.1.and.incy.eq.1)go to20
c
ccodeforunequalincrementsorequalincrementsnotequal
cto1
c
ix=1
iy=1
if(incx.lt.0)ix=(-n+1)*incx+1
if(incy.lt.0)iy=(-n+1)*incy+1
do10i=1,n
dtemp=dx(ix)
dx(ix)=dy(iy)
dy(iy)=dtemp
ix=ix+incx
iy=iy+incy
10continue
return
c
ccodeforbothincrementsequalto1
c
c
cclean-uploop
c
20m=mod(n,3)
if(m.eq.0)go to40
do30i=1,m
dtemp=dx(i)
dx(i)=dy(i)
dy(i)=dtemp
30continue
if(n.lt.3)return
40mp1=m+1
do50i=mp1,n,3
dtemp=dx(i)
dx(i)=dy(i)
dy(i)=dtemp
dtemp=dx(i+1)
dx(i+1)=dy(i+1)
dy(i+1)=dtemp
dtemp=dx(i+2)
dx(i+2)=dy(i+2)
dy(i+2)=dtemp
50continue
return
end