1! program fortran_test 2subroutine f_main 3! A simple test of the fortran wrappers 4 5 use qdmodule 6 implicit none 7 integer*4 old_cw 8 integer i 9 type (qd_real) x, y, z 10 11 call f_fpu_fix_start (old_cw) 12 13 ! Test for read/write 14 z = "3.14159265358979323846264338327950288419716939937510582097494459230" 15 call write_scalar(6, z) 16 17 ! Test for atan/write 18 do i=1,3 19 x = qdreal(dble(i)) 20 call write_scalar(6, x) 21 y = atan(x) 22 call write_scalar(6, y) 23 end do 24 25 call write_scalar(6, nan(x)) 26 call write_scalar(6, qdcomplex(x, y)) 27 28 call f_fpu_fix_end (old_cw) 29end 30 31