orbiter 0.0.2

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 0ff7911d0476677d3acb652adbe1d4723fa61199
+  data.tar.gz: f4a056a95ff34c1ccaf12796c2a719cee4e6a4ea
+SHA512:
+  metadata.gz: 1ab6aa144eb18f4a5a95c1f48d2984e14e6174c41bde0a9691ccc77b84f977275af662e83946895eaa032994e647f8b2ca6ccf8af5caa7506f47b118a596d1a2
+  data.tar.gz: 5fbebd1b4434c79dfd959f9837592f4e4f3ea43b169d8e0a25ca3ce9583434183acbb18ba50dfb6f55492357214b292a42d037813fe40850f32cb30f195376ce

data/lib/free_body.rb

@@ -0,0 +1,17 @@
+module Orbiter	
+	class Free_body
+		attr_accessor :mass, :x, :y, :vel_x, :vel_y, :acc_x, :acc_y, :old_acc_x, :old_acc_y
+	
+		def initialize(attrs = {})
+			@mass = attrs.fetch(:mass, 10)
+			@x = attrs.fetch(:x, 0)
+			@y = attrs.fetch(:y, 0)
+			@vel_x = attrs.fetch(:vel_x, 0)
+			@vel_y = attrs.fetch(:vel_y, 0)
+			@acc_x = attrs.fetch(:acc_x, 0)
+			@acc_y = attrs.fetch(:acc_y, 0)
+			@old_acc_x = attrs.fetch(:old_acc_x, 0)
+			@old_acc_y = attrs.fetch(:old_acc_y, 0)
+		end
+	end
+end

data/lib/orbit_helper.rb

@@ -0,0 +1,2 @@
+require_relative 'free_body.rb'
+require_relative 'orbit_updater.rb'

data/lib/orbit_updater.rb

@@ -0,0 +1,50 @@
+module Orbiter
+	class Orbit_updater
+	
+		#Changing these breaks tests, which expect both to = 1
+		G_CONSTANT = 1
+		TIME_CONSTANT = 1
+	
+		def self.calc_dist(body_a, body_b)
+			Math.sqrt( (body_a.x - body_b.x)**2 + (body_a.y - body_b.y)**2 )
+		end
+	
+		def self.calc_grav_acc(target, other_bodies)
+			target.old_acc_x = target.acc_x
+			target.old_acc_y = target.acc_y
+			target.acc_x = 0.0
+			target.acc_y = 0.0
+			other_bodies.each do |other_body|
+				if (dist = calc_dist(target, other_body)) > 0
+					acc = G_CONSTANT*other_body.mass / (dist ** 2)
+					target.acc_x += acc * Math.cos(calc_angle(target, other_body)) 
+					target.acc_y += acc * Math.sin(calc_angle(target, other_body))
+				end
+			end
+		end
+	
+		def self.calc_angle(body_a, body_b)
+			y_diff = body_b.y - body_a.y
+			x_diff = body_b.x - body_a.x
+			Math.atan2(y_diff, x_diff)
+		end
+	
+		def self.update_velocity(target)
+			#binding.pry
+			target.vel_x += TIME_CONSTANT * (target.old_acc_x + target.acc_x ) / 2
+			target.vel_y += TIME_CONSTANT * (target.old_acc_y + target.acc_y ) / 2
+		end
+	
+		def self.update_position(target)
+			target.x += TIME_CONSTANT * ( target.vel_x + TIME_CONSTANT * target.acc_x / 2 )
+			target.y += TIME_CONSTANT * ( target.vel_y + TIME_CONSTANT * target.acc_y / 2 )
+		end
+	
+		def self.track(file, bodies)
+			bodies.each do |body|
+				file.print("#{body.x}\t#{body.y}\t")
+			end
+			file.puts()
+		end
+	end
+end

data/lib/orbiter.rb

@@ -0,0 +1,27 @@
+require_relative 'orbit_helper.rb'
+
+module Orbiter	
+
+	class Orbit_runner
+
+		def self.run(attrs = {})
+			@bodies = attrs.fetch(:bodies, [Orbiter::Free_body.new(mass: 10, x: 100, y: 100, vel_x: -5, vel_y: -1),
+																			Orbiter::Free_body.new(mass:1000)])
+			@output = attrs[:output] ? File.open(attrs[:output], 'w') : $stdout
+	
+			100.times do 
+				@bodies.each do |body|
+					Orbiter::Orbit_updater.calc_grav_acc(body, @bodies)
+					Orbiter::Orbit_updater.update_position(body)
+					Orbiter::Orbit_updater.calc_grav_acc(body, @bodies)
+					Orbiter::Orbit_updater.update_velocity(body)
+				end
+				Orbiter::Orbit_updater.track(@output, @bodies)
+			end
+	
+			@output.close
+		
+			puts "Thanks for running!"
+		end
+	end
+end

metadata

@@ -0,0 +1,47 @@
+--- !ruby/object:Gem::Specification
+name: orbiter
+version: !ruby/object:Gem::Version
+  version: 0.0.2
+platform: ruby
+authors:
+- Will Weiss
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2014-02-25 00:00:00.000000000 Z
+dependencies: []
+description: 
+email: wweiss00@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- lib/orbiter.rb
+- lib/free_body.rb
+- lib/orbit_helper.rb
+- lib/orbit_updater.rb
+homepage: 
+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.0.7
+signing_key: 
+specification_version: 4
+summary: For tracking the orbits of n free bodies in a 2D space
+test_files: []