Skip to content

Commit 2cb17f9

Browse files
HugoMValeHugoMVale
authored
fix wp->dp
1 parent b22bd51 commit 2cb17f9

File tree

3 files changed

+34
-34
lines changed

3 files changed

+34
-34
lines changed

‎extern/blas.f‎

Lines changed: 32 additions & 32 deletions
Original file line numberDiff line numberDiff line change
@@ -30,7 +30,7 @@ FUNCTION DASUM(N,DX,INCX) RESULT(DASUMR)
3030
C***End Prologue DASUM
3131

3232
C...Used modules
33-
use odrpack_kinds, only: wp
33+
use odrpack_kinds, only: dp
3434

3535
C...Scalar arguments
3636
INTEGER
@@ -56,7 +56,7 @@ FUNCTION DASUM(N,DX,INCX) RESULT(DASUMR)
5656
C***First executable statement DASUM
5757

5858

59-
DASUMR = 0.E0_wp
59+
DASUMR = 0.E0_dp
6060
IF(N.LE.0)RETURN
6161
IF(INCX.EQ.1)GOTO 20
6262

@@ -121,7 +121,7 @@ SUBROUTINE DAXPY(N,DA,DX,INCX,DY,INCY)
121121
C***End Prologue DAXPY
122122

123123
C...Used modules
124-
use odrpack_kinds, only: wp
124+
use odrpack_kinds, only: dp
125125

126126
C...Scalar arguments
127127
REAL(wp)
@@ -145,7 +145,7 @@ SUBROUTINE DAXPY(N,DA,DX,INCX,DY,INCY)
145145
C***First executable statement DAXPY
146146

147147

148-
IF(N.LE.0.OR.DA.EQ.0.E0_wp) RETURN
148+
IF(N.LE.0.OR.DA.EQ.0.E0_dp) RETURN
149149
IF(INCX.EQ.INCY) THEN
150150
IF(INCX-1.LT.0) THEN
151151
GOTO 5
@@ -234,7 +234,7 @@ SUBROUTINE DCOPY(N,DX,INCX,DY,INCY)
234234
C***End Prologue DCOPY
235235

236236
C...Used modules
237-
use odrpack_kinds, only: wp
237+
use odrpack_kinds, only: dp
238238

239239
C...Scalar arguments
240240
INTEGER
@@ -348,7 +348,7 @@ FUNCTION DDOT(N,DX,INCX,DY,INCY) RESULT(DDOTR)
348348
C***End Prologue DDOT
349349

350350
C...Used modules
351-
use odrpack_kinds, only: wp
351+
use odrpack_kinds, only: dp
352352

353353
C...Scalar arguments
354354
INTEGER
@@ -374,7 +374,7 @@ FUNCTION DDOT(N,DX,INCX,DY,INCY) RESULT(DDOTR)
374374
C***First executable statement DDOT
375375

376376

377-
DDOTR = 0.E0_wp
377+
DDOTR = 0.E0_dp
378378
IF(N.LE.0)RETURN
379379
IF(INCX.EQ.INCY) THEN
380380
IF(INCX-1.LT.0) THEN
@@ -474,14 +474,14 @@ FUNCTION DNRM2(N,DX,INCX) RESULT(DNRM2R)
474474
C document the limiting values are as follows..
475475
C CUTLO, S.P. U/EPS = 2**(-102) for Honeywell. Close seconds are
476476
C UNIVAC and DEC at 2**(-103)
477-
C thus CUTLO = 2**(-51) = 4.44089E-16_wp
477+
C thus CUTLO = 2**(-51) = 4.44089E-16_dp
478478
C CUTHI, S.P. V = 2**127 for UNIVAC, Honeywell, and DEC.
479-
C thus CUTHI = 2**(63.5) = 1.30438E19_wp
479+
C thus CUTHI = 2**(63.5) = 1.30438E19_dp
480480
C CUTLO, D.P. U/EPS = 2**(-67) for Honeywell and DEC.
481-
C thus CUTLO = 2**(-33.5) = 8.23181E-11_wp
482-
C CUTHI, D.P. Same as S.P. CUTHI = 1.30438E19_wp
483-
C DATA CUTLO, CUTHI / 8.232E-11_wp, 1.304E19_wp /
484-
C DATA CUTLO, CUTHI / 4.441E-16_wp, 1.304E19_wp /
481+
C thus CUTLO = 2**(-33.5) = 8.23181E-11_dp
482+
C CUTHI, D.P. Same as S.P. CUTHI = 1.30438E19_dp
483+
C DATA CUTLO, CUTHI / 8.232E-11_dp, 1.304E19_dp /
484+
C DATA CUTLO, CUTHI / 4.441E-16_dp, 1.304E19_dp /
485485
C***References Lawson C.L., Hanson R.J., Kincaid D.R., Krogh F.T.,
486486
C *Basic Linear Algebra Subprograms for FORTRAN Usage*,
487487
C Algorithm No. 539, Transactions on Mathematical
@@ -490,7 +490,7 @@ FUNCTION DNRM2(N,DX,INCX) RESULT(DNRM2R)
490490
C***End Prologue DNRM2
491491

492492
C...Used modules
493-
use odrpack_kinds, only: wp
493+
use odrpack_kinds, only: dp
494494

495495
C...Scalar arguments
496496
INTEGER
@@ -516,9 +516,9 @@ FUNCTION DNRM2(N,DX,INCX) RESULT(DNRM2R)
516516

517517
C...Data statements
518518
DATA
519-
& ZERO,ONE/0.0E0_wp,1.0E0_wp/
519+
& ZERO,ONE/0.0E0_dp,1.0E0_dp/
520520
DATA
521-
& CUTLO,CUTHI/8.232E-11_wp,1.304E19_wp/
521+
& CUTLO,CUTHI/8.232E-11_dp,1.304E19_dp/
522522

523523

524524
C***First executable statement DNRM2
@@ -646,7 +646,7 @@ SUBROUTINE DROT(N,DX,INCX,DY,INCY,DC,DS)
646646
C***End Prologue DROT
647647

648648
C...Used modules
649-
use odrpack_kinds, only: wp
649+
use odrpack_kinds, only: dp
650650

651651
C...Scalar arguments
652652
REAL(wp)
@@ -666,7 +666,7 @@ SUBROUTINE DROT(N,DX,INCX,DY,INCY,DC,DS)
666666

667667
C...Data statements
668668
DATA
669-
& ZERO,ONE/0.E0_wp,1.E0_wp/
669+
& ZERO,ONE/0.E0_dp,1.E0_dp/
670670

671671

672672
C***First executable statement DROT
@@ -735,7 +735,7 @@ SUBROUTINE DROTG(DA,DB,DC,DS)
735735
C the quantity R = (+/-)DSQRT(DA**2 + DB**2) overwrites DA in
736736
C storage. The value of DB is overwritten by a value Z which
737737
C allows DC and DS to be recovered by the following algorithm.
738-
C If Z=1 set DC=0.E0_wp and DS=1.E0_wp
738+
C If Z=1 set DC=0.E0_dp and DS=1.E0_dp
739739
C If DABS(Z) .LT. 1 set DC=DSQRT(1-Z**2) and DS=Z
740740
C If DABS(Z) .GT. 1 set DC=1/Z and DS=DSQRT(1-DC**2)
741741
C Normally, the subprogram DROT(N,DX,INCX,DY,INCY,DC,DS) will
@@ -748,7 +748,7 @@ SUBROUTINE DROTG(DA,DB,DC,DS)
748748
C***End Prologue DROTG
749749

750750
C...Used modules
751-
use odrpack_kinds, only: wp
751+
use odrpack_kinds, only: dp
752752

753753
C...Scalar arguments
754754
REAL(wp)
@@ -775,7 +775,7 @@ SUBROUTINE DROTG(DA,DB,DC,DS)
775775

776776
C Note that U and R have the sign of DA
777777

778-
R = DSQRT(.25E0_wp + V**2) * U
778+
R = DSQRT(.25E0_dp + V**2) * U
779779

780780
C Note that DC is positive
781781

@@ -787,29 +787,29 @@ SUBROUTINE DROTG(DA,DB,DC,DS)
787787

788788
C *** Here DABS(DA) .LE. DABS(DB) ***
789789

790-
10 IF (DB .EQ. 0.E0_wp) GO TO 20
790+
10 IF (DB .EQ. 0.E0_dp) GO TO 20
791791
U = DB + DB
792792
V = DA / U
793793

794794
C Note that U and R have the sign of DB
795795
C (R is immediately stored in DA)
796796

797-
DA = DSQRT(.25E0_wp + V**2) * U
797+
DA = DSQRT(.25E0_dp + V**2) * U
798798

799799
C Note that DS is positive
800800

801801
DS = DB / DA
802802
DC = V * (DS + DS)
803-
IF (DC .EQ. 0.E0_wp) GO TO 15
804-
DB = 1.E0_wp / DC
803+
IF (DC .EQ. 0.E0_dp) GO TO 15
804+
DB = 1.E0_dp / DC
805805
RETURN
806-
15 DB = 1.E0_wp
806+
15 DB = 1.E0_dp
807807
RETURN
808808

809-
C *** Here DA = DB = 0.E0_wp ***
809+
C *** Here DA = DB = 0.E0_dp ***
810810

811-
20 DC = 1.E0_wp
812-
DS = 0.E0_wp
811+
20 DC = 1.E0_dp
812+
DS = 0.E0_dp
813813
RETURN
814814

815815
END
@@ -845,7 +845,7 @@ SUBROUTINE DSCAL(N,DA,DX,INCX)
845845
C***End Prologue DSCAL
846846

847847
C...Used modules
848-
use odrpack_kinds, only: wp
848+
use odrpack_kinds, only: dp
849849

850850
C...Scalar arguments
851851
REAL(wp)
@@ -937,7 +937,7 @@ SUBROUTINE DSWAP(N,DX,INCX,DY,INCY)
937937
C***End Prologue DSWAP
938938

939939
C...Used modules
940-
use odrpack_kinds, only: wp
940+
use odrpack_kinds, only: dp
941941

942942
C...Scalar arguments
943943
INTEGER
@@ -1058,7 +1058,7 @@ FUNCTION IDAMAX(N,DX,INCX) RESULT(IDAMAXR)
10581058
C***End Prologue IDAMAX
10591059

10601060
C...Used modules
1061-
use odrpack_kinds, only: wp
1061+
use odrpack_kinds, only: dp
10621062

10631063
C...Scalar arguments
10641064
INTEGER

‎fpm.toml‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,5 @@
11
name = "odrpack95"
2-
version = "3.0.0"
2+
version = "3.0.1"
33
license = "MIT"
44
author = "HugoMVale"
55
maintainer = "57530119+HugoMVale@users.noreply.github.com"

‎meson.build‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
project(
22
'odrpack95',
33
['fortran', 'c'],
4-
version: '3.0.0',
4+
version: '3.0.1',
55
meson_version: '>=1.1',
66
license: 'MIT',
77
default_options: [

0 commit comments

Comments
 (0)