← Back to branch summary

~ubuntu-branches/ubuntu/wily/julia/wily

« back to all changes in this revision

Viewing changes to deps/openlibm/slatec/srot.f

  • Committer: Package Import Robot
  • Author(s): Sébastien Villemot
  • Date: 2013-01-16 12:29:42 UTC
  • mfrom: (1.1.2)
  • Revision ID: package-import@ubuntu.com-20130116122942-x86e42akjq31repw
Tags: 0.0.0+20130107.gitd9656f41-1
* New upstream snashot
* No longer try to rebuild helpdb.jl.
   + debian/rules: remove helpdb.jl from build-arch rule
   + debian/control: move back python-sphinx to Build-Depends-Indep
* debian/copyright: reflect upstream changes
* Add Build-Conflicts on libatlas3-base (makes linalg tests fail)
* debian/rules: replace obsolete USE_DEBIAN makeflag by a list of
  USE_SYSTEM_* flags
* debian/rules: on non-x86 systems, use libm instead of openlibm
* dpkg-buildflags.patch: remove patch, applied upstream
* Refreshed other patches

Show diffs side-by-side

added added

removed removed

Lines of Context:
deps/openlibm/slatec/srot.f
 
1
*DECK SROT
 
2
      SUBROUTINE SROT (N, SX, INCX, SY, INCY, SC, SS)
 
3
C***BEGIN PROLOGUE  SROT
 
4
C***PURPOSE  Apply a plane Givens rotation.
 
5
C***LIBRARY   SLATEC (BLAS)
 
6
C***CATEGORY  D1A8
 
7
C***TYPE      SINGLE PRECISION (SROT-S, DROT-D, CSROT-C)
 
8
C***KEYWORDS  BLAS, GIVENS ROTATION, GIVENS TRANSFORMATION,
 
9
C             LINEAR ALGEBRA, PLANE ROTATION, VECTOR
 
10
C***AUTHOR  Lawson, C. L., (JPL)
 
11
C           Hanson, R. J., (SNLA)
 
12
C           Kincaid, D. R., (U. of Texas)
 
13
C           Krogh, F. T., (JPL)
 
14
C***DESCRIPTION
 
15
C
 
16
C                B L A S  Subprogram
 
17
C    Description of Parameters
 
18
C
 
19
C     --Input--
 
20
C        N  number of elements in input vector(s)
 
21
C       SX  single precision vector with N elements
 
22
C     INCX  storage spacing between elements of SX
 
23
C       SY  single precision vector with N elements
 
24
C     INCY  storage spacing between elements of SY
 
25
C       SC  element of rotation matrix
 
26
C       SS  element of rotation matrix
 
27
C
 
28
C     --Output--
 
29
C       SX  rotated vector SX (unchanged if N .LE. 0)
 
30
C       SY  rotated vector SY (unchanged if N .LE. 0)
 
31
C
 
32
C     Multiply the 2 x 2 matrix  ( SC SS) times the 2 x N matrix (SX**T)
 
33
C                                (-SS SC)                        (SY**T)
 
34
C     where **T indicates transpose.  The elements of SX are in
 
35
C     SX(LX+I*INCX), I = 0 to N-1, where LX = 1 if INCX .GE. 0, else
 
36
C     LX = 1+(1-N)*INCX, and similarly for SY using LY and INCY.
 
37
C
 
38
C***REFERENCES  C. L. Lawson, R. J. Hanson, D. R. Kincaid and F. T.
 
39
C                 Krogh, Basic linear algebra subprograms for Fortran
 
40
C                 usage, Algorithm No. 539, Transactions on Mathematical
 
41
C                 Software 5, 3 (September 1979), pp. 308-323.
 
42
C***ROUTINES CALLED  (NONE)
 
43
C***REVISION HISTORY  (YYMMDD)
 
44
C   791001  DATE WRITTEN
 
45
C   861211  REVISION DATE from Version 3.2
 
46
C   891214  Prologue converted to Version 4.0 format.  (BAB)
 
47
C   920310  Corrected definition of LX in DESCRIPTION.  (WRB)
 
48
C   920501  Reformatted the REFERENCES section.  (WRB)
 
49
C***END PROLOGUE  SROT
 
50
      REAL SX, SY, SC, SS, ZERO, ONE, W, Z
 
51
      DIMENSION SX(*), SY(*)
 
52
      SAVE ZERO, ONE
 
53
      DATA ZERO, ONE /0.0E0, 1.0E0/
 
54
C***FIRST EXECUTABLE STATEMENT  SROT
 
55
      IF (N .LE. 0 .OR. (SS .EQ. ZERO .AND. SC .EQ. ONE)) GO TO 40
 
56
      IF (.NOT. (INCX .EQ. INCY .AND. INCX .GT. 0)) GO TO 20
 
57
C
 
58
C          Code for equal and positive increments.
 
59
C
 
60
           NSTEPS=INCX*N
 
61
           DO 10 I = 1,NSTEPS,INCX
 
62
                W=SX(I)
 
63
                Z=SY(I)
 
64
                SX(I)=SC*W+SS*Z
 
65
                SY(I)=-SS*W+SC*Z
 
66
   10           CONTINUE
 
67
           GO TO 40
 
68
C
 
69
C     Code for unequal or nonpositive increments.
 
70
C
 
71
   20 CONTINUE
 
72
           KX=1
 
73
           KY=1
 
74
C
 
75
           IF (INCX .LT. 0) KX = 1-(N-1)*INCX
 
76
           IF (INCY .LT. 0) KY = 1-(N-1)*INCY
 
77
C
 
78
           DO 30 I = 1,N
 
79
                W=SX(KX)
 
80
                Z=SY(KY)
 
81
                SX(KX)=SC*W+SS*Z
 
82
                SY(KY)=-SS*W+SC*Z
 
83
                KX=KX+INCX
 
84
                KY=KY+INCY
 
85
   30           CONTINUE
 
86
   40 CONTINUE
 
87
C
 
88
      RETURN
 
89
      END