dedicated to W. Van Snyder[rg]
A KISS pure Fortran Library for improving reliability of physical math computations by taking into account units of measure:
- FURY is a pure Fortran (KISS) library for improving reliability of physical math computations;
- FURY is a pure Fortran (KISS) library for symbolic algebra on units symbols;
- FURY is Fortran 2003+ standard compliant;
- FURY is OOP designed;
- FURY is a Free, Open Source Project.
use, intrinsic :: iso_fortran_env, only : real64
use fury
type(uom) :: meter
type(uom) :: second
type(qreal) :: distance_to_arrival
type(qreal) :: time_to_arrival
type(qreal) :: mean_velocity
meter = uom('m = meter = metre [length] {meter}')
second = uom('s = sec = second [time] {second}')
distance_to_arrival = qreal(100._real64, meter)
time_to_arrival = qreal(9.58_real64, second)
mean_velocity = distance_to_arrival / time_to_arrival
print "(A)", 'Bolt''s record speed: '//mean_velocity%stringify(with_dimensions=.true.)
! print
! Bolt's record speed: +0.104384133611691E+002 m.s-1 [length.time-1]
What is FURY? | Main features | Copyrights | Download | Compilation | Documentation | References
FURY is a poor-by-blow of reliability concept
- reliability is expression of the consistency-of-a-set-of-measures
- reliability measures successful probability, i.e. 1 - probability-of-failure
- reliability (absence) played a role on the Mars Climate Orbiter disastermco
ergo reliability counts
FURY project started to respond to two questions:
- is a units-consistency-check facility (UCCF) a desirable feature for a programming language?
- is useful? or
- is dangerous?
- how it can be implemented in Fortran?
- is feasible?
These questions are generated by an inspiring Fortran Language Google Group discussion started by W. Van Snyder[rg]. They are really controversial: someone want UCCF to be a built-in feature in Fortran, others claim that such a feature is not useful at all or even dangerous just because it could add more space for errors than it could reduce (thus reducing the reliability).
Starting from an agnostic point of view, the FURY case study conveys us to think that a UCCF is a desirable feature for the next Fortran standards.
Moreover, the large experience borrowed from other programming languages, see references, proves that UCCF like features can be feasible in modern programming languages and they are of proved usefulness. As a consequence, FURY started to begin a fully featured, self-contained project that is presently usable for production-applications, even if the first comprehensive testing phase (of pre v1.0.0 bug-fix one) is still not completed.
FURY is the result of the tentative to respond the above questions:
FURY is a pure Fortran library providing a new set of types (the quantities) having a built-in units-consistency-check facility that is able to improve the reliability of physical computations. Moreover, FURY can perform symbolic math on units symbols.
FURY can do reliable math computations, preserving quantities dimensionality, checking consistency and propagating multiplicative scaling factors. For example, if you want to know how fast is Bolt you can do it with FURY, in a reliable way 😄
use, intrinsic :: iso_fortran_env, only : real64
use fury
type(uom) :: meter
type(uom) :: second
type(qreal) :: distance_to_arrival
type(qreal) :: time_to_arrival
type(qreal) :: mean_velocity
meter = uom('m = meter = metre [length] {meter}')
second = uom('s = sec = second [time] {second}')
distance_to_arrival = qreal(100._real64, meter)
time_to_arrival = qreal(9.58_real64, second)
mean_velocity = distance_to_arrival / time_to_arrival
print "(A)", 'Bolt''s record speed: '//mean_velocity%stringify(with_dimensions=.true.)
! print
! Bolt's record speed: +0.104384133611691E+002 m.s-1 [length.time-1]
FURY is based on a powerful parser of unit-of-measure (UOM) definitions: once a (set of) UOM is defined its scaling factor, symbol and dimensions are consistently checked, preserved and propagated in all computations. For set of UOMs with defined conversion-formulas FURY can check the consistency of all conversions. Eventually, if a quantity has a defined UOM, any tentative to assign it to another quantity with a different UOM will raise an error as well as for all inconsistent computations.
Go to Top
FURY is inspired by the python great module pint, thus many features are taken from it. Here the main features are listed.
- User-friendly classes to add units of measure to numbers;
- errors trapping for invalid computations/assignments;
- unit parsing: prefixed forms of units are recognized without explicitly defining them, i.e. as the prefix kilo and the unit meter are defined, FURY understands kilometer;
- effortless conversion even for complex derived UOM;
- standalone unit definitions: units definitions are loaded from simple and easy to edit text file;
- advanced string parsing:
- symbolic algebra on units symbols;
- Buckingham Pi Theorem;
- replicate all the useful features of pint;
- Test Driven Developed (TDD);
- collaborative developed;
- well documented;
- free!
Any feature request is welcome.
Go to Top
FURY is an open source project, it is distributed under a multi-licensing system:
- for FOSS projects:
- for closed source/commercial projects:
Anyone is interest to use, to develop or to contribute to FURY is welcome, feel free to select the license that best matches your soul!
More details can be found on wiki.
Go to Top
FURY home is at https://github.com/szaghi/FURY. To download all the source files you can:
- clone recursively this repository:
git clone --recursive https://github.com/szaghi/FURY
- download the latest master-branch archive at https://github.com/szaghi/FURY/archive/master.zip
- download a release archive at https://github.com/szaghi/FURY/releases
Go to Top
FURY is a modern Fortran project thus a modern Fortran compiler is need to compile the project.
The library is modular, namely it exploits Fortran modules. As a consequence, there is compilation-cascade hierarchy to build the library. To correctly build the library the following approaches are supported
- Build by means of FoBiS: full support;
- Build by means of GNU Make: to be implemented.
- Build by means of CMake: to be implemented.
The FoBiS building support is the most complete, as it is the one used for the developing FURY.
A fobos
file is provided to build the library by means of the Fortran Building System FoBiS.
Type
FoBiS.py build
After (a successful) building a directory ./exe
is created containing all the compiled tests that constitute the FURY regression-tests-suite, e.g.
→ FoBiS.py build
Builder options
Directories
Building directory: "exe"
Compiled-objects .o directory: "exe/obj"
Compiled-objects .mod directory: "exe/mod"
Compiler options
Vendor: "gnu"
Compiler command: "gfortran"
Module directory switch: "-J"
Compiling flags: "-c -frealloc-lhs -std=f2008 -fall-intrinsics -O2 -Dr16p"
Linking flags: "-O2"
Preprocessing flags: "-Dr16p"
Coverage: False
Profile: False
PreForM.py used: False
PreForM.py output directory: None
PreForM.py extensions processed: []
Building src/tests/basic_use.f90
Compiling src/lib/penf.F90 serially
Compiling src/lib/string_t.F90 serially
Compiling src/lib/stringifor.F90 serially
Compiling ...
Linking exe/basic_use
Target src/tests/basic_use.f90 has been successfully built
...
→ tree -L 1 exe/
exe/
├── basic_use
...
To be implemented.
To be implemented.
Go to Top
Besides this README file the FURY documentation is contained into its own wiki. Detailed documentation of the API is contained into the GitHub Pages that can also be created locally by means of ford tool.
FURY is designed upon a simple, yet powerful UOM grammar, for example the definition of the Pascal unit could be one the followings:
'Pa'
'Pa [pressure]'
'Pa [pressure] {pascal}'
'kg.m-1.s-2'
'kg.m-1.s-2 = Pa'
- ...
'kg [mass].m-1 [length-1].s-2 [time-2] (Pa[pressure]) {pascal}'
where the first is the most basic definition defining only the main symbol while the others are definitions with increasing complexity where we are defining dimensions, symbol aliases, name, definition with respect other units... The flexibility left to the end user is really high. More in details, the UOM grammar is composed as following:
'kg [mass].m-1 = metre-1 = meter-1 [length-1].s-2 = seconds-2 = sec-2 [time-2] (Pa[pressure]) {pascal}'
where
- the terms
[...]
define the dimensions of each unit and are optional (the white spaces are ignored); moreover the exponents of dimensions can be omitted: in this case they are inferred from the symbol reference exponents; in the case they are explicitly written they must match the corresponding symbol ones or an error is raised; - the term
(...)
defines a main unit alias that is optional and must come always after unit reference definition; - the term
{...}
defines the unit name that is optional and must be always the last term.
Reference units used for the definition of another one are separated by the .
symbol, thus the Pascal definition is kg/(m * s^2)
is defined as kg.m-1.s-2
to mimic the rules of The Unified Code for Units of Measure: each reference unit is intended multiplied the subsequent with the proper exponent with sign. The FURY parser is designed on this grammar, thus each reference unit is tokenized with the .
separator.
Symbol aliases are used also for defining units conversion formulas, e.g.
'km = 1000 * m [length].h-1 = 3600 s-1 [time-1] (km/h[speed]) {km/h}'
This is the base, more concise conversion formula grammar, but it could be unsafe if the conversion implies real-valued scaling factor, as it often happens. In this case, a more verbose but more safe grammar should be used, e.g.
'km< = 1000._real64 * m > [length].h-1< = 3600._real64 s-1> [time-1] (km/h[speed]) {km/h}'
where the conversion formulas are protected by the <>
symbols. Defining a UOM with fury is simple as
use fury
type(uom) :: pascal
pascal = uom('kg [mass].m-1 [length-1].s-2 [time-2] (Pa[pressure]) {pascal}')
Having a unit so defined, it is allowed to perform your reliable computations, namely all the symbolic computations necessary when dealing with units, e.g. computing the force from a pressure on a surface
use fury
type(uom) :: meter
type(uom) :: pascal
type(uom) :: newton
meter = uom(m [length] (meter[length]) {meter}')
pascal = uom('kg [mass].m-1 [length-1].s-2 [time-2] (Pa[pressure]) {pascal}')
newton = pascal * meter**2
print "(A)", ' newton: '//newton%stringify(with_dimensions=.true.)
! output
! newton: kg.m.s-2 [mass.length.time-2]
For make easier the life of snappy coders a far complete SI system of predefined units is provided, as well an abstract class for defining new units systems.
We indicate a UOM to be a reference unit when it is defined by only main symbol having exponent 1 and multiplicative factor 1. With this categorization the base units of SI are all reference units for us. Moreover, the user can define his/her own reference units, e.g. apple [fruit] {not a pc}
is a valid UOM reference definition that can be used to define other (complex) units like GNU< = 1.e20 * apple> [salad_of_fruits] {HPC}
. To define a reference UOM you need to specify just the main symbol, but other data can be optionally provided to cover more complex purposes. For example, defining a reference unit you can provide aliases defining conversion formulas to other UOM, you can provide the dimensions or UOM name useful to query systems of units.
required | alias formula(s), optional | optional | optional | optional | ||||||
optional | = | optional if == 1 | Hz | optional if == 1 | optional | |||||
grammar_element | main_symbol | < | = | multiplicative_factor | symbol | symbol_exponent | > | [dimensions] | (alias[dimensions]) | {unit_name} |
example | kHz | < | = | 1000. * | Hz | 1 | > | [time-1] | (kilohertz[frequency]) | {kilohertz} |
Examples
kHz
, only the main symbol is strictly required;kHz = 1000 * Hz
, one alias (for conversion formula);kHz = 1000 * Hz = 1000 * s-1
, two aliases (for conversion formulas);kHz< = 1000. * Hz>
, one alias protected (due to real factor);kHz< = 1000. * Hz = 1000. * s-1 >
, two aliases protected;kHz [time-1]
, dimensions specification;kHz [time-1] (kilohertz[frequency])
, a main alias specification, useful for units like pressure:kg.m-1.s-2 (Pa[pressure]) {pascal}
kHz [time-1] (kilohertz[frequency]) {my kilohertz}
, unit name specification, useful for queries in units systems usage;kHz< = 1000. * Hz = 1000. * s-1 > [frequency] (kilohertz[frequency]) {kilohertz}
In the case you want to define a UOM with a respect a set of reference UOM, you can do using the .
operator that mimics the definition of The Unified Code for Units of Measure group.
required at least 2 UOM definitions | optional | optional | |||
grammar_element | UOMreference | . | otherUOMreference | (alias[dimensions]) | {unit_name} |
example | m [length] | . | s-1 [time-1] | (m/s [speed]) | {speed SI} |
Examples
kg [mass].m [length].s-2 [time-2] (N[force]) {newton}
kg [mass].m2 [length2].s-3 [time-3].A-2 [current-2]{ohm}
kg [mass].m-1 [length-1].s-2 [time-2] (Pa[pressure]) {pascal}
It worth to note that the user can always define a complex unit as a reference one, let us consider for example two possible definitions of the pressure:
kg [mass].m-1 [length-1].s-2 [time-2] (Pa[pressure]) {pascal_complex}
Pa [pressure] {pascal_simple}
Both are valid, but the first one provides much more informations thus resulting into a more powerful symbolic math. For example consider the following example where the above 2 definitions are used to define the Newton unit:
use fury
type(uom) :: meter
type(uom) :: pascal_complex
type(uom) :: pascal_simple
type(uom) :: newton
meter = uom('m [length] (meter[length]) {meter}')
pascal_complex = uom('kg [mass].m-1 [length-1].s-2 [time-2] (Pa[pressure]) {pascal}')
pascal_simple = uom('Pa[pressure] {pascal}')
newton = pascal_complex * meter**2
print "(A)", ' newton: '//newton%stringify(with_dimensions=.true.)
! output
! newton: kg.m.s-2 [mass.length.time-2]
call newton%unset
newton = pascal_simple * meter**2
print "(A)", ' newton: '//newton%stringify(with_dimensions=.true.)
! output
! newton: Pa.m2 [pressure.length2]
Both results are correct, but the first is more informative if a consistent set of reference units has been defined, as it happens for the SI system.
FURY's quantities are built up upon the above UOM: a quantity is essentially a number (real or integer of any kind) with attached a UOM. When a quantity has a defined UOM all UOM's check are pushed to the quantity's magnitude computations, thus improving the math computations reliability.
Unfortunately, we have not find a simple way to attach UOM to standard Fortran types, thus the end user must do reliable computations on a new set of types, namely qreal
and qinteger
(or their kinded version).
Let us assume to consider the summation of 2 length, the first expressed in meters and the second in kilometers. With FURY this can be implemented as:
use fury
type(system_si) :: SI ! SI system.
type(qreal) :: q ! A physical quantity.
type(qreal) :: q1 ! A physical quantity.
type(qreal) :: q2 ! A physical quantity.
type(uom) :: kilometre ! The kilometre unit based on the metre unit.
! initialize SI units system
call SI%initialize
! define a new unit, i.e. the kilometre, that must
! be compatible with SI's metre
kilometre = uom('km = 1000 * m')
! define 2 quantities
q1 = qreal(magnitude=1._R_P, unit=SI%unit('metre'))
q2 = qreal(magnitude=1._R_P, unit=kilometre)
! try to blindly to sum them
q = q1 + q2
! error: cannot convert from "metre [length]" to "km"
! a more conscious sum
q = q1 + q2%to(SI%unit('metre'))
! no error is raised
print "(A)", '1 m + 1 km = '//q%stringify(format='(F6.1)')
! output "1 m + 1 km = 1001.0 m"
To use FURY the minimal steps necessary are:
use fury
type(qreal) :: q1 ! A real physical quantity.
type(qinteger) :: q2 ! An integer physical quantity.
type(uom) :: centimetre ! The centimetre unit based on the metre unit.
type(uom) :: kilometre ! The kilometre unit based on the metre unit.
type(uom) :: hour ! The hour unit based on the second unit.
centimetre = uom('cm = 0.01 * m')
kilometre = uom('km = 1000.0 * m')
hour = uom('h = 3600.0 * s')
type(system_si) :: SI ! SI system
call SI%initialize
Note that currently on the SI system is provided.
q1 = qreal(magnitude=1._R_P, unit=kilometre)
q2 = qinteger(magnitude=3, unit=hour)
Now you can do your reliable physical computations :party:
Go to Top
A non comprehensive references list of UCCF related contents.
- The Unified Code for Units of Measure
- Meyers, Dimensional analysis. C++ Report
- Petty, Automated computation and consistency checking of physical dimensions and units in scientific programs
- Brown, Introduction to the SI Library of Unit-Based Computation
- International System of Units
- .NET
- Ada
- C#
- C++
- Haskell
- JavaScript
- Julia
- Python
- Ruby
The Mars Climate Orbiter (MCO) Mishap Investigation Board (MIB) has determined that the root cause for the loss of the MCO spacecraft was the failure to use metric units in the coding of a ground software file, “Small Forces” used in trajectory models, specifically:
thruster performance data in English units instead of metric units was used in the software application code titled
SM_FORCES
(small forces).
The units mismatching generates an erroneous trajectory that was underestimated by a factor of 4.45, which is the required conversion factor from force in pounds to Newtons.
if the software application
SM_FORCES
had been developed with a units-consistency-check facility, $300 million of spacecraft could be saved.
Go to Top