Path: | ZolotarevModule/Errorbound.F90 |
Last Update: | Wed Jan 18 16:14:48 +0900 2012 |
Main Program : |
program errorbound use signfunc_mod implicit none integer, parameter :: QP=KIND(1.0q0) integer :: nss,j real(QP) :: kappa,dkappa,eps do j=6,14,2 eps=10.0_QP**(-j) kappa = 1.0_QP/1000 dkappa = exp(1.0_QP/1000) do while(kappa<1.0_QP) nss = get_ns_exact(kappa,eps) write(*,'(2ES24.15,2I3)')eps,kappa,nss kappa=kappa*dkappa enddo write(*,*) write(*,*) enddo stop end program