Skip to content

Update TWOPI and PIOVER2 to 39 significant digits (issue #428). Added make.inc for gfortran in quadruple precision #433

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 1 commit into from
Aug 7, 2020
Merged

Update TWOPI and PIOVER2 to 39 significant digits (issue #428). Added make.inc for gfortran in quadruple precision #433

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
90 changes: 90 additions & 0 deletions INSTALL/make.inc.gfortran-quad
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,90 @@
####################################################################
# LAPACK make include file. #
# LAPACK, Version 3.7.0 #
# November 2017 #
####################################################################

SHELL = /bin/sh

# CC is the C compiler, normally invoked with options CFLAGS.
#
CC = gcc
CFLAGS = -O3

# Modify the FC and FFLAGS definitions to the desired compiler
# and desired compiler options for your machine. NOOPT refers to
# the compiler options desired when NO OPTIMIZATION is selected.
#
# Note: During a regular execution, LAPACK might create NaN and Inf
# and handle these quantities appropriately. As a consequence, one
# should not compile LAPACK with flags such as -ffpe-trap=overflow.
#
# The -freal-4-real-16 -freal-8-real-16 flags instruct gfortran
# to treal REAL, DOUBLE PRECISION, and corresponding COMPLEX
# declarations as 128-bit floating-point types. These flags
# are supported from gfortran 4.7 on. Please note that both
# s*/c* and d*/z* entry points will expect 128-bit floating-point
# type.
#
FC = gfortran
FFLAGS = -O2 -frecursive -freal-4-real-16 -freal-8-real-16
FFLAGS_DRV = $(FFLAGS)
FFLAGS_NOOPT = -O0 -frecursive -freal-4-real-16 -freal-8-real-16

# Define LDFLAGS to the desired linker options for your machine.
#
LDFLAGS =

# The archiver and the flag(s) to use when building an archive
# (library). If your system has no ranlib, set RANLIB = echo.
#
AR = ar
ARFLAGS = cr
RANLIB = ranlib

# Timer for the SECOND and DSECND routines
#
# Default: SECOND and DSECND will use a call to the
# EXTERNAL FUNCTION ETIME
#TIMER = EXT_ETIME
# For RS6K: SECOND and DSECND will use a call to the
# EXTERNAL FUNCTION ETIME_
#TIMER = EXT_ETIME_
# For gfortran compiler: SECOND and DSECND will use a call to the
# INTERNAL FUNCTION ETIME
#TIMER = INT_ETIME
# If your Fortran compiler does not provide etime (like Nag Fortran
# Compiler, etc...) SECOND and DSECND will use a call to the
# INTERNAL FUNCTION CPU_TIME
#TIMER = INT_CPU_TIME
# If none of these work, you can use the NONE value.
# In that case, SECOND and DSECND will always return 0.
TIMER = NONE

# Uncomment the following line to include deprecated routines in
# the LAPACK library.
#
#BUILD_DEPRECATED = Yes

# LAPACKE has the interface to some routines from tmglib.
# If LAPACKE_WITH_TMG is defined, add those routines to LAPACKE.
#
#LAPACKE_WITH_TMG = Yes

# Location of the extended-precision BLAS (XBLAS) Fortran library
# used for building and testing extended-precision routines. The
# relevant routines will be compiled and XBLAS will be linked only
# if USEXBLAS is defined.
#
#USEXBLAS = Yes
#XBLASLIB = -lxblas

# The location of the libraries to which you will link. (The
# machine-specific, optimized BLAS library should be used whenever
# possible.)
#
BLASLIB = $(TOPSRCDIR)/librefblas.a
CBLASLIB = $(TOPSRCDIR)/libcblas.a
LAPACKLIB = $(TOPSRCDIR)/liblapack.a
TMGLIB = $(TOPSRCDIR)/libtmglib.a
LAPACKELIB = $(TOPSRCDIR)/liblapacke.a
7 changes: 4 additions & 3 deletions SRC/cbbcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -354,12 +354,13 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
* .. Parameters ..
INTEGER MAXITR
PARAMETER ( MAXITR = 6 )
REAL HUNDRED, MEIGHTH, ONE, PIOVER2, TEN, ZERO
REAL HUNDRED, MEIGHTH, ONE, TEN, ZERO
PARAMETER ( HUNDRED = 100.0E0, MEIGHTH = -0.125E0,
$ ONE = 1.0E0, PIOVER2 = 1.57079632679489662E0,
$ TEN = 10.0E0, ZERO = 0.0E0 )
$ ONE = 1.0E0, TEN = 10.0E0, ZERO = 0.0E0 )
COMPLEX NEGONECOMPLEX
PARAMETER ( NEGONECOMPLEX = (-1.0E0,0.0E0) )
REAL PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210E0 )
* ..
* .. Local Scalars ..
LOGICAL COLMAJOR, LQUERY, RESTART11, RESTART12,
Expand Down
2 changes: 1 addition & 1 deletion SRC/clarnv.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -120,7 +120,7 @@ SUBROUTINE CLARNV( IDIST, ISEED, N, X )
INTEGER LV
PARAMETER ( LV = 128 )
REAL TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663E+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839E+0 )
* ..
* .. Local Scalars ..
INTEGER I, IL, IV
Expand Down
7 changes: 4 additions & 3 deletions SRC/dbbcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -354,12 +354,13 @@ SUBROUTINE DBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
* .. Parameters ..
INTEGER MAXITR
PARAMETER ( MAXITR = 6 )
DOUBLE PRECISION HUNDRED, MEIGHTH, ONE, PIOVER2, TEN, ZERO
DOUBLE PRECISION HUNDRED, MEIGHTH, ONE, TEN, ZERO
PARAMETER ( HUNDRED = 100.0D0, MEIGHTH = -0.125D0,
$ ONE = 1.0D0, PIOVER2 = 1.57079632679489662D0,
$ TEN = 10.0D0, ZERO = 0.0D0 )
$ ONE = 1.0D0, TEN = 10.0D0, ZERO = 0.0D0 )
DOUBLE PRECISION NEGONE
PARAMETER ( NEGONE = -1.0D0 )
DOUBLE PRECISION PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210D0 )
* ..
* .. Local Scalars ..
LOGICAL COLMAJOR, LQUERY, RESTART11, RESTART12,
Expand Down
2 changes: 1 addition & 1 deletion SRC/dlarnv.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -118,7 +118,7 @@ SUBROUTINE DLARNV( IDIST, ISEED, N, X )
INTEGER LV
PARAMETER ( LV = 128 )
DOUBLE PRECISION TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663D+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839D+0 )
* ..
* .. Local Scalars ..
INTEGER I, IL, IL2, IV
Expand Down
7 changes: 4 additions & 3 deletions SRC/sbbcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -354,12 +354,13 @@ SUBROUTINE SBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
* .. Parameters ..
INTEGER MAXITR
PARAMETER ( MAXITR = 6 )
REAL HUNDRED, MEIGHTH, ONE, PIOVER2, TEN, ZERO
REAL HUNDRED, MEIGHTH, ONE, TEN, ZERO
PARAMETER ( HUNDRED = 100.0E0, MEIGHTH = -0.125E0,
$ ONE = 1.0E0, PIOVER2 = 1.57079632679489662E0,
$ TEN = 10.0E0, ZERO = 0.0E0 )
$ ONE = 1.0E0, TEN = 10.0E0, ZERO = 0.0E0 )
REAL NEGONE
PARAMETER ( NEGONE = -1.0E0 )
REAL PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210E0 )
* ..
* .. Local Scalars ..
LOGICAL COLMAJOR, LQUERY, RESTART11, RESTART12,
Expand Down
2 changes: 1 addition & 1 deletion SRC/slarnv.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -118,7 +118,7 @@ SUBROUTINE SLARNV( IDIST, ISEED, N, X )
INTEGER LV
PARAMETER ( LV = 128 )
REAL TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663E+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839E+0 )
* ..
* .. Local Scalars ..
INTEGER I, IL, IL2, IV
Expand Down
7 changes: 4 additions & 3 deletions SRC/zbbcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -354,12 +354,13 @@ SUBROUTINE ZBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
* .. Parameters ..
INTEGER MAXITR
PARAMETER ( MAXITR = 6 )
DOUBLE PRECISION HUNDRED, MEIGHTH, ONE, PIOVER2, TEN, ZERO
DOUBLE PRECISION HUNDRED, MEIGHTH, ONE, TEN, ZERO
PARAMETER ( HUNDRED = 100.0D0, MEIGHTH = -0.125D0,
$ ONE = 1.0D0, PIOVER2 = 1.57079632679489662D0,
$ TEN = 10.0D0, ZERO = 0.0D0 )
$ ONE = 1.0D0, TEN = 10.0D0, ZERO = 0.0D0 )
COMPLEX*16 NEGONECOMPLEX
PARAMETER ( NEGONECOMPLEX = (-1.0D0,0.0D0) )
DOUBLE PRECISION PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210D0 )
* ..
* .. Local Scalars ..
LOGICAL COLMAJOR, LQUERY, RESTART11, RESTART12,
Expand Down
2 changes: 1 addition & 1 deletion SRC/zlarnv.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -120,7 +120,7 @@ SUBROUTINE ZLARNV( IDIST, ISEED, N, X )
INTEGER LV
PARAMETER ( LV = 128 )
DOUBLE PRECISION TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663D+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839D+0 )
* ..
* .. Local Scalars ..
INTEGER I, IL, IV
Expand Down
5 changes: 3 additions & 2 deletions TESTING/EIG/cckcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -208,13 +208,14 @@ SUBROUTINE CCKCSD( NM, MVAL, PVAL, QVAL, NMATS, ISEED, THRESH,
PARAMETER ( NTESTS = 15 )
INTEGER NTYPES
PARAMETER ( NTYPES = 4 )
REAL GAPDIGIT, ORTH, PIOVER2, REALONE, REALZERO, TEN
REAL GAPDIGIT, ORTH, REALONE, REALZERO, TEN
PARAMETER ( GAPDIGIT = 10.0E0, ORTH = 1.0E-4,
$ PIOVER2 = 1.57079632679489662E0,
$ REALONE = 1.0E0, REALZERO = 0.0E0,
$ TEN = 10.0E0 )
COMPLEX ONE, ZERO
PARAMETER ( ONE = (1.0E0,0.0E0), ZERO = (0.0E0,0.0E0) )
REAL PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210E0 )
* ..
* .. Local Scalars ..
LOGICAL FIRSTT
Expand Down
7 changes: 4 additions & 3 deletions TESTING/EIG/ccsdts.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -248,11 +248,12 @@ SUBROUTINE CCSDTS( M, P, Q, X, XF, LDX, U1, LDU1, U2, LDU2, V1T,
* =====================================================================
*
* .. Parameters ..
REAL PIOVER2, REALONE, REALZERO
PARAMETER ( PIOVER2 = 1.57079632679489662E0,
$ REALONE = 1.0E0, REALZERO = 0.0E0 )
REAL REALONE, REALZERO
PARAMETER ( REALONE = 1.0E0, REALZERO = 0.0E0 )
COMPLEX ZERO, ONE
PARAMETER ( ZERO = (0.0E0,0.0E0), ONE = (1.0E0,0.0E0) )
REAL PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210E0 )
* ..
* .. Local Scalars ..
INTEGER I, INFO, R
Expand Down
5 changes: 3 additions & 2 deletions TESTING/EIG/dckcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -208,11 +208,12 @@ SUBROUTINE DCKCSD( NM, MVAL, PVAL, QVAL, NMATS, ISEED, THRESH,
PARAMETER ( NTESTS = 15 )
INTEGER NTYPES
PARAMETER ( NTYPES = 4 )
DOUBLE PRECISION GAPDIGIT, ONE, ORTH, PIOVER2, TEN, ZERO
DOUBLE PRECISION GAPDIGIT, ONE, ORTH, TEN, ZERO
PARAMETER ( GAPDIGIT = 18.0D0, ONE = 1.0D0,
$ ORTH = 1.0D-12,
$ PIOVER2 = 1.57079632679489662D0,
$ TEN = 10.0D0, ZERO = 0.0D0 )
DOUBLE PRECISION PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210D0 )
* ..
* .. Local Scalars ..
LOGICAL FIRSTT
Expand Down
7 changes: 4 additions & 3 deletions TESTING/EIG/dcsdts.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -248,11 +248,12 @@ SUBROUTINE DCSDTS( M, P, Q, X, XF, LDX, U1, LDU1, U2, LDU2, V1T,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION PIOVER2, REALONE, REALZERO
PARAMETER ( PIOVER2 = 1.57079632679489662D0,
$ REALONE = 1.0D0, REALZERO = 0.0D0 )
DOUBLE PRECISION REALONE, REALZERO
PARAMETER ( REALONE = 1.0D0, REALZERO = 0.0D0 )
DOUBLE PRECISION ZERO, ONE
PARAMETER ( ZERO = 0.0D0, ONE = 1.0D0 )
DOUBLE PRECISION PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210D0 )
* ..
* .. Local Scalars ..
INTEGER I, INFO, R
Expand Down
8 changes: 4 additions & 4 deletions TESTING/EIG/sckcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -208,11 +208,11 @@ SUBROUTINE SCKCSD( NM, MVAL, PVAL, QVAL, NMATS, ISEED, THRESH,
PARAMETER ( NTESTS = 15 )
INTEGER NTYPES
PARAMETER ( NTYPES = 4 )
REAL GAPDIGIT, ONE, ORTH, PIOVER2, TEN, ZERO
REAL GAPDIGIT, ONE, ORTH, TEN, ZERO
PARAMETER ( GAPDIGIT = 10.0E0, ONE = 1.0E0,
$ ORTH = 1.0E-4,
$ PIOVER2 = 1.57079632679489662E0,
$ TEN = 10.0E0, ZERO = 0.0E0 )
$ ORTH = 1.0E-4, TEN = 10.0E0, ZERO = 0.0E0 )
REAL PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210E0 )
* ..
* .. Local Scalars ..
LOGICAL FIRSTT
Expand Down
7 changes: 4 additions & 3 deletions TESTING/EIG/scsdts.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -248,11 +248,12 @@ SUBROUTINE SCSDTS( M, P, Q, X, XF, LDX, U1, LDU1, U2, LDU2, V1T,
* =====================================================================
*
* .. Parameters ..
REAL PIOVER2, REALONE, REALZERO
PARAMETER ( PIOVER2 = 1.57079632679489662E0,
$ REALONE = 1.0E0, REALZERO = 0.0E0 )
REAL REALONE, REALZERO
PARAMETER ( REALONE = 1.0E0, REALZERO = 0.0E0 )
REAL ZERO, ONE
PARAMETER ( ZERO = 0.0E0, ONE = 1.0E0 )
REAL PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210E0 )
* ..
* .. Local Scalars ..
INTEGER I, INFO, R
Expand Down
5 changes: 3 additions & 2 deletions TESTING/EIG/zckcsd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -208,13 +208,14 @@ SUBROUTINE ZCKCSD( NM, MVAL, PVAL, QVAL, NMATS, ISEED, THRESH,
PARAMETER ( NTESTS = 15 )
INTEGER NTYPES
PARAMETER ( NTYPES = 4 )
DOUBLE PRECISION GAPDIGIT, ORTH, PIOVER2, REALONE, REALZERO, TEN
DOUBLE PRECISION GAPDIGIT, ORTH, REALONE, REALZERO, TEN
PARAMETER ( GAPDIGIT = 18.0D0, ORTH = 1.0D-12,
$ PIOVER2 = 1.57079632679489662D0,
$ REALONE = 1.0D0, REALZERO = 0.0D0,
$ TEN = 10.0D0 )
COMPLEX*16 ONE, ZERO
PARAMETER ( ONE = (1.0D0,0.0D0), ZERO = (0.0D0,0.0D0) )
DOUBLE PRECISION PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210D0 )
* ..
* .. Local Scalars ..
LOGICAL FIRSTT
Expand Down
7 changes: 4 additions & 3 deletions TESTING/EIG/zcsdts.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -248,11 +248,12 @@ SUBROUTINE ZCSDTS( M, P, Q, X, XF, LDX, U1, LDU1, U2, LDU2, V1T,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION PIOVER2, REALONE, REALZERO
PARAMETER ( PIOVER2 = 1.57079632679489662D0,
$ REALONE = 1.0D0, REALZERO = 0.0D0 )
DOUBLE PRECISION REALONE, REALZERO
PARAMETER ( REALONE = 1.0D0, REALZERO = 0.0D0 )
COMPLEX*16 ZERO, ONE
PARAMETER ( ZERO = (0.0D0,0.0D0), ONE = (1.0D0,0.0D0) )
DOUBLE PRECISION PIOVER2
PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210D0 )
* ..
* .. Local Scalars ..
INTEGER I, INFO, R
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/clarnd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -93,7 +93,7 @@ COMPLEX FUNCTION CLARND( IDIST, ISEED )
REAL ZERO, ONE, TWO
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0, TWO = 2.0E+0 )
REAL TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663E+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839E+0 )
* ..
* .. Local Scalars ..
REAL T1, T2
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/clatms.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -358,7 +358,7 @@ SUBROUTINE CLATMS( M, N, DIST, ISEED, SYM, D, MODE, COND, DMAX,
COMPLEX CZERO
PARAMETER ( CZERO = ( 0.0E+0, 0.0E+0 ) )
REAL TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663E+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839E+0 )
* ..
* .. Local Scalars ..
LOGICAL CSYM, GIVENS, ILEXTR, ILTEMP, TOPDWN
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/clatmt.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -366,7 +366,7 @@ SUBROUTINE CLATMT( M, N, DIST, ISEED, SYM, D, MODE, COND, DMAX,
COMPLEX CZERO
PARAMETER ( CZERO = ( 0.0E+0, 0.0E+0 ) )
REAL TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663E+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839E+0 )
* ..
* .. Local Scalars ..
COMPLEX C, CT, CTEMP, DUMMY, EXTRA, S, ST
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/dlarnd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -91,7 +91,7 @@ DOUBLE PRECISION FUNCTION DLARND( IDIST, ISEED )
DOUBLE PRECISION ONE, TWO
PARAMETER ( ONE = 1.0D+0, TWO = 2.0D+0 )
DOUBLE PRECISION TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663D+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839D+0 )
* ..
* .. Local Scalars ..
DOUBLE PRECISION T1, T2
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/dlatms.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -344,7 +344,7 @@ SUBROUTINE DLATMS( M, N, DIST, ISEED, SYM, D, MODE, COND, DMAX,
DOUBLE PRECISION ONE
PARAMETER ( ONE = 1.0D0 )
DOUBLE PRECISION TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663D+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839D+0 )
* ..
* .. Local Scalars ..
LOGICAL GIVENS, ILEXTR, ILTEMP, TOPDWN
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/dlatmt.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -354,7 +354,7 @@ SUBROUTINE DLATMT( M, N, DIST, ISEED, SYM, D, MODE, COND, DMAX,
DOUBLE PRECISION ONE
PARAMETER ( ONE = 1.0D0 )
DOUBLE PRECISION TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663D+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839D+0 )
* ..
* .. Local Scalars ..
DOUBLE PRECISION ALPHA, ANGLE, C, DUMMY, EXTRA, S, TEMP
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/slarnd.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -91,7 +91,7 @@ REAL FUNCTION SLARND( IDIST, ISEED )
REAL ONE, TWO
PARAMETER ( ONE = 1.0E+0, TWO = 2.0E+0 )
REAL TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663E+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839E+0 )
* ..
* .. Local Scalars ..
REAL T1, T2
Expand Down
2 changes: 1 addition & 1 deletion TESTING/MATGEN/slatms.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -344,7 +344,7 @@ SUBROUTINE SLATMS( M, N, DIST, ISEED, SYM, D, MODE, COND, DMAX,
REAL ONE
PARAMETER ( ONE = 1.0E0 )
REAL TWOPI
PARAMETER ( TWOPI = 6.2831853071795864769252867663E+0 )
PARAMETER ( TWOPI = 6.28318530717958647692528676655900576839E+0 )
* ..
* .. Local Scalars ..
LOGICAL GIVENS, ILEXTR, ILTEMP, TOPDWN
Expand Down
Loading