eph_bpn 0.1.1 → 0.1.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 0d37575850f490a07120adab907b267cf602104f
-  data.tar.gz: 5ca94950a2bd3d28a6487fc731e8608b4895cd58
+  metadata.gz: 82ade37ba8efd0aa1c294324fddc96d79229a5cf
+  data.tar.gz: b4b1b3d03b1f41edc3f7596270100a1a88cc9f38
 SHA512:
-  metadata.gz: 6dd46b12e43ed9a38f98e70c43b94c0d1ef5c6e54fa0b29cbc4ce37763596063c785587b27f20e2d8e169b11e4f92340fcca0f0111866f32d4b2f8e2eef13ca1
-  data.tar.gz: 03b649e1f06eee7c89ac8f0bfd12dfa7f33f343a2fbeb74f36db01d83a592060d59eb5a18b12a8753d45899e5f39d3ced6e78e1415a71f851e085da2f69fc880
+  metadata.gz: 4c885105a0264105eb4ce7143e65f6a76894ebda64c1666cd0d75298d2fc2bc0bd9ab03581be55937d5a2028c1d6dc291d83ecffdfb6de1e4f48915f0df9a08f
+  data.tar.gz: 9477203afe7ff658c6529fb2c00058cc553942cc746f8e0697c986754930aeb9b8c7ed771ef1b0a91cf44520762b77003135c026b3a99d2e615c92235f01e9db

data/README.md

@@ -24,6 +24,8 @@ Or install it yourself as:
 
 ### Instantiation
 
+    require 'eph_bpn'
+    
     e = EphBpn.new("20160723")
     e = EphBpn.new("20160723123456")
     e = EphBpn.new
@@ -33,18 +35,20 @@ Or install it yourself as:
 
 ### Calculation
 
-    r_i = [-0.50787065, 0.80728228, 0.34996714]  # ICRS Coordinate
-    r_b = e.apply_bias(r_i)  # Apply Bias
-    r_p = e.apply_prec(r_b)  # Apply Precession
-    r_n = e.apply_nut(r_p)   # Apply Nutation
-    puts "TDB: #{e.tdb}"     #=> TDB: 2016-07-23 00:00:00 +0000
-    puts " JD: #{e.jd}"      #=>  JD: 2457592.5
-    puts " JC: #{e.jc}"      #=>  JC: 0.16557152635181382
-    puts "EPS: #{e.eps}"     #=> EPS: 0.40905500411767176
-    p r_i                    #=> [-0.50787065, 0.80728228, 0.34996714]
-    p r_b                    #=> [-0.5078706789483137, 0.8072822556207571, 0.34996715422685276]
-    p r_p                    #=> [-0.5114184398597698, 0.8053953379861463, 0.3491472536549302]
-    p r_n                    #=> [-0.5114184385245731, 0.8053953403876852, 0.3491472500709293]
+    pos_g  = [-0.50787065, 0.80728228, 0.34996714]  # GCRS Coordinate
+    pos_b  = e.apply_bias(pos_g)       # Apply Bias
+    pos_p  = e.apply_prec(pos_b)       # Apply Precession
+    pos_n  = e.apply_nut(pos_p)        # Apply Nutation
+    pos_bp = e.apply_bias_prec(pos_g)  # Apply Bias + Precession
+    puts "TDB: #{e.tdb}"               #=> TDB: 2016-07-23 00:00:00 +0000
+    puts " JD: #{e.jd}"                #=>  JD: 2457592.5
+    puts " JC: #{e.jc}"                #=>  JC: 0.16557152635181382
+    puts "EPS: #{e.eps}"               #=> EPS: 0.40905500411767176
+    p pos_g                            #=> [-0.50787065, 0.80728228, 0.34996714]
+    p pos_b                            #=> [-0.507870315369686, 0.8072824634004779, 0.34996720255697145]
+    p pos_p                            #=> [-0.5114180771311418, 0.8053955471104202, 0.34914730256926324]
+    p pos_n                            #=> [-0.5114180757959449, 0.8053955495119587, 0.34914729898526137]
+    p pos_bp                           #=> [-0.5114184398598124, 0.805395337986056, 0.3491472536550768]
 
 ## Development
 

data/exe/eph_bpn

@@ -3,16 +3,18 @@
 require "eph_bpn"
 
 e = EphBpn.new(ARGV[0])
-r_i = [-0.50787065, 0.80728228, 0.34996714]  # ICRS Coordinate
-r_b = e.apply_bias(r_i)  # Apply Bias
-r_p = e.apply_prec(r_b)  # Apply Precession
-r_n = e.apply_nut(r_p)   # Apply Nutation
-puts "TDB: #{e.tdb}"     #=> TDB: 2016-07-23 00:00:00 +0000
-puts " JD: #{e.jd}"      #=>  JD: 2457592.5
-puts " JC: #{e.jc}"      #=>  JC: 0.16557152635181382
-puts "EPS: #{e.eps}"     #=> EPS: 0.40905500411767176
-p r_i                    #=> [-0.50787065, 0.80728228, 0.34996714]
-p r_b                    #=> [-0.5078706789483137, 0.8072822556207571, 0.34996715422685276]
-p r_p                    #=> [-0.5114184398597698, 0.8053953379861463, 0.3491472536549302]
-p r_n                    #=> [-0.5114184385245731, 0.8053953403876852, 0.3491472500709293]
+pos_g  = [-0.50787065, 0.80728228, 0.34996714]  # GCRS Coordinate
+pos_b  = e.apply_bias(pos_g)       # Apply Bias
+pos_p  = e.apply_prec(pos_b)       # Apply Precession
+pos_n  = e.apply_nut(pos_p)        # Apply Nutation
+pos_bp = e.apply_bias_prec(pos_g)  # Apply Bias + Precession
+puts "TDB: #{e.tdb}"               #=> TDB: 2016-07-23 00:00:00 +0000
+puts " JD: #{e.jd}"                #=>  JD: 2457592.5
+puts " JC: #{e.jc}"                #=>  JC: 0.16557152635181382
+puts "EPS: #{e.eps}"               #=> EPS: 0.40905500411767176
+p pos_g                            #=> [-0.50787065, 0.80728228, 0.34996714]
+p pos_b                            #=> [-0.507870315369686, 0.8072824634004779, 0.34996720255697145]
+p pos_p                            #=> [-0.5114180771311418, 0.8053955471104202, 0.34914730256926324]
+p pos_n                            #=> [-0.5114180757959449, 0.8053955495119587, 0.34914729898526137]
+p pos_bp                           #=> [-0.5114184398598124, 0.805395337986056, 0.3491472536550768]
 

data/lib/eph_bpn/compute.rb

@@ -107,49 +107,49 @@ module EphBpn
     # @param:  <none>
     # @return: r  (回転行列)
     #=========================================================================
-    def r_generic
-      r = r_x( -5.1 * 180 * Const::MAS2R)
-      r = r_y(-17.3 * 180 * Const::MAS2R, r)
-      r = r_z( 78.0 * 180 * Const::MAS2R, r)
+    def r_bias
+      r = r_x( -5.1 * Const::MAS2R)
+      r = r_y(-17.3 * Const::MAS2R, r)
+      r = r_z( 78.0 * Const::MAS2R, r)
       return r
     rescue => e
       raise
     end
 
     #=========================================================================
-    # Bias 変換行列（IAU 2006 (Fukushima-Williams 4-angle formulation) 理論）
+    # Bias + Precession 変換行列
+    #
+    # * IAU 2006 (Fukushima-Williams 4-angle formulation) 理論
     #
     # @param:  <none>
     # @return: r  (変換行列)
     #=========================================================================
     def r_fw_iau_06
-      jc = jd2jc(2451545.0)
-      eps = compute_obliquity(jc)
       gamma = (-0.052928    + \
               (10.556378    + \
               ( 0.4932044   + \
               (-0.00031238  + \
               (-0.000002788 + \
               ( 0.0000000260) \
-              * jc) * jc) * jc) * jc) * jc) * Const::AS2R
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
       phi   = (84381.412819    + \
               (  -46.811016    + \
               (    0.0511268   + \
               (    0.00053289  + \
               (   -0.000000440 + \
               (   -0.0000000176) \
-              * jc) * jc) * jc) * jc) * jc) * Const::AS2R
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
       psi   = (  -0.041775    + \
               (5038.481484    + \
               (   1.5584175   + \
               (  -0.00018522  + \
               (  -0.000026452 + \
               (  -0.0000000148) \
-              * jc) * jc) * jc) * jc) * jc) * Const::AS2R
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
       r = r_z(gamma)
       r = r_x(phi,   r)
       r = r_z(-psi,  r)
-      r = r_x(-eps,  r)
+      r = r_x(-@eps, r)
       return r
     rescue => e
       raise

data/lib/eph_bpn/ephemeris.rb

@@ -17,8 +17,7 @@ module EphBpn
     # @return: pos_b  (適用後直角座標)
     #=========================================================================
     def apply_bias(pos)
-      #return rotate(r_generic, pos)   # 一般的な理論
-      return rotate(r_fw_iau_06, pos)  # IAU 2006 (Fukushima-Williams) 理論
+      return rotate(r_bias, pos)
     end
 
     #=========================================================================
@@ -31,6 +30,16 @@ module EphBpn
       return rotate(r_prec, pos)
     end
 
+    #=========================================================================
+    # Bias + Precession（歳差） 適用
+    #
+    # @param:  pos    (適用前直角座標)
+    # @return: pos_p  (適用後直角座標)
+    #=========================================================================
+    def apply_bias_prec(pos)
+      return rotate(r_fw_iau_06, pos)  # IAU 2006 (Fukushima-Williams) 理論
+    end
+
     #=========================================================================
     # Nutation（章動） 適用
     #

data/lib/eph_bpn/version.rb

@@ -1,3 +1,3 @@
 module EphBpn
-  VERSION = "0.1.1"
+  VERSION = "0.1.2"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: eph_bpn
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 0.1.2
 platform: ruby
 authors:
 - komasaru
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2016-08-24 00:00:00.000000000 Z
+date: 2016-09-01 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler