123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542 |
- // Boost.Units - A C++ library for zero-overhead dimensional analysis and
- // unit/quantity manipulation and conversion
- //
- // Copyright (C) 2003-2008 Matthias Christian Schabel
- // Copyright (C) 2008 Steven Watanabe
- //
- // Distributed under the Boost Software License, Version 1.0. (See
- // accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- /**
- \file
-
- \brief kitchen_sink.cpp
- \details
- More extensive quantity tests.
- Output:
- @verbatim
- //[kitchen_sink_output_1
- S1 : 2
- X1 : 2
- X2 : (4/3)
- U1 : N
- U2 : J
- Q1 : 1 N
- Q2 : 2 J
- //]
- //[kitchen_sink_output_2
- U1*S1 : 2 N
- S1*U1 : 2 N
- U1/S1 : 0.5 N
- S1/U1 : 2 m^-1 kg^-1 s^2
- //]
- //[kitchen_sink_output_3
- U1+U1 : N
- U1-U1 : N
- U1*U1 : m^2 kg^2 s^-4
- U1/U1 : dimensionless
- U1*U2 : m^3 kg^2 s^-4
- U1/U2 : m^-1
- U1^X : m^2 kg^2 s^-4
- X1vU1 : m^(1/2) kg^(1/2) s^-1
- U1^X2 : m^(4/3) kg^(4/3) s^(-8/3)
- X2vU1 : m^(3/4) kg^(3/4) s^(-3/2)
- //]
- //[kitchen_sink_output_4
- Q1*S1 : 2 N
- S1*Q1 : 2 N
- Q1/S1 : 0.5 N
- S1/Q1 : 2 m^-1 kg^-1 s^2
- //]
- //[kitchen_sink_output_5
- U1*Q1 : 1 m^2 kg^2 s^-4
- Q1*U1 : 1 m^2 kg^2 s^-4
- U1/Q1 : 1 dimensionless
- Q1/U1 : 1 dimensionless
- //]
- //[kitchen_sink_output_6
- +Q1 : 1 N
- -Q1 : -1 N
- Q1+Q1 : 2 N
- Q1-Q1 : 0 N
- Q1*Q1 : 1 m^2 kg^2 s^-4
- Q1/Q1 : 1 dimensionless
- Q1*Q2 : 2 m^3 kg^2 s^-4
- Q1/Q2 : 0.5 m^-1
- Q1^X1 : 1 m^2 kg^2 s^-4
- X1vQ1 : 1 m^(1/2) kg^(1/2) s^-1
- Q1^X2 : 1 m^(4/3) kg^(4/3) s^(-8/3)
- X2vQ1 : 1 m^(3/4) kg^(3/4) s^(-3/2)
- //]
- //[kitchen_sink_output_7
- l1 == l2 false
- l1 != l2 true
- l1 <= l2 true
- l1 < l2 true
- l1 >= l2 false
- l1 > l2 false
- //]
- dimless = 1
- //[kitchen_sink_output_8
- v1 = 2 m s^-1
- //]
- //[kitchen_sink_output_9
- F = 1 N
- dx = 1 m
- E = 1 J
- //]
- //[kitchen_sink_output_10
- r = 5e-07 m
- P = 101325 Pa
- V = 5.23599e-19 m^3
- T = 310 K
- n = 2.05835e-17 mol
- R = 8.314472 m^2 kg s^-2 K^-1 mol^-1 (rel. unc. = 1.8e-06)
- //]
- //[kitchen_sink_output_11
- theta = 0.375 rd
- sin(theta) = 0.366273 dimensionless
- asin(sin(theta)) = 0.375 rd
- //]
- //[kitchen_sink_output_12
- V = (12.5,0) V
- I = (3,4) A
- Z = (1.5,-2) Ohm
- I*Z = (12.5,0) V
- //]
- //[kitchen_sink_output_13
- x+y-w = 0.48(+/-0.632772) m
- w*x = 9.04(+/-0.904885) m^2
- x/y = 0.666667(+/-0.149071) dimensionless
- //]
- //[kitchen_sink_output_14
- w*y^2/(u*x)^2 = 10.17(+/-3.52328) m^-1
- w/(u*x)^(1/2) = 3.19612(+/-0.160431) dimensionless
- //]
- //[kitchen_sink_output_15
- I*w = m^2 kg s^-1 rad^-1
- I*w/L = dimensionless
- I*w^2 = J
- //]
- //[kitchen_sink_output_16
- 1 F
- 1 kat
- 1 S
- 1 C
- 1 V
- 1 J
- 1 N
- 1 Hz
- 1 lx
- 1 H
- 1 lm
- 1 Wb
- 1 T
- 1 W
- 1 Pa
- 1 Ohm
- //]
- //[kitchen_sink_output_18
- 1 farad
- 1 katal
- 1 siemen
- 1 coulomb
- 1 volt
- 1 joule
- 1 newton
- 1 hertz
- 1 lux
- 1 henry
- 1 lumen
- 1 weber
- 1 tesla
- 1 watt
- 1 pascal
- 1 ohm
- //]
- @endverbatim
- **/
- #include <cmath>
- #include <complex>
- #include <iostream>
- #include <boost/typeof/std/complex.hpp>
- #include <boost/units/cmath.hpp>
- #include <boost/units/io.hpp>
- #include <boost/units/systems/si.hpp>
- #include <boost/units/systems/si/codata/physico-chemical_constants.hpp>
- #include <boost/units/systems/si/io.hpp>
- #include "measurement.hpp"
- namespace boost {
- namespace units {
- //[kitchen_sink_function_snippet_3
- /// the physical definition of work - computed for an arbitrary unit system
- template<class System,class Y>
- constexpr
- quantity<unit<energy_dimension,System>,Y>
- work(quantity<unit<force_dimension,System>,Y> F,
- quantity<unit<length_dimension,System>,Y> dx)
- {
- return F*dx;
- }
- //]
- //[kitchen_sink_function_snippet_4
- /// the ideal gas law in si units
- template<class Y>
- constexpr
- quantity<si::amount,Y>
- idealGasLaw(const quantity<si::pressure,Y>& P,
- const quantity<si::volume,Y>& V,
- const quantity<si::temperature,Y>& T)
- {
- using namespace boost::units::si;
-
- using namespace constants::codata;
- return (P*V/(R*T));
- }
- //]
- } // namespace units
- } // namespace boost
- int main()
- {
- using namespace boost::units;
- using namespace boost::units::si;
- {
- //[kitchen_sink_snippet_1
- /// scalar
- const double s1 = 2;
-
- const long x1 = 2;
- const static_rational<4,3> x2;
-
- /// define some units
- force u1 = newton;
- energy u2 = joule;
-
- /// define some quantities
- quantity<force> q1(1.0*u1);
- quantity<energy> q2(2.0*u2);
- //]
-
- /// check scalar, unit, and quantity io
- std::cout << "S1 : " << s1 << std::endl
- << "X1 : " << x1 << std::endl
- << "X2 : " << x2 << std::endl
- << "U1 : " << u1 << std::endl
- << "U2 : " << u2 << std::endl
- << "Q1 : " << q1 << std::endl
- << "Q2 : " << q2 << std::endl
- << std::endl;
-
- /// check scalar-unit algebra
- std::cout //<< "U1+S1 : " << u1+s1 << std::endl // illegal
- //<< "S1+U1 : " << s1+u1 << std::endl // illegal
- //<< "U1-S1 : " << u1-s1 << std::endl // illegal
- //<< "S1-U1 : " << s1-u1 << std::endl // illegal
- << "U1*S1 : " << u1*s1 << std::endl
- << "S1*U1 : " << s1*u1 << std::endl
- << "U1/S1 : " << u1/s1 << std::endl
- << "S1/U1 : " << s1/u1 << std::endl
- << std::endl;
-
- /// check unit-unit algebra
- std::cout << "U1+U1 : " << u1+u1 << std::endl
- << "U1-U1 : " << u1-u1 << std::endl
- << "U1*U1 : " << u1*u1 << std::endl
- << "U1/U1 : " << u1/u1 << std::endl
- //<< "U1+U2 : " << u1+u2 << std::endl // illegal
- //<< "U1-U2 : " << u1-u2 << std::endl // illegal
- << "U1*U2 : " << u1*u2 << std::endl
- << "U1/U2 : " << u1/u2 << std::endl
- << "U1^X : " << pow<2>(u1) << std::endl
- << "X1vU1 : " << root<2>(u1) << std::endl
- << "U1^X2 : " << pow<static_rational<4,3> >(u1) << std::endl
- << "X2vU1 : " << root<static_rational<4,3> >(u1) << std::endl
- << std::endl;
-
- /// check scalar-quantity algebra
- std::cout //<< "Q1+S1 : " << q1+s1 << std::endl // illegal
- //<< "S1+Q1 : " << s1+q1 << std::endl // illegal
- //<< "Q1-S1 : " << q1-s1 << std::endl // illegal
- //<< "S1-Q1 : " << s1-q1 << std::endl // illegal
- << "Q1*S1 : " << q1*s1 << std::endl
- << "S1*Q1 : " << s1*q1 << std::endl
- << "Q1/S1 : " << q1/s1 << std::endl
- << "S1/Q1 : " << s1/q1 << std::endl
- << std::endl;
-
- /// check unit-quantity algebra
- std::cout //<< "U1+Q1 : " << u1+q1 << std::endl // illegal
- //<< "Q1+U1 : " << q1+u1 << std::endl // illegal
- //<< "U1-Q1 : " << u1-q1 << std::endl // illegal
- //<< "Q1-U1 : " << q1-u1 << std::endl // illegal
- << "U1*Q1 : " << u1*q1 << std::endl
- << "Q1*U1 : " << q1*u1 << std::endl
- << "U1/Q1 : " << u1/q1 << std::endl
- << "Q1/U1 : " << q1/u1 << std::endl
- << std::endl;
-
- /// check quantity-quantity algebra
- std::cout << "+Q1 : " << +q1 << std::endl
- << "-Q1 : " << -q1 << std::endl
- << "Q1+Q1 : " << q1+q1 << std::endl
- << "Q1-Q1 : " << q1-q1 << std::endl
- << "Q1*Q1 : " << q1*q1 << std::endl
- << "Q1/Q1 : " << q1/q1 << std::endl
- //<< "Q1+Q2 : " << q1+q2 << std::endl // illegal
- //<< "Q1-Q2 : " << q1-q2 << std::endl // illegal
- << "Q1*Q2 : " << q1*q2 << std::endl
- << "Q1/Q2 : " << q1/q2 << std::endl
- << "Q1^X1 : " << pow<2>(q1) << std::endl
- << "X1vQ1 : " << root<2>(q1) << std::endl
- << "Q1^X2 : " << pow<static_rational<4,3> >(q1) << std::endl
- << "X2vQ1 : " << root<static_rational<4,3> >(q1) << std::endl
- << std::endl;
-
- //[kitchen_sink_snippet_2
- /// check comparison tests
- quantity<length> l1(1.0*meter),
- l2(2.0*meters);
- //]
-
- std::cout << std::boolalpha
- << "l1 == l2" << "\t" << (l1 == l2) << std::endl
- << "l1 != l2" << "\t" << (l1 != l2) << std::endl
- << "l1 <= l2" << "\t" << (l1 <= l2) << std::endl
- << "l1 < l2 " << "\t" << (l1 < l2) << std::endl
- << "l1 >= l2" << "\t" << (l1 >= l2) << std::endl
- << "l1 > l2 " << "\t" << (l1 > l2) << std::endl
- << std::endl;
-
- //[kitchen_sink_snippet_3
- /// check implicit unit conversion from dimensionless to value_type
- const double dimless = (q1/q1);
- //]
-
- std::cout << "dimless = " << dimless << std::endl
- << std::endl;
-
- quantity<velocity> v1 = 2.0*meters/second;
-
- std::cout << "v1 = " << v1 << std::endl
- << std::endl;
-
- //[kitchen_sink_snippet_4
- /// test calcuation of work
- quantity<force> F(1.0*newton);
- quantity<length> dx(1.0*meter);
- quantity<energy> E(work(F,dx));
- //]
-
- std::cout << "F = " << F << std::endl
- << "dx = " << dx << std::endl
- << "E = " << E << std::endl
- << std::endl;
-
- {
- //[kitchen_sink_snippet_5
- /// test ideal gas law
- quantity<temperature> T = (273.+37.)*kelvin;
- quantity<pressure> P = 1.01325e5*pascals;
- quantity<length> r = 0.5e-6*meters;
- quantity<volume> V = (4.0/3.0)*3.141592*pow<3>(r);
- quantity<amount> n(idealGasLaw(P,V,T));
- //]
-
- std::cout << "r = " << r << std::endl
- << "P = " << P << std::endl
- << "V = " << V << std::endl
- << "T = " << T << std::endl
- << "n = " << n << std::endl
- #if BOOST_UNITS_HAS_TYPEOF
- << "R = " << constants::codata::R << std::endl
- #else
- << "no typeof" << std::endl
- #endif // BOOST_UNITS_HAS_TYPEOF
- << std::endl;
- }
-
- //[kitchen_sink_snippet_6
- /// test trig stuff
- quantity<plane_angle> theta = 0.375*radians;
- quantity<dimensionless> sin_theta = sin(theta);
- quantity<plane_angle> thetap = asin(sin_theta);
- //]
-
- std::cout << "theta = " << theta << std::endl
- << "sin(theta) = " << sin_theta << std::endl
- << "asin(sin(theta)) = " << thetap << std::endl
- << std::endl;
-
- /// test implicit conversion of dimensionless to value
- double tmp = sin_theta;
-
- tmp = sin_theta;
-
- /// test implicit conversion from value to dimensionless
- quantity<dimensionless> tmpp = tmp;
-
- tmpp = tmp;
-
- /// check complex quantities
- typedef std::complex<double> complex_type;
-
- //[kitchen_sink_snippet_7
- quantity<electric_potential,complex_type> v = complex_type(12.5,0.0)*volts;
- quantity<current,complex_type> i = complex_type(3.0,4.0)*amperes;
- quantity<resistance,complex_type> z = complex_type(1.5,-2.0)*ohms;
- //]
-
- std::cout << "V = " << v << std::endl
- << "I = " << i << std::endl
- << "Z = " << z << std::endl
- << "I*Z = " << i*z << std::endl
- << std::endl;
-
- /// check quantities using user-defined type encapsulating error propagation
- //[kitchen_sink_snippet_8
- quantity<length,measurement<double> >
- u(measurement<double>(1.0,0.0)*meters),
- w(measurement<double>(4.52,0.02)*meters),
- x(measurement<double>(2.0,0.2)*meters),
- y(measurement<double>(3.0,0.6)*meters);
- //]
-
- std::cout << "x+y-w = " << x+y-w << std::endl
- << "w*x = " << w*x << std::endl
- << "x/y = " << x/y << std::endl
- << "w*y^2/(u*x)^2 = " << w*y*y/pow<2>(u*x) << std::endl
- << "w/(u*x)^(1/2) = " << w/pow< static_rational<1,2> >(u*x)
- << std::endl << std::endl;
- }
-
- /// check moment of inertia/angular momentum/rotational energy
-
- //[kitchen_sink_snippet_9
- std::cout << symbol_format
- << "I*w = " << moment_of_inertia()*angular_velocity() << std::endl
- << "I*w/L = " << moment_of_inertia()*angular_velocity()/angular_momentum() << std::endl
- << "I*w^2 = " << moment_of_inertia()*pow<2>(angular_velocity()) << std::endl
- << std::endl;
- //]
-
- //[kitchen_sink_snippet_10
- // std::cout << typename_format
- // << quantity<capacitance>(1.0*farad) << std::endl
- // << quantity<catalytic_activity>(1.0*katal) << std::endl
- // << quantity<conductance>(1.0*siemen) << std::endl
- // << quantity<electric_charge>(1.0*coulomb) << std::endl
- // << quantity<electric_potential>(1.0*volt) << std::endl
- // << quantity<energy>(1.0*joule) << std::endl
- // << quantity<force>(1.0*newton) << std::endl
- // << quantity<frequency>(1.0*hertz) << std::endl
- // << quantity<illuminance>(1.0*lux) << std::endl
- // << quantity<inductance>(1.0*henry) << std::endl
- // << quantity<luminous_flux>(1.0*lumen) << std::endl
- // << quantity<magnetic_flux>(1.0*weber) << std::endl
- // << quantity<magnetic_flux_density>(1.0*tesla) << std::endl
- // << quantity<power>(1.0*watt) << std::endl
- // << quantity<pressure>(1.0*pascals) << std::endl
- // << quantity<resistance>(1.0*ohm) << std::endl
- // << std::endl;
- //]
-
- //[kitchen_sink_snippet_11
- // std::cout << raw_format
- // << quantity<capacitance>(1.0*farad) << std::endl
- // << quantity<catalytic_activity>(1.0*katal) << std::endl
- // << quantity<conductance>(1.0*siemen) << std::endl
- // << quantity<electric_charge>(1.0*coulomb) << std::endl
- // << quantity<electric_potential>(1.0*volt) << std::endl
- // << quantity<energy>(1.0*joule) << std::endl
- // << quantity<force>(1.0*newton) << std::endl
- // << quantity<frequency>(1.0*hertz) << std::endl
- // << quantity<illuminance>(1.0*lux) << std::endl
- // << quantity<inductance>(1.0*henry) << std::endl
- // << quantity<luminous_flux>(1.0*lumen) << std::endl
- // << quantity<magnetic_flux>(1.0*weber) << std::endl
- // << quantity<magnetic_flux_density>(1.0*tesla) << std::endl
- // << quantity<power>(1.0*watt) << std::endl
- // << quantity<pressure>(1.0*pascals) << std::endl
- // << quantity<resistance>(1.0*ohm) << std::endl
- // << std::endl;
- //]
-
- //[kitchen_sink_snippet_12
- std::cout << symbol_format
- << quantity<capacitance>(1.0*farad) << std::endl
- << quantity<catalytic_activity>(1.0*katal) << std::endl
- << quantity<conductance>(1.0*siemen) << std::endl
- << quantity<electric_charge>(1.0*coulomb) << std::endl
- << quantity<electric_potential>(1.0*volt) << std::endl
- << quantity<energy>(1.0*joule) << std::endl
- << quantity<force>(1.0*newton) << std::endl
- << quantity<frequency>(1.0*hertz) << std::endl
- << quantity<illuminance>(1.0*lux) << std::endl
- << quantity<inductance>(1.0*henry) << std::endl
- << quantity<luminous_flux>(1.0*lumen) << std::endl
- << quantity<magnetic_flux>(1.0*weber) << std::endl
- << quantity<magnetic_flux_density>(1.0*tesla) << std::endl
- << quantity<power>(1.0*watt) << std::endl
- << quantity<pressure>(1.0*pascals) << std::endl
- << quantity<resistance>(1.0*ohm) << std::endl
- << std::endl;
- //]
-
- //[kitchen_sink_snippet_13
- std::cout << name_format
- << quantity<capacitance>(1.0*farad) << std::endl
- << quantity<catalytic_activity>(1.0*katal) << std::endl
- << quantity<conductance>(1.0*siemen) << std::endl
- << quantity<electric_charge>(1.0*coulomb) << std::endl
- << quantity<electric_potential>(1.0*volt) << std::endl
- << quantity<energy>(1.0*joule) << std::endl
- << quantity<force>(1.0*newton) << std::endl
- << quantity<frequency>(1.0*hertz) << std::endl
- << quantity<illuminance>(1.0*lux) << std::endl
- << quantity<inductance>(1.0*henry) << std::endl
- << quantity<luminous_flux>(1.0*lumen) << std::endl
- << quantity<magnetic_flux>(1.0*weber) << std::endl
- << quantity<magnetic_flux_density>(1.0*tesla) << std::endl
- << quantity<power>(1.0*watt) << std::endl
- << quantity<pressure>(1.0*pascals) << std::endl
- << quantity<resistance>(1.0*ohm) << std::endl
- << std::endl;
- //]
- return 0;
- }
|