~ubuntu-branches/ubuntu/vivid/atlas/vivid

« back to all changes in this revision

Viewing changes to src/blas/f77reference/dnrm2.f

  • Committer: Package Import Robot
  • Author(s): Sébastien Villemot, Sylvestre Ledru, Sébastien Villemot
  • Date: 2013-06-11 15:58:16 UTC
  • mfrom: (1.1.4) (25 sid)
  • mto: This revision was merged to the branch mainline in revision 26.
  • Revision ID: package-import@ubuntu.com-20130611155816-8xeeiziu1iml040c
Tags: 3.10.1-1
[ Sylvestre Ledru ]
* New upstream release (Closes: #609287)

[ Sébastien Villemot ]
* Provide architectural defaults (i.e. precomputed timings) for all
  release archs (except armel and mips for the time being, due to slow
  porterboxes). This will make the package build much faster and should
  eliminate transient build failures due to excessive variance in the
  timings.
* Move symlinks for lib{cblas,f77blas,atlas,lapack_atlas} out of the
  libblas.so.3 alternative and make them always present, so that
  software relying on these libs do not break when another alternative
  is selected for BLAS
* ATLAS now has improved ARM support with native asm constructs. This required
  the following tunes:
  + armel-is-v4t.diff: new patch, prevents FTBFS on armel; otherwise,
    ATLAS uses asm constructs too recent for the platform (armel is only v4t)
  + debian/rules: on armhf, define the ATL_ARM_HARDFP flag; otherwise the asm
    constructs use the soft-float ABI for passing floating points
  + on armhf, ensure that -mfloat-abi=softfp and -mcpu=vfpv3 flags are never
    used; this is implemented via a patch (armhf.diff) and by the use of fixed
    archdefs
* The generic package is now built without multi-threading, because otherwise
  the package fails to build on some single-processor machines (this required
  the introduction of a patch: fix-non-threaded-build.diff). As a side effect,
  the build of the custom package gracefully handles non-threaded
  builds. (Closes: #602524)
* Add libblas.a as slave in the libblas.so alternative (Closes: #701921)
* Add symlinks for lib{f77blas,atlas}.a in /usr/lib (Closes: #666203)
* Modify shlibs file of libatlas3-base, such that packages using
  libblas/liblapack depend on any BLAS/LAPACK alternative, while packages
  depending on ATLAS-specific libraries (e.g. libatlas.so) depend specifically
  on libatlas3-base.
* corei1.diff: remove patch, applied upstream
* Use my @debian.org email address
* Remove obsolete DM-Upload-Allowed flag
* Switch VCS to git
* Remove Conflicts/Replaces against pre-squeeze packages
* libatlas-base-dev now provides libblas.so, as libblas-dev
* No longer use -Wa,--noexecstack in CFLAGS, it makes the package FTBFS
* Do not use POWER3 arch for powerpcspe port (Closes: #701068)
* Bump to debhelper compat level 9
* README.Debian: mention that devscripts is needed to compile the custom
  package (Closes: #697431)
* Bump Standards-Version to 3.9.4. As a consequence, add Built-Using
  fields because the package embeds stuff from liblapack-pic

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
      DOUBLE PRECISION FUNCTION DNRM2(N,X,INCX)
 
2
*     .. Scalar Arguments ..
 
3
      INTEGER INCX,N
 
4
*     ..
 
5
*     .. Array Arguments ..
 
6
      DOUBLE PRECISION X(*)
 
7
*     ..
 
8
*
 
9
*  Purpose
 
10
*  =======
 
11
*
 
12
*  DNRM2 returns the euclidean norm of a vector via the function
 
13
*  name, so that
 
14
*
 
15
*     DNRM2 := sqrt( x'*x )
 
16
*
 
17
*
 
18
*  -- This version written on 25-October-1982.
 
19
*     Modified on 14-October-1993 to inline the call to DLASSQ.
 
20
*     Sven Hammarling, Nag Ltd.
 
21
*
 
22
*
 
23
*     .. Parameters ..
 
24
      DOUBLE PRECISION ONE,ZERO
 
25
      PARAMETER (ONE=1.0D+0,ZERO=0.0D+0)
 
26
*     ..
 
27
*     .. Local Scalars ..
 
28
      DOUBLE PRECISION ABSXI,NORM,SCALE,SSQ
 
29
      INTEGER IX
 
30
*     ..
 
31
*     .. Intrinsic Functions ..
 
32
      INTRINSIC ABS,SQRT
 
33
*     ..
 
34
      IF (N.LT.1 .OR. INCX.LT.1) THEN
 
35
          NORM = ZERO
 
36
      ELSE IF (N.EQ.1) THEN
 
37
          NORM = ABS(X(1))
 
38
      ELSE
 
39
          SCALE = ZERO
 
40
          SSQ = ONE
 
41
*        The following loop is equivalent to this call to the LAPACK
 
42
*        auxiliary routine:
 
43
*        CALL DLASSQ( N, X, INCX, SCALE, SSQ )
 
44
*
 
45
          DO 10 IX = 1,1 + (N-1)*INCX,INCX
 
46
              IF (X(IX).NE.ZERO) THEN
 
47
                  ABSXI = ABS(X(IX))
 
48
                  IF (SCALE.LT.ABSXI) THEN
 
49
                      SSQ = ONE + SSQ* (SCALE/ABSXI)**2
 
50
                      SCALE = ABSXI
 
51
                  ELSE
 
52
                      SSQ = SSQ + (ABSXI/SCALE)**2
 
53
                  END IF
 
54
              END IF
 
55
   10     CONTINUE
 
56
          NORM = SCALE*SQRT(SSQ)
 
57
      END IF
 
58
*
 
59
      DNRM2 = NORM
 
60
      RETURN
 
61
*
 
62
*     End of DNRM2.
 
63
*
 
64
      END