cf_sim 1.0.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 9378e236733ed66fb70343fa5e33c508cfdba69f
+  data.tar.gz: 603031535dbc93e9dacc2ec9fce5f6415f9e6ba4
+SHA512:
+  metadata.gz: 4a7aca9a5cf5bae58b1e715c8276f9faee4418718f00fd7230a34a877ef59631e911dffc1e03af3bf34796aedc36955b68ffc85a890640ab4f2458a738904683
+  data.tar.gz: aca9a38f0f726005d54b7a40cb3901dabb27551cf6521af330523a5f44233b2339633b1a1a2114b2b6b2425ea8d34f23275b0f8d20a86676534c7b4e6eecfcc9

data/.gitignore

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

data/.travis.yml

@@ -0,0 +1,5 @@
+language: ruby
+rvm:
+  - 2.0
+  - 2.1
+  - 2.2

data/CODE_OF_CONDUCT.md

@@ -0,0 +1,13 @@
+# Contributor Code of Conduct
+
+As contributors and maintainers of this project, we pledge to respect all people who contribute through reporting issues, posting feature requests, updating documentation, submitting pull requests or patches, and other activities.
+
+We are committed to making participation in this project a harassment-free experience for everyone, regardless of level of experience, gender, gender identity and expression, sexual orientation, disability, personal appearance, body size, race, age, or religion.
+
+Examples of unacceptable behavior by participants include the use of sexual language or imagery, derogatory comments or personal attacks, trolling, public or private harassment, insults, or other unprofessional conduct.
+
+Project maintainers have the right and responsibility to remove, edit, or reject comments, commits, code, wiki edits, issues, and other contributions that are not aligned to this Code of Conduct. Project maintainers who do not follow the Code of Conduct may be removed from the project team.
+
+Instances of abusive, harassing, or otherwise unacceptable behavior may be reported by opening an issue or contacting one or more of the project maintainers.
+
+This Code of Conduct is adapted from the [Contributor Covenant](http:contributor-covenant.org), version 1.0.0, available at [http://contributor-covenant.org/version/1/0/0/](http://contributor-covenant.org/version/1/0/0/)

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in cf_sim.gemspec
+gemspec

data/LICENSE

@@ -0,0 +1,22 @@
+The MIT License (MIT)
+
+Copyright (c) 2015 pinzolo
+
+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,41 @@
+# CF Simulator
+[![Build Status](https://travis-ci.org/pinzolo/cf_sim.svg)](http://travis-ci.org/pinzolo/cf_sim)
+[![Coverage Status](https://coveralls.io/repos/pinzolo/cf_sim/badge.png)](https://coveralls.io/r/pinzolo/cf_sim)
+
+Welcome to your new gem! In this directory, you'll find the files you need to be able to package up your Ruby library into a gem. Put your Ruby code in the file `lib/cf_sim`. To experiment with that code, run `bin/console` for an interactive prompt.
+
+TODO: Delete this and the text above, and describe your gem
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'cf_sim'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install cf_sim
+
+## Usage
+
+TODO: Write usage instructions here
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `bin/console` for an interactive prompt that will allow you to experiment. Run `bundle exec cf_sim` to use the code located 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` to create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+1. Fork it ( https://github.com/pinzolo/cf_sim/fork )
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create a new Pull Request

data/Rakefile

@@ -0,0 +1,9 @@
+require "bundler/gem_tasks"
+require 'rake/testtask'
+
+Rake::TestTask.new do |t|
+  t.libs << 'test'
+  t.test_files = FileList['test/cf_sim/*_test.rb']
+  t.verbose = true
+end
+task default: :test

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "cf_sim"
+
+# 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,7 @@
+#!/bin/bash
+set -euo pipefail
+IFS=$'\n\t'
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/cf_sim.gemspec

@@ -0,0 +1,27 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'cf_sim/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "cf_sim"
+  spec.version       = CfSim::VERSION
+  spec.authors       = ["pinzolo"]
+  spec.email         = ["pinzolo@gmail.com"]
+  spec.summary       = %q{Multi layered CF simulator}
+  spec.description   = %q{Propose multi layered control fields pattern by given portals.}
+  spec.homepage      = "https://github.com/pinzolo/cf_sim"
+  spec.license       = "MIT"
+
+  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.8"
+  spec.add_development_dependency "rake"#, "~> 10.0"
+  spec.add_development_dependency 'test-unit'
+  spec.add_development_dependency 'coveralls'
+
+  spec.add_dependency 'thor'
+end

data/exe/cf_sim

@@ -0,0 +1,4 @@
+#!/usr/bin/env ruby
+
+require "cf_sim"
+CfSim::CLI.start

data/lib/cblxy.rb

@@ -0,0 +1,262 @@
+# see https://github.com/KMR-zoar/cblxy
+include Math
+
+
+#定数の定義
+#ベッセル楕円体
+$aB = 6377397.155 #長半径
+$FB = 299.152813  #逆扁平率
+#GRS80楕円体
+$aG = 6378137.0  #長半径
+$FG = 298.257222101 #逆扁平率
+#共通定数
+$sbyS = 0.9999 #原点における縮尺係数
+$a = 0.0 #長半径
+$F = 0.0 #逆扁平率
+$e = 0.0 #第一離心率
+$phi0 = 0.0 #原点緯度
+$lambda0 = 0.0 #原点経度
+$bits = ARGV[2] #出力値(0 = 緯度 or X, 1 = 経度 or Y)
+
+#世界測地系で緯度経度を平面直角座標に変換するメソッド
+def blxy(phi1, lamda1,k)
+   $a = $aG #GRS80楕円体の定数を代入
+   $F = $FG
+   genten(k)
+
+   bl2xy(phi1, lamda1)
+
+end
+
+#世界測地系で平面直角座標を緯度経度に変換するメソッド
+def xybl(phi1,lamda1,k)
+   $a = $aG #GRS80楕円体の定数を代入
+   $F = $FG
+   genten(k)
+
+   xy2bl(phi1, lamda1)
+end
+
+#緯度経度を平面直角座標に変換する計算
+def bl2xy(phi, lamda)
+   $e = sqrt(2.0 * $F - 1.0) / $F
+
+   phi1 = deg2rad(phi)
+   lamda1 = deg2rad(lamda)
+
+   s0 = kocyou($phi0)
+   s1 = kocyou(phi1)
+
+   ut = $a / sqrt(1.0 - $e ** 2.0 * sin(phi1) ** 2.0)
+   conp = cos(phi1)
+   t1 = tan(phi1)
+   dlamda = lamda1 - $lamda0
+   eta2 = ($e ** 2.0 / (1.0 - $e ** 2.0)) * conp ** 2.0
+
+   #X座標値の算出
+   v1 = 5.0 - t1 ** 2.0 + 9.0 * eta2 + 4.0 * eta2 ** 2.0
+   v2 = -61.0 + 58.0 * t1 ** 2.0 - t1 ** 4.0 - 270.0 * eta2 + 330.0 * t1 ** 2.0 * eta2
+   v3 = -1385.0 + 3111.0 * t1 ** 2.0 - 543.0 * t1 ** 4.0 + t1 ** 6.0
+
+   x = ((s1 - s0) + ut * conp ** 2.0 * t1 * dlamda ** 2.0 / 2.0 + ut * conp ** 4.0 * t1 * v1 * dlamda ** 4.0 / 24.0 - ut * conp ** 6.0 * t1 * v2 * dlamda ** 6.0 / 720.0 - ut * conp ** 8.0 * t1 * v3 * dlamda ** 8.0 / 40320.0) * $sbyS
+
+   #Y座標値の算出
+   v1 = -1.0 + t1 ** 2.0 - eta2
+   v2 = -5.0 + 18.0 * t1 ** 2.0 - t1 ** 4.0 - 14.0 * eta2 + 58.0 * t1 ** 2.0 * eta2
+   v3 = -61.0 + 479.0 * t1 ** 2.0 - 179.0 * t1 ** 4.0 + t1 ** 6.0
+
+   y = (ut * conp * dlamda - ut * conp ** 3.0 * v1 * dlamda ** 3.0 / 6.0 - ut * conp ** 5.0 * v2 * dlamda ** 5.0 / 120.0 - ut * conp ** 7.0 * v3 * dlamda ** 7.0 / 5040.0) * $sbyS
+
+   return [x,y]
+
+end
+
+#平面直角座標から緯度経度に変換する計算
+def xy2bl(x,y)
+   $e = sqrt(2.0 * $F - 1.0) / $F
+
+   phi1 = suisen(x)
+
+   ut = $a / sqrt( 1.0 - $e ** 2.0 * sin(phi1) ** 2.0)
+   conp = cos(phi1)
+   t1 = tan(phi1)
+   eta2 = ($e ** 2.0 / (1.0 - $e ** 2.0)) * conp ** 2.0
+
+   #緯度算出
+   yy = y / $sbyS
+   v1 = 1.0 + eta2
+   v2 = 5.0 + 3.0 * t1 ** 2.0 + 6.0 * eta2 - 6.0 * t1 ** 2.0 * eta2 - 3.0 * eta2 ** 2.0 - 9.0 * t1 ** 2.0 * eta2 ** 2.0
+   v3 = 61.0 + 90.0 * t1 ** 2.0 + 45.0 * t1 ** 4.0 + 107.0 * eta2 - 162.0 * t1 ** 2.0 * eta2 - 45.0 * t1 ** 4.0 * eta2
+   v4 = 1385.0 + 3633.0 * t1 ** 2.0 + 4095.0 * t1 ** 4.0 + 1575.0 * t1 ** 6.0
+
+   phir = -(v1 / (2.0 * ut ** 2.0)) * yy ** 2.0
+   phir = phir + (v2 / (24.0 * ut ** 4.0)) * yy ** 4.0
+   phir = phir - (v3 / (720.0 * ut ** 6.0)) * yy ** 6.0
+   phir = phir + (v4 / (40320.0 * ut ** 8.0)) * yy ** 8.0
+   phir = phir * t1
+   phir = phir + phi1
+   phir = rad2deg(phir)
+
+   #経度算出
+   v1 = ut * conp
+   v2 = 1.0 + 2.0 * t1 ** 2.0 + eta2
+   v3 = 5.0 + 28.0 * t1 ** 2.0 + 24.0 * t1 ** 4.0 + 6.0 * eta2 + 8.0 * t1 ** 2.0 * eta2
+   v4 = 61.0 + 662.0 * t1 ** 2.0 + 1320.0 * t1 ** 4.0 + 720.0 * t1 ** 6.0
+
+   lamdar = (1.0 / v1) * yy
+   lamdar = lamdar - (v2 / (6.0 * ut ** 2.0 * v1)) * yy ** 3.0
+   lamdar = lamdar + (v3 / (120.0 * ut ** 4.0 * v1)) * yy ** 5.0
+   lamdar = lamdar - (v4 / (5040.0 * ut ** 6.0 * v1)) * yy ** 7.0
+   lamdar = lamdar + $lamda0
+
+   lamdar = rad2deg(lamdar)
+
+   return [phir,lamdar]
+end
+
+#座標系原点の緯度経度
+def genten(k)
+
+   err = "座標系が正しく指定されていません"
+
+   case k
+   when 1
+      degen = [330000.0,1293000.0]
+   when 2
+      degen = [330000.0,1310000.0]
+   when 3
+      degen = [360000.0,1321000.0]
+   when 4
+      degen = [330000.0,1333000.0]
+   when 5
+      degen = [360000.0,1342000.0]
+   when 6
+      degen = [360000.0,1360000.0]
+   when 7
+      degen = [360000.0,1371000.0]
+   when 8
+      degen = [360000.0,1383000.0]
+   when 9
+      degen = [360000.0,1395000.0]
+   when 10
+      degen = [400000.0,1405000.0]
+   when 11
+      degen = [440000.0,1401500.0]
+   when 12
+      degen = [440000.0,1421500.0]
+   when 13
+      degen = [440000.0,1441500.0]
+   when 14
+      degen = [260000.0,1420000.0]
+   when 15
+      degen = [260000.0,1273000.0]
+   when 16
+      degen = [260000.0,1240000.0]
+   when 17
+      degen = [260000.0,1310000.0]
+   when 18
+      degen = [200000.0,1360000.0]
+   when 19
+      degen = [260000.0,1540000.0]
+   else
+      degen = err
+   end
+
+   unless degen == err
+      $phi0 = (deg2rad degen[0])
+      $lamda0 = (deg2rad degen[1])
+      return [$phi0,$lamda0]
+   else
+      return err
+   end
+end
+
+
+
+#度分秒からラジアンへの変換
+def deg2rad(deg)
+   radbase = PI / 180.0
+   if deg < 0
+      fugou  = -1.0
+   else
+      fugou = 1.0
+   end
+
+   deg = deg.abs
+
+   angle = (deg / 10000.0).to_i
+   minute = ((deg / 100.0).to_i - (angle * 100.0)).to_f
+   second =  deg - (angle * 10000.0) - (minute * 100.0)
+   rad = fugou * (angle + (minute + (second / 60.0)) / 60.0)
+   rad = rad * radbase
+
+end
+
+#ラジアンから度分秒への変換
+def rad2deg(rad)
+
+   degbase = 180.0 / PI
+
+   if rad < 0
+      fugou = -1.0
+   else
+      fugou = 1.0
+   end
+
+   rad = rad.abs
+
+   rad = rad * degbase
+   angle = rad.to_i
+   rad = (rad - angle).to_f * 60.0
+   minute = rad.to_i
+   second = (rad - minute).to_f * 60.0
+
+   deg = fugou * (angle * 10000.0 + minute * 100.0 + second)
+end
+
+#緯度から赤道への子午線弧長を計算する
+def kocyou(ido)
+   e2 = $e ** 2.0
+   e4 = $e ** 4.0
+   e6 = $e ** 6.0
+   e8 = $e ** 8.0
+   e10 = $e ** 10.0
+   e12 = $e ** 12.0
+   e14 = $e ** 14.0
+   e16 = $e ** 16.0
+
+   a = 1.0 + 3.0 / 4.0 * e2 + 45.0 / 64.0 * e4 + 175.0 / 256.0 * e6 + 11025.0 / 16384.0 * e8 + 43659.0 / 65536.0 * e10 + 693693.0 / 1048576.0 * e12 + 19324305.0 / 29360128.0 * e14 + 4927697775.0 / 7516192768.0 * e16
+   b = 3.0 / 4.0 * e2 + 15.0 / 16.0 * e4 + 525.0 / 512.0 * e6 + 2205.0 / 2048.0 * e8 + 72765.0 / 65536.0 * e10 + 297297.0 / 262144.0 * e12 + 135270135.0 / 117440512.0 * e14 + 547521975.0 / 469762048.0 * e16
+   c = 15.0 / 64.0 * e4 + 105.0 / 256.0 * e6 + 2205.0 / 4096.0 * e8 + 10395.0 / 16384.0 * e10 + 1486485.0 / 2097152.0 * e12 + 45090045.0 / 58720256.0 * e14 + 766530765.0 / 939524096.0 * e16
+   d = 35.0 / 512.0 * e6 + 315.0 / 2048.0 * e8 + 31185.0 / 131072.0 * e10 + 165165.0 / 524288.0 * e12 + 45090045.0 / 117440512.0 * e14 + 209053845.0 / 469762048.0 * e16
+   e = 315.0 / 16384.0 * e8 + 3465.0 / 65536.0 * e10 + 99099.0 / 1048576.0 * e12 + 4099095.0 / 29360128.0 * e14 + 348423075.0 / 1879048192.0 * e16
+   f = 693.0 / 131072 * e10 + 9009.0 / 524288.0 * e12 + 4099095.0 / 117440512.0 * e14 + 26801775.0 / 469762048.0 * e16
+   g = 3003 / 2097152.0 * e12 + 315315.0 / 58720256.0 * e14 + 11486475.0 / 939524096.0 * e16
+   h = 45045.0 / 117440512.0 * e14 + 765765.0 / 469762048.0 * e16
+   i = 765765.0 / 7516192768.0 * e16
+
+   sigosen = $a * (1.0 - e2) * (a * ido - b * sin(ido * 2.0) / 2.0 + c * sin(ido * 4.0) / 4.0 - d * sin(ido * 6.0) / 6.0 + e * sin(ido * 8.0) / 8.0 - f * sin(ido * 10.0) / 10.0 + g * sin(ido * 12.0) / 12.0 - h * sin(ido * 14.0) / 14.0 + i * sin(ido * 16.0) / 16.0)
+end
+
+def suisen(x)
+   s0 = kocyou($phi0)
+   m = s0 + x / $sbyS
+   cnt = 0
+   phin = $phi0
+   e2 = $e ** 2.0
+   phi0 = phin
+
+   while
+      cnt = cnt + 1
+      phi0 = phin
+      sn = kocyou(phin)
+      v1 = 2.0 * (sn - m) * ((1.0 - e2 * sin(phin) ** 2.0) ** 1.5)
+      v2 = 3.0 * e2 * (sn - m) * sin(phin) * cos(phin) * sqrt(1.0 - e2 * sin(phin) ** 2.0) - 2.0 * $a * (1.0 - e2)
+      phin = phin +v1 / v2
+      if ((phin - phi0).abs < 0.00000000000001) or cnt > 100
+         break
+      end
+   end
+   return phin
+
+end

data/lib/cf_sim/cli.rb

@@ -0,0 +1,31 @@
+require 'thor'
+
+class CfSim::CLI < Thor
+  desc 'max_area DATA_FILE', 'print max area control fields URL.'
+  def max_area(data_file)
+    finder = load_finder_from(data_file)
+    fields = finder.find_max_area_fields
+    puts @generator.fields_link(fields)
+  end
+
+  desc 'max_count DATA_FILE', 'print max count control fields URL list.'
+  def max_count(data_file)
+    finder = load_finder_from(data_file)
+    fields_list = finder.find_max_count_fields_list
+    max_area = fields_list.first.total_area
+    fields_list.each do |fields|
+      area_rate = (fields.total_area * 100 / max_area).round(3)
+      puts "#{@generator.fields_link(fields)} (Area rate: #{area_rate})"
+    end
+  end
+
+  private
+
+  def load_finder_from(data_file)
+    portals = CfSim::DataFileParser.new(data_file).parse
+    portal_map = CfSim::PortalMap.new(portals)
+    @generator = CfSim::IntelMapUrlGenerator.new(portal_map)
+    CfSim::ControlFieldSetFinder.new(CfSim::PointList.new(portal_map.points).creatable_fields)
+  end
+end
+

data/lib/cf_sim/coexistable_field_map.rb

@@ -0,0 +1,19 @@
+require 'cf_sim/control_field_set'
+
+class CfSim::CoexistableFieldMap
+  include Enumerable
+
+  def initialize(fields)
+    @map = fields.each_with_object({}) do |field, map|
+      map[field] = CfSim::ControlFieldSet.new(fields.reject { |other| field == other || field.intersected?(other) })
+    end
+  end
+
+  def [](field)
+    @map[field]
+  end
+
+  def each
+    @map.each { |field, coexistable_fields| yield field, coexistable_fields }
+  end
+end

data/lib/cf_sim/control_field.rb

@@ -0,0 +1,55 @@
+require 'cf_sim/link'
+
+class CfSim::ControlField
+  attr_reader :point1, :point2, :point3, :link1, :link2, :link3
+
+  def initialize(point1, point2, point3)
+    @point1, @point2, @point3 = point1, point2, point3
+    @link1 = CfSim::Link.new(point1, point2)
+    @link2 = CfSim::Link.new(point2, point3)
+    @link3 = CfSim::Link.new(point3, point1)
+  end
+
+  def ==(other)
+    eql?(other)
+  end
+
+  def eql?(other)
+    (points - other.points).empty?
+  end
+
+  def hash
+    point1.hash + point2.hash + point3.hash
+  end
+
+  def points
+    @points ||= [@point1, @point2, @point3].freeze
+  end
+
+  def links
+    @links ||= [@link1, @link2, @link3].freeze
+  end
+
+  def area
+    @area ||= calculate_area
+  end
+
+  def intersected?(other)
+    links.any? do |link1|
+      other.links.any? do |link2|
+        link1.intersected?(link2)
+      end
+    end
+  end
+
+  def to_s
+    "#{point1.name} -> #{point2.name} -> #{point3.name}"
+  end
+
+  private
+
+  def calculate_area
+    s = links.map(&:length).inject(:+) / 2.0
+    Math.sqrt(s * (s - link1.length) * (s - link2.length) * (s - link3.length))
+  end
+end

data/lib/cf_sim/control_field_set.rb

@@ -0,0 +1,64 @@
+require 'forwardable'
+
+class CfSim::ControlFieldSet
+  extend Forwardable
+  include Enumerable
+
+  attr_reader :control_fields
+  def_delegators :@control_fields, :clear, :each, :empty?, :hash, :size
+
+  def initialize(*fields)
+    @control_fields = fields.flatten.dup
+  end
+
+  def total_area
+    @total_area ||= @control_fields.inject(0) { |area, field| area += field.area }
+  end
+
+  def -(fields)
+    CfSim::ControlFieldSet.new(@control_fields - extract_field_array(fields))
+  end
+
+  def &(fields)
+    CfSim::ControlFieldSet.new(@control_fields & extract_field_array(fields))
+  end
+
+  def +(fields)
+    CfSim::ControlFieldSet.new(@control_fields + extract_field_array(fields))
+  end
+
+  def ==(other)
+    eql?(other)
+  end
+
+  def eql?(other)
+    other.class == CfSim::ControlFieldSet &&
+      @control_fields.size == other.control_fields.size &&
+      (@control_fields - other.control_fields).empty?
+  end
+
+  def dup
+    # initialize 内で dup しているためここでする必要はない
+    CfSim::ControlFieldSet.new(@control_fields)
+  end
+
+  def to_s
+    text = "Fields:\n"
+    @control_fields.each { |field| text << "  #{field}\n" }
+    text << "Total area: #{total_area}"
+    text
+  end
+
+  private
+
+  def extract_field_array(fields)
+    if fields.is_a?(Array)
+      fields
+    elsif fields.is_a?(CfSim::ControlFieldSet)
+      fields.control_fields
+    else
+      raise 'Only for Array or ControlFieldSet'
+    end
+  end
+end
+

data/lib/cf_sim/control_field_set_finder.rb

@@ -0,0 +1,83 @@
+require 'cf_sim/coexistable_field_map'
+require 'cf_sim/control_field_set'
+
+class CfSim::ControlFieldSetFinder
+  def initialize(fields, limit_field_count = 1)
+    @fields = CfSim::ControlFieldSet.new(fields.sort_by(&:area).reverse)
+    @coexistable_field_map = CfSim::CoexistableFieldMap.new(@fields)
+    @limit_field_count = limit_field_count
+  end
+
+  def find_max_area_fields
+    @max_area_control_fields = nil
+    find_recursively(CfSim::ControlFieldSet.new, @fields, :max_area)
+    @max_area_control_fields
+  end
+
+  def find_all_fields_list
+    @control_fields_list = []
+    find_recursively(CfSim::ControlFieldSet.new, @fields, :all)
+    @control_fields_list
+  end
+
+  def find_max_count_fields_list
+    @max_count_control_fields_list = []
+    find_recursively(CfSim::ControlFieldSet.new, @fields, :max_count)
+    @max_count_control_fields_list.sort_by(&:total_area).reverse
+  end
+
+  private
+
+  def find_recursively(fields, coexistable_fields, type)
+    return if ignore?(fields, coexistable_fields, type)
+
+    if coexistable_fields.empty?
+      deploy(fields, type)
+    else
+      coexistable_fields.each do |field|
+        find_recursively(fields + [field], coexistable_fields & (@coexistable_field_map[field].reject { |f| f.area > field.area }) , type)
+      end
+    end
+  end
+
+  def deploy(fields, type)
+    case type
+    when :max_area
+      replace_if_max(fields)
+    when :all
+      append_to_list(fields, @control_fields_list)
+    when :max_count
+      append_or_create_if_max_count(fields)
+    end
+  end
+
+  def replace_if_max(fields)
+    @max_area_control_fields = fields if @max_area_control_fields.nil? || fields.total_area > @max_area_control_fields.total_area
+  end
+
+  def append_to_list(fields, fields_list)
+    fields_list << fields unless fields_list.include?(fields)
+  end
+
+  def append_or_create_if_max_count(fields)
+    if @max_count_control_fields_list.empty? || @max_count_control_fields_list.first.size == fields.size
+      append_to_list(fields, @max_count_control_fields_list)
+    elsif @max_count_control_fields_list.first.size < fields.size
+      @max_count_control_fields_list = [fields]
+    end
+  end
+
+  def ignore?(fields, coexistable_fields, type)
+    field_count = fields.size + coexistable_fields.size
+    case type
+    when :max_area
+      @max_area_control_fields && fields.total_area + coexistable_fields.total_area <= @max_area_control_fields.total_area
+    when :all
+      field_count < @limit_field_count
+    when :max_count
+      field_count < @limit_field_count && @max_count_control_fields_list.any? && field_count < @max_count_control_fields_list.first.size
+    else
+      false
+    end
+  end
+end

data/lib/cf_sim/cordination_system.rb

@@ -0,0 +1,43 @@
+# 本来なら都道府県によって座標系が決定されるべきだが、とりあえず簡易的に最も近いものを対象とする
+class CfSim::CordinationSystem
+  attr_reader :latitude, :longitude
+
+  ORIGIN_MAP = { 1  => [33.00000,129.30000],
+                 2  => [33.00000,131.00000],
+                 3  => [36.00000,132.10000],
+                 4  => [33.00000,133.30000],
+                 5  => [36.00000,134.20000],
+                 6  => [36.00000,136.00000],
+                 7  => [36.00000,137.10000],
+                 8  => [36.00000,138.30000],
+                 9  => [36.00000,139.50000],
+                 10 => [40.00000,140.50000],
+                 11 => [44.00000,140.15000],
+                 12 => [44.00000,142.15000],
+                 13 => [44.00000,144.15000],
+                 14 => [26.00000,142.00000],
+                 15 => [26.00000,127.30000],
+                 16 => [26.00000,124.00000],
+                 17 => [26.00000,131.00000],
+                 18 => [20.00000,136.00000],
+                 19 => [26.00000,154.00000] }
+
+  def initialize(latitude, longitude)
+    @latitude = latitude
+    @longitude = longitude
+  end
+
+  def nearest_system_number
+    @nearest_system_number ||= ORIGIN_MAP.min_by { |_, origin| distance(origin) }.first
+  end
+
+  def nearest_origin
+    @nearest_origin ||= ORIGIN_MAP[nearest_system_number]
+  end
+
+  private
+
+  def distance(origin)
+    Math.sqrt((origin.first - @latitude) ** 2 + (origin.last - @longitude) ** 2)
+  end
+end

data/lib/cf_sim/data_file_parser.rb

@@ -0,0 +1,23 @@
+class CfSim::DataFileParser
+  def initialize(data_file_path)
+    @data_file_path = data_file_path
+  end
+
+  def parse
+    @portals = []
+    File.open(@data_file_path) do |file|
+      file.each_line do |line|
+        next if line.strip.empty?
+        @portals << create_portal_from_line(line)
+      end
+    end
+    @portals
+  end
+
+  private
+
+  def create_portal_from_line(line)
+    latitude, longitude, *rest = line.split(',')
+    CfSim::Portal.new(rest.join(',').strip, latitude.to_f, longitude.to_f)
+  end
+end

data/lib/cf_sim/intel_map_url_generator.rb

@@ -0,0 +1,21 @@
+class CfSim::IntelMapUrlGenerator
+  def initialize(portal_map)
+    @portal_map = portal_map
+  end
+
+  def portal_link(portal)
+    "https://www.ingress.com/intel?ll=#{portal.latitude},#{portal.longitude}&pll=#{portal.latitude},#{portal.longitude}"
+  end
+
+  def fields_link(fields)
+    links = fields.each_with_object([]) { |field, list| list << field.links }.flatten.uniq
+    base_point = @portal_map.find_portal(links.first.source_point)
+    url_base = "https://www.ingress.com/intel?ll=#{base_point.latitude},#{base_point.longitude}&pls="
+    link_url_params = links.map do |link|
+      source_portal = @portal_map.find_portal(link.source_point)
+      destination_portal = @portal_map.find_portal(link.destination_point)
+      "#{source_portal.latitude},#{source_portal.longitude},#{destination_portal.latitude},#{destination_portal.longitude}"
+    end
+    url_base + link_url_params.join('_')
+  end
+end

data/lib/cf_sim/link.rb

@@ -0,0 +1,38 @@
+class CfSim::Link
+  attr_reader :source_point, :destination_point
+
+  def initialize(source_point, destination_point)
+    raise 'Cannot link to same point' if source_point == destination_point
+    @source_point = source_point
+    @destination_point = destination_point
+  end
+
+  def length
+    @length ||= Math.sqrt((source_point.x - destination_point.x) ** 2 + (source_point.y - destination_point.y) ** 2)
+  end
+
+  def ==(other)
+    eql?(other)
+  end
+
+  def eql?(other)
+    (source_point == other.source_point && destination_point == other.destination_point) ||
+      (source_point == other.destination_point && destination_point == other.source_point)
+  end
+
+  def hash
+    source_point.hash + destination_point.hash
+  end
+
+  def intersected?(other)
+    ta = (other.source_point.x - other.destination_point.x) * (source_point.y - other.source_point.y) + (other.source_point.y - other.destination_point.y) * (other.source_point.x - source_point.x)
+    tb = (other.source_point.x - other.destination_point.x) * (destination_point.y - other.source_point.y) + (other.source_point.y - other.destination_point.y) * (other.source_point.x - destination_point.x)
+    tc = (source_point.x - destination_point.x) * (other.source_point.y - source_point.y) + (source_point.y - destination_point.y) * (source_point.x - other.source_point.x)
+    td = (source_point.x - destination_point.x) * (other.destination_point.y - source_point.y) + (source_point.y - destination_point.y) * (source_point.x - other.destination_point.x)
+    tc * td < 0 && ta * tb < 0
+  end
+
+  def to_s
+    "#{source_point}->#{destination_point}"
+  end
+end

data/lib/cf_sim/point.rb

@@ -0,0 +1,24 @@
+class CfSim::Point
+  attr_reader :name, :x, :y
+
+  def initialize(name, x, y)
+    @name = name
+    @x, @y = x, y
+  end
+
+  def ==(other)
+    eql?(other)
+  end
+
+  def eql?(other)
+    @x == other.x && @y == other.y
+  end
+
+  def hash
+    @x + @y
+  end
+
+  def to_s
+    "#{name}(#{x}, #{y})"
+  end
+end

data/lib/cf_sim/point_list.rb

@@ -0,0 +1,36 @@
+require 'forwardable'
+
+require 'cf_sim/control_field'
+require 'cf_sim/control_field_set'
+
+class CfSim::PointList
+  extend Forwardable
+  include Enumerable
+
+  attr_reader :points
+  def_delegators :@points, :clear, :each, :empty?, :hash, :size
+
+  def initialize(*points)
+    @points = points.flatten.freeze
+  end
+
+  def max_field_count
+    @max_field_count ||= 3 * (points.length - 3) + 1
+  end
+
+  def creatable_field_count
+    @creatable_field_count ||= points.combination(3).size
+  end
+
+  def creatable_fields
+    @creatable_fields ||= CfSim::ControlFieldSet.new(points.combination(3).map { |a, b, c| CfSim::ControlField.new(a, b, c) })
+  end
+
+  def eql?(other)
+    other.class == CfSim::PointList && @points == other.points
+  end
+
+  def ==(other)
+    eql?(other)
+  end
+end

data/lib/cf_sim/portal.rb

@@ -0,0 +1,25 @@
+class CfSim::Portal
+  attr_reader :name, :latitude, :longitude
+
+  def initialize(name, latitude, longitude)
+    @name = name
+    @latitude = latitude
+    @longitude = longitude
+  end
+
+  def ==(other)
+    eql?(other)
+  end
+
+  def eql?(other)
+    @name == other.name && @latitude == other.latitude && @longitude == other.longitude
+  end
+
+  def hash
+    @name.hash + @latitude + @longitude
+  end
+
+  def to_s
+    "#{name}(#{latitude}, #{longitude})"
+  end
+end

data/lib/cf_sim/portal_map.rb

@@ -0,0 +1,49 @@
+require 'cblxy'
+require 'cf_sim/cordination_system'
+require 'cf_sim/point'
+
+class CfSim::PortalMap
+
+  attr_reader :portals, :points
+
+  def initialize(*portals)
+    @portals = portals.flatten.dup.freeze
+    @cordination_system = CfSim::CordinationSystem.new(@portals.first.latitude, @portals.first.longitude)
+    @portal_map = {}
+    @point_map = {}
+    convert_to_points
+  end
+
+  def find_portal(point)
+    @portal_map[point]
+  end
+
+  def find_point(portal)
+    @point_map[portal]
+  end
+
+  private
+
+  def convert_to_points
+    @points = []
+    @portals.each do |portal|
+      point = convert_to_point(portal)
+      @points << point
+      @portal_map[point] = portal
+      @point_map[portal] = point
+    end
+    @points.freeze
+  end
+
+  def convert_to_point(portal)
+    x, y = blxy(convert_point_unit(portal.latitude), convert_point_unit(portal.longitude), @cordination_system.nearest_system_number)
+    CfSim::Point.new(portal.name, x, y)
+  end
+
+  # 10進数による座標を、60進数（度分秒）に変換し 10000 倍する
+  def convert_point_unit(point_unit)
+    degree, rest = point_unit.divmod(1)
+    minute, rest = (rest * 60).divmod(1)
+    degree * 10000 + minute * 100 + rest * 60
+  end
+end

data/lib/cf_sim/version.rb

@@ -0,0 +1,3 @@
+module CfSim
+  VERSION = "1.0.0"
+end

data/lib/cf_sim.rb

@@ -0,0 +1,19 @@
+require 'cblxy'
+require 'cf_sim/point'
+require 'cf_sim/portal'
+require 'cf_sim/cordination_system'
+require 'cf_sim/portal_map'
+require 'cf_sim/link'
+require 'cf_sim/control_field'
+require 'cf_sim/control_field_set'
+require 'cf_sim/point_list'
+require 'cf_sim/coexistable_field_map'
+require 'cf_sim/control_field_set_finder'
+require 'cf_sim/intel_map_url_generator'
+require 'cf_sim/data_file_parser'
+require 'cf_sim/version'
+require 'cf_sim/cli'
+
+module CfSim
+  # Your code goes here...
+end

metadata

@@ -0,0 +1,142 @@
+--- !ruby/object:Gem::Specification
+name: cf_sim
+version: !ruby/object:Gem::Version
+  version: 1.0.0
+platform: ruby
+authors:
+- pinzolo
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2015-10-31 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: test-unit
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: coveralls
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: thor
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: Propose multi layered control fields pattern by given portals.
+email:
+- pinzolo@gmail.com
+executables:
+- cf_sim
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".travis.yml"
+- CODE_OF_CONDUCT.md
+- Gemfile
+- LICENSE
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- cf_sim.gemspec
+- exe/cf_sim
+- lib/cblxy.rb
+- lib/cf_sim.rb
+- lib/cf_sim/cli.rb
+- lib/cf_sim/coexistable_field_map.rb
+- lib/cf_sim/control_field.rb
+- lib/cf_sim/control_field_set.rb
+- lib/cf_sim/control_field_set_finder.rb
+- lib/cf_sim/cordination_system.rb
+- lib/cf_sim/data_file_parser.rb
+- lib/cf_sim/intel_map_url_generator.rb
+- lib/cf_sim/link.rb
+- lib/cf_sim/point.rb
+- lib/cf_sim/point_list.rb
+- lib/cf_sim/portal.rb
+- lib/cf_sim/portal_map.rb
+- lib/cf_sim/version.rb
+homepage: https://github.com/pinzolo/cf_sim
+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.4.5
+signing_key: 
+specification_version: 4
+summary: Multi layered CF simulator
+test_files: []