fury_test_qreal_sub_mixed.f90 Source File

!< FURY test of [[qreal]].
program fury_test_qreal_sub_mixed
!-----------------------------------------------------------------------------------------------------------------------------------
!< FURY test of [[qreal]].
!-----------------------------------------------------------------------------------------------------------------------------------
use fury
!-----------------------------------------------------------------------------------------------------------------------------------

!-----------------------------------------------------------------------------------------------------------------------------------
type(qreal32)  :: force1         !< A force.
type(qreal64)  :: force2         !< A force.
type(qreal128) :: force3         !< A force.
type(qreal128) :: force          !< A force.
type(uom64)    :: u_acceleration !< Acceleration unit.
type(uom64)    :: u_force        !< Force unit.
type(uom64)    :: u_length       !< Length unit.
type(uom64)    :: u_mass         !< Mass unit.
type(uom64)    :: u_speed        !< Speed unit.
type(uom64)    :: u_time         !< Time unit.
logical        :: test_passed(6) !< List of passed tests.
!-----------------------------------------------------------------------------------------------------------------------------------

!-----------------------------------------------------------------------------------------------------------------------------------
test_passed = .false.

u_length = uom64(source='m [length]', name='metre')
u_mass = uom64(source='kg [mass]', name='kilogram')
u_time = uom64(source='s [time]', name='second')

u_speed = u_length / u_time
u_acceleration = u_speed / u_time
u_force = u_mass * u_acceleration

force1 = 1 * u_force
force2 = 2 * u_force
force3 = 3 * u_force

call force%unset
force = force1 - force2
test_passed(1) = force%stringify(format='(F4.1)')=='-1.0 kg.m.s-2'
print "(A,L1)", '1 * force(float32) - 2 * force(float64) = '//force%stringify(format='(F4.1)')//', is correct? ', test_passed(1)

call force%unset
force = force2 - force1
test_passed(2) = force%stringify(format='(F3.1)')=='1.0 kg.m.s-2'
print "(A,L1)", '2 * force(float64) - 1 * force(float32) = '//force%stringify(format='(F3.1)')//', is correct? ', test_passed(2)

call force%unset
force = force1 - force3
test_passed(3) = force%stringify(format='(F4.1)')=='-2.0 kg.m.s-2'
print "(A,L1)", '1 * force(float32) - 3 * force(float128) = '//force%stringify(format='(F4.1)')//', is correct? ', test_passed(3)

call force%unset
force = force3 - force1
test_passed(4) = force%stringify(format='(F3.1)')=='2.0 kg.m.s-2'
print "(A,L1)", '3 * force(float128) - 1 * force(float32) = '//force%stringify(format='(F3.1)')//', is correct? ', test_passed(4)

call force%unset
force = force2 - force3
test_passed(5) = force%stringify(format='(F4.1)')=='-1.0 kg.m.s-2'
print "(A,L1)", '2 * force(float64) - 3 * force(float128) = '//force%stringify(format='(F4.1)')//', is correct? ', test_passed(5)

call force%unset
force = force3 - force2
test_passed(6) = force%stringify(format='(F3.1)')=='1.0 kg.m.s-2'
print "(A,L1)", '3 * force(float164) - 2 * force(float32) = '//force%stringify(format='(F3.1)')//', is correct? ', test_passed(6)

print "(A,L1)", new_line('a')//'Are all tests passed? ', all(test_passed)
stop
!-----------------------------------------------------------------------------------------------------------------------------------
endprogram fury_test_qreal_sub_mixed