
##############################################################################
#
# Copyright (c) 2003-2018 by The University of Queensland
# http://www.uq.edu.au
#
# Primary Business: Queensland, Australia
# Licensed under the Apache License, version 2.0
# http://www.apache.org/licenses/LICENSE-2.0
#
# Development until 2012 by Earth Systems Science Computational Center (ESSCC)
# Development 2012-2013 by School of Earth Sciences
# Development from 2014 by Centre for Geoscience Computing (GeoComp)
#
##############################################################################

from __future__ import print_function, division

__copyright__="""Copyright (c) 2003-2018 by The University of Queensland
http://www.uq.edu.au
Primary Business: Queensland, Australia"""
__license__="""Licensed under the Apache License, version 2.0
http://www.apache.org/licenses/LICENSE-2.0"""
__url__="https://launchpad.net/escript-finley"

import esys.escriptcore.utestselect as unittest
from esys.escriptcore.testing import *
import sys
from esys.escript.unitsSI import *
from esys.escript.util import EPSILON

class Test_UnitsSI(unittest.TestCase):
    TOL=EPSILON*100.
    def testUnit(self):
       s=Unit("s","something",1.,4.)
       self.assertTrue(s.getName() == "s", "wrong name")
       s.setName("t")
       self.assertTrue(s.getName() == "t", "wrong reset name")
       self.assertTrue(s.getLongName() == "something", "wrong long name")
       s.setLongName("test")
       self.assertTrue(s.getLongName() == "test", "wrong reset long name")
       self.assertTrue(abs( s(3.)-13.) <= self.TOL*13. , " s(3.) wrong.")
       self.assertTrue(abs( 3.*s-13.) <= self.TOL*13. , " 13.s wrong.")
       self.assertTrue(abs( 10.*s-41.) <= self.TOL*41. , " 10.s wrong.")
       self.assertTrue(abs( 41./s-10.) <= self.TOL*10. , " 10/s wrong.")
       self.assertTrue(abs( 13./s-3.) <= self.TOL*3. , " 13/s wrong.")

       x=Unit("x","X",10.,2.)
       y=Unit("y","Y",0.,4.)
       z=Unit("z","Z",0.,0.5)
       p_m=x*s
       self.assertTrue(p_m.getName() == "xt", "wrong p_m name")
       self.assertTrue(p_m.getLongName() == "X*test", "wrong p_m long name")
       self.assertTrue(abs(p_m(0.)-12.) <= self.TOL*12 , " p_m(0) wrong.")
       self.assertTrue(abs(p_m(1.)-20.) <= self.TOL*20 , " p_m(1) wrong.")

       p_d=y/z
       self.assertTrue(p_d.getName() == "y/z", "wrong p_d name")
       self.assertTrue(p_d.getLongName() == "Y/Z", "wrong p_d long name")
       self.assertTrue(abs(p_d(1.)-8.) <= self.TOL*8. , " p_d(0) wrong.")
       self.assertTrue(abs(p_d(2.)-16.) <= self.TOL*16. , " p_d(1) wrong.")

       p_p=y**3
       self.assertTrue(p_p.getName() == "y^3", "p_p mult name")
       self.assertTrue(p_p.getLongName() == "Y^3", "wrong p_p long name")
       self.assertTrue(abs(p_p(1.)-64.) <= self.TOL*64. , " p_p(0) wrong.")
       self.assertTrue(abs(p_p(.5)-32) <= self.TOL*32. , " p_p(1) wrong.")

        
       x=Unit("x","X",0.,2.)
       y=Unit("y","Y",0.,4.)
       z=Unit("z","Z",0.,0.5)

       self.assertTrue(((x*x)).getName() == "xx", "wrong name")
       self.assertTrue((x*(x*z)).getName() == "xxz", "wrong name")
       self.assertTrue((x*(x/z)).getName() == "xx/z", "wrong name")
       self.assertTrue((x*(y**3)).getName() == "x y^3", "wrong name")

       self.assertTrue(((x*z)*x).getName() == "xzx", "wrong name")
       self.assertTrue(((x*z)*(x*z)).getName() == "xzxz", "wrong name")
       self.assertTrue(((x*z)*(x/z)).getName() == "xzx/z", "wrong name")
       self.assertTrue(((x*z)*(y**3)).getName() == "xz y^3", "wrong name")

       self.assertTrue(((x/z)*x).getName() == "x/z x", "wrong name")
       self.assertTrue(((x/z)*(x*z)).getName() == "x/z xz", "wrong name")
       self.assertTrue(((x/z)*(x/z)).getName() == "x/z x/z", "wrong name")
       self.assertTrue(((x/z)*(y**3)).getName() == "x/z y^3", "wrong name")

       self.assertTrue(((y**3)*x).getName() == "y^3 x", "wrong name")
       self.assertTrue(((y**3)*(x*z)).getName() == "y^3 xz", "wrong name")
       self.assertTrue(((y**3)*(x/z)).getName() == "y^3 x/z", "wrong name")
       self.assertTrue(((y**3)*(y**3)).getName() == "y^3 y^3", "wrong name")

       self.assertTrue((x/x).getName() == "x/x", "wrong name")
       self.assertTrue((x/(x*z)).getName() == "x/(xz)", "wrong name")
       self.assertTrue((x/(x/z)).getName() == "x/(x/z)", "wrong name")
       self.assertTrue((x/(y**3)).getName() == "x/y^3", "wrong name")

       self.assertTrue(((x*z)/x).getName() == "xz/x", "wrong name")
       self.assertTrue(((x*z)/(x*z)).getName() == "xz/(xz)", "wrong name")
       self.assertTrue(((x*z)/(x/z)).getName() == "xz/(x/z)", "wrong name")
       self.assertTrue(((x*z)/(y**3)).getName() == "xz/y^3", "wrong name")

       self.assertTrue(((x/z)/x).getName() == "x/z/x", "wrong name")
       self.assertTrue(((x/z)/(x*z)).getName() == "x/z/(xz)", "wrong name")
       self.assertTrue(((x/z)/(x/z)).getName() == "x/z/(x/z)", "wrong name")
       self.assertTrue(((x/z)/(y**3)).getName() == "x/z/y^3", "wrong name")

       self.assertTrue(((y**3)/x).getName() == "y^3/x", "wrong name")
       self.assertTrue(((y**3)/(x*z)).getName() == "y^3/(xz)", "wrong name")
       self.assertTrue(((y**3)/(x/z)).getName() == "y^3/(x/z)", "wrong name")
       self.assertTrue(((y**3)/(y**3)).getName() == "y^3/y^3", "wrong name")

       self.assertTrue((x**2).getName() == "x^2", "wrong name")
       self.assertTrue(((x*z)**2).getName() == "(xz)^2", "wrong name")
       self.assertTrue(((x/z)**2).getName() == "(x/z)^2", "wrong name")
       self.assertTrue(((y**3)**2).getName() == "(y^3)^2", "wrong name")

       self.assertTrue(((x*x)).getLongName() == "X*X", "wrong long name")
       self.assertTrue((x*(x*z)).getLongName() == "X*X*Z", "wrong long name")
       self.assertTrue((x*(x/z)).getLongName() == "X*X/Z", "wrong long name")
       self.assertTrue((x*(y**3)).getLongName() == "X*Y^3", "wrong long name")

       self.assertTrue(((x*z)*x).getLongName() == "X*Z*X", "wrong long name")
       self.assertTrue(((x*z)*(x*z)).getLongName() == "X*Z*X*Z", "wrong long name")
       self.assertTrue(((x*z)*(x/z)).getLongName() == "X*Z*X/Z", "wrong long name")
       self.assertTrue(((x*z)*(y**3)).getLongName() == "X*Z*Y^3", "wrong long name")

       self.assertTrue(((x/z)*x).getLongName() == "X/Z*X", "wrong long name")
       self.assertTrue(((x/z)*(x*z)).getLongName() == "X/Z*X*Z", "wrong long name")
       self.assertTrue(((x/z)*(x/z)).getLongName() == "X/Z*X/Z", "wrong long name")
       self.assertTrue(((x/z)*(y**3)).getLongName() == "X/Z*Y^3", "wrong long name")

       self.assertTrue(((y**3)*x).getLongName() == "Y^3*X", "wrong long name")
       self.assertTrue(((y**3)*(x*z)).getLongName() == "Y^3*X*Z", "wrong long name")
       self.assertTrue(((y**3)*(x/z)).getLongName() == "Y^3*X/Z", "wrong long name")
       self.assertTrue(((y**3)*(y**3)).getLongName() == "Y^3*Y^3", "wrong long name")

       self.assertTrue((x/x).getLongName() == "X/X", "wrong long name")
       self.assertTrue((x/(x*z)).getLongName() == "X/(X*Z)", "wrong long name")
       self.assertTrue((x/(x/z)).getLongName() == "X/(X/Z)", "wrong long name")
       self.assertTrue((x/(y**3)).getLongName() == "X/Y^3", "wrong long name")

       self.assertTrue(((x*z)/x).getLongName() == "X*Z/X", "wrong long name")
       self.assertTrue(((x*z)/(x*z)).getLongName() == "X*Z/(X*Z)", "wrong long name")
       self.assertTrue(((x*z)/(x/z)).getLongName() == "X*Z/(X/Z)", "wrong long name")
       self.assertTrue(((x*z)/(y**3)).getLongName() == "X*Z/Y^3", "wrong long name")

       self.assertTrue(((x/z)/x).getLongName() == "X/Z/X", "wrong long name")
       self.assertTrue(((x/z)/(x*z)).getLongName() == "X/Z/(X*Z)", "wrong long name")
       self.assertTrue(((x/z)/(x/z)).getLongName() == "X/Z/(X/Z)", "wrong long name")
       self.assertTrue(((x/z)/(y**3)).getLongName() == "X/Z/Y^3", "wrong long name")

       self.assertTrue(((y**3)/x).getLongName() == "Y^3/X", "wrong long name")
       self.assertTrue(((y**3)/(x*z)).getLongName() == "Y^3/(X*Z)", "wrong long name")
       self.assertTrue(((y**3)/(x/z)).getLongName() == "Y^3/(X/Z)", "wrong long name")
       self.assertTrue(((y**3)/(y**3)).getLongName() == "Y^3/Y^3", "wrong long name")

       self.assertTrue((x**2).getLongName() == "X^2", "wrong long name")
       self.assertTrue(((x*z)**2).getLongName() == "(X*Z)^2", "wrong long name")
       self.assertTrue(((x/z)**2).getLongName() == "(X/Z)^2", "wrong long name")
       self.assertTrue(((y**3)**2).getLongName() == "(Y^3)^2", "wrong long name")

    def testPrefix(self):


       self.assertTrue(Yotta.getName() == "Y","Yotta name wrong")
       self.assertTrue(Zetta.getName() == "Z","Zetta name wrong")
       self.assertTrue(Exa.getName() == "E","Exa name wrong")
       self.assertTrue(Peta.getName() == "P","Peta name wrong")
       self.assertTrue(Tera.getName() == "T","Tera name wrong")
       self.assertTrue(Giga.getName() == "G","Giga name wrong")
       self.assertTrue(Mega.getName() == "M","Mega name wrong")
       self.assertTrue(Kilo.getName() == "k","Kilo name wrong")
       self.assertTrue(Hecto.getName() == "h","Hecto name wrong")
       self.assertTrue(Deca.getName() == "da","Deca name wrong")
       self.assertTrue(Deci.getName() == "d","Deci name wrong")
       self.assertTrue(Centi.getName() == "c","Centi name wrong")
       self.assertTrue(Milli.getName() == "m","Milli name wrong")
       self.assertTrue(Micro.getName() == "mu","Micro name wrong")
       self.assertTrue(Nano.getName() == "n","Nano name wrong")
       self.assertTrue(Pico.getName() == "p","Pico name wrong")
       self.assertTrue(Femto.getName() == "f","Femto name wrong")
       self.assertTrue(Atto.getName() == "a","Atto name wrong")
       self.assertTrue(Zepto.getName() == "z","Zepto name wrong")
       self.assertTrue(Yocto.getName() == "y","Yocto name wrong")

       self.assertTrue(Yotta.getLongName() == "Yotta","wrong long name")
       self.assertTrue(Zetta.getLongName() == "Zetta","wrong long name")
       self.assertTrue(Exa.getLongName() == "Exa","wrong long name")
       self.assertTrue(Peta.getLongName() == "Peta","wrong long name")
       self.assertTrue(Tera.getLongName() == "Tera","wrong long name")
       self.assertTrue(Giga.getLongName() == "Giga","wrong long name")
       self.assertTrue(Mega.getLongName() == "Mega","wrong long name")
       self.assertTrue(Kilo.getLongName() == "Kilo","wrong long name")
       self.assertTrue(Hecto.getLongName() == "Hecto","wrong long name")
       self.assertTrue(Deca.getLongName() == "Deca","wrong long name")
       self.assertTrue(Deci.getLongName() == "Deci","wrong long name")
       self.assertTrue(Centi.getLongName() == "Centi","wrong long name")
       self.assertTrue(Milli.getLongName() == "Milli","wrong long name")
       self.assertTrue(Micro.getLongName() == "Micro","wrong long name")
       self.assertTrue(Nano.getLongName() == "Nano","wrong long name")
       self.assertTrue(Pico.getLongName() == "Pico","wrong long name")
       self.assertTrue(Femto.getLongName() == "Femto","wrong long name")
       self.assertTrue(Atto.getLongName() == "Atto","wrong long name")
       self.assertTrue(Zepto.getLongName() == "Zepto","wrong long name")
       self.assertTrue(Yocto.getLongName() == "Yocto","wrong long name")

       self.assertTrue(abs(Yotta(1.)/1.e24 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Zetta(1.)/1.e21 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Exa(1.)/1.e18 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Peta(1.)/1.e15 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Tera(1.)/1.e12 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Giga(1.)/1.e9 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Mega(1.)/1.e6 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Kilo(1.)/1.e3 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Hecto(1.)/1.e2 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Deca(1.)/1.e1 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Deci(1.)/1.e-1 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Centi(1.)/1.e-2 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Milli(1.)/1.e-3 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Micro(1.)/1.e-6 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Nano(1.)/1.e-9 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Pico(1.)/1.e-12 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Femto(1.)/1.e-15 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Atto(1.)/1.e-18 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Zepto(1.)/1.e-21 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Yocto(1.)/1.e-24 -1.)  <= self.TOL , "wrong b value")

       self.assertTrue(abs( Yotta(1.) * Yocto(1.) - 1.) <= self.TOL , " Yotta or Yocto wrong.")
       self.assertTrue(Yotta(0.) == 0. , " Yotta a value wrong.")
       self.assertTrue(Yocto(0.) == 0. , " Yocto a value wrong.")
       self.assertTrue(abs( Zetta(1.) * Zepto(1.) - 1.) <= self.TOL , "  Zetta or Zepto  wrong.")
       self.assertTrue(Zetta(0.) == 0. , "  Zetta a value  wrong.")
       self.assertTrue(Zepto(0.) == 0. , "  Zepto a value  wrong.")
       self.assertTrue(abs( Exa(1.) * Atto(1.) - 1.) <= self.TOL , " Exa or Atto  wrong.")
       self.assertTrue(Exa(0.) == 0. , " Exa a value  wrong.")
       self.assertTrue(Atto(0.) == 0. , "  Atto a value  wrong.")
       self.assertTrue(abs( Peta(1.) *Femto(1.) - 1.) <= self.TOL , " Peta or Femto  wrong.")
       self.assertTrue(Peta(0.) == 0. , " Peta  a value wrong.")
       self.assertTrue(Femto(0.) == 0. , " Femto a value  wrong.")
       self.assertTrue(abs( Tera(1.) * Pico(1.) - 1.) <= self.TOL , " Tera or Pico  wrong.")
       self.assertTrue(Tera(0.) == 0. , " Tera a value wrong.")
       self.assertTrue(Pico(0.) == 0. , " Pico a value  wrong.")
       self.assertTrue(abs( Giga(1.) * Nano(1.) - 1.) <= self.TOL , "  Giga or Nano  wrong.")
       self.assertTrue(Giga(0.) == 0. , "  Giga a value  wrong.")
       self.assertTrue(Nano(0.) == 0. , "  Nano a value  wrong.")
       self.assertTrue(abs( Mega(1.) *Micro(1.) - 1.) <= self.TOL , " Mega or Micro  wrong.")
       self.assertTrue(Mega(0.) == 0. , "Mega a value wrong.")
       self.assertTrue(Micro(0.) == 0. , " Micro a value  wrong.")
       self.assertTrue(abs( Kilo(1.) * Milli(1.) - 1.) <= self.TOL , "Kilo or Milli   wrong.")
       self.assertTrue( Kilo(0.) == 0. , "Kilo a value wrong.")
       self.assertTrue(Milli(0.) == 0. , "Milli a value   wrong.")
       self.assertTrue(abs( Hecto(1.) * Centi(1.)- 1.) <= self.TOL , " Hecto or Centi  wrong.")
       self.assertTrue(Hecto(0.) == 0., " Hecto a value wrong.")
       self.assertTrue(Centi(0.)== 0. , " Centi a value  wrong.")
       self.assertTrue(abs( Deca(1.) * Deci(1.) - 1.) <= self.TOL , "Deca or Deci  wrong.")
       self.assertTrue(Deca(0.) == 0. , "Deca a value wrong.")
       self.assertTrue(Deci(0.) == 0. , "Deci a value  wrong.")
    def testLength(self):
       self.assertTrue(m.getName() == "m", "meter name wrong.")
       self.assertTrue(m.getLongName() == "meter", "meter long name wrong.")
       self.assertTrue(abs( m(1.) - 1.) <= self.TOL , " meter wrong.")
       self.assertTrue(abs( km(1.) - 1000.) <= self.TOL*1000. , " km wrong.")
       self.assertTrue(abs( cm(1.) - 1./100.) <= self.TOL/100. , " cm wrong.")
       self.assertTrue(abs( mm(1.) - 1./1000.) <= self.TOL/1000. , " mm wrong.")
       self.assertTrue(m(0.) == 0., " meter wrong.")
       self.assertTrue(km(0.) == 0., " km wrong.")
       self.assertTrue(cm(0.) == 0. , " cm wrong.")
       self.assertTrue(mm(0.) ==0.," mm wrong.")
    def testTime(self):
       self.assertTrue(sec.getName() == "sec" , "sec name wrong")
       self.assertTrue(sec.getLongName() == "second" , "sec long name wrong")
       self.assertTrue(minute.getName() == "min" , "minute name wrong")
       self.assertTrue(minute.getLongName() == "minute" , "minute long name wrong")
       self.assertTrue(h.getName() == "h" , "hour name wrong")
       self.assertTrue(h.getLongName() == "hour" , "hour long name wrong")
       self.assertTrue(day.getName() == "d" , "day name wrong")
       self.assertTrue(day.getLongName() == "day" , "day long name wrong")
       self.assertTrue(yr.getName() == "yr" , "year name wrong")
       self.assertTrue(yr.getLongName() == "year" , "year long name wrong")

       self.assertTrue(abs( sec(1.) - 1.) <= self.TOL , " second wrong.")
       self.assertTrue(abs( minute(1.) - 60.) <= self.TOL*60 , " minute wrong.")
       self.assertTrue(abs( h(1.) - 3600) <= self.TOL*3600., " h wrong.")
       self.assertTrue(abs( day(1.) - 86400.) <= self.TOL*86400 , " day wrong.")
       self.assertTrue(abs( yr(1.) - 31556952) <= self.TOL*31556952 , " yr wrong.")
       self.assertTrue(abs( Myr(1.) - 31556952.e6) <= self.TOL*31556952.e6 , " Myr wrong.")
       self.assertTrue(abs( Gyr(1.) - 31556952.e9) <= self.TOL*31556952.e9 , " Gyr wrong.")

       self.assertTrue(sec(0.) == 0. , " second wrong.")
       self.assertTrue(minute(0.) == 0. , " minute wrong.")
       self.assertTrue(h(0.) == 0., " h wrong.")
       self.assertTrue(day(0.) == 0. , " day wrong.")
       self.assertTrue(yr(0.) == 0. , " yr wrong.")
       self.assertTrue(Myr(0.) == 0. , " Myr wrong.")
       self.assertTrue(Gyr(0.) == 0., " Gyr wrong.")
    def testMass(self):
       self.assertTrue(kg.getName() == "kg" , " kg wrong.")
       self.assertTrue(lb.getName() == "lb" , " lb wrong.")
       self.assertTrue(kg.getLongName() == "kg" , " kg wrong.")
       self.assertTrue(lb.getLongName() == "pound" , " lb wrong.")

       self.assertTrue(abs( gram(1.) - 1./1000.) <= self.TOL/1000. , " gram wrong.")
       self.assertTrue(abs( kg(1.) - 1.) <= self.TOL , " kg wrong.")
       self.assertTrue(abs( ton(1.) - 1.e3) <= self.TOL*1.e3 , " ton wrong.")
       self.assertTrue(abs( lb(1.) - 0.45359237) <= self.TOL , " lb wrong.")

       self.assertTrue(gram(0.) == 0. , " gram wrong.")
       self.assertTrue(kg(0.) == 0. , " kg wrong.")
       self.assertTrue(ton(0.) == 0. , " ton wrong.")
       self.assertTrue(lb(0.) == 0. , " lb wrong.")

    def testCurrent(self):
       self.assertTrue(A.getName() == "A" , " A wrong.")
       self.assertTrue(A.getLongName() == "Ampere" , " A wrong.")
       self.assertTrue(A(0.) == 0. , " A wrong.")
       self.assertTrue(abs( A(1.)- 1.) <= self.TOL , " A wrong.")
    def testFrequency(self):
       self.assertTrue(Hz(0.) == 0. , " Hz wrong.")
       self.assertTrue(abs( Hz(1)- 1.) <= self.TOL , " Hz wrong.")

    def testForce(self):
       self.assertTrue(N.getName() == "N" , "N wrong.")
       self.assertTrue(N.getLongName() == "Newton" , " N wrong.")
       self.assertTrue(N(0.) == 0. , " N wrong.")
       self.assertTrue(abs( N(1)- 1.) <= self.TOL , " N wrong.")
    def testPressure(self):
       self.assertTrue(Pa.getName() == "Pa" , " Pa wrong.")
       self.assertTrue(Pa.getLongName() == "Pascal" , " Pa wrong.")
       self.assertTrue(Pa(0.) == 0. , " Pa wrong.")
       self.assertTrue(abs( Pa(1)- 1.) <= self.TOL , " Pa wrong.")

       self.assertTrue(atm.getName() == "atm" , " atm wrong.")
       self.assertTrue(atm.getLongName() == "atmosphere" , " atm wrong.")
       self.assertTrue(atm(0.) == 0. , " atm wrong.")
       self.assertTrue(abs( atm(1)- 101325.024) <= self.TOL*101325.024 , " atm wrong.")

    def testEnegry(self):
       self.assertTrue(J.getName() == "J" , " J wrong.")
       self.assertTrue(J.getLongName() == "Joule" , " J wrong.")
       self.assertTrue(J(0.) == 0. , " J wrong.")
       self.assertTrue(abs( J(1)- 1.) <= self.TOL , " J wrong.")

    def testPower(self):
       self.assertTrue(W.getName() == "W" , " W wrong.")
       self.assertTrue(W.getLongName() == "Watt" , " W wrong.")
       self.assertTrue(W(0.) == 0. , " W wrong.")
       self.assertTrue(abs( W(1)- 1.) <= self.TOL , " W wrong.")

    def testPotential(self):
       self.assertTrue(V.getName() == "V" , " V wrong.")
       self.assertTrue(V.getLongName() == "Volt" , " V wrong.")
       self.assertTrue(V(0.) == 0. , " V wrong.")
       self.assertTrue(abs( V(1)- 1.) <= self.TOL , " V wrong.")

    def testCharge(self):
       self.assertTrue(C.getName() == "C" , " C wrong.")
       self.assertTrue(C.getLongName() == "Coulomb" , " C wrong.")
       self.assertTrue(C(0.) == 0. , " C wrong.")
       self.assertTrue(abs( C(1)- 1.) <= self.TOL , " C wrong.")

    def testResistance(self):
       self.assertTrue(Ohm.getName() == "Ohm" , " Ohm wrong.")
       self.assertTrue(Ohm.getLongName() == "Ohm" , " Ohm wrong.")
       self.assertTrue(Ohm(0.) == 0. , " Ohm wrong.")
       self.assertTrue(abs( Ohm(1)- 1.) <= self.TOL , " Ohm wrong.")
    def testTemperature(self):
       self.assertTrue(K.getName() == "K" , " Kelvin wrong.")
       self.assertTrue(K.getLongName() == "Kelvin" , " Kelvin wrong.")
       self.assertTrue(K(0.) == 0. , " Kelvin wrong.")
       self.assertTrue(abs( K(1.)- 1.) <= self.TOL , " Kelvin wrong.")

       self.assertTrue(Celsius.getName() == "C" , " Celsius wrong.")
       self.assertTrue(Celsius.getLongName() == "Celsius" , " Celsius wrong.")
       self.assertTrue(abs(Celsius(20.) - 293.15) <= self.TOL*293.15 , " Celsius 20 conversion wrong.")
       self.assertTrue(abs(Celsius(-2.)- 271.15) <= self.TOL*271.15 , " Celsius conversion wrong.")
       self.assertTrue(abs(Celsius(100.)- 373.15) <= self.TOL*373.15 , " Celsius conversion wrong.")

       self.assertTrue(Fahrenheit.getName() == "F" , " Fahrenheit wrong.")
       self.assertTrue(Fahrenheit.getLongName() == "Fahrenheit" , " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(0.)- 255.3722222) <= 1.e-8*255.3722222, " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(20.)- 266.4833333) <=  1.e-8*266.4833333 , " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(100.)- 310.9277778 ) <= 1.e-8 * 310.9277778 , " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(-1.)- 254.8166667) <= 1.e-8 *254.8166667, "Fahrenheit wrong.")


    def testEarth(self):
       R=6367444.65
       self.assertTrue(abs( R_Earth_equator-6378137) <= self.TOL*R , " R_Earth_equator wrong.")
       self.assertTrue(abs( R_Earth_poles-6356752.3) <= self.TOL*R , " R_Earth_poles wrong.")
       self.assertTrue(abs( R_Earth-R) <= self.TOL*R , "R_Earth wrong.")
    def testVlight(self):
       self.assertTrue(abs(v_light - 299792458.) <= self.TOL*299792458. , " speed of light is wrong.")

if __name__ == '__main__':
    run_tests(__name__, exit_on_failure=True)


##############################################################################
#
# Copyright (c) 2003-2018 by The University of Queensland
# http://www.uq.edu.au
#
# Primary Business: Queensland, Australia
# Licensed under the Apache License, version 2.0
# http://www.apache.org/licenses/LICENSE-2.0
#
# Development until 2012 by Earth Systems Science Computational Center (ESSCC)
# Development 2012-2013 by School of Earth Sciences
# Development from 2014 by Centre for Geoscience Computing (GeoComp)
#
##############################################################################

from __future__ import print_function, division

__copyright__="""Copyright (c) 2003-2018 by The University of Queensland
http://www.uq.edu.au
Primary Business: Queensland, Australia"""
__license__="""Licensed under the Apache License, version 2.0
http://www.apache.org/licenses/LICENSE-2.0"""
__url__="https://launchpad.net/escript-finley"

import esys.escriptcore.utestselect as unittest
from esys.escriptcore.testing import *
import sys
from esys.escript.unitsSI import *
from esys.escript.util import EPSILON

class Test_UnitsSI(unittest.TestCase):
    TOL=EPSILON*100.
    def testUnit(self):
       s=Unit("s","something",1.,4.)
       self.assertTrue(s.getName() == "s", "wrong name")
       s.setName("t")
       self.assertTrue(s.getName() == "t", "wrong reset name")
       self.assertTrue(s.getLongName() == "something", "wrong long name")
       s.setLongName("test")
       self.assertTrue(s.getLongName() == "test", "wrong reset long name")
       self.assertTrue(abs( s(3.)-13.) <= self.TOL*13. , " s(3.) wrong.")
       self.assertTrue(abs( 3.*s-13.) <= self.TOL*13. , " 13.s wrong.")
       self.assertTrue(abs( 10.*s-41.) <= self.TOL*41. , " 10.s wrong.")
       self.assertTrue(abs( 41./s-10.) <= self.TOL*10. , " 10/s wrong.")
       self.assertTrue(abs( 13./s-3.) <= self.TOL*3. , " 13/s wrong.")

       x=Unit("x","X",10.,2.)
       y=Unit("y","Y",0.,4.)
       z=Unit("z","Z",0.,0.5)
       p_m=x*s
       self.assertTrue(p_m.getName() == "xt", "wrong p_m name")
       self.assertTrue(p_m.getLongName() == "X*test", "wrong p_m long name")
       self.assertTrue(abs(p_m(0.)-12.) <= self.TOL*12 , " p_m(0) wrong.")
       self.assertTrue(abs(p_m(1.)-20.) <= self.TOL*20 , " p_m(1) wrong.")

       p_d=y/z
       self.assertTrue(p_d.getName() == "y/z", "wrong p_d name")
       self.assertTrue(p_d.getLongName() == "Y/Z", "wrong p_d long name")
       self.assertTrue(abs(p_d(1.)-8.) <= self.TOL*8. , " p_d(0) wrong.")
       self.assertTrue(abs(p_d(2.)-16.) <= self.TOL*16. , " p_d(1) wrong.")

       p_p=y**3
       self.assertTrue(p_p.getName() == "y^3", "p_p mult name")
       self.assertTrue(p_p.getLongName() == "Y^3", "wrong p_p long name")
       self.assertTrue(abs(p_p(1.)-64.) <= self.TOL*64. , " p_p(0) wrong.")
       self.assertTrue(abs(p_p(.5)-32) <= self.TOL*32. , " p_p(1) wrong.")

        
       x=Unit("x","X",0.,2.)
       y=Unit("y","Y",0.,4.)
       z=Unit("z","Z",0.,0.5)

       self.assertTrue(((x*x)).getName() == "xx", "wrong name")
       self.assertTrue((x*(x*z)).getName() == "xxz", "wrong name")
       self.assertTrue((x*(x/z)).getName() == "xx/z", "wrong name")
       self.assertTrue((x*(y**3)).getName() == "x y^3", "wrong name")

       self.assertTrue(((x*z)*x).getName() == "xzx", "wrong name")
       self.assertTrue(((x*z)*(x*z)).getName() == "xzxz", "wrong name")
       self.assertTrue(((x*z)*(x/z)).getName() == "xzx/z", "wrong name")
       self.assertTrue(((x*z)*(y**3)).getName() == "xz y^3", "wrong name")

       self.assertTrue(((x/z)*x).getName() == "x/z x", "wrong name")
       self.assertTrue(((x/z)*(x*z)).getName() == "x/z xz", "wrong name")
       self.assertTrue(((x/z)*(x/z)).getName() == "x/z x/z", "wrong name")
       self.assertTrue(((x/z)*(y**3)).getName() == "x/z y^3", "wrong name")

       self.assertTrue(((y**3)*x).getName() == "y^3 x", "wrong name")
       self.assertTrue(((y**3)*(x*z)).getName() == "y^3 xz", "wrong name")
       self.assertTrue(((y**3)*(x/z)).getName() == "y^3 x/z", "wrong name")
       self.assertTrue(((y**3)*(y**3)).getName() == "y^3 y^3", "wrong name")

       self.assertTrue((x/x).getName() == "x/x", "wrong name")
       self.assertTrue((x/(x*z)).getName() == "x/(xz)", "wrong name")
       self.assertTrue((x/(x/z)).getName() == "x/(x/z)", "wrong name")
       self.assertTrue((x/(y**3)).getName() == "x/y^3", "wrong name")

       self.assertTrue(((x*z)/x).getName() == "xz/x", "wrong name")
       self.assertTrue(((x*z)/(x*z)).getName() == "xz/(xz)", "wrong name")
       self.assertTrue(((x*z)/(x/z)).getName() == "xz/(x/z)", "wrong name")
       self.assertTrue(((x*z)/(y**3)).getName() == "xz/y^3", "wrong name")

       self.assertTrue(((x/z)/x).getName() == "x/z/x", "wrong name")
       self.assertTrue(((x/z)/(x*z)).getName() == "x/z/(xz)", "wrong name")
       self.assertTrue(((x/z)/(x/z)).getName() == "x/z/(x/z)", "wrong name")
       self.assertTrue(((x/z)/(y**3)).getName() == "x/z/y^3", "wrong name")

       self.assertTrue(((y**3)/x).getName() == "y^3/x", "wrong name")
       self.assertTrue(((y**3)/(x*z)).getName() == "y^3/(xz)", "wrong name")
       self.assertTrue(((y**3)/(x/z)).getName() == "y^3/(x/z)", "wrong name")
       self.assertTrue(((y**3)/(y**3)).getName() == "y^3/y^3", "wrong name")

       self.assertTrue((x**2).getName() == "x^2", "wrong name")
       self.assertTrue(((x*z)**2).getName() == "(xz)^2", "wrong name")
       self.assertTrue(((x/z)**2).getName() == "(x/z)^2", "wrong name")
       self.assertTrue(((y**3)**2).getName() == "(y^3)^2", "wrong name")

       self.assertTrue(((x*x)).getLongName() == "X*X", "wrong long name")
       self.assertTrue((x*(x*z)).getLongName() == "X*X*Z", "wrong long name")
       self.assertTrue((x*(x/z)).getLongName() == "X*X/Z", "wrong long name")
       self.assertTrue((x*(y**3)).getLongName() == "X*Y^3", "wrong long name")

       self.assertTrue(((x*z)*x).getLongName() == "X*Z*X", "wrong long name")
       self.assertTrue(((x*z)*(x*z)).getLongName() == "X*Z*X*Z", "wrong long name")
       self.assertTrue(((x*z)*(x/z)).getLongName() == "X*Z*X/Z", "wrong long name")
       self.assertTrue(((x*z)*(y**3)).getLongName() == "X*Z*Y^3", "wrong long name")

       self.assertTrue(((x/z)*x).getLongName() == "X/Z*X", "wrong long name")
       self.assertTrue(((x/z)*(x*z)).getLongName() == "X/Z*X*Z", "wrong long name")
       self.assertTrue(((x/z)*(x/z)).getLongName() == "X/Z*X/Z", "wrong long name")
       self.assertTrue(((x/z)*(y**3)).getLongName() == "X/Z*Y^3", "wrong long name")

       self.assertTrue(((y**3)*x).getLongName() == "Y^3*X", "wrong long name")
       self.assertTrue(((y**3)*(x*z)).getLongName() == "Y^3*X*Z", "wrong long name")
       self.assertTrue(((y**3)*(x/z)).getLongName() == "Y^3*X/Z", "wrong long name")
       self.assertTrue(((y**3)*(y**3)).getLongName() == "Y^3*Y^3", "wrong long name")

       self.assertTrue((x/x).getLongName() == "X/X", "wrong long name")
       self.assertTrue((x/(x*z)).getLongName() == "X/(X*Z)", "wrong long name")
       self.assertTrue((x/(x/z)).getLongName() == "X/(X/Z)", "wrong long name")
       self.assertTrue((x/(y**3)).getLongName() == "X/Y^3", "wrong long name")

       self.assertTrue(((x*z)/x).getLongName() == "X*Z/X", "wrong long name")
       self.assertTrue(((x*z)/(x*z)).getLongName() == "X*Z/(X*Z)", "wrong long name")
       self.assertTrue(((x*z)/(x/z)).getLongName() == "X*Z/(X/Z)", "wrong long name")
       self.assertTrue(((x*z)/(y**3)).getLongName() == "X*Z/Y^3", "wrong long name")

       self.assertTrue(((x/z)/x).getLongName() == "X/Z/X", "wrong long name")
       self.assertTrue(((x/z)/(x*z)).getLongName() == "X/Z/(X*Z)", "wrong long name")
       self.assertTrue(((x/z)/(x/z)).getLongName() == "X/Z/(X/Z)", "wrong long name")
       self.assertTrue(((x/z)/(y**3)).getLongName() == "X/Z/Y^3", "wrong long name")

       self.assertTrue(((y**3)/x).getLongName() == "Y^3/X", "wrong long name")
       self.assertTrue(((y**3)/(x*z)).getLongName() == "Y^3/(X*Z)", "wrong long name")
       self.assertTrue(((y**3)/(x/z)).getLongName() == "Y^3/(X/Z)", "wrong long name")
       self.assertTrue(((y**3)/(y**3)).getLongName() == "Y^3/Y^3", "wrong long name")

       self.assertTrue((x**2).getLongName() == "X^2", "wrong long name")
       self.assertTrue(((x*z)**2).getLongName() == "(X*Z)^2", "wrong long name")
       self.assertTrue(((x/z)**2).getLongName() == "(X/Z)^2", "wrong long name")
       self.assertTrue(((y**3)**2).getLongName() == "(Y^3)^2", "wrong long name")

    def testPrefix(self):


       self.assertTrue(Yotta.getName() == "Y","Yotta name wrong")
       self.assertTrue(Zetta.getName() == "Z","Zetta name wrong")
       self.assertTrue(Exa.getName() == "E","Exa name wrong")
       self.assertTrue(Peta.getName() == "P","Peta name wrong")
       self.assertTrue(Tera.getName() == "T","Tera name wrong")
       self.assertTrue(Giga.getName() == "G","Giga name wrong")
       self.assertTrue(Mega.getName() == "M","Mega name wrong")
       self.assertTrue(Kilo.getName() == "k","Kilo name wrong")
       self.assertTrue(Hecto.getName() == "h","Hecto name wrong")
       self.assertTrue(Deca.getName() == "da","Deca name wrong")
       self.assertTrue(Deci.getName() == "d","Deci name wrong")
       self.assertTrue(Centi.getName() == "c","Centi name wrong")
       self.assertTrue(Milli.getName() == "m","Milli name wrong")
       self.assertTrue(Micro.getName() == "mu","Micro name wrong")
       self.assertTrue(Nano.getName() == "n","Nano name wrong")
       self.assertTrue(Pico.getName() == "p","Pico name wrong")
       self.assertTrue(Femto.getName() == "f","Femto name wrong")
       self.assertTrue(Atto.getName() == "a","Atto name wrong")
       self.assertTrue(Zepto.getName() == "z","Zepto name wrong")
       self.assertTrue(Yocto.getName() == "y","Yocto name wrong")

       self.assertTrue(Yotta.getLongName() == "Yotta","wrong long name")
       self.assertTrue(Zetta.getLongName() == "Zetta","wrong long name")
       self.assertTrue(Exa.getLongName() == "Exa","wrong long name")
       self.assertTrue(Peta.getLongName() == "Peta","wrong long name")
       self.assertTrue(Tera.getLongName() == "Tera","wrong long name")
       self.assertTrue(Giga.getLongName() == "Giga","wrong long name")
       self.assertTrue(Mega.getLongName() == "Mega","wrong long name")
       self.assertTrue(Kilo.getLongName() == "Kilo","wrong long name")
       self.assertTrue(Hecto.getLongName() == "Hecto","wrong long name")
       self.assertTrue(Deca.getLongName() == "Deca","wrong long name")
       self.assertTrue(Deci.getLongName() == "Deci","wrong long name")
       self.assertTrue(Centi.getLongName() == "Centi","wrong long name")
       self.assertTrue(Milli.getLongName() == "Milli","wrong long name")
       self.assertTrue(Micro.getLongName() == "Micro","wrong long name")
       self.assertTrue(Nano.getLongName() == "Nano","wrong long name")
       self.assertTrue(Pico.getLongName() == "Pico","wrong long name")
       self.assertTrue(Femto.getLongName() == "Femto","wrong long name")
       self.assertTrue(Atto.getLongName() == "Atto","wrong long name")
       self.assertTrue(Zepto.getLongName() == "Zepto","wrong long name")
       self.assertTrue(Yocto.getLongName() == "Yocto","wrong long name")

       self.assertTrue(abs(Yotta(1.)/1.e24 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Zetta(1.)/1.e21 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Exa(1.)/1.e18 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Peta(1.)/1.e15 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Tera(1.)/1.e12 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Giga(1.)/1.e9 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Mega(1.)/1.e6 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Kilo(1.)/1.e3 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Hecto(1.)/1.e2 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Deca(1.)/1.e1 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Deci(1.)/1.e-1 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Centi(1.)/1.e-2 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Milli(1.)/1.e-3 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Micro(1.)/1.e-6 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Nano(1.)/1.e-9 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Pico(1.)/1.e-12 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Femto(1.)/1.e-15 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Atto(1.)/1.e-18 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Zepto(1.)/1.e-21 -1.)  <= self.TOL , "wrong b value")
       self.assertTrue(abs(Yocto(1.)/1.e-24 -1.)  <= self.TOL , "wrong b value")

       self.assertTrue(abs( Yotta(1.) * Yocto(1.) - 1.) <= self.TOL , " Yotta or Yocto wrong.")
       self.assertTrue(Yotta(0.) == 0. , " Yotta a value wrong.")
       self.assertTrue(Yocto(0.) == 0. , " Yocto a value wrong.")
       self.assertTrue(abs( Zetta(1.) * Zepto(1.) - 1.) <= self.TOL , "  Zetta or Zepto  wrong.")
       self.assertTrue(Zetta(0.) == 0. , "  Zetta a value  wrong.")
       self.assertTrue(Zepto(0.) == 0. , "  Zepto a value  wrong.")
       self.assertTrue(abs( Exa(1.) * Atto(1.) - 1.) <= self.TOL , " Exa or Atto  wrong.")
       self.assertTrue(Exa(0.) == 0. , " Exa a value  wrong.")
       self.assertTrue(Atto(0.) == 0. , "  Atto a value  wrong.")
       self.assertTrue(abs( Peta(1.) *Femto(1.) - 1.) <= self.TOL , " Peta or Femto  wrong.")
       self.assertTrue(Peta(0.) == 0. , " Peta  a value wrong.")
       self.assertTrue(Femto(0.) == 0. , " Femto a value  wrong.")
       self.assertTrue(abs( Tera(1.) * Pico(1.) - 1.) <= self.TOL , " Tera or Pico  wrong.")
       self.assertTrue(Tera(0.) == 0. , " Tera a value wrong.")
       self.assertTrue(Pico(0.) == 0. , " Pico a value  wrong.")
       self.assertTrue(abs( Giga(1.) * Nano(1.) - 1.) <= self.TOL , "  Giga or Nano  wrong.")
       self.assertTrue(Giga(0.) == 0. , "  Giga a value  wrong.")
       self.assertTrue(Nano(0.) == 0. , "  Nano a value  wrong.")
       self.assertTrue(abs( Mega(1.) *Micro(1.) - 1.) <= self.TOL , " Mega or Micro  wrong.")
       self.assertTrue(Mega(0.) == 0. , "Mega a value wrong.")
       self.assertTrue(Micro(0.) == 0. , " Micro a value  wrong.")
       self.assertTrue(abs( Kilo(1.) * Milli(1.) - 1.) <= self.TOL , "Kilo or Milli   wrong.")
       self.assertTrue( Kilo(0.) == 0. , "Kilo a value wrong.")
       self.assertTrue(Milli(0.) == 0. , "Milli a value   wrong.")
       self.assertTrue(abs( Hecto(1.) * Centi(1.)- 1.) <= self.TOL , " Hecto or Centi  wrong.")
       self.assertTrue(Hecto(0.) == 0., " Hecto a value wrong.")
       self.assertTrue(Centi(0.)== 0. , " Centi a value  wrong.")
       self.assertTrue(abs( Deca(1.) * Deci(1.) - 1.) <= self.TOL , "Deca or Deci  wrong.")
       self.assertTrue(Deca(0.) == 0. , "Deca a value wrong.")
       self.assertTrue(Deci(0.) == 0. , "Deci a value  wrong.")
    def testLength(self):
       self.assertTrue(m.getName() == "m", "meter name wrong.")
       self.assertTrue(m.getLongName() == "meter", "meter long name wrong.")
       self.assertTrue(abs( m(1.) - 1.) <= self.TOL , " meter wrong.")
       self.assertTrue(abs( km(1.) - 1000.) <= self.TOL*1000. , " km wrong.")
       self.assertTrue(abs( cm(1.) - 1./100.) <= self.TOL/100. , " cm wrong.")
       self.assertTrue(abs( mm(1.) - 1./1000.) <= self.TOL/1000. , " mm wrong.")
       self.assertTrue(m(0.) == 0., " meter wrong.")
       self.assertTrue(km(0.) == 0., " km wrong.")
       self.assertTrue(cm(0.) == 0. , " cm wrong.")
       self.assertTrue(mm(0.) ==0.," mm wrong.")
    def testTime(self):
       self.assertTrue(sec.getName() == "sec" , "sec name wrong")
       self.assertTrue(sec.getLongName() == "second" , "sec long name wrong")
       self.assertTrue(minute.getName() == "min" , "minute name wrong")
       self.assertTrue(minute.getLongName() == "minute" , "minute long name wrong")
       self.assertTrue(h.getName() == "h" , "hour name wrong")
       self.assertTrue(h.getLongName() == "hour" , "hour long name wrong")
       self.assertTrue(day.getName() == "d" , "day name wrong")
       self.assertTrue(day.getLongName() == "day" , "day long name wrong")
       self.assertTrue(yr.getName() == "yr" , "year name wrong")
       self.assertTrue(yr.getLongName() == "year" , "year long name wrong")

       self.assertTrue(abs( sec(1.) - 1.) <= self.TOL , " second wrong.")
       self.assertTrue(abs( minute(1.) - 60.) <= self.TOL*60 , " minute wrong.")
       self.assertTrue(abs( h(1.) - 3600) <= self.TOL*3600., " h wrong.")
       self.assertTrue(abs( day(1.) - 86400.) <= self.TOL*86400 , " day wrong.")
       self.assertTrue(abs( yr(1.) - 31556952) <= self.TOL*31556952 , " yr wrong.")
       self.assertTrue(abs( Myr(1.) - 31556952.e6) <= self.TOL*31556952.e6 , " Myr wrong.")
       self.assertTrue(abs( Gyr(1.) - 31556952.e9) <= self.TOL*31556952.e9 , " Gyr wrong.")

       self.assertTrue(sec(0.) == 0. , " second wrong.")
       self.assertTrue(minute(0.) == 0. , " minute wrong.")
       self.assertTrue(h(0.) == 0., " h wrong.")
       self.assertTrue(day(0.) == 0. , " day wrong.")
       self.assertTrue(yr(0.) == 0. , " yr wrong.")
       self.assertTrue(Myr(0.) == 0. , " Myr wrong.")
       self.assertTrue(Gyr(0.) == 0., " Gyr wrong.")
    def testMass(self):
       self.assertTrue(kg.getName() == "kg" , " kg wrong.")
       self.assertTrue(lb.getName() == "lb" , " lb wrong.")
       self.assertTrue(kg.getLongName() == "kg" , " kg wrong.")
       self.assertTrue(lb.getLongName() == "pound" , " lb wrong.")

       self.assertTrue(abs( gram(1.) - 1./1000.) <= self.TOL/1000. , " gram wrong.")
       self.assertTrue(abs( kg(1.) - 1.) <= self.TOL , " kg wrong.")
       self.assertTrue(abs( ton(1.) - 1.e3) <= self.TOL*1.e3 , " ton wrong.")
       self.assertTrue(abs( lb(1.) - 0.45359237) <= self.TOL , " lb wrong.")

       self.assertTrue(gram(0.) == 0. , " gram wrong.")
       self.assertTrue(kg(0.) == 0. , " kg wrong.")
       self.assertTrue(ton(0.) == 0. , " ton wrong.")
       self.assertTrue(lb(0.) == 0. , " lb wrong.")

    def testCurrent(self):
       self.assertTrue(A.getName() == "A" , " A wrong.")
       self.assertTrue(A.getLongName() == "Ampere" , " A wrong.")
       self.assertTrue(A(0.) == 0. , " A wrong.")
       self.assertTrue(abs( A(1.)- 1.) <= self.TOL , " A wrong.")
    def testFrequency(self):
       self.assertTrue(Hz(0.) == 0. , " Hz wrong.")
       self.assertTrue(abs( Hz(1)- 1.) <= self.TOL , " Hz wrong.")

    def testForce(self):
       self.assertTrue(N.getName() == "N" , "N wrong.")
       self.assertTrue(N.getLongName() == "Newton" , " N wrong.")
       self.assertTrue(N(0.) == 0. , " N wrong.")
       self.assertTrue(abs( N(1)- 1.) <= self.TOL , " N wrong.")
    def testPressure(self):
       self.assertTrue(Pa.getName() == "Pa" , " Pa wrong.")
       self.assertTrue(Pa.getLongName() == "Pascal" , " Pa wrong.")
       self.assertTrue(Pa(0.) == 0. , " Pa wrong.")
       self.assertTrue(abs( Pa(1)- 1.) <= self.TOL , " Pa wrong.")

       self.assertTrue(atm.getName() == "atm" , " atm wrong.")
       self.assertTrue(atm.getLongName() == "atmosphere" , " atm wrong.")
       self.assertTrue(atm(0.) == 0. , " atm wrong.")
       self.assertTrue(abs( atm(1)- 101325.024) <= self.TOL*101325.024 , " atm wrong.")

    def testEnegry(self):
       self.assertTrue(J.getName() == "J" , " J wrong.")
       self.assertTrue(J.getLongName() == "Joule" , " J wrong.")
       self.assertTrue(J(0.) == 0. , " J wrong.")
       self.assertTrue(abs( J(1)- 1.) <= self.TOL , " J wrong.")

    def testPower(self):
       self.assertTrue(W.getName() == "W" , " W wrong.")
       self.assertTrue(W.getLongName() == "Watt" , " W wrong.")
       self.assertTrue(W(0.) == 0. , " W wrong.")
       self.assertTrue(abs( W(1)- 1.) <= self.TOL , " W wrong.")

    def testPotential(self):
       self.assertTrue(V.getName() == "V" , " V wrong.")
       self.assertTrue(V.getLongName() == "Volt" , " V wrong.")
       self.assertTrue(V(0.) == 0. , " V wrong.")
       self.assertTrue(abs( V(1)- 1.) <= self.TOL , " V wrong.")

    def testCharge(self):
       self.assertTrue(C.getName() == "C" , " C wrong.")
       self.assertTrue(C.getLongName() == "Coulomb" , " C wrong.")
       self.assertTrue(C(0.) == 0. , " C wrong.")
       self.assertTrue(abs( C(1)- 1.) <= self.TOL , " C wrong.")

    def testResistance(self):
       self.assertTrue(Ohm.getName() == "Ohm" , " Ohm wrong.")
       self.assertTrue(Ohm.getLongName() == "Ohm" , " Ohm wrong.")
       self.assertTrue(Ohm(0.) == 0. , " Ohm wrong.")
       self.assertTrue(abs( Ohm(1)- 1.) <= self.TOL , " Ohm wrong.")
    def testTemperature(self):
       self.assertTrue(K.getName() == "K" , " Kelvin wrong.")
       self.assertTrue(K.getLongName() == "Kelvin" , " Kelvin wrong.")
       self.assertTrue(K(0.) == 0. , " Kelvin wrong.")
       self.assertTrue(abs( K(1.)- 1.) <= self.TOL , " Kelvin wrong.")

       self.assertTrue(Celsius.getName() == "C" , " Celsius wrong.")
       self.assertTrue(Celsius.getLongName() == "Celsius" , " Celsius wrong.")
       self.assertTrue(abs(Celsius(20.) - 293.15) <= self.TOL*293.15 , " Celsius 20 conversion wrong.")
       self.assertTrue(abs(Celsius(-2.)- 271.15) <= self.TOL*271.15 , " Celsius conversion wrong.")
       self.assertTrue(abs(Celsius(100.)- 373.15) <= self.TOL*373.15 , " Celsius conversion wrong.")

       self.assertTrue(Fahrenheit.getName() == "F" , " Fahrenheit wrong.")
       self.assertTrue(Fahrenheit.getLongName() == "Fahrenheit" , " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(0.)- 255.3722222) <= 1.e-8*255.3722222, " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(20.)- 266.4833333) <=  1.e-8*266.4833333 , " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(100.)- 310.9277778 ) <= 1.e-8 * 310.9277778 , " Fahrenheit wrong.")
       self.assertTrue(abs( Fahrenheit(-1.)- 254.8166667) <= 1.e-8 *254.8166667, "Fahrenheit wrong.")


    def testEarth(self):
       R=6367444.65
       self.assertTrue(abs( R_Earth_equator-6378137) <= self.TOL*R , " R_Earth_equator wrong.")
       self.assertTrue(abs( R_Earth_poles-6356752.3) <= self.TOL*R , " R_Earth_poles wrong.")
       self.assertTrue(abs( R_Earth-R) <= self.TOL*R , "R_Earth wrong.")
    def testVlight(self):
       self.assertTrue(abs(v_light - 299792458.) <= self.TOL*299792458. , " speed of light is wrong.")

if __name__ == '__main__':
    run_tests(__name__, exit_on_failure=True)

############################################################################## # # Copyright (c) 2003-2018 by The University of Queensland # http://www.uq.edu.au # # Primary Business: Queensland, Australia # Licensed under the Apache License, version 2.0 # http://www.apache.org/licenses/LICENSE-2.0 # # Development until 2012 by Earth Systems Science Computational Center (ESSCC) # Development 2012-2013 by School of Earth Sciences # Development from 2014 by Centre for Geoscience Computing (GeoComp) # ############################################################################## Copyright (c) 2003-2018 by The University of Queensland http://www.uq.edu.au Primary Business: Queensland, Australia Licensed under the Apache License, version 2.0 http://www.apache.org/licenses/LICENSE-2.0