|
|
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.