Mercurial > octave-nkf

comparison libcruft/odessa/odessa_intdy.f @ 3987:258c1d15ad78

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2002-07-11 19:33:35 by tenny] changed real to double precision
author tenny
date Thu, 11 Jul 2002 19:33:35 +0000
parents 7a37caf6ed43
children 3995501ce1d9
comparison
equal deleted inserted replaced
3986:7cb85d5c7aad 3987:258c1d15ad78
35 IC = 1 35 IC = 1
36 IF (K .EQ. 0) GO TO 15 36 IF (K .EQ. 0) GO TO 15
37 JJ1 = L - K 37 JJ1 = L - K
38 DO 10 JJ = JJ1,NQ 38 DO 10 JJ = JJ1,NQ
39 10 IC = IC*JJ 39 10 IC = IC*JJ
40 15 C = REAL(IC) 40 15 C = DBLE(IC)
41 DO 20 I = 1,NYH 41 DO 20 I = 1,NYH
42 20 DKY(I) = C*YH(I,L) 42 20 DKY(I) = C*YH(I,L)
43 IF (K .EQ. NQ) GO TO 55 43 IF (K .EQ. NQ) GO TO 55
44 JB2 = NQ - K 44 JB2 = NQ - K
45 DO 50 JB = 1,JB2 45 DO 50 JB = 1,JB2
48 IC = 1 48 IC = 1
49 IF (K .EQ. 0) GO TO 35 49 IF (K .EQ. 0) GO TO 35
50 JJ1 = JP1 - K 50 JJ1 = JP1 - K
51 DO 30 JJ = JJ1,J 51 DO 30 JJ = JJ1,J
52 30 IC = IC*JJ 52 30 IC = IC*JJ
53 35 C = REAL(IC) 53 35 C = DBLE(IC)
54 DO 40 I = 1,NYH 54 DO 40 I = 1,NYH
55 40 DKY(I) = C*YH(I,JP1) + S*DKY(I) 55 40 DKY(I) = C*YH(I,JP1) + S*DKY(I)
56 50 CONTINUE 56 50 CONTINUE
57 IF (K .EQ. 0) RETURN 57 IF (K .EQ. 0) RETURN
58 55 R = H**(-K) 58 55 R = H**(-K)