phys-units 0.9.1 → 0.9.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 8fd589f3ac5dded9010b0d83ba9f31abccc40a4f
-  data.tar.gz: 17427a7a2bcb4631adea755b4fed53e9f0adccdd
+  metadata.gz: 562362594d0457c4492bc6a2467e2ee271c81dab
+  data.tar.gz: 86d22f70ec1f1962631fde49b63fcf72e408650b
 SHA512:
-  metadata.gz: a5018dbf704ba64a03afcfa92718389636f9ac058ee633eebe81ef5e27353178bb7eb9e9c88b127508258d12584071ab32a37df3827e66f26e46fe5909266cff
-  data.tar.gz: 02d92bee85d27cfc47555a3b14bad2b0c2df683786700e73747d0b1a2c65737b9241f081bb03d2ff569d76af54ad1561ad711342ac3c80db7eb5b840e0dd86ac
+  metadata.gz: 1124141150ac160dcd3568ba3507e0cf86d1521ba80e49a7c6307b7130710943fa1e7452b3f0df54fe2bc885f67cd690f3f433ecdd7f28605897d3736cd37060
+  data.tar.gz: d92bfe06366359ff027c0eb599d2bd8c78c7795a8c4743789c017bdb23b6d9725e0c2a1d46ca30815f0f854f817c012e6ed546db680b4bc7fa75ed97c136ca4f

data/README.md

@@ -5,7 +5,17 @@ Former name is [Quanty](http://narray.rubyforge.org/quanty/quanty-en.html),
 the first Ruby units library released in 2001.
 
 ## Phys::Quantity
-is a primary class of Phys-Units library, to be manipulated by users.
+is the primary class of Phys-Units library, to be manipulated by users.
+It contains:
+
+* *Value* :
+  +value+ must be a class instance having arithmetic methods,
+  but it is not necessary to be a Numeric.
+  This is a duck typing way. 
+* *Unit* :
+  +unit+ is a instance of Phys::Unit class
+  obtained by parsing a +expr+ string.
+
 See Documentation at [Rubygems Site](https://rubygems.org/gems/phys-units)
 
 ## Installation
@@ -36,7 +46,7 @@ Phys-Units library is discriminated from many other units libraies for Ruby,
 by the following features:
 
 * Compatible with GNU Units except nonlinear units.
-* Provides 2415 units, 85 prefixes, including UTF-8 unit names.
+* Provides 2331 units, 85 prefixes, including UTF-8 unit names.
 * All units are defined in a unit data file from GNU Units
   and not defined as a Ruby codes, except temperature definitions.
 * No addition or modification to Ruby standard classes by default,
@@ -47,8 +57,18 @@ by the following features:
   in the decimal form such as `1.0e10'.
 * PI number has a dimension.
 
+## Platforms tested
+
+* ruby 2.0.0p0 (2013-02-24 revision 39474) [x86_64-linux]
+
 ## Copying License
-GPL3
+
+This program is free software.
+You can distribute/modify this program
+under the same terms as GPL3.
+See "COPYING" file.
+NO WARRANTY.
 
 ## Author
+
 Masahiro TANAKA

data/lib/phys/units/load_units.rb

@@ -927,8 +927,8 @@ wc                      water            # water column
 mach                    331.46 m/s       # speed of sound in dry air at STP
 standardtemp            273.15 K         # standard temperature
 stdtemp                 standardtemp
-normaltemp              tempF(70)        # for gas density, from NIST
-normtemp                normaltemp       # Handbook 44
+#normaltemp              tempF(70)        # for gas density, from NIST
+#normtemp                normaltemp       # Handbook 44
 
 # Weight of mercury and water at different temperatures using the standard
 # force of gravity.
@@ -3446,62 +3446,6 @@ megalerg                megaerg    # 'L' added to make it pronounceable [18].
 
 unitedstatesdollar      US$
 $                       dollar
-mark                    germanymark
-bolivar                 venezuelanbolivarfuerte
-bolivarfuerte           bolivar        # The currency was revalued by
-oldbolivar              1|1000 bolivar # a factor of 1000.
-peseta                  spainpeseta
-rand                    southafricarand
-escudo                  portugalescudo
-guilder                 netherlandsguilder
-hollandguilder          netherlandsguilder
-peso                    mexicopeso
-yen                     japanyen
-lira                    italylira
-rupee                   indiarupee
-drachma                 greecedrachma
-franc                   francefranc
-markka                  finlandmarkka
-britainpound            unitedkingdompound
-greatbritainpound       unitedkingdompound
-poundsterling           britainpound
-yuan                    chinayuan
-
-# Some European currencies have permanent fixed exchange rates with
-# the Euro.  These rates were taken from the EC's web site:
-# http://ec.europa.eu/economy_finance/euro/adoption/conversion/index_en.htm
-
-austriaschilling        1|13.7603 euro
-belgiumfranc            1|40.3399 euro
-estoniakroon            1|15.6466 euro # Equal to 1|8 germanymark
-finlandmarkka           1|5.94573 euro
-francefranc             1|6.55957 euro
-germanymark             1|1.95583 euro
-greecedrachma           1|340.75 euro
-irelandpunt             1|0.787564 euro
-italylira               1|1936.27 euro
-luxembourgfranc         1|40.3399 euro
-netherlandsguilder      1|2.20371 euro
-portugalescudo          1|200.482 euro
-spainpeseta             1|166.386 euro
-cypruspound             1|0.585274 euro
-maltalira               1|0.429300 euro
-sloveniatolar           1|239.640 euro
-slovakiakoruna          1|30.1260 euro
-
-UKP                     GBP        # Not an ISO code, but looks like one, and
-                                   # sometimes used on usenet.
-VEB                     1|1000 VEF # old venezuelan bolivar
-
-!include currency.units
-
-# Money on the gold standard, used in the late 19th century and early
-# 20th century.
-
-olddollargold           23.22 grains goldprice  # Used until 1934
-newdollargold           96|7 grains goldprice   # After Jan 31, 1934
-dollargold              newdollargold
-poundgold               113 grains goldprice
 
 # Nominal masses of US coins.  Note that dimes, quarters and half dollars
 # have weight proportional to value.  Before 1965 it was $40 / kg.
@@ -3513,45 +3457,6 @@ USquarterweight         US$ 0.25 / (20 US$ / lb)   # Since 1965
 UShalfdollarweight      US$ 0.50 / (20 US$ / lb)   # Since 1971
 USdollarmass            8.1 grams
 
-# British currency
-
-quid                    britainpound        # Slang names
-fiver                   5 quid
-tenner                  10 quid
-monkey                  500 quid
-brgrand                 1000 quid
-bob                     shilling
-
-shilling                1|20 britainpound   # Before decimalisation, there
-oldpence                1|12 shilling       # were 20 shillings to a pound,
-farthing                1|4 oldpence        # each of twelve old pence
-guinea                  21 shilling         # Still used in horse racing
-crown                   5 shilling
-florin                  2 shilling
-groat                   4 oldpence
-tanner                  6 oldpence
-brpenny                 0.01 britainpound
-pence                   brpenny
-tuppence                2 pence
-tuppenny                tuppence
-ha'penny                halfbrpenny
-hapenny                 ha'penny
-oldpenny                oldpence
-oldtuppence             2 oldpence
-oldtuppenny             oldtuppence
-threepence              3 oldpence    # threepence never refers to new money
-threepenny              threepence
-oldthreepence           threepence
-oldthreepenny           threepence
-oldhalfpenny            halfoldpenny
-oldha'penny             oldhalfpenny
-oldhapenny              oldha'penny
-brpony                  25 britainpound
-
-# Canadian currency
-
-loony                   1 canadadollar    # This coin depicts a loon
-toony                   2 canadadollar
 
 #
 # Units used for measuring volume of wood
@@ -5205,7 +5110,7 @@ hogshead                brhogshead
 # U+2150-               1|7  For some reason these characters are getting
 # U+2151-               1|9  flagged as invalid UTF8.
 # U+2152-               1|10
-ℯ                       exp(1)      # U+212F, base of natural log
+#ℯ                       exp(1)      # U+212F, base of natural log
 
 µ-                      micro       # micro sign U+00B5
 μ-                      micro       # small mu U+03BC
@@ -5223,13 +5128,6 @@ röntgen                 roentgen
 K                       K          # Kelvin symbol, U+212A
 ℓ                       liter      # unofficial abbreviation used in some places
 ¢                       cent
-£                       britainpound
-¥                       japanyen
-€                       euro
-₩                       southkoreawon
-₪                       israelnewshekel
-₤                       lira
-₨                       rupee
 
 Ω                       ohm       # Ohm symbol U+2126
 Ω                       ohm       # Greek capital omega U+03A9
@@ -5403,4 +5301,6 @@ module Phys
   zero_degc = Unit["stdtemp"].conversion_factor
   OffsetUnit.define( "tempC", Unit["K"], zero_degc )
   OffsetUnit.define( "tempF", Unit["K"]*Rational(5,9), zero_degc-32*Rational(5,9) )
+
+  #puts "#{Phys::Unit::LIST.size} units, #{Phys::Unit::PREFIX.size} prefixes"
 end

data/lib/phys/units/parse.rb

@@ -481,7 +481,6 @@ module_eval(<<'...end parse.y/module_eval...', 'parse.y', 63)
 
   def parse( str )
     return Unit.new(str) if str.empty?
-    #p str
     @q = []
 
     c = Unit.unit_chars
@@ -490,15 +489,24 @@ module_eval(<<'...end parse.y/module_eval...', 'parse.y', 63)
       case str
       when /\A[\s]+/o
       when /\A(\d+)(?:(?:\.(\d*))?(?:[eE]([+-]?\d+))?)/o
-	@q.push [:NUMBER, build_num($1,$2,$3)]
+        @q.push [:NUMBER, build_num($1,$2,$3)]
       when /\A(sin|cos|tan|log|ln|log2)\b/o
         @q.push [:UFUNC, $&]
       when /\A\//o
         @q.push [:DIV, $&]
       when /\Aper\b/o
         @q.push [:DIV, $&]
-      when /\A[^#{c}0-9,.-]+([^#{c}$-]*[^#{c}1-9,.])?/o
-        @q.push [:WORD, $&]
+        when /\A([^#{c}0-9,.-]([^#{c}$-]*[^#{c}0-9,.])?)(\d)?/o
+        a = $&
+        x = $1
+        n = $3
+        if n && n.to_i >= 2
+            @q.push [:WORD, x]
+            @q.push [:POW, '']
+            @q.push [:NUMBER, n.to_i]
+        else
+          @q.push [:WORD, a]
+        end
       when /\A[%'"]'?/o
         @q.push [:WORD, $&]
       when /\A\^|\A\*\*/o

data/lib/phys/units/parse.y

@@ -41,7 +41,7 @@ rule
 
  list: numexpr
      | pexpr
-     | WORD	               { result = Unit.word(val[0]) }
+     | WORD                    { result = Unit.word(val[0]) }
      | list list    = MULTIPLY { result = val[0] * val[1] }
      | list POW list           { result = val[0]** val[2] }
      | list POW '-' list = POW { result = val[0]**(-val[3]) }
@@ -81,7 +81,6 @@ class Unit
 
   def parse( str )
     return Unit.new(str) if str.empty?
-    #p str
     @q = []
 
     c = Unit.unit_chars
@@ -90,15 +89,24 @@ class Unit
       case str
       when /\A[\s]+/o
       when /\A(\d+)(?:(?:\.(\d*))?(?:[eE]([+-]?\d+))?)/o
-	@q.push [:NUMBER, build_num($1,$2,$3)]
+        @q.push [:NUMBER, build_num($1,$2,$3)]
       when /\A(sin|cos|tan|log|ln|log2)\b/o
         @q.push [:UFUNC, $&]
       when /\A\//o
         @q.push [:DIV, $&]
       when /\Aper\b/o
         @q.push [:DIV, $&]
-      when /\A[^#{c}0-9,.-]+([^#{c}$-]*[^#{c}1-9,.])?/o
-        @q.push [:WORD, $&]
+        when /\A([^#{c}0-9,.-]([^#{c}$-]*[^#{c}0-9,.])?)(\d)?/o
+        a = $&
+        x = $1
+        n = $3
+        if n && n.to_i >= 2
+            @q.push [:WORD, x]
+            @q.push [:POW, '']
+            @q.push [:NUMBER, n.to_i]
+        else
+          @q.push [:WORD, a]
+        end
       when /\A[%'"]'?/o
         @q.push [:WORD, $&]
       when /\A\^|\A\*\*/o

data/lib/phys/units/quantity.rb

@@ -9,14 +9,21 @@
 
 module Phys
 
+  # alias to Phys::Quantity.new
   def Quantity(*a)
     Quantity.new(*a)
   end
 
   # Phys::Quantity is a class to represent physical quantities
   # with unit of measure.
-  # It contains a value in Numeric or similar class,
-  # and unit in Phys::Unit class.
+  # It contains:
+  # * *Value* of the quantity.
+  #   +value+ must be a class instance having arithmetic methods,
+  #   but it is not necessary to be a Numeric.
+  #   This is a duck typing way. 
+  # * *Unit* of the quantity.
+  #   +unit+ is a instance of Phys::Unit class
+  #   obtained by parsing a +expr+ string.
   #== Usage
   #   require 'phys/units'
   #   Q=Phys::Quantity
@@ -33,7 +40,7 @@ module Phys
       # Alias to Phys::Quantity.new.
       # @param  [Object] value
       #         Value of quantity.
-      # @param  [String] expr  unit expression.
+      # @param  [String] expr  a string of unit expression.
       #         If +expr+ is not supplied, it becomes dimeinsionless.
       # @return [Phys::Quantity]
       # @raise  [Phys::UnitConversionError] if unit conversion is failed.
@@ -45,9 +52,7 @@ module Phys
     # Initialize a new quantity.
     # @param  [Object] value
     #         Value of quantity.
-    #         +value+ must be a class instance having arithmetic methods
-    #         in a duck typing way.
-    # @param  [String] expr  unit expression.
+    # @param  [String] expr  a string of unit expression.
     #         If +expr+ is not supplied, it becomes dimeinsionless.
     # @param  [Phys::Unit] unit  (optional)
     # @raise  [Phys::UnitConversionError] if unit conversion is failed.
@@ -86,7 +91,7 @@ module Phys
     end
     alias convert want
 
-    # Addition of two quantities.
+    # Addition.
     # Before the operation, it converts +other+ to the unit of +self+.
     # * If the +other+ param is Phys::Quantity,
     #   +other+ is converted to the unit of +self+.
@@ -101,7 +106,7 @@ module Phys
       self.class.new( val, @expr, @unit )
     end
 
-    # Subtraction of two quantities.
+    # Subtraction.
     # Before the operation, it converts +other+ to the unit of +self+.
     # * If the +other+ param is Phys::Quantity,
     #   +other+ is converted to the unit of +self+.
@@ -116,32 +121,32 @@ module Phys
       self.class.new( val, @expr, @unit )
     end
 
-    # @return [Phys::Quantity] abs quantity in the unit of +self+.
+    # @return [Phys::Quantity] Abs. in the unit of +self+.
     def abs
       self.class.new( @val.abs, @expr, @unit )
     end
 
-    # @return [Phys::Quantity] abs2 quantity in the squared unit of +self+.
+    # @return [Phys::Quantity] Abs2. in the squared unit of +self+.
     def abs2
       self**2
     end
 
-    # @return [Phys::Quantity] ceil quantity in the unit of +self+.
+    # @return [Phys::Quantity] Ceil. in the unit of +self+.
     def ceil
       self.class.new( @val.ceil, @expr, @unit )
     end
 
-    # @return [Phys::Quantity] round quantity in the unit of +self+.
+    # @return [Phys::Quantity] Round. in the unit of +self+.
     def round
       self.class.new( @val.round, @expr, @unit )
     end
 
-    # @return [Phys::Quantity] floor quantity in the unit of +self+.
+    # @return [Phys::Quantity] Floor. in the unit of +self+.
     def floor
       self.class.new( @val.floor, @expr, @unit )
     end
 
-    # @return [Phys::Quantity] truncate quantity in the unit of +self+.
+    # @return [Phys::Quantity] Truncate. in the unit of +self+.
     def truncate
       self.class.new( @val.truncate, @expr, @unit )
     end
@@ -153,58 +158,81 @@ module Phys
     end
 
     # Unary Minus.
-    # @return [Phys::Quantity] a quantity in the unit of +self+.
+    # @return [Phys::Quantity] in the unit of +self+.
     def -@
       self.class.new( -@val, @expr, @unit )
     end
 
     # Comparison of quantities.
     # Before the comparison, it converts +other+ to the unit of +self+.
+    # @param  [Phys::Quantity] other
     # @return [Integer]
+    # @raise  [Phys::UnitConversionError] if unit conversion is failed.
     def <=> (other)
       @val <=> @unit.convert(other)
     end
 
-    # Comparison of quantities.
+    # Comparison. Returns +true+ if +self+ has the same value as +other+. 
     # Before the comparison, it converts +other+ to the unit of +self+.
+    # @param  [Phys::Quantity] other
     # @return [Boolean]
+    # @raise  [Phys::UnitConversionError] if unit conversion is failed.
     def  == (other)
       @val  == @unit.convert(other)
     end
 
-    # Comparison of quantities.
+    # Comparison. Returns +true+ if +self+ is greather-than or equal-to +other+.
     # Before the comparison, it converts +other+ to the unit of +self+.
+    # @param  [Phys::Quantity] other
     # @return [Boolean]
+    # @raise  [Phys::UnitConversionError] if unit conversion is failed.
     def  >= (other)
       @val  >= @unit.convert(other)
     end
 
-    # Comparison of quantities.
+    # Comparison. Returns +true+ if +self+ is less-than or equal-to +other+.
     # Before the comparison, it converts +other+ to the unit of +self+.
+    # @param  [Phys::Quantity] other
     # @return [Boolean]
+    # @raise  [Phys::UnitConversionError] if unit conversion is failed.
     def  <= (other)
       @val  <= @unit.convert(other)
     end
 
-    # Comparison of quantities.
+    # Comparison. Returns +true+ if +self+ is less than +other+.
     # Before the comparison, it converts +other+ to the unit of +self+.
+    # @param  [Phys::Quantity] other
     # @return [Boolean]
+    # @raise  [Phys::UnitConversionError] if unit conversion is failed.
     def  <  (other)
       @val  <  @unit.convert(other)
     end
 
-    # Comparison of quantities.
+    # Comparison. Returns +true+ if +self+ is greater than +other+.
     # Before the comparison, it converts +other+ to the unit of +self+.
+    # @param  [Phys::Quantity] other
     # @return [Boolean]
+    # @raise  [Phys::UnitConversionError] if unit conversion is failed.
     def  >  (other)
       @val  >  @unit.convert(other)
     end
 
+    # Closeness. Returns +true+ if +self+ and +other+ is close to each other
+    # within +epsilon+ relative to the absolute values.
+    # Before the comparison, it converts +other+ to the unit of +self+.
+    # @param  [Phys::Quantity] other
+    # @param  [Numeric] epsilon
+    # @return [Boolean]
+    # @raise  [Phys::UnitConversionError] if unit conversion is failed.
+    def close_to(other,epsilon=Float::EPSILON)
+      other_value = @unit.convert(other)
+      (@val-other_value).abs/(@val.abs+other_value.abs) <= epsilon
+    end
+
     # Power of a quantity.
     # @param  [Numeric] n
     # @return [Phys::Quantity] a quantity in the +n+ -powered unit of +self+.
     # @raise  [Phys::UnitConversionError] if unit conversion is failed.
-    #
     def **(n)
       if @expr.nil?
         expr = nil
@@ -240,7 +268,7 @@ module Phys
       end
     end
 
-    # Multiplication of two quantities.
+    # Multiplication.
     # If the +other+ param is *not* Phys::Quantity,
     # +other+ is regarded as a dimensionless value.
     # @param  [Object] other
@@ -258,7 +286,7 @@ module Phys
       end
     end
 
-    # Division of two quantities.
+    # Division.
     # If the +other+ param is *not* Phys::Quantity,
     # +other+ is regarded as a dimensionless value.
     # @param  [Object] other
@@ -276,7 +304,7 @@ module Phys
       end
     end
 
-    # Division of two quantities.
+    # Division more correctly.
     # If the +other+ param is *not* Phys::Quantity,
     # +other+ is regarded as a dimensionless value.
     # @param  [Object] other
@@ -295,7 +323,7 @@ module Phys
     end
     alias fdiv quo
 
-    # Division of two quantities without Remainder.
+    # Division without Modulo.
     # * If the +other+ param is Phys::Quantity,
     #   +other+ is converted to the unit of +self+,
     #   and returns +div+ of values.
@@ -314,7 +342,7 @@ module Phys
       end
     end
 
-    # Remainder of two quantities.
+    # Remainder.
     # * If the +other+ param is Phys::Quantity,
     #   +other+ is converted to the unit of +self+,
     #   and returns +remainder+ of values.
@@ -333,7 +361,7 @@ module Phys
       end
     end
 
-    # Modulo of two quantities.
+    # Modulo.
     # * If the +other+ param is Phys::Quantity,
     #   +other+ is converted to the unit of +self+,
     #   and returns +modulo+ of values.
@@ -353,7 +381,7 @@ module Phys
     end
     alias modulo %
 
-    # Division and Modulo of two quantities.
+    # Division with Modulo.
     # * If the +other+ param is Phys::Quantity,
     #   +other+ is converted to the unit of +self+,
     #   and returns +divmod+ of values.
@@ -383,7 +411,6 @@ module Phys
     # Returns the quantity converted to a base unit.
     # @return [Phys::Quantity] a quantity in the base unit.
     # @raise  [Phys::UnitConversionError] if unit conversion is failed.
-    #
     def to_base_unit
       unit = @unit.base_unit
       val  = unit.convert(self)
@@ -396,7 +423,6 @@ module Phys
     # Conversion to Numeric.
     # @return [Numeric]
     # @raise  [Phys::UnitConversionError] if the unit is non-dminensionless.
-    #
     def to_numeric
       @unit.convert_to_numeric(@val)
     end
@@ -405,7 +431,6 @@ module Phys
     # Conversion to Float.
     # @return [Float]
     # @raise  [Phys::UnitConversionError] if the unit is non-dminensionless.
-    #
     def to_f
       to_numeric.to_f
     end
@@ -414,7 +439,6 @@ module Phys
     # Conversion to Integer.
     # @return [Integer]
     # @raise  [Phys::UnitConversionError] if the unit is non-dminensionless.
-    #
     def to_i
       to_numeric.to_i
     end
@@ -424,7 +448,6 @@ module Phys
     # Conversion to Rational.
     # @return [Rational]
     # @raise  [Phys::UnitConversionError] if the unit is non-dminensionless.
-    #
     def to_r
       to_numeric.to_r
     end
@@ -432,7 +455,6 @@ module Phys
 
     # Conversion to String.
     # @return [String]
-    #
     def to_s
       if @expr
         expr = ",'" +@expr+"'"
@@ -444,7 +466,6 @@ module Phys
 
     # Inspect String.
     # @return [String]
-    #
     def inspect
       if @expr
         expr = "," +@expr.inspect

data/lib/phys/units/unit.rb

@@ -212,7 +212,7 @@ module Phys
     end
 
     def base_unit
-      Unit.new(1,dim)
+      Unit.new(1,dimensionless_deleted)
     end
 
     # Unit operation
@@ -294,40 +294,40 @@ module Phys
     end
 
     def *(x)
-      y = Unit.cast(x)
+      x = Unit.cast(x)
       if scalar?
-        return y
-      elsif y.scalar?
+        return x
+      elsif x.scalar?
         return self
       end
-      check_operable2(y)
-      dims = dimension_binop(y){|a,b| a+b}
-      factor = self.factor * y.factor
+      check_operable2(x)
+      dims = dimension_binop(x){|a,b| a+b}
+      factor = self.factor * x.factor
       Unit.new(factor,dims)
     end
 
     def /(x)
-      y = Unit.cast(x)
+      x = Unit.cast(x)
       if scalar?
-        return y.inv
-      elsif y.scalar?
+        return x.inv
+      elsif x.scalar?
         return self
       end
-      check_operable2(y)
-      dims = dimension_binop(y){|a,b| a-b}
-      factor = self.factor / y.factor
+      check_operable2(x)
+      dims = dimension_binop(x){|a,b| a-b}
+      factor = self.factor / x.factor
       Unit.new(factor,dims)
     end
 
     def rdiv(x)
-      y = Unit.cast(x)
+      x = Unit.cast(x)
       if scalar?
-        return y.inv
-      elsif y.scalar?
+        return x.inv
+      elsif x.scalar?
         return self
       end
-      check_operable2(y)
-      dims = dimension_binop(y){|a,b| a-b}
+      check_operable2(x)
+      dims = dimension_binop(x){|a,b| a-b}
       factor = Rational(self.factor,x.factor)
       Unit.new(factor,dims)
     end
@@ -342,6 +342,10 @@ module Phys
       Unit.new(Rational(1,self.factor), dims)
     end
 
+    def self.inv(x)
+      Unit.cast(x).inv
+    end
+
     def **(x)
       check_operable
       m = Utils.as_numeric(x)