calc_sun 0.1.1 → 1.2.6

data/ext/side_time/extconf.rb

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

data/ext/side_time/side_time.c

@@ -0,0 +1,117 @@
+#include <ruby.h>
+#include <math.h>
+#include <time.h>
+/* if PI's not defined, define it */
+#ifndef PI
+#define PI 3.14159265358979323846
+#endif
+#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.0L
+# define RND12 1000000000000.0
+
+/*
+ * call-seq:
+ *  initialize()
+ *
+ * Create CalcSun class Ruby object.
+ *
+ */
+static VALUE t_init(VALUE self){
+  return self;
+}
+/*
+ * call-seq:
+ *  date('yyyy-mm-dd')
+ *
+ * convert input string to Date object.
+ *
+ */
+static VALUE func_set_datetime(VALUE self, VALUE vdate){
+  VALUE cDate = rb_const_get(rb_cObject, rb_intern("DateTime"));
+  VALUE day = rb_funcall(cDate, rb_intern("parse"), 1, vdate);
+  return day;
+}
+/*
+ * call-seq:
+ *  jd(date)
+ *
+ * convert Date or DateTime object to JD number.
+ *
+ */
+static VALUE func_get_jd(VALUE self, VALUE vdate){
+  double jd = NUM2DBL(rb_funcall(vdate, rb_intern("jd"), 0));
+  return DBL2NUM(jd);
+}
+/*
+ * call-seq:
+ *  ajd(date)
+ *
+ * convert Date or DateTime object to AJD number.
+ *
+ */
+static VALUE func_get_ajd(VALUE self, VALUE vdate){
+  double ajd = NUM2DBL(rb_funcall(vdate, rb_intern("ajd"), 0));
+  return DBL2NUM(ajd);
+}
+/*
+ * call-seq:
+ *  gmst(date)
+ *
+ * calculate Greenwhich Mean Sidereal Time.
+ *
+ */
+static VALUE func_mean_sidetime(VALUE self, VALUE vdate){
+  double ajd = NUM2DBL(func_get_ajd(self, vdate));
+  long double sidereal;
+  long double t;
+  t = (ajd - 2451545.0) / 36525.0;
+  /* calc mean angle */
+  sidereal =
+  280.46061837 +
+  (360.98564736629 * (ajd - 2451545.0)) +
+  (0.000387933 * t * t) -
+  (t * t * t / 38710000.0);
+  sidereal = fmod(sidereal, 360.0);
+  /* change to hours */
+  sidereal *= 24.0 / 360.0;
+  return DBL2NUM(roundf(sidereal * RND12) / RND12);
+}
+/*
+ * call-seq:
+ *  lmst(date)
+ *
+ * calculate Local Mean Sidereal Time.
+ *
+ */
+static VALUE func_local_sidetime(VALUE self, VALUE vdate, VALUE vlon){
+  double sidereal;
+  double ls;
+  double lon = NUM2DBL(vlon);
+  sidereal = NUM2DBL(func_mean_sidetime(self, vdate));
+  sidereal *= 15.0;
+  ls = sidereal + lon;
+  ls *= 24.0 / 360.0;
+  return DBL2NUM(roundf(ls * RND12) / RND12);
+}
+
+
+void Init_side_time(void){
+  VALUE cSideTime = rb_define_class("SideTime", rb_cObject);
+  rb_require("date");
+  rb_define_method(cSideTime, "initialize", t_init, 0);
+  rb_define_attr(cSideTime, "date", 1, 1);
+  rb_define_method(cSideTime, "ajd", func_get_ajd, 1);
+  rb_define_method(cSideTime, "s_datetime", func_set_datetime, 1);
+  rb_define_method(cSideTime, "jd", func_get_jd, 1);
+  //rb_define_method(cSideTime, "jd2000_dif", func_jd_from_2000, 1);
+  //rb_define_method(cSideTime, "jd2000_dif_lon", func_days_from_2000, 2);
+  rb_define_method(cSideTime, "lmst", func_local_sidetime, 2);
+  rb_define_method(cSideTime, "gmst", func_mean_sidetime, 1);
+}

data/lib/calc_sun.rb

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

data/lib/calc_sun/version.rb

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

data/lib/side_time/version.rb

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

data/lib/sidereal_time.rb

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

metadata

@@ -1,10 +1,10 @@
 --- !ruby/object:Gem::Specification
 name: calc_sun
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 1.2.6
 platform: ruby
 authors:
-- kb9agt
+- Douglas Allen
 autorequire: 
 bindir: bin
 cert_chain:
@@ -30,101 +30,76 @@ cert_chain:
   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'
+date: 2017-04-06 00:00:00.000000000 Z
+dependencies: []
 description: |-
-  Calculate Sunrise and Sunset. Now uses native C extension.
-  An example is supplied to show gem usage.
+  supply the date yyyy-mm-dd, latitude decimal, and
+                          longitude decimal in the calls to each method
+                          rise(date, lat, lon),
+                          noon(date, lat, lon), and
+                          set(date, lat, lon)
 email:
 - kb9agt@gmail.com
 executables: []
 extensions:
 - ext/calc_sun/extconf.rb
-extra_rdoc_files:
-- History.rdoc
-- LICENSE.txt
-- Manifest.txt
-- README.rdoc
+extra_rdoc_files: []
 files:
+- .autotest
+- .gitignore
+- .source_index
+- .travis.yml
+- CODE_OF_CONDUCT.md
 - Gemfile
+- Guardfile
 - History.rdoc
 - LICENSE.txt
 - Manifest.txt
 - README.rdoc
 - Rakefile
 - calc_sun.gemspec
+- dev.md
+- example/Almanac.html
+- example/Makefile
+- example/adt_helper.rb
+- example/askgeo_query.rb
+- example/eq_center.md
+- example/sunriset.c
+- example/sunriset.exe
 - example/sunriset.rb
+- example/test_my_times.rb
 - ext/calc_sun/calc_sun.c
 - ext/calc_sun/extconf.rb
+- ext/side_time/extconf.rb
+- ext/side_time/side_time.c
 - lib/calc_sun.rb
+- lib/calc_sun/calc_sun.so
 - lib/calc_sun/version.rb
-- test/calc_sun_test.rb
-homepage: https://rubygems.org/gems/calc_sun
+- lib/side_time/side_time.so
+- lib/side_time/version.rb
+- lib/sidereal_time.rb
+homepage: https://github.com/DouglasAllen/calc_sun
 licenses:
 - MIT
 metadata: {}
 post_install_message: 
-rdoc_options:
-- "--main"
-- README.rdoc
+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.8
+rubygems_version: 2.6.10
 signing_key: 
 specification_version: 4
-summary: Calculate Sunrise and Sunset
+summary: Calculates Sun Times ea. rise, noon, set
 test_files: []

data/test/calc_sun_test.rb

@@ -1,230 +0,0 @@
-
-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