Go to the documentation of this file.00001
00003
00016 module mod_maths
00017 #include "fml_constants.h"
00018
00019 implicit none
00020
00023 character(len=len_what_exception) :: math_what_exception
00024 CONTAINS
00025
00032 function math_machine_eps() result(machine_eps)
00033
00034 type_precision :: machine_eps, eps_test
00035 machine_eps=1.0
00036 eps_test = 1.0 + machine_eps
00037 do while(1.0.ne.eps_test)
00038 machine_eps=machine_eps/2.0
00039 eps_test = 1.0 + machine_eps
00040 end do
00041 end function math_machine_eps
00042
00043
00044
00045 end module mod_maths