[883ebf]: src / daskr.diff  Maximize  Restore  History

Download this file

65 lines (64 with data), 1.9 kB

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
--- daskr.orig/ddaskr.f 2008-03-14 06:34:48.000000000 -0400
+++ daskr/ddaskr.f 2014-03-02 16:10:47.319479666 -0500
@@ -1405,6 +1405,7 @@
DIMENSION RWORK(LRW),IWORK(LIW)
DIMENSION RTOL(*),ATOL(*)
DIMENSION RPAR(*),IPAR(*)
+ DIMENSION JROOT(*)
CHARACTER MSG*80
EXTERNAL RES, JAC, PSOL, RT, DDASID, DDASIK, DNEDD, DNEDK
C
@@ -2515,6 +2516,7 @@
* RWORK, RPAR
DIMENSION Y(*), YP(*), PHI(NEQ,*), PSI(*),
* R0(*), R1(*), RX(*), JROOT(*), RWORK(*), IWORK(*)
+ DIMENSION RPAR(*), IPAR(*)
INTEGER I, JFLAG
DOUBLE PRECISION H
DOUBLE PRECISION HMINR, T1, TEMP1, TEMP2, X, ZERO
--- daskr.orig/dlinpk.f 2008-03-14 06:34:42.000000000 -0400
+++ daskr/dlinpk.f 2014-03-02 16:09:47.580384128 -0500
@@ -760,14 +760,24 @@
dnrm2 = zero
go to 300
c
- 10 assign 30 to next
+ 10 next = 30
sum = zero
nn = n * incx
c begin main loop
i = 1
- 20 go to next,(30, 50, 70, 110)
+ 20 if (next .eq. 30) then
+ go to 30
+ else if (next .eq. 50) then
+ go to 50
+ else if (next .eq. 70) then
+ go to 70
+ else if (next .eq. 110) then
+ go to 110
+ else
+ continue
+ end if
30 if( dabs(dx(i)) .gt. cutlo) go to 85
- assign 50 to next
+ next = 50
xmax = zero
c
c phase 1. sum is zero
@@ -776,13 +786,13 @@
if( dabs(dx(i)) .gt. cutlo) go to 85
c
c prepare for phase 2.
- assign 70 to next
+ next = 70
go to 105
c
c prepare for phase 4.
c
100 i = j
- assign 110 to next
+ next = 110
sum = (sum / dx(i)) / dx(i)
105 xmax = dabs(dx(i))
go to 115

Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:





No, thanks