summaryrefslogtreecommitdiff
path: root/src/lib/blas/drotg.f
diff options
context:
space:
mode:
authorjofret2009-04-28 07:17:27 +0000
committerjofret2009-04-28 07:17:27 +0000
commitb4291825aaa210d693f4f858e1c8c25f514a3f71 (patch)
tree4ba928deac68d7098a7043cd286acb23351825d8 /src/lib/blas/drotg.f
parent8c8d2f518968ce7057eec6aa5cd5aec8faab861a (diff)
downloadscilab2c-b4291825aaa210d693f4f858e1c8c25f514a3f71.tar.gz
scilab2c-b4291825aaa210d693f4f858e1c8c25f514a3f71.tar.bz2
scilab2c-b4291825aaa210d693f4f858e1c8c25f514a3f71.zip
Moving blas to right place
Diffstat (limited to 'src/lib/blas/drotg.f')
-rw-r--r--src/lib/blas/drotg.f27
1 files changed, 0 insertions, 27 deletions
diff --git a/src/lib/blas/drotg.f b/src/lib/blas/drotg.f
deleted file mode 100644
index 67838e2c..00000000
--- a/src/lib/blas/drotg.f
+++ /dev/null
@@ -1,27 +0,0 @@
- subroutine drotg(da,db,c,s)
-c
-c construct givens plane rotation.
-c jack dongarra, linpack, 3/11/78.
-c
- double precision da,db,c,s,roe,scale,r,z
-c
- roe = db
- if( dabs(da) .gt. dabs(db) ) roe = da
- scale = dabs(da) + dabs(db)
- if( scale .ne. 0.0d0 ) go to 10
- c = 1.0d0
- s = 0.0d0
- r = 0.0d0
- z = 0.0d0
- go to 20
- 10 r = scale*dsqrt((da/scale)**2 + (db/scale)**2)
- r = dsign(1.0d0,roe)*r
- c = da/r
- s = db/r
- z = 1.0d0
- if( dabs(da) .gt. dabs(db) ) z = s
- if( dabs(db) .ge. dabs(da) .and. c .ne. 0.0d0 ) z = 1.0d0/c
- 20 da = r
- db = z
- return
- end