eph_bpn 0.1.2 → 0.1.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 82ade37ba8efd0aa1c294324fddc96d79229a5cf
-  data.tar.gz: b4b1b3d03b1f41edc3f7596270100a1a88cc9f38
+  metadata.gz: 255f966f8773fae9d08e41c638c378c4c0c6ae28
+  data.tar.gz: 4dbae44ebd99bcac4fd65ca9f9bcb1c82ecba6e5
 SHA512:
-  metadata.gz: 4c885105a0264105eb4ce7143e65f6a76894ebda64c1666cd0d75298d2fc2bc0bd9ab03581be55937d5a2028c1d6dc291d83ecffdfb6de1e4f48915f0df9a08f
-  data.tar.gz: 9477203afe7ff658c6529fb2c00058cc553942cc746f8e0697c986754930aeb9b8c7ed771ef1b0a91cf44520762b77003135c026b3a99d2e615c92235f01e9db
+  metadata.gz: 7ed2c8c2db252ccf2f905dd1056ab81b8c84ae0fcfcf82b3a1c441dd5fa6bf0b85c0e569024f24af3453b50f6b50de2cc49f180f9e7f4af89792d4e36b6fc63e
+  data.tar.gz: 9f38ffdbd97be6b5af71e3c8110420a37c0b4f2ad2fec4234d2aef89387228a7c72b769b13cec4aa8718becddac4ced33367d61d62b3591c9c805ec386db84a0

data/exe/eph_bpn

@@ -12,6 +12,10 @@ 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 e.r_bias                         # Rotation matrix for Bias
+p e.r_prec                         # Rotation matrix for Precession
+p e.r_nut                          # Rotation matrix for Nutation
+p e.r_bias_prec                    # Rotation matrix for Bias & precession
 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]

data/lib/eph_bpn/compute.rb

@@ -107,7 +107,7 @@ module EphBpn
     # @param:  <none>
     # @return: r  (回転行列)
     #=========================================================================
-    def r_bias
+    def comp_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)
@@ -116,45 +116,6 @@ module EphBpn
       raise
     end
 
-    #=========================================================================
-    # Bias + Precession 変換行列
-    #
-    # * IAU 2006 (Fukushima-Williams 4-angle formulation) 理論
-    #
-    # @param:  <none>
-    # @return: r  (変換行列)
-    #=========================================================================
-    def r_fw_iau_06
-      gamma = (-0.052928    + \
-              (10.556378    + \
-              ( 0.4932044   + \
-              (-0.00031238  + \
-              (-0.000002788 + \
-              ( 0.0000000260) \
-              * @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
-      psi   = (  -0.041775    + \
-              (5038.481484    + \
-              (   1.5584175   + \
-              (  -0.00018522  + \
-              (  -0.000026452 + \
-              (  -0.0000000148) \
-              * @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)
-      return r
-    rescue => e
-      raise
-    end
-
     #=========================================================================
     # precession（歳差）変換行列（J2000.0 用）
     #
@@ -180,7 +141,7 @@ module EphBpn
     # @param:  <none>
     # @return: r  (変換行列)
     #=========================================================================
-    def r_prec
+    def comp_r_prec
       gamma = ((10.556403    + \
                (0.4932044    + \
                (-0.00031238  + \
@@ -209,13 +170,52 @@ module EphBpn
       raise
     end
 
+    #=========================================================================
+    # Bias + Precession 変換行列
+    #
+    # * IAU 2006 (Fukushima-Williams 4-angle formulation) 理論
+    #
+    # @param:  <none>
+    # @return: r  (変換行列)
+    #=========================================================================
+    def comp_r_bias_prec
+      gamma = (-0.052928    + \
+              (10.556378    + \
+              ( 0.4932044   + \
+              (-0.00031238  + \
+              (-0.000002788 + \
+              ( 0.0000000260) \
+              * @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
+      psi   = (  -0.041775    + \
+              (5038.481484    + \
+              (   1.5584175   + \
+              (  -0.00018522  + \
+              (  -0.000026452 + \
+              (  -0.0000000148) \
+              * @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)
+      return r
+    rescue => e
+      raise
+    end
+
     #=========================================================================
     # nutation（章動）変換行列
     #
     # @param:  <none>
     # @return: r  (変換行列)
     #=========================================================================
-    def r_nut
+    def comp_r_nut
       dpsi_ls, deps_ls = compute_lunisolar
       dpsi_pl, deps_pl = compute_planetary
       dpsi, deps = dpsi_ls + dpsi_pl, deps_ls + deps_pl

data/lib/eph_bpn/ephemeris.rb

@@ -1,6 +1,7 @@
 module EphBpn
   class Ephemeris
-    attr_reader :tdb, :jd, :jc, :eps
+    attr_reader :tdb, :jd, :jc, :eps,
+                :r_bias, :r_prec, :r_nut, :r_bias_prec
     include Compute
 
     def initialize(tdb)
@@ -8,6 +9,10 @@ module EphBpn
       @jd  = gc2jd(@tdb)             # TDB -> JD（ユリウス日）
       @jc  = jd2jc(@jd)              # ユリウス世紀
       @eps = compute_obliquity(@jc)  # 平均黄道傾斜角
+      @r_bias      = comp_r_bias       # 回転行列（バイアス）
+      @r_prec      = comp_r_prec       # 回転行列（歳差）
+      @r_nut       = comp_r_nut        # 回転行列（章動）
+      @r_bias_prec = comp_r_bias_prec  # 回転行列（バイアス＆歳差）
     end
 
     #=========================================================================
@@ -17,7 +22,7 @@ module EphBpn
     # @return: pos_b  (適用後直角座標)
     #=========================================================================
     def apply_bias(pos)
-      return rotate(r_bias, pos)
+      return rotate(@r_bias, pos)
     end
 
     #=========================================================================
@@ -27,27 +32,27 @@ module EphBpn
     # @return: pos_p  (適用後直角座標)
     #=========================================================================
     def apply_prec(pos)
-      return rotate(r_prec, pos)
+      return rotate(@r_prec, pos)
     end
 
     #=========================================================================
-    # Bias + Precession（歳差） 適用
+    # Nutation（章動） 適用
     #
     # @param:  pos    (適用前直角座標)
-    # @return: pos_p  (適用後直角座標)
+    # @return: pos_n  (適用後直角座標)
     #=========================================================================
-    def apply_bias_prec(pos)
-      return rotate(r_fw_iau_06, pos)  # IAU 2006 (Fukushima-Williams) 理論
+    def apply_nut(pos)
+      return rotate(@r_nut, pos)
     end
 
     #=========================================================================
-    # Nutation（章動） 適用
+    # Bias + Precession（歳差） 適用
     #
     # @param:  pos    (適用前直角座標)
-    # @return: pos_n  (適用後直角座標)
+    # @return: pos_p  (適用後直角座標)
     #=========================================================================
-    def apply_nut(pos)
-      return rotate(r_nut, pos)
+    def apply_bias_prec(pos)
+      return rotate(@r_bias_prec, pos)  # IAU 2006 (Fukushima-Williams) 理論
     end
   end
 end

data/lib/eph_bpn/version.rb

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

metadata

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