test_eps Subroutine

subroutine test_eps()

logs=.true. logs=.false.

Arguments

None

Contents

Source Code


Source Code

subroutine test_eps ()
   !!logs=.true.
   if(logs)call lala( 'diary(''eps.log'');')
   call lala( [ character(len=256) :: &
     & ' display(ones(80,1)''*''='');help eps;display(ones(80,1)''*''_'')              ', &
     & 'tally=[0];                                                               ', &
     & ' // find the eps for this programming environment by brute force         ', &
     & ' myeps = 1;                                                              ', &
     & ' while 1 + myeps > 1, myeps = myeps/2;                                   ', &
     & ' myeps = 2*myeps                                                         ', &
     & '// compare it to the eps used by LALA                                   ', &
     & 'if myeps=eps, ...                                                        ', &
     & '   tally=[tally,0];display(''eps matches expected value''); ...            ', &
     & 'else, ...                                                                ', &
     & '   tally=[tally,1];display(''eps is NOT the expected value''), ...         ', &
     & 'end;                                                                     ', &
     & 'if sum(tally)=0,display(''eps PASSED'');else,display(''eps FAILED'');tally   ', &
     & ''])
     !!logs=.false.
end subroutine test_eps