eph_bpn 0.1.7 → 0.1.12

checksums.yaml

@@ -1,7 +1,7 @@
 ---
-SHA1:
-  metadata.gz: 3414ccc9da1aeda053afc02428f844734bc1ec93
-  data.tar.gz: f6de30dccbbd50590cc5d7b775f6ff56522f5188
+SHA256:
+  metadata.gz: 5c7d3bf9f89cb4c208fbd82ff4038958663638a9bd42d736c1006a8350b3e0c3
+  data.tar.gz: 43759361d503231634a6b901ac32d20b878a6b1c9fa90d9cd27797c372562b42
 SHA512:
-  metadata.gz: 12f6e701345507c6a27eba4ceb9619032e9fb18762595468b21a65884cdf5d6ee3fb5effabe4494405bbf5ab73b22cfeceafdd84209e626976430889def756fb
-  data.tar.gz: 9f683b67b2bc9bd6a87395fb8e944a552ccddb9ec04946c949ee3629f6aac228579cdbbff862101c8e3cbae1028dd5f9e9fac49af385d5f33d19bc11b046ff38
+  metadata.gz: a0691ed6d1065fc8e1916cc5aa004d4333c6e71491909fe942d771533a054428e47c67f496ffd722d5bb540eda34941ad0ad8d8ca49ea1dabad582003c4af723
+  data.tar.gz: 16161e02f0d135206a6d439da79d3ad90b7ab57cf63dda2f1d50312d83f07ad573dc07b0bb8416c8f8facf0f1be9522e85408dfdfd77a7b7de39ea3ba52c0866

data/README.md

@@ -35,20 +35,39 @@ Or install it yourself as:
 
 ### Calculation
 
-    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]
+    pos_g   = [-0.50787065, 0.80728228, 0.34996714]  # GCRS Coordinate
+    pos_b   = e.apply_bias(pos_g)           # Apply Bias
+    pos_bp  = e.apply_bias_prec(pos_g)      # Apply Bias & Precession
+    pos_bpn = e.apply_bias_prec_nut(pos_g)  # Apply Bias & Precession & Nutation
+    pos_p   = e.apply_prec(pos_b)           # Apply Precession
+    pos_pn  = e.apply_prec_nut(pos_b)       # Apply Precession & Nutation
+    pos_n   = e.apply_nut(pos_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 e.r_bias                              # Rotation matrix for Bias
+    p e.r_bias_prec                         # Rotation matrix for Bias & Precession
+    p e.r_bias_prec_nut                     # Rotation matrix for Bias & Precession & Nutation
+    p e.r_prec                              # Rotation matrix for Precession
+    p e.r_prec_nut                          # Rotation matrix for Precession & Nutation
+    p e.r_nut                               # Rotation matrix for Nutation
+    puts "SRC: #{pos_g}"                    #=> [-0.50787065, 0.80728228, 0.34996714]
+    puts "  B: #{pos_b}"                    #=> [-0.507870315369686, 0.8072824634004779, 0.34996720255697145]
+    puts " BP: #{pos_bp}"                   #=> [-0.5114184398598124, 0.805395337986056, 0.3491472536550768]
+    puts "BPN: #{pos_bpn}"                  #=> [-0.5114038678325663, 0.8054186649898355, 0.3491147860877318]
+    puts "  P: #{pos_p}"                    #=> [-0.5114180771311418, 0.8053955471104202, 0.34914730256926324]
+    puts " PN: #{pos_pn}"                   #=> [-0.5114035051003876, 0.8054188741108519, 0.34911483499021656]
+    puts "  N: #{pos_n}"                    #=> [-0.5114035051003876, 0.8054188741108514, 0.34911483499021656]
+
+About alias:
+
+* `apply_bias` => apply_b`
+* `apply_bias_prec` => apply_bp`
+* `apply_bias_prec_nut` => apply_bpn`
+* `apply_prec` => apply_p`
+* `apply_prec_nut` => apply_pn`
+* `apply_nut` => apply_n`
 
 ## Development
 

data/eph_bpn.gemspec

@@ -27,7 +27,9 @@ Gem::Specification.new do |spec|
   spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
   spec.require_paths = ["lib"]
 
-  spec.add_development_dependency "bundler", "~> 1.12"
-  spec.add_development_dependency "rake", "~> 10.0"
+  #spec.add_development_dependency "bundler", "~> 1.12"
+  spec.add_development_dependency "bundler", ">= 2.2.10"
+  #spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rake", ">= 12.3.3"
   spec.add_development_dependency "rspec", "~> 3.0"
 end

data/exe/eph_bpn

@@ -3,22 +3,29 @@
 require "eph_bpn"
 
 e = EphBpn.new(ARGV[0])
-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 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]
-p pos_n                            #=> [-0.5114180757959449, 0.8053955495119587, 0.34914729898526137]
-p pos_bp                           #=> [-0.5114184398598124, 0.805395337986056, 0.3491472536550768]
+#pos_g   = [-0.50787065, 0.80728228, 0.34996714]  # GCRS Coordinate
+pos_g   = [-1.0020195, 0.0660430, 0.0286337]
+pos_b   = e.apply_bias(pos_g)           # Apply Bias
+pos_bp  = e.apply_bias_prec(pos_g)      # Apply Bias & Precession
+pos_bpn = e.apply_bias_prec_nut(pos_g)  # Apply Bias & Precession & Nutation
+pos_p   = e.apply_prec(pos_b)           # Apply Precession
+pos_pn  = e.apply_prec_nut(pos_b)       # Apply Precession & Nutation
+pos_n   = e.apply_nut(pos_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 e.r_bias                              # Rotation matrix for Bias
+p e.r_bias_prec                         # Rotation matrix for Bias & Precession
+p e.r_bias_prec_nut                     # Rotation matrix for Bias & Precession & Nutation
+p e.r_prec                              # Rotation matrix for Precession
+p e.r_prec_nut                          # Rotation matrix for Precession & Nutation
+p e.r_nut                               # Rotation matrix for Nutation
+puts "SRC: #{pos_g}"                    #=> [-0.50787065, 0.80728228, 0.34996714]
+puts "  B: #{pos_b}"                    #=> [-0.507870315369686, 0.8072824634004779, 0.34996720255697145]
+puts " BP: #{pos_bp}"                   #=> [-0.5114184398598124, 0.805395337986056, 0.3491472536550768]
+puts "BPN: #{pos_bpn}"                  #=> [-0.5114038678325663, 0.8054186649898355, 0.3491147860877318]
+puts "  P: #{pos_p}"                    #=> [-0.5114180771311418, 0.8053955471104202, 0.34914730256926324]
+puts " PN: #{pos_pn}"                   #=> [-0.5114035051003876, 0.8054188741108519, 0.34911483499021656]
+puts "  N: #{pos_n}"                    #=> [-0.5114035051003876, 0.8054188741108514, 0.34911483499021656]
 

data/lib/eph_bpn/compute.rb

@@ -116,6 +116,54 @@ module EphBpn
       raise
     end
 
+    #=========================================================================
+    # Bias + Precession 変換行列
+    #
+    # * IAU 2006 (Fukushima-Williams 4-angle formulation) 理論
+    #
+    # @param:  <none>
+    # @return: r  (変換行列)
+    #=========================================================================
+    def comp_r_bias_prec
+      gamma = comp_gamma_bp
+      phi   = comp_phi_bp
+      psi   = comp_psi_bp
+      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
+
+    #=========================================================================
+    # Bias + Precession + Nutation 変換行列
+    #
+    # * IAU 2006 (Fukushima-Williams 4-angle formulation) 理論
+    #
+    # @param:  <none>
+    # @return: r  (変換行列)
+    #=========================================================================
+    def comp_r_bias_prec_nut
+      gamma = comp_gamma_bp
+      phi   = comp_phi_bp
+      psi   = comp_psi_bp
+      fj2 = -2.7774e-6 * @jc
+      dpsi_ls, deps_ls = compute_lunisolar
+      dpsi_pl, deps_pl = compute_planetary
+      dpsi, deps = dpsi_ls + dpsi_pl, deps_ls + deps_pl
+      dpsi += dpsi * (0.4697e-6 + fj2)
+      deps += deps * fj2
+      r = r_z(gamma)
+      r = r_x(phi,          r)
+      r = r_z(-psi - dpsi,  r)
+      r = r_x(-@eps - deps, r)
+      return r
+    rescue => e
+      raise
+    end
+
     #=========================================================================
     # precession（歳差）変換行列（J2000.0 用）
     #
@@ -142,25 +190,9 @@ module EphBpn
     # @return: r  (変換行列)
     #=========================================================================
     def comp_r_prec
-      gamma = ((10.556403    + \
-               (0.4932044    + \
-               (-0.00031238  + \
-               (-0.000002788 + \
-               (0.0000000260)  \
-               * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
-      phi   =  (84381.406000    + \
-               (  -46.811015    + \
-               (    0.0511269   + \
-               (    0.00053289  + \
-               (   -0.000000440 + \
-               (   -0.0000000176) \
-               * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
-      psi   = (( 5038.481507    + \
-               (    1.5584176   + \
-               (   -0.00018522  + \
-               (   -0.000026452 + \
-               (   -0.0000000148) \
-               * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
+      gamma = comp_gamma_p
+      phi   = comp_phi_p
+      psi   = comp_psi_p
       r = r_z(gamma)
       r = r_x(phi,   r)
       r = r_z(-psi,  r)
@@ -171,39 +203,27 @@ module EphBpn
     end
 
     #=========================================================================
-    # Bias + Precession 変換行列
+    # Precession + Nutation 変換行列
     #
-    # * IAU 2006 (Fukushima-Williams 4-angle formulation) 理論
+    # * IAU 2000A nutation with adjustments to match the IAU 2006 precession.
     #
     # @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
+    def comp_r_prec_nut
+      gamma = comp_gamma_p
+      phi   = comp_phi_p
+      psi   = comp_psi_p
+      fj2 = -2.7774e-6 * @jc
+      dpsi_ls, deps_ls = compute_lunisolar
+      dpsi_pl, deps_pl = compute_planetary
+      dpsi, deps = dpsi_ls + dpsi_pl, deps_ls + deps_pl
+      dpsi += dpsi * (0.4697e-6 + fj2)
+      deps += deps * fj2
       r = r_z(gamma)
-      r = r_x(phi,   r)
-      r = r_z(-psi,  r)
-      r = r_x(-@eps, r)
+      r = r_x(phi,          r)
+      r = r_z(-psi - dpsi,  r)
+      r = r_x(-@eps - deps, r)
       return r
     rescue => e
       raise
@@ -224,7 +244,6 @@ module EphBpn
       dpsi, deps = dpsi_ls + dpsi_pl, deps_ls + deps_pl
       dpsi += dpsi * (0.4697e-6 + fj2)
       deps += deps * fj2
-      p dpsi, deps
       r = r_x(@eps)
       r = r_z(-dpsi, r)
       r = r_x(-@eps-deps, r)
@@ -233,6 +252,112 @@ module EphBpn
       raise
     end
 
+    #=========================================================================
+    # 歳差変換行列用 gamma 計算
+    #
+    # @param:  <none>
+    # @return: gamma
+    #=========================================================================
+    def comp_gamma_p
+      return ((10.556403    + \
+              ( 0.4932044   + \
+              (-0.00031238  + \
+              (-0.000002788 + \
+              ( 0.0000000260) \
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # 歳差変換行列用 phi 計算
+    #
+    # @param:  <none>
+    # @return: phi
+    #=========================================================================
+    def comp_phi_p
+      return (84381.406000    + \
+             (  -46.811015    + \
+             (    0.0511269   + \
+             (    0.00053289  + \
+             (   -0.000000440 + \
+             (   -0.0000000176) \
+             * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # 歳差変換行列用 psi 計算
+    #
+    # @param:  <none>
+    # @return: psi
+    #=========================================================================
+    def comp_psi_p
+      return (( 5038.481507    + \
+              (    1.5584176   + \
+              (   -0.00018522  + \
+              (   -0.000026452 + \
+              (   -0.0000000148) \
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # バイアス＆歳差変換行列用 gamma 計算
+    #
+    # @param:  <none>
+    # @return: gamma
+    #=========================================================================
+    def comp_gamma_bp
+      gamma = (-0.052928    + \
+              (10.556378    + \
+              ( 0.4932044   + \
+              (-0.00031238  + \
+              (-0.000002788 + \
+              ( 0.0000000260) \
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # バイアス＆歳差変換行列用 phi 計算
+    #
+    # @param:  <none>
+    # @return: phi
+    #=========================================================================
+    def comp_phi_bp
+      phi   = (84381.412819    + \
+              (  -46.811016    + \
+              (    0.0511268   + \
+              (    0.00053289  + \
+              (   -0.000000440 + \
+              (   -0.0000000176) \
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # バイアス＆歳差変換行列用 psi 計算
+    #
+    # @param:  <none>
+    # @return: psi
+    #=========================================================================
+    def comp_psi_bp
+      psi   = (  -0.041775    + \
+              (5038.481484    + \
+              (   1.5584175   + \
+              (  -0.00018522  + \
+              (  -0.000026452 + \
+              (  -0.0000000148) \
+              * @jc) * @jc) * @jc) * @jc) * @jc) * Const::AS2R
+    rescue => e
+      raise
+    end
+
     #=========================================================================
     # 日月章動(luni-solar nutation)の計算
     #

data/lib/eph_bpn/ephemeris.rb

@@ -1,18 +1,22 @@
 module EphBpn
+
   class Ephemeris
     attr_reader :tdb, :jd, :jc, :eps,
-                :r_bias, :r_prec, :r_nut, :r_bias_prec
+                :r_bias, :r_bias_prec, :r_bias_prec_nut,
+                :r_prec, :r_prec_nut, :r_nut
     include Compute
 
     def initialize(tdb)
-      @tdb = tdb                     # 太陽系力学時
-      @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  # 回転行列（バイアス＆歳差）
+      @tdb = tdb                               # 太陽系力学時
+      @jd  = gc2jd(@tdb)                       # TDB -> JD（ユリウス日）
+      @jc  = jd2jc(@jd)                        # ユリウス世紀
+      @eps = compute_obliquity(@jc)            # 平均黄道傾斜角
+      @r_bias          = comp_r_bias           # 回転行列（バイアス）
+      @r_bias_prec     = comp_r_bias_prec      # 回転行列（バイアス＆歳差）
+      @r_bias_prec_nut = comp_r_bias_prec_nut  # 回転行列（バイアス＆歳差＆章動）
+      @r_prec          = comp_r_prec           # 回転行列（歳差）
+      @r_prec_nut      = comp_r_prec_nut       # 回転行列（歳差＆章動）
+      @r_nut           = comp_r_nut            # 回転行列（章動）
     end
 
     #=========================================================================
@@ -24,6 +28,29 @@ module EphBpn
     def apply_bias(pos)
       return rotate(@r_bias, pos)
     end
+    alias :apply_b :apply_bias
+
+    #=========================================================================
+    # Bias + Precession（歳差） 適用
+    #
+    # @param:  pos    (適用前直角座標)
+    # @return: pos_p  (適用後直角座標)
+    #=========================================================================
+    def apply_bias_prec(pos)
+      return rotate(@r_bias_prec, pos)
+    end
+    alias :apply_bp :apply_bias_prec
+
+    #=========================================================================
+    # Bias + Precession（歳差） + Nutation（章動） 適用
+    #
+    # @param:  pos    (適用前直角座標)
+    # @return: pos_p  (適用後直角座標)
+    #=========================================================================
+    def apply_bias_prec_nut(pos)
+      return rotate(@r_bias_prec_nut, pos)
+    end
+    alias :apply_bpn :apply_bias_prec_nut
 
     #=========================================================================
     # Precession（歳差） 適用
@@ -34,26 +61,29 @@ module EphBpn
     def apply_prec(pos)
       return rotate(@r_prec, pos)
     end
+    alias :apply_p :apply_prec
 
     #=========================================================================
-    # Nutation（章動） 適用
+    # Precession（歳差） + Nutation（章動） 適用
     #
     # @param:  pos    (適用前直角座標)
-    # @return: pos_n  (適用後直角座標)
+    # @return: pos_p  (適用後直角座標)
     #=========================================================================
-    def apply_nut(pos)
-      return rotate(@r_nut, pos)
+    def apply_prec_nut(pos)
+      return rotate(@r_prec_nut, pos)
     end
+    alias :apply_pn :apply_prec_nut
 
     #=========================================================================
-    # Bias + Precession（歳差） 適用
+    # Nutation（章動） 適用
     #
     # @param:  pos    (適用前直角座標)
-    # @return: pos_p  (適用後直角座標)
+    # @return: pos_n  (適用後直角座標)
     #=========================================================================
-    def apply_bias_prec(pos)
-      return rotate(@r_bias_prec, pos)  # IAU 2006 (Fukushima-Williams) 理論
+    def apply_nut(pos)
+      return rotate(@r_nut, pos)
     end
+    alias :apply_n :apply_nut
   end
 end
 

data/lib/eph_bpn/version.rb

@@ -1,3 +1,3 @@
 module EphBpn
-  VERSION = "0.1.7"
+  VERSION = "0.1.12"
 end

metadata

@@ -1,43 +1,43 @@
 --- !ruby/object:Gem::Specification
 name: eph_bpn
 version: !ruby/object:Gem::Version
-  version: 0.1.7
+  version: 0.1.12
 platform: ruby
 authors:
 - komasaru
-autorequire: 
+autorequire:
 bindir: exe
 cert_chain: []
-date: 2016-09-06 00:00:00.000000000 Z
+date: 2021-05-25 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
-        version: '1.12'
+        version: 2.2.10
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
-        version: '1.12'
+        version: 2.2.10
 - !ruby/object:Gem::Dependency
   name: rake
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
-        version: '10.0'
+        version: 12.3.3
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
-        version: '10.0'
+        version: 12.3.3
 - !ruby/object:Gem::Dependency
   name: rspec
   requirement: !ruby/object:Gem::Requirement
@@ -82,7 +82,7 @@ homepage: https://github.com/komasaru/eph_bpn
 licenses:
 - MIT
 metadata: {}
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -97,9 +97,8 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.6.6
-signing_key: 
+rubygems_version: 3.2.9
+signing_key:
 specification_version: 4
 summary: Ephemeris(Bias, Precession, Nutation) calculation tool.
 test_files: []