calc_sun 0.1.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 29357270031ecabae17fb6934c30a86b3aa49e6b
+  data.tar.gz: 9b95397a0a518191709c2ff26019e756c95eafdf
+SHA512:
+  metadata.gz: c0e3bf0a00afd5b82de2ee91a342109ddaf64634907ac27d078be1b758a9e37d1f818fe04cb79adaff9a30b7ed201eddacfe3d96b1fe6be6d3e382a3fe9256f8
+  data.tar.gz: f337cbfd61716c4b273a1143c08598a3be4b5f3f58e3dd87f7c2a0ba513a067e5b8674081a8c595b13b582cf0388c41720b1dc59a2fe304340b2d1c7f3d5b422

data/Gemfile

@@ -0,0 +1,30 @@
+source 'https://rubygems.org'
+
+gemspec
+
+group :development, :test do
+  gem 'hoe'
+  gem 'rake-compiler', '~> 0.9.3'
+  gem 'mustache'
+  gem 'rubocop'
+  gem 'development'
+  gem 'thor'
+  gem 'guard'
+  gem 'guard-rspec', require: false
+  gem 'minitest'
+  gem 'guard-minitest'
+  gem 'ZenTest'
+  gem 'rspec'
+  gem 'yard'
+  gem 'rdoc'
+  gem 'ffi'
+  gem 'simplecov'
+  gem 'test-unit'
+  gem 'shoulda'
+  gem 'pry'
+  gem 'guard-pytest'
+  gem 'nenv'
+  gem 'lumberjack'
+  gem 'listen'
+  gem 'ruby_dep', '~> 1.2.0'
+end

data/History.rdoc

@@ -0,0 +1,6 @@
+=== 0.1.1 / 2016-12-30
+
+* 1 major enhancement
+
+  * Birthday!
+

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2016 Douglas Allen
+
+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/Manifest.txt

@@ -0,0 +1,13 @@
+History.rdoc
+LICENSE.txt
+Manifest.txt
+Rakefile
+README.rdoc
+calc_sun.gemspec
+ext/calc_sun/calc_sun.c
+ext/calc_sun/extconf.rb
+lib/calc_sun.rb
+lib/calc_sun/version.rb
+test/calc_sun_test.rb
+example/sunriset.rb
+Gemfile

data/README.rdoc

@@ -0,0 +1,89 @@
+
+= calc_sun
+
+Home    :: https://rubygems.org/gems/calc_sun
+Code    :: https://github.com/DouglasAllen/calc_sun
+RDoc    :: http://http://www.rubydoc.info/gems/calc_sun/
+Issues  :: https://github.com/DouglasAllen/calc_sun/issues
+
+
+== DESCRIPTION:
+
+Calculate Sunrise and Sunset. Now uses native C extension.
+An example is supplied to show gem usage.
+
+== FEATURES/PROBLEMS:
+
+  n/a
+
+== SYNOPSIS:
+
+=== start irb
+
+  $ irb --simple-prompt
+
+=== enter Code
+
+        require 'calc_sun'
+        cs = CalcSun.new
+        require 'date'
+        # not included yet for determining your mean time
+        DJ00 = DateTime.parse('2000-01-01T12:00:00').jd
+
+        lat = 51.4770228
+        lon = -0.0001147
+        day = Date.parse('2016-12-25')
+        jd = day.jd - DJ00 - lon / 360.0
+        cs.t_rise(jd, lon, lat)
+        cs.t_mid_day(jd, lon, lat)
+        cs.t_set(jd, lon, lat)require 'date'
+
+
+== INSTALL:
+
+  $> bundle install
+
+  $> bundle exec rake install
+
+
+=== Build with DEVKIT
+
+   cd to devkit dir and run devkitvars.bat
+
+   or make a batchfile to start your envirnment with 'run devkitvars'
+   cd to calc_sun dir
+   Rake -T
+
+== DEVELOPERS:
+
+After checking out the source, run:
+
+  $ bundle exec rake newb
+
+This task will install any missing dependencies, run the tests/specs,
+and generate the RDoc.
+
+== LICENSE:
+
+(The MIT License)
+
+Copyright (c) 2016
+
+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/Rakefile

@@ -0,0 +1,33 @@
+#!/usr/bin/env ruby
+
+require 'hoe'
+require 'rake/testtask'
+require 'rake/extensiontask'
+
+
+Hoe.plugin :manifest
+
+Hoe.spec 'calc_sun' do
+  self.version = '0.1.1'
+  developer('kb9agt', 'kb9agt@gmail.com')
+  self.readme_file   = 'README.rdoc'
+  self.history_file  = 'HISTORY.rdoc'
+
+  extra_dev_deps << ['rake-compiler', '~> 0.9.3']
+  self.spec_extras = { extensions: ['ext/calc_sun/extconf.rb'] }
+
+  Rake::ExtensionTask.new('calc_sun', spec) do |ext|
+    ext.lib_dir = File.join('lib', 'calc_sun')
+  end
+end
+
+Rake::Task[:test].prerequisites << :compile
+# vim: syntax=ruby
+
+Rake::TestTask.new(:test) do |t|
+  t.libs << 'test'
+  t.libs << 'lib'
+  t.test_files = FileList['test/*_test.rb']
+end
+
+task default: :test

data/calc_sun.gemspec

@@ -0,0 +1,33 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'calc_sun/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = 'calc_sun'
+  spec.version       = CalcSun::VERSION
+  spec.author        = 'Douglas Allen'
+  spec.email         = ['kb9agt@gmail.com']
+
+  spec.summary       = 'Calculates Sun Times ea. rise, noon, set'
+  spec.description   = 'supply the jd date, latitude, and longitude in the calls
+                        to each method t_rise(jd, lon, lat),
+                        t_mid_day(jd, lon, lat), and t_set(jd, lon, lat)'
+  spec.homepage      = 'https://github.com/DouglasAllen/calc_sun'
+  spec.license       = 'MIT'
+  require 'rake'
+  spec.files = FileList['lib     .rb',
+                        '[A-Z]*',
+                        'test/   *'].to_a
+  spec.extensions << 'ext/eot/extconf.rb'
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ['lib']
+  spec.platform = Gem::Platform::RUBY
+  spec.add_development_dependency 'bundler', '~> 1.13'
+  spec.add_development_dependency 'rake', '~> 10.0'
+  spec.add_development_dependency 'rake-compiler', '~> 0.9.3'
+  spec.add_development_dependency 'hoe', '~> 3.12.0'
+  spec.add_development_dependency 'guard', '~> 2.6.1'
+  spec.add_development_dependency 'guard-minitest', '~> 2.3.2'
+  spec.add_development_dependency 'minitest', '~> 5.4.1'
+end

data/example/sunriset.rb

@@ -0,0 +1,35 @@
+require 'benchmark'
+lib = File.expand_path('../../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'calc_sun'
+cs = CalcSun.new
+
+require 'date'
+DJ00 = DateTime.parse('2000-01-01T12:00:00').jd
+
+lat = 41.95
+lon = -88.75
+day = Date.parse('2016-12-25')
+jd = day.jd - DJ00 - lon / 360.0
+rise = cs.t_rise(jd, lon, lat)
+midday = cs.t_mid_day(jd, lon, lat)
+set = cs.t_set(jd, lon, lat)
+
+printf("\n")
+
+printf("\tSun rises \t\t\t : %2.0f:%02.0f UTC\n",
+       rise.floor, (rise % 1 * 60.0).floor)
+
+printf("\tSun midday \t\t\t : %2.0f:%02.0f UTC\n",
+       midday.floor, ((midday % 1.0) * 60).floor)
+
+printf("\tSun sets \t\t\t : %2.0f:%02.0f UTC\n",
+       set.floor, (set % 1 * 60.0).floor)
+
+n = 1_000_000
+puts 'running three method calls 1_000_000 times'
+Benchmark.bm(7) do |x|
+  x.report('rise:') { n.times { ; cs.t_rise(jd, lon, lat); } }
+  x.report('midday:') { n.times { ; cs.t_mid_day(jd, lon, lat); } }
+  x.report('set:') { n.times { ; cs.t_set(jd, lon, lat); } }
+end

data/ext/calc_sun/calc_sun.c

@@ -0,0 +1,365 @@
+#include <ruby.h>
+
+#ifndef DBL2NUM
+# define DBL2NUM(dbl) rb_float_new(dbl)
+#endif
+
+#define R2D 57.295779513082320876798154814105
+#define D2R 0.017453292519943295769236907684886
+#define M2PI M_PI * 2.0
+#define INV24 1.0 / 24.0
+#define INV360 1.0 / 360.0
+#define DJ00 2451545.0
+
+// static ID id_status;
+
+static VALUE
+t_init(VALUE self)
+{
+  return self;
+}
+
+static VALUE
+func_rev12(VALUE self, VALUE vx)
+{
+  double x = NUM2DBL(vx);
+  return DBL2NUM(x - 24.0 * floor(x * INV24 + 0.5));
+}
+
+static VALUE
+func_mean_anomally(VALUE self, VALUE vd)
+{
+  double d = NUM2DBL(vd);
+  double vma =
+  fmod(
+    (357.5291  +
+     0.98560028 * d
+    ) * D2R, M2PI);
+  return DBL2NUM(vma);
+}
+
+static VALUE
+func_eccentricity(VALUE self, VALUE vd)
+{
+  double d = NUM2DBL(vd);
+  double ve =
+  0.016709 -
+  1.151e-9 * d;
+  return DBL2NUM(ve);
+}
+
+static VALUE
+func_equation_of_center(VALUE self, VALUE vd)
+{
+  double vma =
+  NUM2DBL(func_mean_anomally(self, vd));
+  double ve =
+  NUM2DBL(func_eccentricity(self, vd));
+  double ve2 = ve * 2.0;
+  double vesqr = ve * ve;
+  double vesqr54 = 5.0 / 4.0 * vesqr;
+  double vecube12 = (vesqr * ve) / 12.0;
+  double veoc =
+  ve2 * sin(vma) +
+  vesqr54 * sin(2 * vma) +
+  vecube12 * (13.0 * sin(3 * vma) - 3.0 * sin(vma));
+  return DBL2NUM(veoc);
+}
+
+static VALUE
+func_true_anomally(VALUE self, VALUE vd)
+{
+  double vma =
+  NUM2DBL(func_mean_anomally(self, vd));
+  double veoc =
+  NUM2DBL(func_equation_of_center(self, vd));
+  double vta = vma + veoc;
+  return DBL2NUM(vta);
+}
+
+static VALUE
+func_mean_longitude(VALUE self, VALUE vd)
+{
+  double d = NUM2DBL(vd);
+  double vml =
+  fmod(
+    (280.4664567 +
+     0.9856473601037645 * d
+    ) * D2R, M2PI);
+  return DBL2NUM(vml);
+}
+
+static VALUE
+func_eccentric_anomally(VALUE self, VALUE vd)
+{
+  double ve =
+  NUM2DBL(func_eccentricity(self, vd));
+  double vml =
+  NUM2DBL(func_mean_longitude(self, vd));
+  double vea =
+  vml + ve * sin(vml) * (1.0 + ve * cos(vml));
+  return DBL2NUM(vea);
+}
+
+static VALUE
+func_obliquity_of_ecliptic(VALUE self, VALUE vd)
+{
+  double d = NUM2DBL(vd);
+  double vooe =
+  (23.439291 - 3.563E-7 * d) * D2R;
+  return DBL2NUM(vooe);
+}
+
+static VALUE
+func_longitude_of_perihelion(VALUE self, VALUE vd)
+{
+  double d = NUM2DBL(vd);
+  double vlop =
+  fmod(
+    (282.9404 +
+     4.70935e-05 * d
+    ) * D2R, M2PI);
+  return DBL2NUM(vlop);
+}
+
+static VALUE
+func_xv(VALUE self, VALUE vd)
+{
+  double vea =
+  NUM2DBL(func_eccentric_anomally(self, vd));
+  double ve =
+  NUM2DBL(func_eccentricity(self, vd));
+  double vxv = cos(vea) - ve;
+  return DBL2NUM(vxv);
+}
+
+static VALUE
+func_yv(VALUE self, VALUE vd)
+{
+  double vea =
+  NUM2DBL(func_eccentric_anomally(self, vd));
+  double ve =
+  NUM2DBL(func_eccentricity(self, vd));
+  double vyv =
+  sqrt(1.0 - ve * ve) * sin(vea);
+  return DBL2NUM(vyv);
+}
+
+static VALUE
+func_true_longitude(VALUE self, VALUE vd)
+{
+  double vta =
+  NUM2DBL(func_true_anomally(self, vd));
+  double vlop =
+  NUM2DBL(func_longitude_of_perihelion(self, vd));
+  double vtl =
+  fmod(vta + vlop, M2PI);
+  return DBL2NUM(vtl);
+}
+
+static VALUE
+func_rv(VALUE self, VALUE vd)
+{
+  double vxv =
+  NUM2DBL(func_xv(self, vd));
+  double vyv =
+  NUM2DBL(func_yv(self, vd));
+  double vrv =
+  sqrt(vxv * vxv + vyv * vyv);
+  return DBL2NUM(vrv);
+}
+
+static VALUE
+func_ecliptic_x(VALUE self, VALUE vd)
+{
+  double vrv =
+  NUM2DBL(func_rv(self, vd));
+  double vtl =
+  NUM2DBL(func_true_longitude(self, vd));
+  double vex = vrv * cos(vtl);
+  return DBL2NUM(vex);
+}
+
+static VALUE
+func_ecliptic_y(VALUE self, VALUE vd)
+{
+  double vrv =
+  NUM2DBL(func_rv(self, vd));
+  double vtl =
+  NUM2DBL(func_true_longitude(self, vd));
+  double vey = vrv * sin(vtl);
+  return DBL2NUM(vey);
+}
+
+static VALUE
+func_right_ascension(VALUE self, VALUE vd)
+{
+  double vey =
+  NUM2DBL(func_ecliptic_y(self, vd));
+  double vooe =
+  NUM2DBL(func_obliquity_of_ecliptic(self, vd));
+  double vex =
+  NUM2DBL(func_ecliptic_x(self, vd));
+  double vra =
+  fmod(atan2(vey * cos(vooe), vex) + M2PI, M2PI);
+  return DBL2NUM(vra * R2D / 15.0);
+}
+
+static VALUE
+func_declination(VALUE self, VALUE vd)
+{
+  double vex =
+  NUM2DBL(func_ecliptic_x(self, vd));
+  double vey =
+  NUM2DBL(func_ecliptic_y(self, vd));
+  double vooe =
+  NUM2DBL(func_obliquity_of_ecliptic(self, vd));
+  double ver = sqrt(vex * vex + vey * vey);
+  double vz = vey * sin(vooe);
+  double vdec = atan2(vz, ver);
+  return DBL2NUM(vdec);
+}
+
+static VALUE
+func_sidetime(VALUE self, VALUE vjd)
+{
+  double vd = NUM2DBL(vjd);
+  double vst =
+  fmod(
+    (180 + 357.52911 + 282.9404) +
+    (0.985600281725 + 4.70935E-5) * vd, 360.0);
+  return DBL2NUM(vst / 15.0);
+}
+
+static VALUE
+func_local_sidetime(VALUE self, VALUE vjd, VALUE vlon)
+{
+  double vst = NUM2DBL(func_sidetime(self, vjd));
+  double vlst = NUM2DBL(vlon) / 15.0 + 12.0 + vst;
+  return DBL2NUM(fmod(vlst, 24.0));
+}
+
+static VALUE
+func_dlt(VALUE self, VALUE vd, VALUE vlat)
+{
+  double vsin_alt = sin(-0.8333 * D2R);
+  double vlat_r = NUM2DBL(vlat) * D2R;
+  double vcos_lat = cos(vlat_r);
+  double vsin_lat = sin(vlat_r);
+  double vooe =
+  NUM2DBL(func_obliquity_of_ecliptic(self, vd));
+  double vtl =
+  NUM2DBL(func_true_longitude(self, vd));
+  double vsin_dec = sin(vooe) * sin(vtl);
+  double vcos_dec =
+  sqrt( 1.0 - vsin_dec * vsin_dec );
+  double vdl =
+  acos(
+    (vsin_alt - vsin_dec * vsin_lat) /
+    (vcos_dec * vcos_lat));
+  double vdla = vdl * R2D;
+  double vdlt = vdla / 15.0 * 2.0;
+  return DBL2NUM(vdlt);
+}
+
+static VALUE
+func_diurnal_arc(VALUE self, VALUE vjd, VALUE vlat)
+{
+  double dlt = NUM2DBL(func_dlt(self, vjd, vlat));
+  double da = dlt / 2.0;
+  return DBL2NUM(da);
+}
+
+static VALUE
+func_t_south(VALUE self, VALUE vjd, VALUE vlon)
+{
+  double lst = NUM2DBL(func_local_sidetime(self, vjd, vlon));
+  double ra = NUM2DBL(func_right_ascension(self, vjd));
+  double vx = lst - ra;
+  double vt = vx - 24.0 * floor(vx * INV24 + 0.5);
+  return DBL2NUM(12 - vt);
+}
+
+static VALUE
+func_t_rise(VALUE self, VALUE vjd, VALUE vlon, VALUE vlat)
+{
+  double ts = NUM2DBL(func_t_south(self, vjd, vlon));
+  double da = NUM2DBL(func_diurnal_arc(self, vjd, vlat));
+  return DBL2NUM(ts - da);
+}
+
+static VALUE
+func_t_set(VALUE self, VALUE vjd, VALUE vlon, VALUE vlat)
+{
+  double ts = NUM2DBL(func_t_south(self, vjd, vlon));
+  double da = NUM2DBL(func_diurnal_arc(self, vjd, vlat));
+  return DBL2NUM(ts + da);
+}
+
+static VALUE
+func_t_mid_day(VALUE self, VALUE vjd, VALUE vlon, VALUE vlat)
+{
+  double tr = NUM2DBL(func_t_rise(self, vjd, vlon, vlat));
+  double ts = NUM2DBL(func_t_set(self, vjd, vlon, vlat));
+  return DBL2NUM((tr + ts) / 2.0);
+}
+
+void Init_calc_sun(void)
+{
+  VALUE cCalcSun =
+  rb_define_class("CalcSun", rb_cObject);
+  rb_define_method(cCalcSun, "initialize", t_init, 0);
+  rb_define_method(cCalcSun,
+  "reverse_12", func_rev12, 1);
+  rb_define_method(cCalcSun,
+  "mean_anomally", func_mean_anomally, 1);
+  rb_define_method(cCalcSun,
+  "eccentricity", func_eccentricity, 1);
+  rb_define_method(cCalcSun,
+  "equation_of_center", func_equation_of_center, 1);
+  rb_define_method(cCalcSun,
+  "true_anomally", func_true_anomally, 1);
+  rb_define_method(cCalcSun,
+  "mean_longitude", func_mean_longitude, 1);
+  rb_define_method(cCalcSun,
+  "eccentric_anomally",
+  func_eccentric_anomally, 1);
+  rb_define_method(cCalcSun,
+  "obliquity_of_ecliptic",
+  func_obliquity_of_ecliptic, 1);
+  rb_define_method(cCalcSun,
+  "longitude_of_perihelion",
+  func_longitude_of_perihelion, 1);
+  // rb_define_method(cCalcSun,
+  // "xv", func_xv, 1);
+  // rb_define_method(cCalcSun,
+  // "yv", func_yv, 1);
+  rb_define_method(cCalcSun,
+  "true_longitude", func_true_longitude, 1);
+  // rb_define_method(cCalcSun,
+  // "rv", func_rv, 1);
+  // rb_define_method(cCalcSun,
+  // "ecliptic_x", func_ecliptic_x, 1);
+  // rb_define_method(cCalcSun,
+  // "ecliptic_y", func_ecliptic_y, 1);
+  rb_define_method(cCalcSun,
+  "right_ascension", func_right_ascension, 1);
+  rb_define_method(cCalcSun,
+  "declination", func_declination, 1);
+  rb_define_method(cCalcSun,
+  "sidereal_time", func_sidetime, 1);
+  rb_define_method(cCalcSun,
+  "local_sidereal_time", func_local_sidetime, 2);
+  rb_define_method(cCalcSun,
+  "dlt", func_dlt, 2);
+  rb_define_method(cCalcSun,
+  "diurnal_arc", func_diurnal_arc, 2);
+  rb_define_method(cCalcSun,
+  "t_south", func_t_south, 2);
+  rb_define_method(cCalcSun,
+  "t_rise", func_t_rise, 3);
+  rb_define_method(cCalcSun,
+  "t_set", func_t_set, 3);
+  rb_define_method(cCalcSun,
+  "t_mid_day", func_t_mid_day, 3);
+}

data/ext/calc_sun/extconf.rb

@@ -0,0 +1,4 @@
+require 'mkmf'
+extension_name = 'calc_sun/calc_sun'
+dir_config(extension_name)
+create_makefile(extension_name)

data/lib/calc_sun/version.rb

@@ -0,0 +1,6 @@
+#!/usr/bin/env ruby
+# class file CalcSun
+class CalcSun
+  # version 0.1.0
+  VERSION = '0.1.1'
+end

data/lib/calc_sun.rb

@@ -0,0 +1,3 @@
+#!/usr/bin/env ruby
+require 'calc_sun/version'
+require 'calc_sun/calc_sun'

data/test/calc_sun_test.rb

@@ -0,0 +1,230 @@
+
+gem 'minitest'
+require 'minitest/autorun'
+require 'test/unit'
+lib = File.expand_path('../../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'calc_sun'
+#
+class CalcSunTest < Minitest::Test
+  def test_that_it_has_a_version_number
+    refute_nil ::CalcSun::VERSION
+  end
+
+  def test_it_does_something_useful
+    assert false
+  end
+end
+require 'date'
+#
+# class TestCalcSun < MiniTest::Test
+class TestCalcSun100 < Test::Unit::TestCase
+  def setup
+    @t = CalcSun.new
+    @t_ajd = 0.0
+    @t_lat = 0.0
+    @t_lon = 0.0
+  end
+
+  def test_mean_anomaly
+    assert_equal(
+      6.240059966692,
+      @t.mean_anomaly(@t_ajd).round(12)
+    )
+  end
+
+  def test_eccentricity
+    assert_equal(
+      0.016709,
+      @t.eccentricity(@t_ajd)
+    )
+  end
+
+  def test_equation_of_center
+    assert_equal(
+      -0.001471380867,
+      @t.equation_of_center(@t_ajd).round(12)
+    )
+  end
+
+  def test_true_anomaly
+    assert_equal(
+      6.238588585825,
+      @t.true_anomaly(@t_ajd).round(12)
+    )
+  end
+
+  def test_mean_longitude
+    assert_equal(
+      4.895063110817,
+      @t.mean_longitude(@t_ajd).round(12)
+    )
+  end
+
+  def test_eccentric_anomaly
+    assert_equal(
+      4.878582250862,
+      @t.eccentric_anomaly(@t_ajd).round(12)
+    )
+  end
+
+  def test_obliquity_of_ecliptic
+    assert_equal(
+      0.409092802283,
+      @t.obliquity_of_ecliptic(@t_ajd).round(12)
+    )
+  end
+
+  def test_xv
+    assert_equal(
+      0.148720277673,
+      @t.xv(@t_ajd).round(12)
+    )
+  end
+
+  def test_yv
+    assert_equal(
+      -0.986083974099,
+      @t.yv(@t_ajd).round(12)
+    )
+  end
+
+  def test_longitude_of_perihelion
+    assert_equal(
+      4.93824156691,
+      @t.longitude_of_perihelion(@t_ajd).round(12)
+    )
+  end
+
+  def test_true_longitude
+    assert_equal(
+      4.893644845555,
+      @t.true_longitude(@t_ajd).round(12)
+    )
+  end
+
+  def test_rv
+    assert_equal(
+      0.997235842199,
+      @t.rv(@t_ajd).round(12)
+    )
+  end
+end
+
+#
+class TestCalcSun200 < Test::Unit::TestCase
+  def setup
+    @t = CalcSun.new
+    @t_ajd = 0.0
+    @t_lat = 0.0
+    @t_lon = 0.0
+  end
+
+  def test_ecliptic_x
+    assert_equal(
+      0.17976672602,
+      @t.ecliptic_x(@t_ajd).round(12)
+    )
+  end
+
+  def test_ecliptic_y
+    assert_equal(
+      -0.980899204395,
+      @t.ecliptic_y(@t_ajd).round(12)
+    )
+  end
+
+  def test_right_ascension
+    assert_equal(
+      18.753078192426,
+      @t.right_ascension(@t_ajd).round(12)
+    )
+  end
+
+  def test_declination
+    assert_equal(
+      -0.372949956542,
+      @t.declination(@t_ajd).round(12)
+    )
+  end
+
+  def test_sidereal_time
+    assert_equal(
+      6.697967333333,
+      @t.sidereal_time(@t_ajd).round(12)
+    )
+  end
+
+  def test_local_sidereal_time
+    assert_equal(
+      18.697967333333,
+      @t.local_sidereal_time(@t_ajd, @t_lon).round(12)
+    )
+  end
+
+  def test_dlt
+    assert_equal(
+      12.120732161881,
+      @t.dlt(@t_ajd, @t_lat).round(12)
+    )
+  end
+
+  def test_diurnal_arc
+    assert_equal(
+      6.06036608094,
+      @t.diurnal_arc(@t_ajd, @t_lat).round(12)
+    )
+  end
+
+  def test_t_south
+    assert_equal(
+      12.055110859092,
+      @t.t_south(@t_ajd, @t_lon).round(12)
+    )
+  end
+
+  def test_t_rise
+    assert_equal(
+      5.994744778152,
+      @t.t_rise(@t_ajd, @t_lon, @t_lat).round(12)
+    )
+  end
+
+  def test_t_mid_day
+    assert_equal(
+      12.055110859092,
+      @t.t_mid_day(@t_ajd, @t_lon, @t_lat).round(12)
+    )
+  end
+
+  def test_t_set
+    assert_equal(
+      18.115476940033,
+      @t.t_set(@t_ajd, @t_lon, @t_lat).round(12)
+    )
+  end
+
+  def test_rise_time
+    rise = @t.t_rise(@t_ajd, @t_lon, @t_lat).round(12)
+    assert_equal(
+      "Sun rises \t\t\t : 5:59 UTC",
+      "Sun rises \t\t\t : #{rise.floor}:#{(rise % 1 * 60.0).floor} UTC"
+    )
+  end
+
+  def test_midday_time
+    midday = @t.t_mid_day(@t_ajd, @t_lon, @t_lat).round(12)
+    assert_equal(
+      "Sun mid day \t\t\t : 12:3 UTC",
+      "Sun mid day \t\t\t : #{midday.floor}:#{((midday % 1.0) * 60).floor} UTC"
+    )
+  end
+
+  def test_set_time
+    set = @t.t_set(@t_ajd, @t_lon, @t_lat).round(12)
+    assert_equal(
+      "Sun sets \t\t\t : 18:6 UTC",
+      "Sun sets \t\t\t : #{set.floor}:#{(set % 1 * 60.0).floor} UTC"
+    )
+  end
+end

metadata

@@ -0,0 +1,130 @@
+--- !ruby/object:Gem::Specification
+name: calc_sun
+version: !ruby/object:Gem::Version
+  version: 0.1.1
+platform: ruby
+authors:
+- kb9agt
+autorequire: 
+bindir: bin
+cert_chain:
+- |
+  -----BEGIN CERTIFICATE-----
+  MIIDaDCCAlCgAwIBAgIBATANBgkqhkiG9w0BAQUFADA9MQ8wDQYDVQQDDAZrYjlh
+  Z3QxFTATBgoJkiaJk/IsZAEZFgVnbWFpbDETMBEGCgmSJomT8ixkARkWA2NvbTAe
+  Fw0xNjEyMzAwNDMxMDhaFw0xNzEyMzAwNDMxMDhaMD0xDzANBgNVBAMMBmtiOWFn
+  dDEVMBMGCgmSJomT8ixkARkWBWdtYWlsMRMwEQYKCZImiZPyLGQBGRYDY29tMIIB
+  IjANBgkqhkiG9w0BAQEFAAOCAQ8AMIIBCgKCAQEAwUJbvvbjVVV2ftkL3R7RcRSO
+  6Mosw9c6vKViXLyKIXyzkXIwtxoygdNj8/WE2Ld/bxxp/tnSeOno1gMWpU5NqpKn
+  FS7+Rd/ihkHa0JUh5of6o+p9MmHVc+e3AZPVRrdq24P9465mTkQRx2j50c3mPmGO
+  AW9gT93A/BJtp91KSpExQmJ6VN0I7AxGYodlegRDJRw4o8Ap5WS8GHjOdpZSKbyq
+  ipD303lPvWEnp9Yra+VzqUXsIdSGVG+As78npWxIn8fZDX2OsYT2qyctYvJiSsqB
+  xMfaBX8VjWPDf3LP3LgLBpUzCGL9KV4iLtcuZaz0OgE9BahA7Kp2T1OjP4fbcQID
+  AQABo3MwcTAJBgNVHRMEAjAAMAsGA1UdDwQEAwIEsDAdBgNVHQ4EFgQUAEuv016V
+  8m/v3zAssXTNrfv76QAwGwYDVR0RBBQwEoEQa2I5YWd0QGdtYWlsLmNvbTAbBgNV
+  HRIEFDASgRBrYjlhZ3RAZ21haWwuY29tMA0GCSqGSIb3DQEBBQUAA4IBAQAXLnar
+  /NKG0Wf2KjYnjPFu3+B+2O6z83VYLp6WqtH6YHpGmYYaYIEpIOOVzwZnG3LVJacj
+  XR+KyspHMyq/7wGmCFzcaUAb77zknqOkpSYVKlRZRdgtJ6oDMIwFFZN9Xwdg7mMj
+  yMevwJb6U/d+XQC8FFxMJsL2u1a5jl1jXo3Sv6ZUxEhlNJSFOsFFQqTfR8Sq8kMm
+  9tEcHmnS7RQS2IQzljk++ha4q0G+qqZz81bjEX0mE6wGYz7mda6lJ40mUKyfA1KH
+  J4L8DCnzZFU8ARwINSdPvhk8WzXpVgErPQezTziVg7gAtBjqdJO5qVlAiOp+z7m2
+  gT57pae1qMGtqvMp
+  -----END CERTIFICATE-----
+date: 2016-12-30 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: rake-compiler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.9.3
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.9.3
+- !ruby/object:Gem::Dependency
+  name: rdoc
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '4.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '4.0'
+- !ruby/object:Gem::Dependency
+  name: hoe
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.15'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.15'
+description: |-
+  Calculate Sunrise and Sunset. Now uses native C extension.
+  An example is supplied to show gem usage.
+email:
+- kb9agt@gmail.com
+executables: []
+extensions:
+- ext/calc_sun/extconf.rb
+extra_rdoc_files:
+- History.rdoc
+- LICENSE.txt
+- Manifest.txt
+- README.rdoc
+files:
+- Gemfile
+- History.rdoc
+- LICENSE.txt
+- Manifest.txt
+- README.rdoc
+- Rakefile
+- calc_sun.gemspec
+- example/sunriset.rb
+- ext/calc_sun/calc_sun.c
+- ext/calc_sun/extconf.rb
+- lib/calc_sun.rb
+- lib/calc_sun/version.rb
+- test/calc_sun_test.rb
+homepage: https://rubygems.org/gems/calc_sun
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options:
+- "--main"
+- README.rdoc
+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.8
+signing_key: 
+specification_version: 4
+summary: Calculate Sunrise and Sunset
+test_files: []