numeric_with_unit 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 48d05ab9f9aa01a815b412324ef7a04b198041db
+  data.tar.gz: 5b55204a0a0b340dfe941ea141e29f6977197a34
+SHA512:
+  metadata.gz: d52e9021767bdbfee29dde88feaa8a64790475db07c0cee5df91c79f546e88604d08097f95a4730088df620336760d74d68f7f1f8ac9bf583f87ab7505b71a18
+  data.tar.gz: e3d277dee997ad68c2b38c8a4c06d6404b60bc1fb8d0328661055b5f75904afc834a8a6513d07c0d0e56688ae0577f6c35352c3fc1020080a92f8d8e0fea20ff

data/lib/numeric_with_unit.rb

@@ -0,0 +1,193 @@
+# coding: utf-8
+
+require 'numeric_with_unit/unit'
+
+class NumericWithUnit
+  include Comparable
+  
+  attr_reader :value, :unit
+  
+  def initialize(value, unit)
+    @value = value
+    @unit = unit.is_a?(Unit) ? unit : Unit[unit]
+  end
+  
+  def inspect
+    "#{@value.inspect} [#{@unit.symbol}] #{unit.dimension.inspect}"
+  end
+  
+  def to_s
+    "#{@value.to_s} #{@unit.symbol}"
+  end
+  
+  # otherがNumericWithUnitで次元が同じだったらsi単位に変換して比較、そうでなければ比較できない(nil)
+  def <=>(other)
+    if other.is_a?(self.class) and @unit.dimension_equal? other.unit
+      @unit.to_si(@value) <=> other.unit.to_si(other.value)
+    end
+  end
+  
+  def succ
+    self.class.new(@value.succ, @unit)
+  end
+  
+  def to_i
+    @value.to_i
+  end
+  
+  def to_f
+    @value.to_f
+  end
+  
+  def to_nwu(unit)
+    new_unit = unit.is_a?(Unit) ? unit : Unit[unit]
+    
+    unless @unit.dimension_equal? new_unit
+      raise DimensionError, "Dimensions are different between #{@unit.symbol}#{@unit.dimension} #{new_unit.symbol}#{new_unit.dimension}"
+    end
+    
+    new_value = new_unit.from_si(@unit.to_si(@value))
+    self.class.new(new_value, new_unit)
+  end
+  alias :[] :to_nwu
+  
+  def +@
+    self
+  end
+  
+  def -@
+    self.class.new(-@value, @unit)
+  end
+  
+  def +(other)
+    nwu = if other.is_a? self.class
+      other
+    else
+      self.class.new(other, Unit.new)
+    end
+    add_with_other_unit(nwu)
+  end
+  
+  def -(other)
+    self + (-other)
+  end
+  
+  def *(other)
+    case other
+    when self.class
+      multiply_with_other_unit(other)
+    else
+      self.class.new(@value*other, @unit)
+    end
+  end
+  
+  def /(other)
+    case other
+    when self.class
+      devide_with_other_unit(other)
+    else
+      self.class.new(@value/other, @unit)
+    end
+  end
+  
+  def coerce(other)
+    if other.is_a?(self.class)
+      [other, self]
+    else
+      [self.class.new(other, Unit.new), self]
+    end
+  end
+  
+  def **(num)
+    # Dimension Check
+    @unit.derivation.each do |k,v|
+      res = v * num
+      raise DimensionError, "Dimension of #{k.symbol}(#{v}*#{num}) must be Integer" unless res.to_i == res # 判定方法見なおせ
+    end
+    
+    self.class.new(@value**num, @unit**num)
+  end
+  
+  def root(num)
+    self**(Rational(1,num))
+  end
+  def sqrt; root(2) end # 平方根
+  def cbrt; root(3) end # 立方根
+  
+  private
+  
+  def add_with_other_unit(other)
+    if @unit.dimension_equal? other.unit
+      v1 = @unit.to_si(@value)
+      v2 = other.unit.to_si(other.value)
+      vr = @unit.from_si(v1+v2)
+      self.class.new(vr, @unit)
+    else
+      raise DimensionError, "Dimensions are different between #{@unit.dimension} #{other.unit.dimension}"
+    end
+  end
+  
+  def multiply_with_other_unit(other)
+    onwu = adjust_other_unit(other)
+    self.class.new(@value * onwu.value, @unit * onwu.unit)
+  end
+  
+  def devide_with_other_unit(other)
+    onwu = adjust_other_unit(other)
+    self.class.new(@value / onwu.value, @unit / onwu.unit)
+  end
+
+  # なるべくselfと同じ単位を使用するようにotherを変換します。
+  def adjust_other_unit(other)
+    if @unit.derivation.any?{|k,v| k == other.unit} # [L/min]*[min]などのケース
+      other
+    elsif h = @unit.derivation.find{|k,v| k.dimension_equal? other.unit} # [L/min]*[s]などのケース
+      other[ h.first ]
+    elsif @unit.dimension_equal? other.unit # [mm]*[cm]などのケース
+      other[@unit]
+    else
+      other
+    end
+  end
+  
+end
+
+
+
+class NumericWithUnit
+  class DimensionError < StandardError; end
+end
+
+
+
+class Fixnum
+  def to_nwu(unit)
+    NumericWithUnit.new(self, unit)
+  end
+end
+
+class Bignum
+  def to_nwu(unit)
+    NumericWithUnit.new(self, unit)
+  end
+end
+
+class Numeric
+  def to_nwu(unit)
+    NumericWithUnit.new(self, unit)
+  end
+end
+
+class String
+  def to_nwu(mthd=:to_r)
+    m = self.match /(?<value>.+) (?<unit>.+)/ # 適当
+    NumericWithUnit[m[:value].__send__(mthd), m[:unit]]
+  end
+end
+
+
+
+# unit definition
+require 'numeric_with_unit/base_unit'
+require 'numeric_with_unit/common_unit'
+

data/lib/numeric_with_unit/base_unit.rb

@@ -0,0 +1,179 @@
+# coding: utf-8
+
+#
+# SI base units & SI derived units & Units in use with SI
+# 
+
+require 'numeric_with_unit/unit'
+
+class NumericWithUnit
+  # Dimensionless
+  Unit << Unit.new do |conf|
+    conf.symbol = ''
+  end
+
+  #
+  # SI base units
+  #
+
+  # Length
+  Unit << Unit.new do |conf|
+    conf.symbol = 'm'
+    conf.dimension[:L] = 1
+    conf.si = true
+  end
+
+  # Mass
+  Unit << Unit.new do |conf|
+    conf.symbol = 'kg'
+    conf.dimension[:M] = 1
+    conf.si = true
+  end
+  Unit['g'] = "1/1000".to_r, 'kg' # for compatible
+
+  # Time
+  Unit << Unit.new do |conf|
+    conf.symbol = 's'
+    conf.dimension[:T] = 1
+    conf.si = true
+  end
+
+  # Electric Current
+  Unit << Unit.new do |conf|
+    conf.symbol = 'A'
+    conf.dimension[:I] = 1
+    conf.si = true
+  end
+
+  # Thermodynamic Temperature
+  Unit << Unit.new do |conf|
+    conf.symbol = 'K'
+    conf.dimension[:Θ] = 1
+    conf.si = true
+  end
+
+  # Amout of Substance
+  Unit << Unit.new do |conf|
+    conf.symbol = 'mol'
+    conf.dimension[:N] = 1
+    conf.si = true
+  end
+
+  # Luminous Intensity
+  Unit << Unit.new do |conf|
+    conf.symbol = 'cd'
+    conf.dimension[:J] = 1
+    conf.si = true
+  end
+
+
+
+  # 
+  # SI derived units
+  # 
+
+  # Frequency
+  Unit['Hz'] = '/s'
+
+  # Angle
+  Unit['rad'] = 'm/m'
+  Unit['°'] = Math::PI/180, 'rad'
+  Unit['′'] = "1/60".to_r, '°'
+  Unit['″'] = "1/60".to_r, '′'
+
+  # Solid Angle
+  Unit['sr'] = 'm2/m2'
+
+  # Force
+  Unit['N'] = 'kg.m/s2'
+
+  # Pressure
+  Unit['Pa'] = 'N/m2'
+
+  # Energy
+  Unit['J'] = 'N.m'
+
+  # Power
+  Unit['W'] = 'J/s'
+
+  # Electric Charge
+  Unit['C'] = 's.A'
+
+  # Voltage
+  Unit['V'] = 'W/A'
+
+  # Electriccal Capacitance
+  Unit['F'] = 'C/V'
+
+  # Electriccal Resistance
+  Unit['Ω'] = 'V/A'
+  Unit['ohm'] = 'Ω'
+
+  # Electriccal Conductance
+  Unit['S'] = 'A/V'
+
+  # Magnetic Flux
+  Unit['Wb'] = 'J/A'
+
+  # Magnetic Field Strength
+  Unit['T'] = 'Wb/m2'
+
+  # Inductance
+  Unit['H'] = 'V.s/A'
+
+  # Temperature
+  Unit << Unit.new do |conf|
+    k = Unit['K']
+    intercept = "273.15".to_r
+    
+    conf.symbol = '℃'
+    conf.dimension = k.dimension
+    conf.from_si{|x| k.from_si(x)-intercept}
+    conf.to_si{|x| k.to_si(x+intercept)}
+  end
+  Unit['degC'] = '℃'
+
+  # Luminouse flux
+  Unit['lx'] = 'cd.sr'
+
+  # Radioactivity
+  Unit['Bq'] = '/s'
+
+  # Absorbed Dose
+  Unit['Gy'] = 'J/kg'
+
+  # Equivalent Dose
+  Unit['Sv'] = 'J/kg'
+
+  # Catalytic Activity
+  Unit['kat'] = 'mol/s'
+
+ 
+
+  #
+  # Units in use with SI
+  #
+
+  # Time
+  Unit['min'] = 60, 's'
+  Unit['hr'] = 60, 'min'
+
+  # Area
+  Unit['a'] = 100, 'm2'
+
+  # Volume
+  Unit['L'] = 'dm3'
+
+  # Mass
+  Unit['t'] = 1000, 'kg'
+
+  # Energy
+  Unit['eV'] = "1.6021765314e-19".to_r, 'J'
+
+  # Mass
+  Unit['u'] = "1.6605388628e-27".to_r, 'kg'
+  Unit['Da'] = 'u'
+
+  # Length
+  Unit['ua'] = "1.495978706916e11".to_r, 'm'
+end

data/lib/numeric_with_unit/cgs_unit.rb

@@ -0,0 +1,19 @@
+# coding: utf-8
+
+require 'numeric_with_unit/base_unit'
+
+Unit['cal'] = "4.184".to_r, 'J'
+
+Unit['Gal'] = 'cm/s2'
+
+Unit['dyn'] = 'g.cm/s2'
+
+Unit['erg'] = 'g.cm2/s2'
+
+Unit['Ba'] = 'g/(cm.s2)'
+
+Unit['P'] = 'g/(cm.s)'
+Unit['poise'] = 'g/(cm.s)'
+
+Unit['St'] = 'cm2/s'
+

data/lib/numeric_with_unit/common_unit.rb

@@ -0,0 +1,112 @@
+# coding: utf-8
+
+# 
+# 独断と偏見による一般的な単位
+# 
+
+require 'numeric_with_unit/unit'
+require 'numeric_with_unit/base_unit'
+
+class NumericWithUnit
+  # Dimensionless
+  Unit['-'] = ''
+  Unit['1'] = ''
+
+  # Time
+  Unit['day'] = 24, 'hr'
+  Unit['week'] = 7, 'day'
+  Unit['month'] = 30, 'day' # 30日固定
+  Unit['year'] = 12, 'month' # 360日固定
+
+  # Mass
+  Unit['ton'] = 1000, 'kg'
+  Unit['oz'] = "28.349523125".to_r, 'g'
+  Unit['lb'] = 16, 'oz'
+
+  # Temperature
+  Unit['degC'] = '℃'
+  Unit['degR'] = "5/9".to_r, 'K'
+  Unit << Unit.new do |conf|
+    degr = Unit['degR']
+    intercept = "459.67".to_r
+    
+    conf.symbol = 'degF'
+    conf.dimension = degr.dimension
+    conf.from_si{|x| degr.from_si(x)-intercept}
+    conf.to_si{|x| degr.to_si(x+intercept)}
+  end
+
+  # Length
+  Unit['Å'] = 10**-10, 'm'
+  Unit['yd'] =  "0.9144".to_r, 'm'
+  Unit['ft'] = "1/3".to_r, 'yd'
+  Unit['in'] = "1/12".to_r, 'ft'
+  Unit['mi'] = 5280, 'ft'
+
+  # Volume
+  Unit['cc'] = 'cm3'
+  Unit['bbl'] = "0.158987294928".to_r, 'm3'
+
+  # Force
+  Unit['kgf'] = "9.80665".to_r, 'N'
+  Unit['lbf'] = "4.4482216152605".to_r, 'N'
+
+  # Power
+  Unit['PS'] = 75, 'kgf.m/s' # 仏馬力。小文字[ps]だとpico secondと区別がつかないため大文字で定義
+  Unit['HP'] = 550, 'lbf.ft/s' # 英馬力
+
+  # Pressure
+  Unit['bar'] = 1e5, 'Pa'
+  Unit['atm'] = 101325, 'Pa'
+  Unit['Torr'] = "101325/760".to_r, 'Pa'
+  Unit['mmHg'] = "101325/760".to_r, 'Pa'
+  Unit['mHg'] = "1/1000".to_r, 'mmHg' # for compatible
+  Unit['mH2O'] = "9806.65".to_r, 'Pa'
+  Unit['mAq'] = 'mH2O'
+  Unit['psi'] = "6894.76".to_r, 'Pa'
+
+  # Guage圧の扱いどうしようか？
+  Unit << Unit.new do |conf|
+    pa = Unit['Pa']
+    atm = 101325
+
+    conf.symbol = 'PaG'
+    conf.dimension = pa.dimension
+    conf.from_si{|x| pa.from_si(x)-atm}
+    conf.to_si{|x| pa.to_si(x+atm)}
+  end
+
+  # Speed
+  Unit['kph'] = 'km/hr'
+
+  # Flowrate
+  Unit['lpm'] = 'L/min'
+
+  # Viscosity
+  Unit['P'] = "1/10".to_r, 'Pa.s'
+
+  # Kinetic Viscosity
+  Unit['St'] = 'cm2/s'
+
+  # Energy
+  Unit['cal'] = "4.184".to_r, 'J' # 熱力学カロリー
+  Unit['MMkcal'] = 10**6 * 10**3, 'cal'
+  Unit['Btu'] = "1055.05585262".to_r, 'J' # 国際蒸気表(IT)Btu
+  Unit['MMBtu'] = 10**6, 'Btu'
+
+  # Ratio
+  Unit['%'] = 10**-2, ''
+  Unit['‰'] = 10**-3, ''
+  Unit['ppm'] = 10**-6, ''
+  Unit['ppb'] = 10**-9, ''
+
+
+  # Infomation
+  Unit << Unit.new do |conf|
+    conf.symbol = 'bit'
+  #  conf.dimension[:] = # 情報量の次元って何？
+  end
+  Unit['B'] = 8, 'bit'
+  Unit['bps'] = 'bit/s'
+  Unit['Bps'] = 'B/s'
+end

data/lib/numeric_with_unit/imperial_unit.rb

@@ -0,0 +1,43 @@
+# coding: utf-8
+
+require 'numeric_with_unit/base_unit'
+
+Unit['yd'] =  "0.9144".to_r , 'm'
+Unit['ft'] = "1/3".to_r, 'yd'
+Unit['th'] = "1/12000".to_r, 'ft'
+Unit['in'] = "1/12".to_r, 'ft'
+Unit['ch'] = 66, 'ft'
+Unit['fur'] = 660, 'ft'
+Unit['mi'] = 5280, 'ft'
+Unit['lea'] = 15840, 'ft'
+
+
+Unit['floz'] = "2.84130625e-05".to_r, 'm3'
+Unit['gi'] = 5, 'floz'
+Unit['pt'] = 20, 'floz'
+Unit['qt'] = 40, 'floz'
+Unit['gal'] = 160, 'floz'
+
+Unit['bbl'] = "0.158987294928".to_r, 'm3'
+
+
+Unit['oz'] = "0.45359237".to_r, 'kg'
+Unit['gr'] = "1/7000".to_r, 'oz'
+Unit['dr'] = "1/256".to_r, 'oz'
+Unit['lb'] = "1/16".to_r, 'oz'
+Unit['st'] = 14, 'oz'
+Unit['qtr'] = 28, 'oz'
+Unit['cwt'] = 112, 'oz'
+Unit['t'] = 2240, 'oz'
+
+
+Unit['degR'] = "9/5".to_r, 'K'
+
+Unit << Unit.new do |conf|
+  deg_r = Unit['degR']
+  conf.symbol = '℉'
+  conf.dimension = deg_r.dimension
+  conf.from_si = ->(x){(deg_r.from_si(x)) - 459.67}
+  conf.to_si = ->(x){deg_r.to_si(x + 459.67)}
+end
+Unit['degF'] = '℉'

data/lib/numeric_with_unit/natural_unit.rb

@@ -0,0 +1,14 @@
+# coding: utf-8
+
+require 'numeric_with_unit/base_unit'
+
+# Naural Units
+# Speed
+Unit['c0'] = "299792458.0".to_r, 'm/s'
+# Action
+Unit['ħ'] = "1.0545716818e−34".to_r, 'J.s'
+Unit['h'] = 'ħ'  # alias
+# Mass
+Unit['me'] = "9.109382616e−31".to_r, 'kg'
+# Time
+Unit['ħ/(me.(c0)2)'] = "1.288088667786e-21".to_r, 's'  # なんかきもい

data/lib/numeric_with_unit/unit.rb

@@ -0,0 +1,353 @@
+# coding: utf-8
+
+
+class NumericWithUnit
+  class Unit
+    class Config
+      attr_reader :symbol, :dimension, :derivation
+      attr_reader :si, :proportional
+      
+      def initialize(parent=nil)
+        @symbol = nil
+        @dimension = Hash.new(0)
+        @from_si = nil
+        @to_si = nil
+        @derivation = Hash.new(0)
+        @si = false
+        
+        @parent = parent
+      end
+      
+      def compile
+        @dimension.delete_if{|k,v| v.zero?}
+        @derivation.delete_if{|k,v| v.zero?}
+        @derivation.delete_if{|k,v| k.symbol.nil?}
+        
+        if @derivation.empty?
+          @from_si ||= ->(x){x}
+          @to_si ||= ->(x){x}
+          @derivation[@parent] += 1 unless @parent.nil?
+        else # configにderivationが与えられた時は、derivationをもとに@dimension,@symbol,@to_si,@from_siを設定
+          h = @derivation.sort_by{|u,v| u.symbol}.sort_by{|u,v| v} # ←どうしよう
+          
+          s1 = h.select{|u,v| v > 0}.map{|u,v| u.symbol + ((v.abs>1) ? v.abs.to_s : '')}.join('.')
+          s2 = h.select{|u,v| v < 0}.map{|u,v| u.symbol + ((v.abs>1) ? v.abs.to_s : '')}.join('.')
+          @symbol = s1 + (s2.empty? ? '' : "/(#{s2})")
+          
+          @derivation.each do |u,v|
+            u.dimension.each do |d,i|
+              @dimension[d] += i*v
+            end
+          end
+          
+          @from_si = @derivation.map{|u,v|
+            prc = if v > 0
+              ->(x){u.from_si(x)}
+            else
+              ->(x){x.quo(u.from_si(1)-u.from_si(0))} # ℃とKの変換のような場合に、変換式の切片を消すため。変換式が線形じゃないケースは想定していない
+            end
+            [prc, v.abs]
+          }.map{|prc,v|
+            ->(x){ v.times{x = prc[x]}; x }
+          }.reduce{|memo, prc|
+            ->(x){memo[prc[x]]}
+          }
+          
+          @to_si = @derivation.map{|u,v|
+            prc = if v > 0
+              ->(x){u.to_si(x)}
+            else
+              ->(x){x.quo(u.to_si(1)-u.to_si(0))} # ℃とKの変換のような場合に、変換式の切片を消すため。変換式が線形じゃないケースは想定していない
+            end
+            [prc, v.abs]
+          }.map{|prc,v|
+            ->(x){ v.times{x = prc[x]}; x }
+          }.reduce{|memo, prc|
+            ->(x){memo[prc[x]]}
+          }
+        end
+        
+        self
+      end
+      
+      def symbol=(arg)
+        @symbol = arg.to_s
+      end
+      
+      def dimension=(arg)
+        raise unless arg.is_a?(Hash)
+        @dimension = arg
+      end
+      
+      def from_si=(arg)
+        raise unless arg.is_a?(Proc)
+        @from_si = arg
+      end
+      
+      def from_si(&block)
+        if block_given?
+          @from_si = block
+        else
+          @from_si
+        end
+      end
+      
+      def to_si=(arg)
+        raise unless arg.is_a?(Proc)
+        @to_si = arg
+      end
+      
+      def to_si(&block)
+        if block_given?
+          @to_si = block
+        else
+          @to_si
+        end
+      end
+      
+      def si=(arg)
+        raise unless [TrueClass, FalseClass].any?{|klass|arg.is_a?(klass)}
+        @si = arg
+      end
+      
+      def derivation=(arg)
+        raise unless arg.is_a?(Hash)
+        @derivation = arg
+      end
+    end
+  end
+end
+
+
+
+class NumericWithUnit
+  class Unit
+    @@list = []
+    @@prefix = {
+      'Y' => 10**24, 
+      'Z' => 10**21, 
+      'E' => 10**18, 
+      'P' => 10**15, 
+      'T' => 10**12, 
+      'G' => 10**9, 
+      'M' => 10**6, 
+      'k' => 10**3, 
+      'h' => 10**2, 
+      'da' => 10**1, 
+      'd' => 10**-1, 
+      'c' => 10**-2, 
+      'm' => 10**-3, 
+      'μ' => 10**-6, 
+      'u' => 10**-6, 
+      'n' => 10**-9, 
+      'p' => 10**-12, 
+      'f' => 10**-15, 
+      'a' => 10**-18, 
+      'z' => 10**-21, 
+      'y' => 10**-24,
+      'Ki' => 2**10,
+      'Mi' => 2**20,
+      'Gi' => 2**30,
+      'Ti' => 2**40,
+      'Pi' => 2**50,
+      'Ei' => 2**60,
+      'Zi' => 2**70,
+      'Yi' => 2**80
+    }
+    
+    # class methods
+    
+    def self.[](arg)
+      self.parse(arg.to_s)
+    end
+    
+    def self.[]=(key, arg)
+      if arg.is_a?(Array) and arg.size == 2
+        a = [key, arg.first]
+        u = arg.last
+      else
+        a = [key]
+        u = arg
+      end
+      @@list << (u.is_a?(self) ? u : self[u]).cast(*a)
+    end
+    
+    def self.list
+      @@list.map(&:symbol)
+    end
+    
+    def self.<<(arg)
+      if arg.is_a?(self)
+        @@list << arg
+      else
+        @@list << Unit[arg]
+      end
+    end
+    
+    def self.delete(unit_symbol)
+      @@list.delete_if{|unit| unit.symbol == unit_symbol}
+    end
+    
+    def self.assign
+      @@list << self.new{|config| yield(config)}
+    end
+    
+    def self.parse(unit_str)
+      a = parse_1st(unit_str)
+      parse_2nd([a])
+    end
+    
+
+    # 文字列を配列にパース
+    # ex. 'J/(kg.K)' -> [#<Unit:J>, ['/', #<Unit:kg>, '.', #<Unit:K>]]
+    # とても手続き的な書き方で禿げる
+    def self.parse_1st(unit_str) 
+      i = @@list.rindex{|u| u.symbol == unit_str}
+      return @@list[i] if i
+      
+      return unit_str if unit_str =~ /^[\.\/]$/
+      
+      # 再帰で呼び出す用
+      rec = ->(arg){__send__(__method__, arg)}
+      
+      a = unit_str.scan(/(?<=\().*(?=\))|[\.\/]|[^\(\)\.\/]+/)
+      return a.map{|elem| rec[elem]} if a.size > 1
+      
+      m = unit_str.match(/-?\d+$/)
+      return m.to_s if m and m.pre_match.empty?
+      return [rec[m.pre_match], m.to_s] if m
+      
+      m = unit_str.match(/^(?<prefix>#{@@prefix.keys.join('|')})(?<unit>#{list.join('|')})$/)
+      return rec[m[:unit]].cast(unit_str, @@prefix[m[:prefix]]) if m
+      
+      raise NoUnitError, "\"#{unit_str}\" is not assigned"
+    end
+    private_class_method :parse_1st
+
+    # 配列を組立単位に変換
+    # derivationにそのまま突っ込んだほうがすっきりする気がする
+    def self.parse_2nd(unit_array)
+      # 再帰で呼び出す用
+      rec = ->(arg){__send__(__method__, arg)}
+      
+      buff_ary = []
+      buff_unit = ''
+      buff_sign = 1
+      
+      unit_array.each do |elem|
+        case elem
+        when self
+          buff_ary << elem ** buff_sign
+        when '.'
+          buff_sign = 1
+        when '/'
+          buff_sign = -1
+        when Array
+          buff_ary << rec[elem] ** buff_sign
+        when /^-?\d+$/
+          buff_ary[-1] **= elem.to_i
+        end
+      end
+      
+      buff_ary.reduce(:*)
+    end
+    private_class_method :parse_2nd
+    
+  end
+
+
+
+  class Unit
+    
+    # Instance Methods
+    
+  #  attr_accessor :symbol
+    attr_reader :symbol
+    attr_reader :dimension, :derivation
+    
+    def initialize
+      
+      # Unit::Configとinitializeの役割が分離できていないので見なおせ
+      config = Config.new(self)
+      yield(config) if block_given?
+      config.compile
+      
+      @symbol = config.symbol
+      @dimension = config.dimension
+      @from_si = config.from_si
+      @to_si = config.to_si
+      
+      @derivation = config.derivation
+    end
+    
+    def cast(new_symbol, factor = 1)
+      self.class.new do |conf|
+        conf.symbol = new_symbol
+        conf.dimension = @dimension
+        conf.from_si = ->(x){from_si(x.quo(factor))}
+        conf.to_si = ->(x){to_si(x * factor)}
+      end
+    end
+    
+    def to_s
+      @symbol
+    end
+    
+    def to_si(value)
+      @to_si[value]
+    end
+    
+    def from_si(value)
+      @from_si[value]
+    end
+    
+    
+    def dimensionless?
+      @dimension.all?{|k,v| v.zero?}
+    end
+    
+    def dimension_equal?(other_unit)
+      (@dimension.keys | other_unit.dimension.keys).all?{|k|
+        @dimension[k] == other_unit.dimension[k]
+      }
+    end
+    
+    def ==(other)
+      if other.is_a?(self.class)
+        symbol == other.symbol and dimension == other.dimension
+      else
+        super
+      end
+    end
+    
+    def *(other_unit)
+      self.class.new do |conf|
+        @derivation.each{|k, v| conf.derivation[k] += v}
+        other_unit.derivation.each{|k, v| conf.derivation[k] += v}
+      end
+    end
+    
+    def /(other_unit)
+      self.class.new do |conf|
+        @derivation.each{|k, v| conf.derivation[k] += v}
+        other_unit.derivation.each{|k, v| conf.derivation[k] -= v}
+      end
+    end
+    
+    def **(num)
+      if num.zero?
+        self.class.new
+      else
+        self.class.new do |conf|
+          # ここto_iでOKか？v*numが整数じゃなければraiseすべき？→すべき→NumericWithUnitでやるべき？
+          # Unitでは整数じゃない次数の単位は許容すべきか否か→していい気がする
+          @derivation.each{|k, v| conf.derivation[k] = (v*num).to_i}
+        end
+      end
+    end
+  end
+
+  class Unit
+    class NoUnitError < StandardError; end
+  end
+end

data/lib/numeric_with_unit/util.rb

@@ -0,0 +1,27 @@
+# coding: utf-8
+
+# 1234['m3/kg']のように書けるようにします。
+# Numeric#[]、Fixnum#[]、Bignum#[]をオーバーライドします。
+
+require 'numeric_with_unit'
+
+class NumericWithUnit
+  module NumUtil
+    def [](unit)
+      NumericWithUnit.new(self.rationalize, unit) # ratoinalizeする？
+    end
+  end
+end
+
+class Fixnum
+  prepend NumericWithUnit::NumUtil
+end
+
+class Bignum
+  prepend NumericWithUnit::NumUtil
+end
+
+class Numeric
+  prepend NumericWithUnit::NumUtil
+end
+

data/lib/numeric_with_unit/util2.rb

@@ -0,0 +1,66 @@
+# coding: utf-8
+
+# 123.m2.K_W のように書けるようにします。
+# method_missing をオーバーライドします。
+
+require 'numeric_with_unit'
+
+class NumericWithUnit
+  def method_missing(*args)
+    if args.size == 1
+      unit_str = args.first.to_s.gsub('_', '/')
+      unit_chain_util(Unit[unit_str])
+    else
+      raise Unit::NoUnitError
+    end
+  rescue Unit::NoUnitError
+    super
+  end
+  
+  attr_writer :unit_chain
+  
+  private
+  def unit_chain_util(unit)
+    ucs = @unit_chain || []
+    ucs.map!{|nwu, u| [nwu, u * unit]}
+    ucs << [self, unit]
+    
+    if i = ucs.index{|nwu, u| nwu.unit.dimension_equal? u}
+      nwu, nu = *ucs[i]
+      nwu[nu]
+    else
+      nnwu = self.class.new(@value, @unit*unit)
+      nnwu.unit_chain = ucs
+      nnwu
+    end
+  end
+end
+
+
+class NumericWithUnit
+  module NumUtil
+    def method_missing(*args)
+      if args.size == 1
+        unit_str = args.first.to_s.gsub('_', '/')
+        self.rationalize.to_nwu(unit_str) # util2は利便性優先なのでratoinalizeしてしまいます
+      else
+        raise Unit::NoUnitError
+      end
+    rescue Unit::NoUnitError
+      super
+    end
+  end
+end
+
+
+class Fixnum
+  prepend NumericWithUnit::NumUtil
+end
+
+class Bignum
+  prepend NumericWithUnit::NumUtil
+end
+
+class Numeric
+  prepend NumericWithUnit::NumUtil
+end

metadata

@@ -0,0 +1,52 @@
+--- !ruby/object:Gem::Specification
+name: numeric_with_unit
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- diaphragm
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2014-11-18 00:00:00.000000000 Z
+dependencies: []
+description: calclation numeric with unit!!!!
+email: 
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- lib/numeric_with_unit.rb
+- lib/numeric_with_unit/base_unit.rb
+- lib/numeric_with_unit/cgs_unit.rb
+- lib/numeric_with_unit/common_unit.rb
+- lib/numeric_with_unit/imperial_unit.rb
+- lib/numeric_with_unit/natural_unit.rb
+- lib/numeric_with_unit/unit.rb
+- lib/numeric_with_unit/util.rb
+- lib/numeric_with_unit/util2.rb
+homepage: https://github.com/diaphragm/ruby-numeric-with-unit
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.2.2
+signing_key: 
+specification_version: 4
+summary: calclation numeric with unit
+test_files: []