Return to Test_float CVS log

Up to [CSRG BSD Unix] / 43BSD / usr.lib / libF77

Annotation of 43BSD/usr.lib/libF77/Test_float, revision 1.1.1.1

1.1       root        1: #!     /bin/csh -f 
                      2: #
                      3: #      this tests if C rounds the return value to single precision in
                      4: #      functions of type float. (see sec. 9.10 of C ref. manual).
                      5: #
                      6: cat << 'EOT' >! /tmp/test_fltc.c
                      7: float temp;
                      8: 
                      9: float f1_(arg1,arg2)
                     10: float *arg1, *arg2;
                     11: {
                     12:        /* force float by storing in global */
                     13:        temp =  *arg1 / *arg2;
                     14:        return temp;
                     15: }
                     16: 
                     17: float f2_(arg1,arg2)
                     18: float *arg1, *arg2;
                     19: {
                     20:        /* should round since function is type float */
                     21:        return  *arg1 / *arg2;
                     22: }
                     23: 
                     24: 
                     25: float f3_(arg1,arg2)
                     26: float *arg1, *arg2;
                     27: {
                     28:        /* use a cast to try to force rounding */
                     29:        return ((float) (*arg1 / *arg2));
                     30: }
                     31: 
                     32: 'EOT'
                     33: cat << 'EOT' >! /tmp/test_fltf.f
                     34:        integer f2ok, f3ok
                     35:        data f2ok/0/, f3ok/0/
                     36: 
                     37:        do 20 i = 1,10
                     38:        do 10 j = 1,10
                     39:        x = 0.1d0*i
                     40:        y = 0.1d0*j
                     41:        temp = f1(x,y)
                     42:        if( f2(x,y).eq.temp) f2ok = f2ok + 1
                     43:        if( f3(x,y).eq.temp) f3ok = f3ok + 1
                     44: 10     continue
                     45: 20     continue
                     46:        print *, "out of 100 tries, f2 was ok", f2ok, "times"
                     47:        print *, "out of 100 tries, f3 was ok", f3ok, "times"
                     48:        end
                     49: 'EOT'
                     50: pushd /tmp
                     51: f77 test_fltc.c test_fltf.f -o test_flt
                     52: test_flt

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.