seis_ruby 0.1.2 → 0.2.0

data/lib/seis_ruby/data/sac.rb

@@ -1,6 +1,102 @@
-module SeisRuby
-  module Data
-    module Sac
+# @see http://www.iris.edu/software/sac/manual/file_format.html SAC User Manual/SAC Data Format
+class ::SeisRuby::Data::Sac < ::SeisRuby::Data
+  require 'seis_ruby/data/sac/head'
+  require 'seis_ruby/data/sac/body'
+  require 'seis_ruby/data/sac/binary'
+  require 'seis_ruby/data/sac/ascii'
+
+  def self.load_file(file)
+    raw_data = File.read_uri(file)
+    type = type_from_file(file)
+    meta_data = {file: file}
+
+    load(raw_data, type, meta_data)
+  end
+
+  def self.load(raw_data, type = :binary, meta_data = {})
+    new(raw_data, type, meta_data)
+  end
+
+  def self.type_from_file(file)
+    if ::SeisRuby::Data::Sac::Binary.uri_for_self?(file)
+      :binary
+    elsif ::SeisRuby::Data::Sac::Ascii.uri_for_self?(file)
+      :ascii
+    else
+      raise ArgumentError, "Unsupported extension: #{ext}"
+    end
+  end
+
+  def self.uri_for_self?(uri)
+    ::SeisRuby::Data::Sac::Binary.uri_for_self?(uri) || ::SeisRuby::Data::Sac::Ascii.uri_for_self?(uri)
+  end
+
+  def initialize(raw_data, type = :binary, meta_data = {})
+    @raw_data = raw_data
+    @meta_data = meta_data
+    parse!(type)
+  end
+
+  attr_accessor :head           # Header
+  attr_accessor :body           # Wave
+
+  ::SeisRuby::Data::Sac::Head::NAMES.each{|name|
+    define_method(name){@head[name]}
+    define_method("#{name}="){|val| @head[name] = val}}
+
+  def dump_file(file)
+    open(file, 'wb'){|io|
+      io.write(module_from_type(type_from_file(file)).dump(@head, @body))
+      io.flush
+    }
+  end
+
+  def dump(type = :binary)
+    module_from_type(type).dump(@head, @body)
+  end
+
+  def eql?(other)
+    self.class == other.class\
+    && @meta_data.eql?(other.meta_data)\
+    && @raw_data.eql?(other.raw_data)\
+    && @head.eql?(other.head)\
+    && @body.eql?(other.body)
+  end
+
+  def hash
+    [self.class, @meta_data, @raw_data, @head, @body].hash
+  end
+
+  def ==(other)
+    @head == other.head\
+    && @body == other.body
+  end
+
+  private
+
+  def type_from_file(file)
+    self.class.type_from_file(file)
+  end
+
+  def module_from_type(type)
+    case type
+    when :binary
+      ::SeisRuby::Data::Sac::Binary
+    when :ascii
+      ::SeisRuby::Data::Sac::Ascii
+    else
+      raise ArgumentError, "Unsupported type: #{type}"
     end
   end
+
+  def parse!(type)
+    @head, @body = module_from_type(type).parse(@raw_data)
+    check_consistency
+  end
+
+  private
+
+  def check_consistency
+
+  end
 end

data/lib/seis_ruby/data.rb

@@ -1,6 +1,7 @@
-module SeisRuby
-  module Data
-    require 'ruby_patch'
-    extend ::RubyPatch::AutoLoad
-  end
+class ::SeisRuby::Data
+  require 'seis_ruby/data/sac'
+  require 'seis_ruby/data/cmtsolution'
+
+  attr_accessor :meta_data      # Meta data about raw data.
+  attr_accessor :raw_data       # Raw data.
 end

data/lib/seis_ruby/database/global_cmt_catalog_search/cmtsolution_format.rb

@@ -0,0 +1,30 @@
+module ::SeisRuby::Database::GlobalCmtCatalogSearch::CmtsolutionFormat
+  require 'mechanize'
+
+  class << self
+    def load_file(uri)
+      agent = Mechanize.new{|a|
+        a.pluggable_parser.html = CustomHtmlParser
+      }
+      agent.get(uri)
+
+      events\
+      = agent.page.search('pre')\
+        .find{|e| e.text =~ /event name:/}\
+        .text.split("\n\n")\
+        .delete_if{|t| t.strip.empty?}\
+        .map{|text| ::SeisRuby::Data::Cmtsolution.load(text, {file: uri})}
+
+      unless next_link = agent.page.links.find{|elem| elem.text == 'More solutions'}
+        return events
+      end
+      events.concat(load_file(next_link.uri))
+    end
+  end
+
+  class CustomHtmlParser < ::Mechanize::Page
+    def initialize(uri = nil, response = nil, body = nil, code = nil, mech = nil)
+      super(uri, response, body.gsub(/<=/, '&lt;='), code, mech)
+    end
+  end
+end

data/lib/seis_ruby/database/global_cmt_catalog_search.rb

@@ -0,0 +1,44 @@
+module ::SeisRuby::Database::GlobalCmtCatalogSearch
+  require 'seis_ruby/database/global_cmt_catalog_search/cmtsolution_format'
+
+  # Examples of target URI
+  #   http://www.globalcmt.org/cgi-bin/globalcmt-cgi-bin/CMT4/form?itype=ymd&yr=1976&mo=1&day=1&oyr=1976&omo=1&oday=1&jyr=1976&jday=1&ojyr=1&ojday=1&otype=nd&nday=400&lmw=0&umw=10&lms=0&ums=10&lmb=0&umb=10&llat=-90&ulat=90&llon=-180&ulon=180&lhd=0&uhd=1000&lts=-9999&uts=9999&lpe1=0&upe1=90&lpe2=0&upe2=90&list=4
+  #   http://www.ldeo.columbia.edu/cgi-bin/globalcmt-cgi-bin/CMT4/form?itype=ymd&yr=1976&mo=1&day=1&oyr=1976&omo=1&oday=1&jyr=1976&jday=1&ojyr=1&ojday=1&otype=nd&nday=400&lmw=0&umw=10&lms=0&ums=10&lmb=0&umb=10&llat=-90&ulat=90&llon=-180&ulon=180&lhd=0&uhd=1000&lts=-9999&uts=9999&lpe1=0&upe1=90&lpe2=0&upe2=90&list=4&start=99
+  URI_FOR_SELF_REG = /\Ahttp:\/\/www\.(?:globalcmt\.org|ldeo\.columbia\.edu)\/cgi-bin\/globalcmt-cgi-bin\/CMT4\/form\?.*list=([0-5]+).*/
+
+  class << self
+    # @param [String] uri
+    # @return [Boolean]
+    def uri_for_self?(uri)
+      uri =~ URI_FOR_SELF_REG
+    end
+
+    # @param [String] uri URI of a first page of Global CMT Catalog Search result.
+    # @return [Array<Hash>] {::SeisRuby::Data::Cmtsolution#parse}
+    def load_file(uri)
+      raise ArgumentError, "Not a URI of #{self}: #{uri}" unless uri =~ URI_FOR_SELF_REG
+      module_from_list_number(Regexp.last_match(1).to_i).load_file(uri)
+    end
+
+    private
+
+    def module_from_list_number(list_number)
+      case list_number
+      when 0
+        raise NotImplementedError, "Standard"
+      when 1
+        raise NotImplementedError, "List of event names"
+      when 2
+        raise NotImplementedError, "GMT psvelomeca input"
+      when 3
+        raise NotImplementedError, "GMT psmeca input"
+      when 4
+        ::SeisRuby::Database::GlobalCmtCatalogSearch::CmtsolutionFormat
+      when 5
+        raise NotImplementedError, "Full format"
+      else
+        raise ArgumentError, "Unsupported list number: #{list_number}"
+      end
+    end
+  end
+end

data/lib/seis_ruby/database.rb

@@ -0,0 +1,3 @@
+module ::SeisRuby::Database
+  require 'seis_ruby/database/global_cmt_catalog_search'
+end

data/lib/seis_ruby/version.rb

@@ -1,3 +1,3 @@
 module SeisRuby
-  VERSION = '0.1.2'
+  VERSION = '0.2.0'
 end

data/lib/seis_ruby.rb

@@ -1,5 +1,14 @@
-module SeisRuby
+# @see ::SeisRuby::Command Module functions.
+module ::SeisRuby
+  AUTHOR = 'Kasahara Amato'
+  LICENSE = 'GNU GENERAL PUBLIC LICENSE Version 3'
+
+  require 'seis_ruby/core_ext'
   require 'seis_ruby/version'
-  require 'ruby_patch'
-  extend ::RubyPatch::AutoLoad
+  require 'seis_ruby/data'
+  require 'seis_ruby/database'
+  require 'seis_ruby/command'
+  require 'seis_ruby/application'
+
+  extend ::SeisRuby::Command
 end

data/seis_ruby.gemspec

@@ -1,5 +1,4 @@
 gem_name = File.basename(Dir.pwd)
-require 'rainbow'
 require 'ruby_patch'
 
 dir = File.join(__DIR__, 'lib')
@@ -11,29 +10,24 @@ Gem::Specification.new do |s|
   s.name = gem_name
   s.summary = "Ruby library for earthquake science."
   s.version = SeisRuby::VERSION
-  s.add_development_dependency 'rspec', '~> 2.10'
-  s.add_development_dependency 'simplecov', '~> 0.6'
-  s.add_runtime_dependency 'rainbow', '~> 1.1'
+  s.add_development_dependency 'watchr', '~> 0.7'
   s.add_runtime_dependency 'thor', '~> 0.15'
   s.add_runtime_dependency 'mechanize', '~> 2.5'
-  s.add_runtime_dependency 'ruby_patch', '~> 1.1'
+  s.add_development_dependency 'pry', '>= 0.9.10'
+  s.add_development_dependency 'ruby_patch', '>= 1.1.0'
   s.author = 'kshramt'
-  s.description = "Ruby library for earthquake science. Provides Parser for various data formats used in seismology."
+  s.description = "Ruby library for earthquake science."
   s.executables << 'seis_ruby'
-  s.post_install_message = (<<-EOS
+  s.post_install_message = <<-EOS
 
 
-SeisRuby #{::SeisRuby::VERSION}
-Thank you for installing.
 Please execute following command:
 
-#{::SeisRuby::Bin::SeisRubyRunner::COMMAND_NAME} generate_completion
-
-Enjoy!
+#{::SeisRuby::Application::COMMAND} generate_completion
 
 
   EOS
-  ).color(:green)
-  s.required_ruby_version = '~> 1.9'
+  s.required_ruby_version = '>= 1.9.3'
   s.test_files.concat `git ls-files spec`.split.select{|path| path =~ /_spec\.rb/}
+  s.test_files.concat `git ls-files test`.split.select{|path| path =~ /_test\.rb/}
 end

data/test/data/000000001.SAC_ascii

@@ -0,0 +1,32 @@
+     0.02500000      -12345.00      -12345.00   1.000000e+10      -12345.00
+   0.0005380000       725.6500      -292.0595      -12345.00       2.000000
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+      -12345.00       0.000000       90.00000       1000.000       2.000000
+       0.000000       0.000000      -12345.00       10000.00       7.000000
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+       10018.76       90.00000       270.0000       90.00000      -12345.00
+      -12345.00      -12345.00       0.000000       0.000000      -12345.00
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+      -12345.00      -12345.00      -12345.00      -12345.00      -12345.00
+      2012       100         1         2         3
+         4         6         0         0         9
+    -12345    -12345    -12345    -12345    -12345
+         1         5    -12345    -12345    -12345
+    -12345    -12345    -12345    -12345    -12345
+    -12345    -12345    -12345    -12345    -12345
+    -12345    -12345    -12345    -12345    -12345
+         1         0         0         1         0
+TEST    -12345          
+10      -12345  -12345  
+-12345  -12345  -12345  
+-12345  -12345  -12345  
+-12345  -12345  -12345  
+-12345  -12345  -12345  
+-12345  -12345  BHZ     
+XX      -12345  -12345  
+      -24370.00      -29946.00      -24963.00      -26414.00      -27634.00
+      -23068.00      -20123.00      -24239.00      -21714.00               

data/test/seis_ruby/command_test.rb

@@ -0,0 +1,19 @@
+require 'test_helper'
+
+class ::CommandTest < ::MiniTest::Unit::TestCase
+  SAC_ASCII_FILE = "test/data/000000001.SAC_ascii"
+  SAC_BINARY_FILE = "test/data/000000001.SAC"
+
+  def test_load_file_sac_binary
+    assert(::SeisRuby.load_file(SAC_BINARY_FILE))
+  end
+
+  def test_load_file_sac_ascii
+    assert(::SeisRuby.load_file(SAC_ASCII_FILE))
+  end
+
+  def test_load_file_global_cmt_catalog_search
+    assert(::SeisRuby.load_file('http://www.globalcmt.org/cgi-bin/globalcmt-cgi-bin/CMT4/form?itype=ymd&yr=1976&mo=1&day=1&oyr=1976&omo=1&oday=1&jyr=1976&jday=1&ojyr=1&ojday=1&otype=nd&nday=400&lmw=0&umw=10&lms=0&ums=10&lmb=0&umb=10&llat=-90&ulat=90&llon=-180&ulon=180&lhd=0&uhd=1000&lts=-9999&uts=9999&lpe1=0&upe1=90&lpe2=0&upe2=90&list=4'))
+    assert(::SeisRuby.load_file('http://www.ldeo.columbia.edu/cgi-bin/globalcmt-cgi-bin/CMT4/form?itype=ymd&yr=1976&mo=1&day=1&oyr=1976&omo=1&oday=1&jyr=1976&jday=1&ojyr=1&ojday=1&otype=nd&nday=400&lmw=0&umw=10&lms=0&ums=10&lmb=0&umb=10&llat=-90&ulat=90&llon=-180&ulon=180&lhd=0&uhd=1000&lts=-9999&uts=9999&lpe1=0&upe1=90&lpe2=0&upe2=90&list=4&start=99'))
+  end
+end

data/test/seis_ruby/core_ext/file_test.rb

@@ -0,0 +1,7 @@
+require 'test_helper'
+
+class FileTest_ < ::MiniTest::Unit::TestCase
+  def test_read_uri
+    assert(::File.read_uri('http://www.google.com'))
+  end
+end

data/test/seis_ruby/data/cmtsolution_test.rb

@@ -0,0 +1,63 @@
+require 'test_helper'
+
+class ::CmtsolutionTest < ::MiniTest::Unit::TestCase
+  RAW_DATA = <<-EOS
+
+   	
+
+ MLI 1976  1  2  3 29 39.60 -28.6100 -177.6400  59.0 6.2 0.0 KERMADEC ISLANDS REGION                 
+event name:     010176A        
+time shift:     13.8000
+half duration:   9.4000
+latitude:      -29.2500
+longitude:    -176.9600
+depth:          47.8000
+Mrr:       7.680000e+26
+Mtt:       9.000000e+24
+Mpp:      -7.770000e+26
+Mrt:       1.390000e+26
+Mrp:       4.520000e+26
+Mtp:      -3.260000e+26
+
+
+  EOS
+
+  def test_load
+    cmtsolution = ::SeisRuby::Data::Cmtsolution.load(RAW_DATA, {file: '(test)'})
+    hypocenter_table = {
+      data_source: 'MLI',
+      year: 1976,
+      month: 1,
+      day: 2,
+      hour: 3,
+      minute: 29,
+      second: 39.6,
+      latitude: -28.61,
+      longitude: -177.64,
+      depth: 59.0,
+      mb: 6.2,
+      ms: 0.0,
+      region_name: 'KERMADEC ISLANDS REGION',
+    }
+    centroid_table = {
+      event_name: '010176A',
+      time_shift: 13.8,
+      half_duration: 9.4,
+      latitude: -29.25,
+      longitude: -176.96,
+      depth: 47.8,
+      mrr: 7.68e+26,
+      mtt: 9.0e+24,
+      mpp: -7.77e+26,
+      mrt: 1.39e+26,
+      mrp: 4.52e+26,
+      mtp: -3.26e+26,
+    }
+    assert_equal(hypocenter_table, cmtsolution.hypocenter.table)
+    assert_equal(centroid_table, cmtsolution.centroid.table)
+  end
+
+  def test_time
+    assert_equal(Time.gm(1976, 1, 2, 3, 29, 39, 600000), ::SeisRuby::Data::Cmtsolution.load(RAW_DATA, {file: '(test)'}).hypocenter.origin_time)
+  end
+end

data/test/seis_ruby/data/sac/ascii_test.rb

@@ -0,0 +1,15 @@
+require 'test_helper'
+
+class ::AsciiTest < ::MiniTest::Unit::TestCase
+  RAW_DATA = File.read "test/data/000000001.SAC_ascii"
+
+  def test_parse
+    assert_equal(
+                 [Hash, Array],
+                 ::SeisRuby::Data::Sac::Ascii.parse(RAW_DATA).map(&:class))
+  end
+
+  def test_dump
+    assert(::SeisRuby::Data::Sac::Binary.dump({iftype: :itime}, []))
+  end
+end

data/test/seis_ruby/data/sac/binary_test.rb

@@ -0,0 +1,15 @@
+require 'test_helper'
+
+class ::BinaryTest < ::MiniTest::Unit::TestCase
+  RAW_DATA = File.read("test/data/000000001.SAC")
+
+  def test_parse
+    assert_equal(
+                 [Hash, Array],
+                 ::SeisRuby::Data::Sac::Binary.parse(RAW_DATA).map(&:class))
+  end
+
+  def test_dump
+    assert(::SeisRuby::Data::Sac::Binary.dump({iftype: :itime}, []))
+  end
+end

data/test/seis_ruby/data/sac/body_test.rb

@@ -0,0 +1,27 @@
+require 'test_helper'
+
+class ::BodyTest < ::MiniTest::Unit::TestCase
+  T = ::SeisRuby::Data::Sac::Body
+
+  def test_shape_body
+    [
+     [[1, 2, 3, 4], {iftype: :itime}, [1, 2, 3, 4]],
+     [[[3, 1], [4, 2]], {iftype: :ixy}, [1, 2, 3, 4]],
+     [[[3, 1], [4, 2]], {iftype: :iamph}, [1, 2, 3, 4]],
+     [[Complex(1, 3), Complex(2, 4)], {iftype: :irlim}, [1, 2, 3, 4]],
+     [[[1, 3, 5], [2, 4, 6]], {iftype: :ixyz, nxsize: 2}, [1, 2, 3, 4, 5, 6]],
+    ].each{|expect, head, body|
+      assert_equal(expect, T.shape_body(body, head))}
+  end
+
+  def test_body_for_dump
+    [
+     [[1, 2, 3, 4], {iftype: :itime}, [1, 2, 3, 4]],
+     [[1, 2, 3, 4], {iftype: :ixy}, [[3, 1], [4, 2]]],
+     [[1, 2, 3, 4], {iftype: :iamph}, [[3, 1], [4, 2]]],
+     [[1, 2, 3, 4], {iftype: :irlim}, [Complex(1, 3), Complex(2, 4)]],
+     [[1, 2, 3, 4, 5, 6], {iftype: :ixyz, nxsize: 2}, [[1, 3, 5], [2, 4, 6]]],
+    ].each{|expect, head, body|
+      assert_equal(expect, T.body_for_dump(body, head))}
+  end
+end

data/test/seis_ruby/data/sac/head_test.rb

@@ -0,0 +1,78 @@
+require 'test_helper'
+
+class HeadTest < ::MiniTest::Unit::TestCase
+  require 'ostruct'
+
+  #   undefined, valid, invalid
+  # :f -12345.0, 42.0, :sym
+  # :n -12345, 42, :sym
+  # :i -12445, :sym, (-1, 0, 87, 88)
+  # :l -12345, (1, 0), /
+  # :k "-12345", "42", /
+
+  T = ::SeisRuby::Data::Sac::Head
+
+  F_FROM = T::CONVERTER_FROM_HEAD_FLOAT
+  N_FROM = T::CONVERTER_FROM_HEAD_INTEGER
+  I_FROM = T::CONVERTER_FROM_HEAD_ENUMERATED_VALUE
+  L_FROM = T::CONVERTER_FROM_HEAD_LOGICAL
+  K_FROM = T::CONVERTER_FROM_HEAD_STRING
+
+  def test_convert_from_head_undefined
+    assert_nil(F_FROM.call(-12345.0))
+    assert_nil(N_FROM.call(-12345))
+    assert_nil(I_FROM.call(-12345))
+    assert_nil(K_FROM.call("-12345"))
+  end
+
+  def test_convert_from_head_valid
+    assert_equal(42.0, F_FROM.call(42.0))
+    assert_equal(42, N_FROM.call(42))
+    assert_equal(:itime, I_FROM.call(1))
+    assert_equal(:ipostq, I_FROM.call(42))
+    assert_equal(:iu, I_FROM.call(86))
+    assert_equal(true, L_FROM.call(1))
+    assert_equal(false, L_FROM.call(0))
+    assert_equal("42", K_FROM.call("42"))
+  end
+
+  def test_convert_from_head_invalid
+    assert_raises(ArgumentError, NoMethodError){F_FROM.call(:'1.0')}
+    assert_raises(ArgumentError, NoMethodError){N_FROM.call(:'1')}
+    assert_raises(ArgumentError, NoMethodError){I_FROM.call(-1)}
+    assert_raises(ArgumentError, NoMethodError){I_FROM.call(0)}
+    assert_raises(ArgumentError, NoMethodError){I_FROM.call(87)}
+    assert_raises(ArgumentError, NoMethodError){I_FROM.call(88)}
+  end
+
+  F_TO = T::CONVERTER_TO_HEAD_FLOAT
+  N_TO = T::CONVERTER_TO_HEAD_INTEGER
+  I_TO = T::CONVERTER_TO_HEAD_ENUMERATED_VALUE
+  L_TO = T::CONVERTER_TO_HEAD_LOGICAL
+  K_TO = T::CONVERTER_TO_HEAD_STRING
+
+  def test_convert_to_head_undefined
+    assert_equal(-12345.0, F_TO.call(nil))
+    assert_equal(-12345, N_TO.call(nil))
+    assert_equal(-12345, I_TO.call(nil))
+    assert_equal(0, L_TO.call(nil))
+    assert_equal("-12345", K_TO.call(nil))
+  end
+
+  def test_convert_to_head_valid
+    assert_equal(42.0, F_TO.call(42.0))
+    assert_equal(42, N_TO.call(42))
+    assert_equal(1, I_TO.call(:itime))
+    assert_equal(42, I_TO.call(:ipostq))
+    assert_equal(86, I_TO.call(:iu))
+    assert_equal(1, L_TO.call(true))
+    assert_equal(0, L_TO.call(false))
+    assert_equal("42", K_TO.call("42"))
+  end
+
+  def test_convert_to_head_invalid
+    assert_raises(ArgumentError, NoMethodError){F_TO.call(:"42.0")}
+    assert_raises(ArgumentError, NoMethodError){N_TO.call(:"42")}
+    assert_raises(ArgumentError, NoMethodError){I_TO.call(:this_is_not_exist)}
+  end
+end

data/test/seis_ruby/data/sac_test.rb

@@ -0,0 +1,52 @@
+require 'test_helper'
+
+class ::SacTest < ::MiniTest::Unit::TestCase
+  ASCII_FILE = "test/data/000000001.SAC_ascii"
+  BINARY_FILE = "test/data/000000001.SAC"
+  def test_new_binary
+    assert(::SeisRuby::Data::Sac.new(File.read(BINARY_FILE), :binary))
+  end
+
+  def test_new_ascii
+    assert(::SeisRuby::Data::Sac.new(File.read(ASCII_FILE), :ascii))
+  end
+
+  def test_load
+    assert(::SeisRuby::Data::Sac.load(File.read(BINARY_FILE), :binary))
+  end
+
+  def test_parse_binary
+    sac = ::SeisRuby::Data::Sac.load(File.read(BINARY_FILE), :binary)
+    assert(sac.head)
+    assert_equal(Hash, sac.head.class)
+    assert(sac.body)
+    assert_equal(sac.body.class, Array)
+  end
+
+  def test_load_file
+    assert(::SeisRuby::Data::Sac.load_file(BINARY_FILE))
+  end
+
+  def test_dump
+    sac = ::SeisRuby::Data::Sac.load_file(BINARY_FILE)
+    assert(sac.dump(:binary))
+  end
+
+  def test_dump_file
+    assert(::SeisRuby::Data::Sac.load_file(BINARY_FILE).dump_file(BINARY_FILE + '.sac'))
+  end
+
+  def test_binary_not_break_file
+    sac_binary_original = ::SeisRuby::Data::Sac.load_file(BINARY_FILE)
+    sac_binary_original.dump_file(BINARY_FILE + '.sac')
+    sac_binary_dumped = ::SeisRuby::Data::Sac.load_file(BINARY_FILE + '.sac')
+    assert_equal(sac_binary_original, sac_binary_dumped)
+  end
+
+  def test_ascii_not_break_file
+    sac_ascii_original = ::SeisRuby::Data::Sac.load_file(ASCII_FILE)
+    sac_ascii_original.dump_file(ASCII_FILE + '.sac_ascii')
+    sac_ascii_dumped = ::SeisRuby::Data::Sac.load_file(ASCII_FILE + '.sac_ascii')
+    assert_equal(sac_ascii_original, sac_ascii_dumped)
+  end
+end