eph_jpl 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 9b03b2d3373fd9bc0620fbf5872f3a6723b296af
+  data.tar.gz: c9ac3a034cea4f40a52f7da183f7db0a039a05f2
+SHA512:
+  metadata.gz: b8581d2d5933ca5912540f400fbd799eaf7eee22e3dc59cf4f1e1e8a99b543424f4371215202840626bb62f3b099280d2b9872a98b069a71911d433174876368
+  data.tar.gz: 3f8521f73b9a180ae73f307bd00e80649c25c100c00f3e577b43da1cce8ffe709e26c4f61491595771a49645dec5cbd4976f93fb1ccd1a8d596600574b74410d

data/.gitignore

@@ -0,0 +1,9 @@
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/

data/.rspec

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

data/.travis.yml

@@ -0,0 +1,4 @@
+language: ruby
+rvm:
+  - 2.3.1
+before_install: gem install bundler -v 1.11.2

data/Gemfile

@@ -0,0 +1,10 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in eph_jpl.gemspec
+gemspec
+
+group :development do
+  gem "guard"
+  gem "guard-rspec", "~> 4.7.0"
+end
+

data/Guardfile

@@ -0,0 +1,70 @@
+# A sample Guardfile
+# More info at https://github.com/guard/guard#readme
+
+## Uncomment and set this to only include directories you want to watch
+# directories %w(app lib config test spec features) \
+#  .select{|d| Dir.exists?(d) ? d : UI.warning("Directory #{d} does not exist")}
+
+## Note: if you are using the `directories` clause above and you are not
+## watching the project directory ('.'), then you will want to move
+## the Guardfile to a watched dir and symlink it back, e.g.
+#
+#  $ mkdir config
+#  $ mv Guardfile config/
+#  $ ln -s config/Guardfile .
+#
+# and, you'll have to watch "config/Guardfile" instead of "Guardfile"
+
+# Note: The cmd option is now required due to the increasing number of ways
+#       rspec may be run, below are examples of the most common uses.
+#  * bundler: 'bundle exec rspec'
+#  * bundler binstubs: 'bin/rspec'
+#  * spring: 'bin/rspec' (This will use spring if running and you have
+#                          installed the spring binstubs per the docs)
+#  * zeus: 'zeus rspec' (requires the server to be started separately)
+#  * 'just' rspec: 'rspec'
+
+guard :rspec, cmd: "bundle exec rspec" do
+  require "guard/rspec/dsl"
+  dsl = Guard::RSpec::Dsl.new(self)
+
+  # Feel free to open issues for suggestions and improvements
+
+  # RSpec files
+  rspec = dsl.rspec
+  watch(rspec.spec_helper) { rspec.spec_dir }
+  watch(rspec.spec_support) { rspec.spec_dir }
+  watch(rspec.spec_files)
+
+  # Ruby files
+  ruby = dsl.ruby
+  dsl.watch_spec_files_for(ruby.lib_files)
+
+  ## Rails files
+  #rails = dsl.rails(view_extensions: %w(erb haml slim))
+  #dsl.watch_spec_files_for(rails.app_files)
+  #dsl.watch_spec_files_for(rails.views)
+
+  #watch(rails.controllers) do |m|
+  #  [
+  #    rspec.spec.call("routing/#{m[1]}_routing"),
+  #    rspec.spec.call("controllers/#{m[1]}_controller"),
+  #    rspec.spec.call("acceptance/#{m[1]}")
+  #  ]
+  #end
+
+  ## Rails config changes
+  #watch(rails.spec_helper)     { rspec.spec_dir }
+  #watch(rails.routes)          { "#{rspec.spec_dir}/routing" }
+  #watch(rails.app_controller)  { "#{rspec.spec_dir}/controllers" }
+
+  ## Capybara features specs
+  #watch(rails.view_dirs)     { |m| rspec.spec.call("features/#{m[1]}") }
+  #watch(rails.layouts)       { |m| rspec.spec.call("features/#{m[1]}") }
+
+  ## Turnip features and steps
+  #watch(%r{^spec/acceptance/(.+)\.feature$})
+  #watch(%r{^spec/acceptance/steps/(.+)_steps\.rb$}) do |m|
+  #  Dir[File.join("**/#{m[1]}.feature")][0] || "spec/acceptance"
+  #end
+end

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2016 Masaru Koizumi
+
+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,103 @@
+# EphJpl
+
+## Introduction
+
+This is the gem library which calculates ephemeris datas by JPL(NASA Jet Propulsion Laboratory) method.
+
+This library calculates rectangular coordinates and velocities of a target astronomical body which made an astronomical body the center on a Julian Day.(Coordinate system is ICRS)
+
+### What's JPL?
+
+Please refer the following link.
+
+* [NASA Jet Propulsion Laboratory (JPL) - Space Mission and Science News, Videos and Images](http://www.jpl.nasa.gov/)
+
+Please refer the following links about DE datas.
+
+* [HORIZONS System](http://ssd.jpl.nasa.gov/?horizons)
+* [ftp://ssd.jpl.nasa.gov/pub/eph/planets/](ftp://ssd.jpl.nasa.gov/pub/eph/planets/)
+
+To understand this library's in-depth specification, you need to comprehend contents of the following link well.
+
+* [testeph.f - JPL](ftp://ssd.jpl.nasa.gov/pub/eph/planets/fortran/testeph.f)
+
+### Supported DE Ver.
+
+This library supports only DE430, now.
+
+### Settable astronomical bodies as target and center numbers.
+
+1. Mercury
+2. Venus
+3. Earth
+4. Mars
+5. Jupiter
+6. Saturn
+7. Uranus
+8. Neptune
+9. Pluto
+10. Moon
+11. Sun
+12. Solar system Barycenter
+13. Earth-Moon barycenter
+14. Earth Nutations
+15. Lunar mantle Librations
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'eph_jpl'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install eph_jpl
+
+## Usage
+
+``` ruby
+obj = EphJpl.new("/path/to/JPLEPH", 11, 3, 2457465.5)
+
+p obj.target       #=> Target atronomical body No. (Integer)
+p obj.target_name  #=> Target atronomical body name (String)
+p obj.center       #=> Center atronomical body No. (Integer)
+p obj.center_name  #=> Center atronomical body name (String)
+p obj.jd           #=> Julian Day (Float)
+p obj.km           #=> KM flag (true: km unit, false: AU unit) (Boolean)
+p obj.unit         #=> UNIT of positions and velocities (String)
+p obj.bin          #=> Acquired data from binary file (Hash)
+p obj.calc         #=> [x, y, z-position, x, y, z-velocity]
+```
+
+About arguments for target and center astronomical bodies.
+
+* You can set a integer betweetn 1 and 15 as a target astronomical body.
+* You can set a integer betweetn 0 and 13 as a center astronomical body.
+
+About units.
+
+* If a target astronomical body number < 14 and KM flag = false, units of position and velocity are au, au/day.
+* If a target astronomical body number < 14 and KM flag = true, units of position and velocity are km, km/sec.
+* If a target astronomical body number = 14 or 15, units of position and velocity are rad, rad/day.
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment. Run `bundle exec eph_jpl` to use the gem in this directory, ignoring other installed copies of this gem.
+
+To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release`, which will create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/komasaru/eph_jpl.
+
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](http://opensource.org/licenses/MIT).
+

data/Rakefile

@@ -0,0 +1,6 @@
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+RSpec::Core::RakeTask.new(:spec)
+
+task :default => :spec

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "eph_jpl"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start

data/bin/setup

@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/eph_jpl.gemspec

@@ -0,0 +1,33 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'eph_jpl/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "eph_jpl"
+  spec.version       = EphJpl::VERSION
+  spec.authors       = ["komasaru"]
+  spec.email         = ["masaru@mk-mode.com"]
+
+  spec.summary       = %q{Ephemeris calculation tool by JPL method.}
+  spec.description   = %q{EphJcg is a ephemeris calculation tool by JPL method.}
+  spec.homepage      = "https://github.com/komasaru/eph_jpl"
+  spec.license       = "MIT"
+
+  # Prevent pushing this gem to RubyGems.org by setting 'allowed_push_host', or
+  # delete this section to allow pushing this gem to any host.
+  #if spec.respond_to?(:metadata)
+  #  spec.metadata['allowed_push_host'] = "TODO: Set to 'http://mygemserver.com'"
+  #else
+  #  raise "RubyGems 2.0 or newer is required to protect against public gem pushes."
+  #end
+
+  spec.files         = `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", "~> 1.11"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rspec", "~> 3.0"
+end

data/exe/eph_jpl

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "eph_jpl"
+
+# e.g. `bundle exec exe/eph_jpl ~/src/ephemeris_jpl/JPLEPH 11 3 2457465.5`  # 2016-03-18
+o = EphJpl.new(*ARGV)
+str =  "TARGET = #{o.target} (#{o.target_name})\n"
+str << "CENTER = #{o.center} (#{o.center_name})\n"
+str << "    JD = #{o.jd}\n"
+str << "  UNIT = #{o.unit}"
+puts str
+#p o.bin
+p o.calc
+

data/lib/eph_jpl/argument.rb

@@ -0,0 +1,87 @@
+module EphJpl
+  class Argument
+    def initialize(*args)
+      @args = *args
+    end
+
+    #=========================================================================
+    # 引数取得
+    #
+    # @return: [BIN_PATH, TARGET, CENTER, JD, KM]
+    #=========================================================================
+    def get_args
+      bin_path = get_binpath
+      target   = get_target
+      center   = get_center
+      jd       = get_jd
+      km       = get_km
+      check_bin_path(bin_path)
+      check_target_center(target, center)
+      return [bin_path, target, center, jd, km]
+    rescue => e
+      raise
+    end
+
+    def get_binpath
+      raise unless bin_path = @args.shift
+      return bin_path
+    rescue => e
+      raise Const::MSG_ERR_1
+    end
+
+    def get_target
+      raise unless target = @args.shift
+      raise unless target.to_s =~ /^\d+$/
+      raise if target.to_i < 1 || 15 < target.to_i
+      return target.to_i
+    rescue => e
+      raise Const::MSG_ERR_3
+    end
+
+    def get_center
+      raise unless center = @args.shift
+      raise unless center.to_s =~ /^\d+$/
+      raise if center.to_i < 0 || 13 < center.to_i
+      return center.to_i
+    rescue => e
+      raise Const::MSG_ERR_4
+    end
+
+    def get_jd
+      raise unless jd = @args.shift
+      if jd.to_s !~ /^[\d\.]+$/ || \
+         jd.to_f < Const::EPOCH_PERIOD[0] || \
+         Const::EPOCH_PERIOD[1] < jd.to_f
+        raise
+      end
+      return jd.to_f
+    rescue => e
+      raise Const::MSG_ERR_7
+    end
+
+    def get_km
+      km = @args.shift
+      km ||= Const::KM
+      raise unless km.to_s =~ /^true|false|[01]$/
+      km = km.to_s =~ /0|false/ ? false : true
+      return km
+    rescue => e
+      raise Const::MSG_ERR_8
+    end
+
+    def check_target_center(target, center)
+      case
+      when target == center
+        raise Const::MSG_ERR_5
+      when target < 14 && center == 0,
+           target > 13 && center != 0
+        raise Const::MSG_ERR_6
+      end
+    end
+
+    def check_bin_path(bin_path)
+      raise Const::MSG_ERR_2 unless File.exist?(bin_path)
+    end
+  end
+end
+

data/lib/eph_jpl/binary.rb

@@ -0,0 +1,270 @@
+module EphJpl
+  class Binary
+    def initialize(*args)
+      @bin_path, @target, @center, @jd = *args
+      @pos = 0
+    end
+
+    def get_binary
+      begin
+        ttl    = get_ttl          # TTL
+        cnams  = get_cnams        # CNAM
+        sss    = get_sss          # SS
+        ncon   = get_ncon         # NCON
+        au     = get_au           # AU
+        emrat  = get_emrat        # EMRAT
+        ipts   = get_ipts         # IPT
+        numde  = get_numde        # NUMDE
+        ipts  << get_ipts_13      # IPT(Month's libration)
+        cvals  = get_cvals(ncon)  # CVAL（定数値）
+        jdepoc = cvals[4]         # JDEPOC
+        coeffs, jds_cheb = get_coeffs(sss, ipts)  # Coefficient, JDs(for Chebyshev polynomial)
+        return {
+          ttl: ttl, cnams: cnams, sss: sss, ncon: ncon, au: au, emrat: emrat,
+          numde: numde, ipts: ipts, cvals: cvals, jdepoc: jdepoc,
+          coeffs: coeffs, jds_cheb: jds_cheb
+        }
+      rescue => e
+        raise
+      end
+    end
+
+    private
+
+    #=========================================================================
+    # TTL
+    #
+    # @param:  <none>
+    # @return: TTL
+    #=========================================================================
+    def get_ttl
+      recl = 84
+
+      begin
+        ttl = (0..2).map do |i|
+          File.binread(@bin_path, recl, @pos + recl * i).unpack("A*")[0]
+        end.join("\n")
+        @pos += recl * 3
+        return ttl
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # CNAM
+    #
+    # @param:  <none>
+    # @return: Array of CNAM
+    #=========================================================================
+    def get_cnams
+      recl = 6
+
+      begin
+        cnams = (0..399).map do |i|
+          File.binread(@bin_path, recl, @pos + recl * i).unpack("A*")[0]
+        end
+        @pos += recl * 400
+        return cnams
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # SS
+    #
+    # @param:  <none>
+    # @return: Array of SS
+    #=========================================================================
+    def get_sss
+      recl = 8
+
+      begin
+        sss = (0..2).map do |i|
+          File.binread(@bin_path, recl, @pos + recl * i).unpack("d*")[0]
+        end
+        @pos += recl * 3
+        return sss
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # NCON
+    #
+    # @param:  <none>
+    # @return: NCON
+    #=========================================================================
+    def get_ncon
+      recl = 4
+
+      begin
+        ncon = File.binread(@bin_path, recl, @pos).unpack("I*")[0]
+        @pos += recl
+        return ncon
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # AU
+    #
+    # @param:  <none>
+    # @return: AU
+    #=========================================================================
+    def get_au
+      recl = 8
+
+      begin
+        au = File.binread(@bin_path, recl, @pos).unpack("d*")[0]
+        @pos += recl
+        return au
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # EMRAT
+    #
+    # @param:  <none>
+    # @return: <none>
+    #=========================================================================
+    def get_emrat
+      recl = 8
+
+      begin
+        emrat = File.binread(@bin_path, recl, @pos).unpack("d*")[0]
+        @pos += recl
+        return emrat
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # IPT
+    #
+    # @param:  <none>
+    # @return: Array of IPT
+    #=========================================================================
+    def get_ipts
+      recl = 4
+
+      begin
+        ipts = (0..11).map do |i|
+          ary = (0..2).map do |j|
+            File.binread(@bin_path, recl, @pos + recl * j).unpack("I*")[0]
+          end
+          @pos += recl * 3
+          ary
+        end
+        return ipts
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # NUMDE
+    #
+    # * If NUMDE != 430, raise error!
+    #
+    # @param:  <none>
+    # @return: NUMDE
+    #=========================================================================
+    def get_numde
+      recl = 4
+
+      begin
+        numde = File.binread(@bin_path, recl, @pos).unpack("I*")[0]
+        raise Const::MSG_ERR_8 unless numde == 430
+        @pos += recl
+        return numde
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # IPT_13(Month's libration)
+    #
+    # @param:  <none>
+    # @return: Array of IPT
+    #=========================================================================
+    def get_ipts_13
+      recl = 4
+
+      begin
+        ipts = (0..2).map do |i|
+          File.binread(@bin_path, recl, @pos + recl * i).unpack("I*")[0]
+        end
+        @pos += recl * 3
+        return ipts
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # CVAL
+    #
+    # @param:  NCON
+    # @return: Array of CVAL
+    #=========================================================================
+    def get_cvals(ncon)
+      pos = Const::KSIZE * Const::RECL
+      recl = 8
+
+      begin
+        return (0..ncon - 1).map do |i|
+          File.binread(@bin_path, recl, pos + recl * i).unpack("d*")[0]
+        end
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # COEFF
+    #
+    # * Set JD(start, end) for Chebyshev polynomial to the array @jd_cheb
+    #
+    # @param: Array of SS
+    # @param: Array of IPT
+    # @return: <none>
+    #=========================================================================
+    def get_coeffs(sss, ipts)
+      idx = ((@jd - sss[0]) / sss[2]).floor  # レコードインデックス
+      pos = Const::KSIZE * Const::RECL * (2 + idx)
+      recl = 8
+      coeffs = Array.new
+
+      begin
+        items = (0..(Const::KSIZE / 2) - 1).map do |i|
+          File.binread(@bin_path, recl, pos + recl * i).unpack("d*")[0]
+        end
+        jd_cheb = [items.shift, items.shift]
+        ipts.each_with_index do |ipt, i|
+          n = i == 11 ? 2 : 3  # 要素数
+          ary_1 = Array.new
+          ipt[2].times do |j|
+            ary_0 = Array.new
+            n.times do |k|
+              ary_0 << items.shift(ipt[1])
+            end
+            ary_1 << ary_0
+          end
+          coeffs << ary_1
+        end
+        return [coeffs, jd_cheb]
+      rescue => e
+        raise
+      end
+    end
+  end
+end
+

data/lib/eph_jpl/const.rb

@@ -0,0 +1,36 @@
+module EphJpl
+  module Const
+    USAGE        = <<-EOS
+[USAGE] EphJpl.new(BIN_PATH, TARGET, CENTER, JD)
+  [ASTRO NO] (TARGET: 1 - 15, CENTER: 0 - 13)
+  1: Mercury, 2: Venus, 3: Earth, 4: Mars, 5: Jupiter,
+  6: Saturn, 7: Uranus, 8: Neptune, 9: Pluto, 10: Moon,
+  11: Sun, 12: Solar system Barycenter, 13: Earth-Moon barycenter,
+  14: Earth Nutations, 15: Lunar mantle Librations
+  * If TARGET = 14 or 15, CENTER = 0
+  * TARGET != CENTER
+  * 2287184.5 <= JD <= 2688976.5
+EOS
+    MSG_ERR_1    = "Binary file path is invalid!"
+    MSG_ERR_2    = "Binary file is not found!"
+    MSG_ERR_3    = "TARGET is invalid!"
+    MSG_ERR_4    = "CENTER is invalid!"
+    MSG_ERR_5    = "TARGET == CENTER ?"
+    MSG_ERR_6    = "TARGET or CENTER is invalid!"
+    MSG_ERR_7    = "JD is invalid!"
+    MSG_ERR_8    = "KM flag is invalid!"
+    MSG_ERR_9    = "This library is supporting only DE430!"
+    EPOCH_PERIOD = [2287184.5, 2688976.5]
+    KSIZE        = 2036
+    RECL         = 4
+    ASTRS        = [
+      "Mercury", "Venus", "Earth", "Mars", "Jupiter", "Saturn", "Uranus",
+      "Neptune", "Pluto", "Moon", "Sun", "Solar system Barycenter",
+      "Earth-Moon barycenter", "Earth Nutations", "Lunar mantle Librations"
+    ]
+    KIND         = 2
+    BARY         = true
+    KM           = false
+  end
+end
+

data/lib/eph_jpl/ephemeris.rb

@@ -0,0 +1,195 @@
+module EphJpl
+  class Ephemeris
+    attr_reader :target, :center, :jd, :bin, :km,
+                :target_name, :center_name, :unit
+
+    def initialize(target, center, jd, bin, km = false)
+      @target, @center, @jd, @km = target, center, jd, km
+      @target_name = Const::ASTRS[@target - 1]
+      @center_name = Const::ASTRS[@center - 1]
+      @unit = @target > 13 ? "rad, rad/day" : @km ? "km, km/sec" : "AU, AU/day"
+      @bin = bin
+      @list = get_list
+    end
+
+    def calc
+      pvs     = Array.new(11).map { |a| Array.new(6, 0.0) }  # for position, velocity
+      pvs_tmp = Array.new(13)      # Temporary array for position, velocity of target and center
+      rrds    = Array.new(6, 0.0)  # for calculated data (difference between target and center)
+
+      begin
+        # Interpolate (11:Sun)
+        pv_sun = interpolate(11)
+        # Interpolate (1:Mercury - 10:Moon)
+        0.upto(9) do |i|
+          next if @list[i] == 0
+          pvs[i] = interpolate(i + 1)
+          next if i > 8
+          next if Const::BARY
+          pvs[i] = pvs[i].map.with_index do |pv, j|
+            pv - pv_sun[j]
+          end
+        end
+        # Interpolate (14:Nutation)
+        if @list[10] > 0 && @bin[:ipts][11][1] > 0
+          p_nut = interpolate(14)
+        end
+        # Interpolate (15:Libration）
+        if @list[11] > 0 && @bin[:ipts][12][1] > 0
+          pvs[10] = interpolate(15)
+        end
+
+        # Difference between target and center
+        case
+        when @target == 14
+          rrds = p_nut if @bin[:ipts][11][1] > 0
+        when @target == 15
+          rrds = pvs[10] if @bin[:ipts][12][1] > 0
+        else
+          0.upto(9) { |i| pvs_tmp[i] = pvs[i] }
+          pvs_tmp[10] = pv_sun            if [@target, @center].include?(11)
+          pvs_tmp[11] = Array.new(6, 0.0) if [@target, @center].include?(12)
+          pvs_tmp[12] = pvs[2]            if [@target, @center].include?(13)
+          if @target * @center == 30 || @target + @center == 13
+            pvs_tmp[2] = Array.new(6, 0.0)
+          else
+            pvs_tmp[2] = pvs[2].map.with_index do |pv, i|
+              pv - pvs[9][i] / (1.0 + @bin[:emrat])
+            end unless @list[2] == 0
+            pvs_tmp[9] = pvs_tmp[2].map.with_index do |pv, i|
+              pv + pvs[9][i]
+            end unless @list[9] == 0
+          end
+          0.upto(5) do |i|
+            rrds[i] = pvs_tmp[@target - 1][i] - pvs_tmp[@center - 1][i]
+          end
+        end
+        return rrds
+      rescue => e
+        raise
+      end
+    end
+
+    private
+
+    #=========================================================================
+    # Computation target list
+    #
+    # @param:  <none>
+    # @return: Array
+    #=========================================================================
+    def get_list
+      list = Array.new(12, 0)
+
+      begin
+        if @target == 14
+          list[10] = Const::KIND if @bin[:ipts][11][1] > 0
+          return list
+        end
+        if @target == 15
+          list[11] = Const::KIND if @bin[:ipts][12][1] > 0
+          return list
+        end
+        [@target, @center].each do |k|
+          list[k - 1] = Const::KIND if k <= 10
+          list[2]     = Const::KIND if k == 10
+          list[9]     = Const::KIND if k ==  3
+          list[2]     = Const::KIND if k == 13
+        end
+        return list
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # Interpolate by Chebyshev's Polynomial
+    #
+    # * Case astro-no
+    #     1 ... 13: Position and velocity of x, y, z (6 items)
+    #           14: Angular position and volocity of delta Psi, delta Epsilon (4 items)
+    #           15: Angular position and volocity of Phi, Theta, Psi (6 items)
+    # * If astro-no = 12, then location and volocity of x, y, z are all 0.0.
+    #
+    # @param:  astr (= astronomical no)
+    # @return: pvs = [
+    #   x-position, y-position, z-position,
+    #   x-velocity, y-velocity, z-velocity
+    # ]
+    #   Case: 14:Nutation
+    #     pvs = [
+    #       delta-psi-angular-position, delta-epsilon-angular-position,
+    #       delta-psi-angular-velocity, delta-epsilon-angular-velocity
+    #     ]
+    #   Case: 15:Libration
+    #     pvs = [
+    #       phi-angular-position,  theta-angular-position, psi-angular-position,
+    #       phi-angular-velocity,  theta-angular-velocity, psi-angular-velocity
+    #     ]
+    #=========================================================================
+    def interpolate(astr)
+      pvs = Array.new
+
+      begin
+        tc, idx_sub = norm_time(astr)
+        n_item = astr == 14 ? 2 : 3  # 要素数
+        i_ipt  = astr > 13 ? astr - 3 : astr - 1
+        i_coef = astr > 13 ? astr - 3 : astr - 1
+
+        # 位置
+        ary_p = [1, tc]
+        2.upto(@bin[:ipts][i_ipt][1] - 1) do |i|
+          ary_p << 2 * tc * ary_p[-1] - ary_p[-2]
+        end  # 各項
+        0.upto(n_item - 1) do |i|
+          val = 0
+          0.upto(@bin[:ipts][i_ipt][1] - 1) do |j|
+            val += @bin[:coeffs][i_coef][idx_sub][i][j] * ary_p[j]
+          end
+          val /= @bin[:au] if !@km && astr < 14
+          pvs << val
+        end  # 値
+
+        # 速度
+        ary_v = [0, 1, 2 * 2 * tc]
+        3.upto(@bin[:ipts][i_ipt][1] - 1) do |i|
+          ary_v << 2 * tc * ary_v[-1] + 2 * ary_p[i - 1] - ary_v[-2]
+        end  # 各項
+        0.upto(n_item - 1) do |i|
+          val = 0
+          0.upto(@bin[:ipts][i_ipt][1] - 1) do |j|
+            val += @bin[:coeffs][i_coef][idx_sub][i][j] * ary_v[j] * 2 * @bin[:ipts][i_ipt][2] / @bin[:sss][2].to_f
+          end
+          if astr < 14
+            val /= @bin[:au] unless @km
+            val /= 86400.0 if @km
+          end
+          pvs << val
+        end  # 値
+
+        return pvs
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # Time normalization, sub-period's index calculation for Chebyshev's Polynomial
+    #
+    # @param:  astr (= atronomical no)
+    # @return: [chebyshev-time, sub-index]
+    #=========================================================================
+    def norm_time(astr)
+      idx = astr > 13 ? astr - 2 : astr
+      jd_start = @bin[:jds_cheb][0]
+      tc = (@jd - jd_start) / @bin[:sss][2].to_f
+      temp = tc * @bin[:ipts][idx - 1][2]
+      idx = (temp - tc.floor).floor         # サブ区間のインデックス
+      tc = (temp % 1.0 + tc.floor) * 2 - 1  # チェビシェフ時間
+      return [tc, idx]
+    rescue => e
+      raise
+    end
+  end
+end
+

data/lib/eph_jpl/version.rb

@@ -0,0 +1,3 @@
+module EphJpl
+  VERSION = "0.1.0"
+end

data/lib/eph_jpl.rb

@@ -0,0 +1,14 @@
+require "eph_jpl/version"
+require "eph_jpl/argument"
+require "eph_jpl/binary"
+require "eph_jpl/const"
+require "eph_jpl/ephemeris"
+
+module EphJpl
+  def self.new(*args)
+    bin_path, target, center, jd, km = Argument.new(*args).get_args
+    bin = Binary.new(bin_path, target, center, jd).get_binary
+    return EphJpl::Ephemeris.new(target, center, jd, bin, km)
+  end
+end
+

metadata

@@ -0,0 +1,105 @@
+--- !ruby/object:Gem::Specification
+name: eph_jpl
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- komasaru
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2016-06-18 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.11'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.11'
+- !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: '3.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+description: EphJcg is a ephemeris calculation tool by JPL method.
+email:
+- masaru@mk-mode.com
+executables:
+- eph_jpl
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- Guardfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- eph_jpl.gemspec
+- exe/eph_jpl
+- lib/eph_jpl.rb
+- lib/eph_jpl/argument.rb
+- lib/eph_jpl/binary.rb
+- lib/eph_jpl/const.rb
+- lib/eph_jpl/ephemeris.rb
+- lib/eph_jpl/version.rb
+homepage: https://github.com/komasaru/eph_jpl
+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.6.4
+signing_key: 
+specification_version: 4
+summary: Ephemeris calculation tool by JPL method.
+test_files: []