|
|
1.1 root 1: c subroutine tst3
2: parameter(ntest=3)
3: c Real arithmetic
4: parameter(R1=1.0, R1H=1.5, R2=2.0, R3=3.0, R4=4.0, RM1=-1.0)
5: parameter(D1=1.0d0, D1H=1.5d0, D2=2.0d0, D3=3.0d0, D4=4.0d0, DM1=-1.0d0)
6: c These constants are machine dependent
7: parameter(eps = 2.0e-7)
8: parameter(dps = 2.0d-15)
9: c common /t3/ x,y,z,dx,dy,dz
10: double precision dx,dy,dz
11:
12: if ( abs(1.0+1.0-2.0) .gt. eps ) call e(1,ntest)
13: if ( abs(1e10-1e10) .gt. eps ) call e(2,ntest)
14: if ( abs(R1+R1H+R2+R3+R4+RM1-10.5) .gt. eps ) call e(3,ntest)
15: if ( abs(1.0e-1 * 1.0e1 - 100e-2) .gt. eps ) call e(4,ntest)
16: if ( abs(10.0/3.0*3.0/10.0-100e-2) .gt. eps ) call e(5,ntest)
17: if ( 0.0e0 .ne. 0 ) call e(6,ntest)
18: if ( abs(32767.0-32767.0) .gt. eps ) call e(7,ntest)
19: if ( abs(1.0+2+5+3.0e0+5.0e+0+140e-1-30.000)/100 .gt. eps ) call e(8,ntest)
20: if ( abs(-1+(-1)+(-1.0)+(-1e0)+(-1e+0)+(-1e-0) + ((((6)))) ) .gt. eps ) call e(9,ntest)
21:
22: x=1.50
23: y=3.00
24: z= 0.10
25: if ( abs(5*y*z-x) .gt. eps ) call e(10,ntest)
26: if ( abs(y*y*y/z*x-405) .gt. eps ) call e(11,ntest)
27: x=1.1
28: y= 1.2
29: if ( y.lt.x ) call e(12,ntest)
30: if ( y .le. x ) call e(13,ntest)
31: if ( y .eq. x ) call e(14,ntest)
32: if ( x .ne. x ) call e(15,ntest)
33: if ( x .ge. y ) call e(16,ntest)
34: if ( x .gt.y ) call e(17,ntest)
35:
36: if ( abs(1.0+1.0-2.0) .gt. dps ) call e(18,ntest)
37: if ( abs(1d10-1d10) .gt. dps ) call e(19,ntest)
38: if ( abs(D1+D1H+D2+D3+D4+DM1-10.5) .gt. dps ) call e(20,ntest)
39: if ( abs(1.0d-1 * 1.0d1 - 100d-2) .gt. dps ) call e(21,ntest)
40: if ( abs(10.0/3.0*3.0/10.0-100d-2) .gt. dps ) call e(22,ntest)
41: if ( 0.0d0 .ne. 0 ) call e(23,ntest)
42: if ( abs(3.2767d4-32767.d0) .gt. dps ) call e(24,ntest)
43: if ( abs(1.0+2+5+3.0d0+5.0d+0+140d-1-30.000)/100 .gt. dps ) call e(25,ntest)
44: if ( abs(-1+(-1)+(-1.0)+(-1d0)+(-1d+0)+(-1d-0) + ((((6)))) ) .gt. dps ) call e(26,ntest)
45:
46: dx=1.50d0
47: dy=3.00d0
48: dz= 0.10d0
49: if ( abs(5*dy*dz-dx) .gt. dps ) call e(27,ntest)
50: if ( abs(dy*dy*dy/dz*dx-405) .gt. dps ) call e(28,ntest)
51: dx=1.1d0
52: dy= 1.2d0
53: if ( dy.lt.dx ) call e(29,ntest)
54: if ( dy .le. dx ) call e(30,ntest)
55: if ( dy .eq. dx ) call e(31,ntest)
56: if ( dx .ne. dx ) call e(32,ntest)
57: if ( dx .ge. dy ) call e(33,ntest)
58: if ( dx .gt.dy ) call e(34,ntest)
59: end
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.