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