The Meson Build System
http://mesonbuild.com/
You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
34 lines
648 B
34 lines
648 B
5 years ago
|
use, intrinsic :: iso_fortran_env, only : error_unit, sp=>real32, dp=>real64
|
||
|
|
implicit none
|
||
|
|
|
||
|
|
real(sp) :: a32
|
||
|
|
real(dp) :: a64
|
||
|
|
|
||
|
|
real(sp), parameter :: pi32 = 4*atan(1._sp)
|
||
|
|
real(dp), parameter :: pi64 = 4*atan(1._dp)
|
||
|
|
|
||
|
|
if (pi32 == pi64) error stop 'real32 values generally do not exactly equal real64 values'
|
||
|
|
|
||
|
|
call timestwo(a32)
|
||
|
|
call timestwo(a64)
|
||
|
|
|
||
|
|
contains
|
||
|
|
|
||
|
|
elemental subroutine timestwo(a)
|
||
|
|
|
||
|
|
class(*), intent(inout) :: a
|
||
|
|
|
||
|
|
select type (a)
|
||
|
|
type is (real(sp))
|
||
|
|
a = 2*a
|
||
|
|
type is (real(dp))
|
||
|
|
a = 2*a
|
||
|
|
type is (integer)
|
||
|
|
a = 2*a
|
||
|
|
class default
|
||
|
|
error stop 'requires real32, real64 or integer'
|
||
|
|
end select
|
||
|
|
|
||
|
|
end subroutine timestwo
|
||
|
|
|
||
|
|
end program
|