tuscan 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 2543badcdb3312b2f5b8137766cddc323887e10f
+  data.tar.gz: 5dc2685148190c63c30277e384fc5b272631bd91
+SHA512:
+  metadata.gz: 3d3dd9a95da487edd7a192cc273df640d9599ef068ec26ad7214e6c0abfb6369a96815cc196f736dbc07343db35d6fb10e68c78bf84930ed865fb75198701665
+  data.tar.gz: 159021a9b23eaa493f9462fde5b363819181364840a878a77f66fb2b8d92ca93d05ad3271545efdfea30a4519b80c642ef3226ce77bd45615593da49d09f0062

data/.gitignore

@@ -0,0 +1,17 @@
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/
+*.bundle
+*.so
+*.o
+*.a
+mkmf.log
+*.swp
+*.swo
+*~

data/.rspec

@@ -0,0 +1,2 @@
+--color
+--format documentation

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in tuscan.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2015 Luis Bacelar
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/README.md

@@ -0,0 +1,74 @@
+# **TUSCAn**
+### **T**emperat**U**re **S**ensor **CA**lculator
+
+Makes available functions for:
+
+* ITS-90 compatible temperature sensors (Standard Platinum Resistance Thermometers), allowing:
+  * Retrieval of reference resistance ratio _Wr_ from _t90r_ reference temperature
+  * Retrieval of reference temperature _t90r_ from _Wr_ reference resistance ratio
+  * Given sensor coefficient's: _sub-range_, _Rtpw_, _a_, _b_, _c_, _d_, _W660_, _c1_, _c2_, _c3_, _c4_, _c5_
+    * Retrieval of resistance ratio deviation  _Wdev_ from _t90_, for each of the 11 sub-ranges of ITS-90
+    * Retrieval of  _t90_ temperature from sensor's resistance _r_. 
+* IEC 60751 compatible temperature sensors (Platinum Resistance Thermometers), allowing:
+  * Given sensor coefficient's: _R0_, _A_, _B_, _C_
+    * Retrieval of sensor's resistance _r_ from _t90_ temperature
+    * Retrieval of temperature _t90_ from sensor's resistance _r_
+* IEC 60584 compatible temperature sensors (Thermocouples), allowing, for a given _type_:
+  * Retrieval of reference voltage _emfr_ from _t90r_ reference temperature
+  * Retrieval of reference temperature _t90r_ from reference voltage _emfr_
+
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'tuscan'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install tuscan
+
+## Usage examples
+
+**ITS-90**
+```ruby
+Tuscan.wr(:its90, 23.456)  # t90 in ºC
+Tuscan.t90r(:its90, 1.15)  # wr in Ohm / Ohm
+
+Tuscan.wdev(:its90, 23.456, subrange: 4, a: -1.2579994e-04, b: 1.0678395e-05)
+
+Tuscan.t90(:its90, r: 28.3433, rtpw: 25.319871, subrange: 7, a: -1.2134e-04, b: -9.9190e-06)
+```
+
+**IEC 60751**
+```ruby
+Tuscan.t90(:iec60751, r: 110)  # t90 in ºC, r in Ohm
+Tuscan.res(:iec60751, t: 10, r0: 99.876, a:3.9083e-03, b:-5.7750e-07, c:-4.1830e-12)
+```
+
+**IEC 60584**
+```ruby
+Tuscan.t90r(:iec60584, emf: 1.234, type: :k)  # emf im mV
+Tuscan.emfir(:iec60584, t: 10.123, type: :k)  # t90 in ºC
+```
+
+
+### Notes
+* `t90` aliased to `t`, `temperature`
+* `res` aliased to `r`, `resistance`
+* `emf` aliased to `v`, `voltage`
+
+
+## Contributing
+
+1. Fork it ( https://github.com/[my-github-username]/tuscan/fork )
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create a new Pull Request

data/Rakefile

@@ -0,0 +1,3 @@
+require 'bundler/gem_tasks'
+
+Dir.glob('tasks/**/*.rake').each(&method(:import))

data/lib/tuscan/iec60584.rb

@@ -0,0 +1,383 @@
+module Tuscan
+  module Iec60584
+    include Math
+    extend self
+
+
+    def emfr t90, type, range_check = true
+      a = emf_a_coeffs t90, type
+      result  = emf_c_coeffs(t90, type).each_with_index.map{ |c,i| c*t90**i }.reduce(:+)
+      result += a[0] * exp(a[1] * (t90 - a[2])**2)
+    end
+
+    def t90r emf, type, err_limit: 0.1, num_iter: 1000
+      t90_coeffs(emf, type).each_with_index.map{ |d,i| d*emf**i }.reduce(:+)
+    end
+
+
+  private
+    def emf_c_coeffs t90, type
+      case type
+      when :b then emf_c_coeffs_b t90
+      when :c then emf_c_coeffs_c t90
+      when :e then emf_c_coeffs_e t90
+      when :j then emf_c_coeffs_j t90
+      when :k then emf_c_coeffs_k t90
+      when :n then emf_c_coeffs_n t90
+      when :r then emf_c_coeffs_r t90
+      when :s then emf_c_coeffs_s t90
+      when :t then emf_c_coeffs_t t90
+      end
+    end
+
+    def emf_c_coeffs_b t90
+      #   0.000..630.615
+      # 630.615..1820.00
+      if t90 < 630.615
+        [  0.000000000000e+00, -0.246508183460e-03,  0.590404211710e-05,
+          -0.132579316360e-08,  0.156682919010e-11, -0.169445292400e-14,
+           0.629903470940e-18                                           ]
+      else
+        [ -0.389381686210e+01,  0.285717474700e-01, -0.848851047850e-04,
+           0.157852801640e-06, -0.168353448640e-09,  0.111097940130e-12,
+          -0.445154310330e-16,  0.989756408210e-20, -0.937913302890e-24 ]
+      end
+    end
+
+    def emf_c_coeffs_c t90
+      # 0.0..2320.0
+      [ -3.109077870000e-04,  1.338547130000e-02,  1.226236040000e-05,
+        -1.050537530000e-08,  3.613274640000e-12, -4.990804550000e-16,
+         6.434651840000e-22                                           ]
+    end
+
+    def emf_c_coeffs_e t90
+      # -270.0..0.0
+      #    0.0..1000.0
+      if t90 < 0.0
+        [  0.000000000000e+00,  0.586655087080e-01,  0.454109771240e-04,
+          -0.779980486860e-06, -0.258001608430e-07, -0.594525830570e-09,
+          -0.932140586670e-11, -0.102876055340e-12, -0.803701236210e-15,
+          -0.439794973910e-17, -0.164147763550e-19, -0.396736195160e-22,
+          -0.558273287210e-25, -0.346578420130e-28                      ]
+      else
+        [  0.000000000000e+00,  0.586655087100e-01,  0.450322755820e-04,
+           0.289084072120e-07, -0.330568966520e-09,  0.650244032700e-12,
+          -0.191974955040e-15, -0.125366004970e-17,  0.214892175690e-20,
+          -0.143880417820e-23,  0.359608994810e-27                      ]
+      end
+    end
+
+    def emf_c_coeffs_j t90
+      # -210.0..760.0
+      #  760.0..1200.0
+      if t90 < 760.0
+        [  0.000000000000e+00,  0.503811878150e-01,  0.304758369300e-04,
+          -0.856810657200e-07,  0.132281952950e-09, -0.170529583370e-12,
+           0.209480906970e-15, -0.125383953360e-18,  0.156317256970e-22 ]
+      else
+        [  0.296456256810e+03, -0.149761277860e+01,  0.317871039240e-02,
+          -0.318476867010e-05,  0.157208190040e-08, -0.306913690560e-12 ]
+      end
+    end
+
+    def emf_c_coeffs_k t90
+      # -270.0..0.0
+      #    0.0..1372.0
+      if t90 < 0.0
+        [  0.000000000000e+00,  0.394501280250e-01,  0.236223735980e-04,
+          -0.328589067840e-06, -0.499048287770e-08, -0.675090591730e-10,
+          -0.574103274280e-12, -0.310888728940e-14, -0.104516093650e-16,
+          -0.198892668780e-19, -0.163226974860e-22                      ]
+      else
+        [ -0.176004136860e-01,  0.389212049750e-01,  0.185587700320e-04,
+          -0.994575928740e-07,  0.318409457190e-09, -0.560728448890e-12,
+           0.560750590590e-15, -0.320207200030e-18,  0.971511471520e-22,
+          -0.121047212750e-25                                           ]
+      end
+    end
+
+    def emf_c_coeffs_n t90
+      # -270.0..0.0
+      #    0.0..1300.0
+      if t90 < 0.0
+        [  0.000000000000e+00,  0.261591059620e-01,  0.109574842280e-04,
+          -0.938411115540e-07, -0.464120397590e-10, -0.263033577160e-11,
+          -0.226534380030e-13, -0.760893007910e-16, -0.934196678350e-19 ]
+      else
+        [  0.000000000000e+00,  0.259293946010e-01,  0.157101418800e-04,
+           0.438256272370e-07, -0.252611697940e-09,  0.643118193390e-12,
+          -0.100634715190e-14,  0.997453389920e-18, -0.608632456070e-21,
+           0.208492293390e-24, -0.306821961510e-28                      ]
+      end
+    end
+
+    def emf_c_coeffs_r t90
+      #  -50.00..1064.18
+      # 1064.18..1664.50
+      # 1664.50..1768.10
+      if t90 < 1064.18
+        [  0.000000000000e+00,  0.528961729765e-02,  0.139166589782e-04,
+          -0.238855693017e-07,  0.356916001063e-10, -0.462347666298e-13,
+           0.500777441034e-16, -0.373105886191e-19,  0.157716482367e-22,
+          -0.281038625251e-26                                           ]
+      elsif t90 < 1664.50
+        [  0.295157925316e+01, -0.252061251332e-02,  0.159564501865e-04,
+          -0.764085947576e-08,  0.205305291024e-11, -0.293359668173e-15 ]
+      else
+        [  0.152232118209e+03, -0.268819888545e+00,  0.171280280471e-03,
+          -0.345895706453e-07, -0.934633971046e-14                      ]
+      end
+    end
+
+    def emf_c_coeffs_s t90
+      #  -50.00..1064.18
+      # 1064.18..1664.50
+      # 1664.50..1768.10
+      if t90 < 1064.18
+        [  0.000000000000e+00,  0.540313308631e-02,  0.125934289740e-04,
+          -0.232477968689e-07,  0.322028823036e-10, -0.331465196389e-13,
+           0.255744251786e-16, -0.125068871393e-19,  0.271443176145e-23 ]
+      elsif t90 < 1664.50
+        [  0.132900444085e+01,  0.334509311344e-02,  0.654805192818e-05,
+          -0.164856259209e-08,  0.129989605174e-13                      ]
+      else
+        [  0.146628232636e+03, -0.258430516752e+00,  0.163693574641e-03,
+          -0.330439046987e-07, -0.943223690612e-14                      ]
+      end
+    end
+
+    def emf_c_coeffs_t t90
+      # -270.0..0.0
+      #    0.0..400.0
+      if t90 < 0.0
+        [  0.000000000000e+00,  0.387481063640e-01,  0.441944343470e-04,
+           0.118443231050e-06,  0.200329735540e-07,  0.901380195590e-09,
+           0.226511565930e-10,  0.360711542050e-12,  0.384939398830e-14,
+           0.282135219250e-16,  0.142515947790e-18,  0.487686622860e-21,
+           0.107955392700e-23,  0.139450270620e-26,  0.797951539270e-30 ]
+      else
+        [  0.000000000000e+00,  0.387481063640e-01,  0.332922278800e-04,
+           0.206182434040e-06, -0.218822568460e-08,  0.109968809280e-10,
+          -0.308157587720e-13,  0.454791352900e-16, -0.275129016730e-19 ]
+      end
+    end
+
+    def emf_a_coeffs t90, type
+      case type
+      when :b then emf_a_coeffs_b t90
+      when :c then emf_a_coeffs_c t90
+      when :e then emf_a_coeffs_e t90
+      when :j then emf_a_coeffs_j t90
+      when :k then emf_a_coeffs_k t90
+      when :n then emf_a_coeffs_n t90
+      when :r then emf_a_coeffs_r t90
+      when :s then emf_a_coeffs_s t90
+      when :t then emf_a_coeffs_t t90
+      end
+    end
+
+    def emf_a_coeffs_b t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+    def emf_a_coeffs_c t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+    def emf_a_coeffs_e t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+    def emf_a_coeffs_j t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+    def emf_a_coeffs_k t90
+      if t90 < 0.0
+        [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+      else
+        [ 0.1185976e+00, -0.1183432e-03,  0.1269686e+03 ]
+      end
+    end
+
+    def emf_a_coeffs_n t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+    def emf_a_coeffs_r t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+    def emf_a_coeffs_s t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+    def emf_a_coeffs_t t90
+      [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00 ]
+    end
+
+
+    def t90_coeffs emf, type
+      case type
+      when :b then t90_coeffs_b emf
+      when :e then t90_coeffs_e emf
+      when :j then t90_coeffs_j emf
+      when :k then t90_coeffs_k emf
+      when :n then t90_coeffs_n emf
+      when :r then t90_coeffs_r emf
+      when :s then t90_coeffs_s emf
+      when :t then t90_coeffs_t emf
+      end
+    end
+
+    def t90_coeffs_b emf
+      # 0.290..2.431
+      # 2.431..13.821
+      if emf < 2.431
+        [  9.8423321e+01,  6.9971500e+02, -8.4765304e+02,
+           1.0052644e+03, -8.3345952e+02,  4.5508542e+02,
+          -1.5523037e+02,  2.9886750e+01, -2.4742860e+00 ]
+      else
+        [  2.1315071e+02,  2.8510504e+02, -5.2742887e+01,
+           9.9160804e+00, -1.2965303e+00,  1.1195870e-01,
+          -6.0625199e-03,  1.8661696e-04, -2.4878585e-06 ]
+      end
+    end
+
+    def t90_coeffs_e emf
+      # -8.825..0.000
+      #  0.000..76.373
+      if emf < 0.0
+        [  0.0000000e+00,  1.6977288e+01, -4.3514970e-01,
+          -1.5859697e-01, -9.2502871e-02, -2.6084314e-02,
+          -4.1360199e-03, -3.4034030e-04, -1.1564890e-05,
+           0.0000000e+00                                 ]
+      else
+        [  0.0000000e+00,  1.7057035e+01, -2.3301759e-01,
+           6.5435585e-03, -7.3562749e-05, -1.7896001e-06,
+           8.4036165e-08, -1.3735879e-09,  1.0629823e-11,
+          -3.2447087e-14                                 ]
+      end
+    end
+
+    def t90_coeffs_j emf
+      # -8.096..0.000
+      #  0.000..42.919
+      # 42.919..69.554
+      if emf < 0.0
+        [  0.0000000e+00,   1.9528268e+01,  -1.2286185e+00,
+          -1.0752178e+00,  -5.9086933e-01,  -1.7256713e-01,
+          -2.8131513e-02,  -2.3963370e-03,  -8.3823321e-05  ]
+      elsif emf < 42.919
+        [  0.000000e+00,    1.978425e+01,   -2.001204e-01,
+           1.036969e-02,   -2.549687e-04,    3.585153e-06,
+          -5.344285e-08,    5.099890e-10                    ]
+      else
+        [ -3.11358187e+03,  3.00543684e+02, -9.94773230e+00,
+           1.70276630e-01, -1.43033468e-03,  4.73886084e-06 ]
+      end
+    end
+
+    def t90_coeffs_k emf
+      # -5.892..0.000
+      #  0.000..20.644
+      # 20.644..54.887
+      if emf < 0.0
+        [  0.0000000e+00,  2.5173462e+01, -1.1662878e+00,
+          -1.0833638e+00, -8.9773540e-01, -3.7342377e-01,
+          -8.6632643e-02, -1.0450598e-02, -5.1920577e-04  ]
+      elsif emf < 20.644
+        [  0.000000e+00,   2.508355e+01,   7.860106e-02,
+          -2.503131e-01,   8.315270e-02,  -1.228034e-02,
+           9.804036e-04,  -4.413030e-05,   1.057734e-06,
+          -1.052755e-08                                   ]
+      else
+        [ -1.318058e+02,   4.830222e+01,  -1.646031e+00,
+           5.464731e-02,  -9.650715e-04,   8.802193e-06,
+          -3.110810e-08                                   ]
+      end
+    end
+
+    def t90_coeffs_n emf
+      # -3.991..0.000
+      #  0.000..20.613
+      # 20.613..47.513
+      if emf < 0.0
+        [  0.0000000e+00,  3.8436847e+01,  1.1010485e+00,
+           5.2229312e+00,  7.2060525e+00,  5.8488586e+00,
+           2.7754916e+00,  7.7075166e-01,  1.1582665e-01,
+           7.3138868e-03                                  ]
+      elsif emf < 20.613
+        [  0.00000e+00,    3.86896e+01,   -1.08267e+00,
+           4.70205e-02,   -2.12169e-06,   -1.17272e-04,
+           5.39280e-06,   -7.98156e-08                    ]
+      else
+        [  1.972485e+01,   3.300943e+01,  -3.915159e-01,
+           9.855391e-03,  -1.274371e-04,   7.767022e-07   ]
+      end
+    end
+
+    def t90_coeffs_r emf
+      # -0.227..1.923
+      #  1.923..13.228
+      # 11.361..19.739
+      # 19.739..21.103
+      if emf < 1.923
+        [  0.0000000E+00,    1.8891380E+02,   -9.3835290E+01,
+           1.3068619E+02,   -2.2703580E+02,    3.5145659E+02,
+          -3.8953900E+02,    2.8239471E+02,   -1.2607281E+02,
+           3.1353611E+01,   -3.3187769E+00                     ]
+      elsif emf < 13.228
+        [  1.334584505E+01,  1.472644573E+02, -1.844024844E+01,
+           4.031129726E+00, -6.249428360E-01,  6.468412046E-02,
+          -4.458750426E-03,  1.994710149E-04, -5.313401790E-06,
+           6.481976217E-08                                     ]
+      elsif emf < 19.739
+        [ -8.199599416E+01,  1.553962042E+02, -8.342197663E+00,
+           4.279433549E-01, -1.191577910E-02,  1.492290091E-04 ]
+      else
+        [  3.406177836E+04, -7.023729171E+03,  5.582903813E+02,
+          -1.952394635E+01,  2.560740231E-01                   ]
+      end
+    end
+
+    def t90_coeffs_s emf
+      # -0.236..1.874
+      #  1.874..11.950
+      # 10.332..17.536
+      # 17.536..18.693
+      if emf < 1.874
+        [  0.00000000E+00,   1.84949460E+02,  -8.00504062E+01,
+           1.02237430E+02,  -1.52248592E+02,   1.88821343E+02,
+          -1.59085941E+02,   8.23027880E+01,  -2.34181944E+01,
+           2.79786260E+00                                      ]
+      elsif emf < 11.950
+        [  1.291507177E+01,  1.466298863E+02, -1.534713402E+01,
+           3.145945973E+00, -4.163257839E-01,  3.187963771E-02,
+          -1.291637500E-03,  2.183475087E-05, -1.447379511E-07,
+           8.211272125E-09                                     ]
+      elsif emf < 17.536
+        [ -8.087801117E+01,  1.621573104E+02, -8.536869453E+00,
+           4.719686976E-01, -1.441693666E-02,  2.081618890E-04 ]
+      else
+        [  5.333875126E+04, -1.235892298E+04,  1.092657613E+03,
+          -4.265693686E+01,  6.247205420E-01                   ]
+      end
+    end
+
+    def t90_coeffs_t emf
+      # -5.603..0.000
+      #  0.000..20.872
+      if emf < 0.0
+        [  0.0000000E+00,  2.5949192E+01, -2.1316967E-01,
+           7.9018692E-01,  4.2527777E-01,  1.3304473E-01,
+           2.0241446E-02,  1.2668171E-03                 ]
+      else
+        [  0.000000E+00,   2.592800E+01,  -7.602961E-01,
+           4.637791E-02,  -2.165394E-03,   6.048144E-05,
+          -7.293422E-07                                  ]
+      end
+    end
+  end
+end

data/lib/tuscan/iec60751.rb

@@ -0,0 +1,35 @@
+module Tuscan
+  module Iec60751
+    extend self
+
+    def r t90, r0: 100.0, a:3.9083e-03, b:-5.7750e-07, c:-4.1830e-12
+      t90 >= 0 ?
+        r0*(1 + a*t90 + b*t90**2) :
+        r0*(1 + a*t90 + b*t90**2 - 100*c*t90**3 + c*t90**4)
+    end
+    alias_method :res, :r
+    alias_method :resistance, :r
+
+    def t90 r, r0: 100.0, a:3.9083e-03, b:-5.7750e-07, c:-4.1830e-12, err_limit: 1e-4, num_iter: 10
+      return 0 if r == r0
+      t = t90_approximation r
+      return Float::NAN if t.is_a? Complex
+      num_iter.times do
+        slope = (r - r0) / t
+        r_calc = r(t)
+        break if (r_calc - r).abs < slope * err_limit
+        t -= (r_calc - r) / slope
+      end
+      t
+    end
+    alias_method :t, :t90
+    alias_method :temperature, :t90
+
+  private
+    def t90_approximation r, r0: 100.0, a:3.9083e-03, b:-5.7750e-07, c:-4.1830e-12
+      r >= r0 ?
+        (-a + (a**2 - 4 * b * (1 - r / r0))**(0.5)) / (2 * b) :
+        ((r / r0) - 1) / (a + 100 * b)
+    end
+  end
+end

data/lib/tuscan/its90.rb

@@ -0,0 +1,107 @@
+module Tuscan
+  module Its90
+    include Math
+    extend self
+
+    A = [ -2.13534729, 3.1832472,  -1.80143597,
+           0.71727204, 0.50344027, -0.61899395,
+          -0.05332322, 0.28021362,  0.10715224,
+          -0.29302865, 0.04459872,  0.11868632,
+          -0.05248134 ]
+
+    B = [  0.183324722,  0.240975303,  0.209108771,
+           0.190439972,  0.142648498,  0.077993465,
+           0.012475611, -0.032267127, -0.075291522,
+          -0.05647067,   0.076201285,  0.123893204,
+          -0.029201193, -0.091173542,  0.001317696,
+           0.026025526 ]
+
+    C = [  2.78157254,  1.64650916, -0.1371439,
+          -0.00649767, -0.00234444,  0.00511868,
+           0.00187982, -0.00204472, -0.00046122,
+           0.00045724 ]
+
+    D = [ 439.932854, 472.41802, 37.684494,
+            7.472018,   2.920828, 0.005184,
+           -0.963864,  -0.188732, 0.191203,
+            0.049025 ]
+
+    def t90 r, rtpw:, **args
+      w = r / rtpw
+      Its90.t90r w - wdev(Its90.t90r(w), args)
+    end
+    alias_method :t, :t90
+    alias_method :temperature, :t90
+
+    # SUBRANGES = 1..11
+
+    # T90_RANGE  = -259.4467..961.88
+    # WR_RANGE   =  0.00119007..4.28642053
+
+    #
+    # VALID TEMPERATURE RANGES
+    #
+    # def range
+    #   WDEV_EQUATIONS[subrange - 1][:valid]
+    # end
+
+
+    #
+    # ITS-90 REFERENCE FUNCTIONS
+    #
+    def self.wr t90
+      # return nil unless T90_RANGE.include? t90
+      if t90< 0.01
+        exp(A.each_with_index.collect{ |a, i| a*( (log((t90+273.15)/273.16) + 1.5)/1.5 )**i }.reduce(:+))
+      else
+        C.each_with_index.collect{ |c, i| c*((t90-481)/481)**i }.reduce(:+)
+      end
+    end
+
+    def self.t90r wr
+      # return nil unless WR_RANGE.include? wr
+      if wr < 1.0
+        B.each_with_index.collect{ |b, i| b*(((wr)**(1.0/6)-0.65)/0.35)**i }.reduce(:+) * 273.16 - 273.15
+      else
+        D.each_with_index.collect{ |d, i| d*((wr-2.64)/1.64)**i }.reduce(:+)
+      end
+    end
+
+    #
+    # ITS-90 DEVIATION FUNCTIONS
+    #
+    WDEV_EQUATIONS = [ { valid: (-259.4467..0.01),   k: %w(c1 c2 c3 c4 c5), n: 2 }, 
+                       { valid: (-248.5939..0.01),   k: %w(c1 c2 c3      ), n: 0 },
+                       { valid: (-218.7916..0.01),   k: %w(c1            ), n: 1 },
+                       { valid: (-189.3442..0.01),   k: %w(              )       },
+                       { valid: (0.0..961.78),       k: %w(a  b  c  d    )       },
+                       { valid: (0.0..660.323),      k: %w(a  b  c       )       },
+                       { valid: (0.0..419.527),      k: %w(a  b          )       },
+                       { valid: (0.0..231.928),      k: %w(a  b          )       },
+                       { valid: (0.0..156.5985),     k: %w(a             )       },
+                       { valid: (0.0..29.7646),      k: %w(a             )       },
+                       { valid: (-38.8344..29.7646), k: %w(a  b          )       } ]
+
+    def wdev t90, subrange: ,
+                  a: 0.0, b: 0.0, c: 0.0, d: 0.0, w660: 0.0,
+                  c1: 0.0, c2: 0.0, c3: 0.0, c4: 0.0, c5: 0.0
+      equation = WDEV_EQUATIONS[subrange - 1]
+      # return nil unless equation[:valid].include? t90
+      wr_t90 = Its90.wr t90
+      case subrange
+      when 1..4
+        wdev  = a * (wr_t90 - 1)
+        subrange == 4 ?
+          wdev += b * (wr_t90 - 1) * log(wr_t90) :
+          wdev += b * (wr_t90 - 1)**2
+        wdev += equation[:k].each_with_index.
+                  collect{ |k, i| eval("#{k}") * log(wr_t90)**(i + equation[:n]) }.reduce(:+) || 0
+      when 5..11
+        wdev   = d * (wr_t90 - w660)**2 if equation[:k].delete('d')
+        wdev ||= 0
+        wdev  += equation[:k].each_with_index.
+                    collect{ |k, i| eval("#{k}") * (wr_t90 - 1)**(i+1) }.reduce(:+)
+      end
+    end
+  end
+end

data/lib/tuscan/version.rb

@@ -0,0 +1,3 @@
+module Tuscan
+  VERSION = "0.0.1"
+end

data/lib/tuscan.rb

@@ -0,0 +1,14 @@
+require 'tuscan/version'
+require 'tuscan/its90'
+require 'tuscan/iec60751'
+require 'tuscan/iec60584'
+
+module Tuscan
+  extend self
+
+private
+  def method_missing method, *args, &block 
+    subject, *remaining_args = args
+    const_get(subject.to_s.capitalize).send method, *remaining_args, &block
+  end
+end

data/spec/spec_helper.rb

@@ -0,0 +1 @@
+require 'tuscan'

data/spec/tuscan/iec60584_spec.rb

@@ -0,0 +1,116 @@
+require 'spec_helper'
+
+#
+# TODO:
+#   1. compute inverse of emfr by iterative process (will pass R, B examples)
+#   2. get aproximate t90 valid fro -270.0..-200.0 (will pass E, K, N, T examples)
+#   3. implement deviation function (polynomial)
+#
+module Tuscan
+  describe Iec60584 do
+    context 'reference functions' do
+      examples = {
+        b: [
+          # { emf:  0.000, t90:   0.0 },
+          { emf:  0.291, t90: 250.0 },
+          { emf:  1.975, t90: 630.0 },
+          { emf:  1.981, t90: 631.0 },
+          { emf:  5.780, t90:1100.0 },
+          { emf: 13.820, t90:1820.0 }
+        ],
+        c: [
+          { emf:  0.000, t90:   0.0 },
+          { emf:  1.451, t90: 100.0 },
+          { emf: 11.778, t90: 660.0 },
+          { emf: 20.066, t90:1100.0 },
+          { emf: 29.402, t90:1680.0 },
+          { emf: 37.107, t90:2320.0 }
+        ],
+        e: [
+          # { emf: -9.835, t90:-270.0 },
+          { emf: -8.825, t90:-200.0 },
+          { emf: -5.237, t90:-100.0 },
+          { emf:  0.000, t90:   0.0 },
+          { emf:  6.319, t90: 100.0 },
+          { emf: 49.917, t90: 660.0 },
+          { emf: 76.373, t90:1000.0 }
+        ],
+        j: [
+          { emf: -8.095, t90:-210.0 },
+          { emf: -4.633, t90:-100.0 },
+          { emf:  0.000, t90:   0.0 },
+          { emf: 36.675, t90: 660.0 },
+          { emf: 43.559, t90: 770.0 },
+          { emf: 69.553, t90:1200.0 }
+        ],
+        k: [
+          # { emf: -6.458, t90:-270.0 },
+          { emf: -5.891, t90:-200.0 },
+          { emf: -3.554, t90:-100.0 },
+          { emf:  0.000, t90:   0.0 },
+          { emf:  5.206, t90: 127.0 },
+          { emf: 27.447, t90: 660.0 },
+          { emf: 54.886, t90:1372.0 }
+        ],
+        n: [
+          # { emf: -4.345, t90:-270.0 },
+          { emf: -3.990, t90:-200.0 },
+          { emf: -2.407, t90:-100.0 },
+          { emf:  0.000, t90:   0.0 },
+          { emf:  2.774, t90: 100.0 },
+          { emf: 22.958, t90: 660.0 },
+          { emf: 47.513, t90:1300.0 }
+        ],
+        r: [
+          # { emf: -0.226, t90: -50.0 },
+          { emf: -0.123, t90: -25.0 },
+          { emf:  0.000, t90:   0.0 },
+          { emf:  6.273, t90: 660.0 },
+          { emf: 11.850, t90:1100.0 },
+          { emf: 21.101, t90:1768.0 }
+        ],
+        s: [
+          { emf: -0.236, t90: -50.0 },
+          { emf: -0.127, t90: -25.0 },
+          { emf:  0.000, t90:   0.0 },
+          { emf:  5.857, t90: 660.0 },
+          { emf: 10.757, t90:1100.0 },
+          { emf: 18.693, t90:1768.0 }
+        ],
+        t: [
+          # { emf: -6.258, t90:-270.0 },
+          { emf: -5.603, t90:-200.0 },
+          { emf: -3.379, t90:-100.0 },
+          { emf:  0.000, t90:   0.0 },
+          { emf:  4.279, t90: 100.0 },
+          { emf:  9.288, t90: 200.0 },
+          { emf: 20.872, t90: 400.0 }
+        ]
+      }
+
+      context 'emfr computation' do
+        %i{ b c e j k n r s t }.each do |type|
+          context "on a type #{type.upcase} thermocouple" do
+            examples[type].each do |example|
+              it "yields #{example[:emf]} mV when t90 equals #{example[:t90]} ºC" do
+                expect(Iec60584.emfr(example[:t90], type)).to be_within(1e-3).of(example[:emf])
+              end
+            end
+          end
+        end
+      end
+
+      context 't90r computation' do
+        %i{ b e j k n r s t }.each do |type|
+          context "on a type #{type.upcase} thermocouple" do
+            examples[type].each do |example|
+              it "yields #{example[:t90]} ºC when emf equals #{example[:emf]} mV" do
+                expect(Iec60584.t90r(example[:emf], type)).to be_within(0.1).of(example[:t90])
+              end
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/spec/tuscan/iec60751_spec.rb

@@ -0,0 +1,47 @@
+require 'spec_helper'
+
+module Tuscan
+  describe Iec60751 do
+    examples = [ { r:  18.5201, t90: -200.0000 },
+                 { r:  60.2558, t90: -100.0000 },
+                 { r: 100.0000, t90:    0.0000 },
+                 { r: 247.0920, t90:  400.0000 },
+                 { r: 390.4811, t90:  850.0000 } ]
+
+    context 'resistance computation' do
+      context 'standard coefficients' do
+        examples.each do |example|
+          it "yields #{example[:r]} Ohm when t90 equals #{example[:t90]} Celsius" do
+            expect(Iec60751.r example[:t90]).to be_within(1e-4).of(example[:r])
+          end
+        end
+      end
+
+      context 'non-standard coefficients' do
+        it 'yields 101.0 Ohm when temperature equals 0.0 Celsius for r0 = 101.0 Ohm' do
+          expect(Iec60751.r 0, r0: 101.0).to be_within(1e-4).of(101)
+        end
+      end
+    end
+
+    context 'temperature computation' do
+      it 'yields complex temperatures as NaN' do
+        expect(Iec60751.t90 1000).to be Float::NAN
+      end
+
+      context 'standard coefficients' do
+        examples.each do |example|
+          it "yields #{example[:t90]} Celius when resistance equals #{example[:r]} Ohm" do
+            expect(Iec60751.t90 example[:r]).to be_within(1e-4).of(example[:t90])
+          end
+        end
+      end
+
+      context 'non-standard coefficients' do
+        it 'yields 0.0 Celsius when resistance equals 101.0 Ohm on a FUNCTION with r0 = 101.0 Ohm' do
+          expect(Iec60751.t90 101, r0: 101.0).to be_within(1e-4).of(0)
+        end
+      end
+    end
+  end
+end

data/spec/tuscan/its90_spec.rb

@@ -0,0 +1,91 @@
+require 'spec_helper'
+
+module Tuscan
+  describe Its90 do
+    context 'reference funtions' do
+      examples = [
+        { fixed_point: 'e-H2 Tp', t90: -259.3467, wr: 0.00119007 },
+        { fixed_point: 'Ne Tp',   t90: -248.5939, wr: 0.00844974 },
+        { fixed_point: 'O2 Tp',   t90: -218.7916, wr: 0.09171804 },
+        { fixed_point: 'Ar Tp',   t90: -189.3442, wr: 0.21585975 },
+        { fixed_point: 'Hg Tp',   t90:  -38.8344, wr: 0.84414211 },
+        { fixed_point: 'H2O Tp',  t90:    0.01,   wr: 1.00000000 },
+        { fixed_point: 'Ga Mp',   t90:   29.7646, wr: 1.11813889 },
+        { fixed_point: 'In Fp',   t90:  156.5985, wr: 1.60980185 },
+        { fixed_point: 'Sn Fp',   t90:  231.9280, wr: 1.89279768 },
+        { fixed_point: 'Zn Fp',   t90:  419.5270, wr: 2.56891730 },
+        { fixed_point: 'Al Fp',   t90:  660.3230, wr: 3.37600860 },
+        { fixed_point: 'Ag Fp',   t90:  961.7800, wr: 4.28642053 }
+      ]
+
+      context 'wr computation' do
+        examples.each do |example|
+          it "complies on #{example[:fixed_point]} fixed point" do
+            expect(Its90.wr example[:t90]).to be_within(1e-8).of(example[:wr])
+          end
+        end
+      end
+
+      context 't90r computation' do
+        examples.each do |example|
+          it "complies on #{example[:fixed_point]} fixed point" do
+            expect(Its90.t90r example[:wr]).to be_within(0.00013).of(example[:t90])
+          end
+        end
+      end
+    end
+
+    context 'deviation functions' do
+      context 'wdev computation' do
+        context 'range 4' do
+          # Examples valid for sub-range = 4, a = -1.2579994e-04, b = 1.0678395e-05 as per NIST SP250-81
+          sprt = { subrange: 4, a: -1.2579994e-04, b: 1.0678395e-05 }
+          examples = [
+            { t90: -189.3442, wdev: 0.000111482 },
+            { t90: -100.0000, wdev: 0.000053258 },
+            { t90:  -38.8344, wdev: 0.000019889 },
+            { t90:    0.0000, wdev: 0.000000005 }
+          ]
+
+          examples.each do |example|
+            it "complies with NIST SP250-81 example on range 4, #{example[:t90]} Celsius" do
+              expect(Its90.wdev example[:t90], sprt).to be_within(1e-8).of(example[:wdev])
+            end
+          end
+        end
+
+        context 'range 6' do
+          sprt = { subrange: 6, a: -1.6462789e-04, b: -8.4598339e-06, c:  1.8898584e-06 }
+          examples = [
+            { t90:   0.000, wdev:  0.000000007},
+            { t90: 100.000, wdev: -0.000065852},
+            { t90: 231.928, wdev: -0.000152378},
+            { t90: 419.527, wdev: -0.000271813},
+            { t90: 660.323, wdev: -0.000413567}
+          ]
+
+          examples.each do |example|
+            it "complies with NIST SP250-81 example on range 6, #{example[:t90]} Celsius" do
+              expect(Its90.wdev example[:t90], sprt).to be_within(1e-8).of(example[:wdev])
+            end
+          end
+        end
+      end
+    end
+
+    context 't90 function' do
+      sprt = { rtpw: 25.319871, subrange: 7, a: -1.2134e-04, b: -9.9190e-06 }
+      examples = [
+        { r: 25.319871, t90:   0.010 },
+        { r: 47.922451, t90: 231.928 },
+        { r: 65.039218, t90: 419.527 }
+      ]
+
+      examples.each do |example|
+        it "complies with IPQ cert. 501.20/1241312 range 7, #{example[:t90]} Celsius" do
+          expect(Its90.t90 example[:r], sprt).to be_within(0.0001).of(example[:t90])
+        end
+      end
+    end
+  end
+end

data/spec/tuscan_spec.rb

@@ -0,0 +1,8 @@
+require 'spec_helper'
+
+describe Tuscan do
+  it 'forwards methods to child modules' do
+    expect(Tuscan::Iec60751).to receive(:t).with(r: 101, r0: 100)
+    Tuscan.t :iec60751, r: 101, r0: 100
+  end
+end

data/tasks/rspec.rake

@@ -0,0 +1,3 @@
+require 'rspec/core/rake_task'
+
+RSpec::Core::RakeTask.new(:spec)

data/tuscan.gemspec

@@ -0,0 +1,25 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'tuscan/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "tuscan"
+  spec.version       = Tuscan::VERSION
+  spec.authors       = ["Luis Bacelar"]
+  spec.email         = ["lmbacelar@gmail.com"]
+  spec.summary       = %q{Converts to and from temperature for standard (S)PRT's and thermocouples}
+  spec.description   = %q{Converts from electrical signal to temperature and vice-versa for ITS-90 SPRT's, IEC 60751 PRT's and IEC 60584 thermocouples}
+  spec.homepage      = ""
+  spec.license       = "MIT"
+
+  spec.files         = `git ls-files -z`.split("\x0")
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", "~> 1.7"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency 'rspec'
+  spec.add_development_dependency 'byebug'
+end

metadata

@@ -0,0 +1,124 @@
+--- !ruby/object:Gem::Specification
+name: tuscan
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Luis Bacelar
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2015-01-13 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.7'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.7'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: byebug
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: Converts from electrical signal to temperature and vice-versa for ITS-90
+  SPRT's, IEC 60751 PRT's and IEC 60584 thermocouples
+email:
+- lmbacelar@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- lib/tuscan.rb
+- lib/tuscan/iec60584.rb
+- lib/tuscan/iec60751.rb
+- lib/tuscan/its90.rb
+- lib/tuscan/version.rb
+- spec/spec_helper.rb
+- spec/tuscan/iec60584_spec.rb
+- spec/tuscan/iec60751_spec.rb
+- spec/tuscan/its90_spec.rb
+- spec/tuscan_spec.rb
+- tasks/rspec.rake
+- tuscan.gemspec
+homepage: ''
+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.4.5
+signing_key: 
+specification_version: 4
+summary: Converts to and from temperature for standard (S)PRT's and thermocouples
+test_files:
+- spec/spec_helper.rb
+- spec/tuscan/iec60584_spec.rb
+- spec/tuscan/iec60751_spec.rb
+- spec/tuscan/its90_spec.rb
+- spec/tuscan_spec.rb