vasputils 0.0.0

data/.document

@@ -0,0 +1,5 @@
+lib/**/*.rb
+bin/*
+- 
+features/**/*.feature
+LICENSE.txt

data/Gemfile

@@ -0,0 +1,18 @@
+source "http://rubygems.org"
+# Add dependencies required to use your gem here.
+# Example:
+#   gem "activesupport", ">= 2.3.5"
+
+# Add dependencies to develop your gem here.
+# Include everything needed to run rake, tests, features, etc.
+group :development do
+  gem "rdoc", "~> 3.12"
+  gem "bundler", "~> 1.1.3"
+  gem "jeweler", "~> 1.8.3"
+  gem "simplecov", ">= 0"
+  gem "crystalcell", ">= 0.0.0"
+  gem "mageo", ">= 0.0.0"
+  gem "malge", ">= 0.0.1"
+  gem "maset", ">= 0.0.0"
+  gem "builtinextension", ">= 0.0.3"
+end

data/LICENSE.txt

@@ -0,0 +1,20 @@
+Copyright (c) 2012 ippei94da
+
+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.rdoc

@@ -0,0 +1,19 @@
+= vasputils
+
+Description goes here.
+
+== Contributing to vasputils
+ 
+* Check out the latest master to make sure the feature hasn't been implemented or the bug hasn't been fixed yet.
+* Check out the issue tracker to make sure someone already hasn't requested it and/or contributed it.
+* Fork the project.
+* Start a feature/bugfix branch.
+* Commit and push until you are happy with your contribution.
+* Make sure to add tests for it. This is important so I don't break it in a future version unintentionally.
+* Please try not to mess with the Rakefile, version, or history. If you want to have your own version, or is otherwise necessary, that is fine, but please isolate to its own commit so I can cherry-pick around it.
+
+== Copyright
+
+Copyright (c) 2012 ippei94da. See LICENSE.txt for
+further details.
+

data/Rakefile

@@ -0,0 +1,54 @@
+# encoding: utf-8
+
+require 'rubygems'
+require 'bundler'
+begin
+  Bundler.setup(:default, :development)
+rescue Bundler::BundlerError => e
+  $stderr.puts e.message
+  $stderr.puts "Run `bundle install` to install missing gems"
+  exit e.status_code
+end
+require 'rake'
+
+require 'jeweler'
+Jeweler::Tasks.new do |gem|
+  # gem is a Gem::Specification... see http://docs.rubygems.org/read/chapter/20 for more options
+  gem.name = "vasputils"
+  gem.homepage = "http://github.com/ippei94da/vasputils"
+  gem.license = "MIT"
+  gem.summary = %Q{Utilities for VASP, first-principles calculation code.}
+  gem.description = %Q{This gem provides parsers for some VASP input and output files.
+    This will provide support command for computations.}
+  gem.email = "ippei94da@gmail.com"
+  gem.authors = ["ippei94da"]
+  # dependencies defined in Gemfile
+end
+Jeweler::RubygemsDotOrgTasks.new
+
+require 'rake/testtask'
+Rake::TestTask.new(:test) do |test|
+  test.libs << 'lib' << 'test'
+  test.pattern = 'test/**/test_*.rb'
+  test.verbose = true
+end
+
+#require 'rcov/rcovtask'
+#Rcov::RcovTask.new do |test|
+#  test.libs << 'test'
+#  test.pattern = 'test/**/test_*.rb'
+#  test.verbose = true
+#  test.rcov_opts << '--exclude "gems/*"'
+#end
+
+task :default => :test
+
+require 'rdoc/task'
+Rake::RDocTask.new do |rdoc|
+  version = File.exist?('VERSION') ? File.read('VERSION') : ""
+
+  rdoc.rdoc_dir = 'rdoc'
+  rdoc.title = "vasputils #{version}"
+  rdoc.rdoc_files.include('README*')
+  rdoc.rdoc_files.include('lib/**/*.rb')
+end

data/VERSION

@@ -0,0 +1 @@
+0.0.0

data/bin/addVolumeToten.rb

@@ -0,0 +1,19 @@
+#! /usr/bin/ruby
+
+#標準入力から文字列リストを受け取り、
+#ファイル名と思しきものがあれば、それに最終的な volume と TOTEN を追加する。
+#CONTCAR も指定できるが、同じディレクトリに OUTCAR が存在する必要がある。
+
+require "vasputils/outcarparser.rb"
+
+STDIN.each do |line|
+	file = line.strip.sub( 'CONTCAR', 'OUTCAR' )
+	if File.exist?( file )
+		outcar = ParseOutcar.new( file )
+		volume = outcar.volumes[-1]
+		toten = outcar.totens[-1]
+		printf( "%s, %6.2f, %f\n", line.chomp, volume, toten )
+	else
+		puts line
+	end
+end

data/bin/latticeconstants

@@ -0,0 +1,14 @@
+#! /usr/bin/env ruby
+
+require "vasputils/poscar.rb"
+require "crystal/latticeaxes.rb"
+require "crystal/cell2.rb"
+
+puts "       a,        b,        c,  alpha,   beta,  gamma,     volume, file"
+ARGV.each do |file|
+	axes = Poscar.load_file(file).axes
+	printf("%8.5f, %8.5f, %8.5f, %6.2f, %6.2f, %6.2f, ",
+		*(axes.get_lattice_constants))
+	printf("%10.5f, ", Cell2.new(axes).calc_volume) # show volume
+	puts file
+end

data/bin/lsvaspdir

@@ -0,0 +1,60 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+# vasp の計算ディレクトリの状態を表示する。
+# 
+#  - YET       未計算(lock がない)
+#  - STARTED   計算中 or 計算が終了したが、何らかのエラーが生じた。
+#              (lock があるが、normal_ended ではない)
+#              計算が中断されたか、計算中かの判断は原理上つけられない。
+#  - NEXT      計算が正常に終了し、収束していない。次の計算が生成された。
+#  - FINISHED  計算が正常に終了し、収束した。
+#  
+#  引数でディレクトリを指定できる。
+#  引数なしだとカレントディレクトリの全てのディレクトリ。
+#
+#  オプションをつけると、指定の状態のディレクトリを出力。
+#    -y : YET
+#    -s : STARTED
+#    -n : NEXT
+#    -f : FINISHED
+#   -c -n のように複数指定できる。
+#  
+#  デフォルトでは引数のディレクトリまたは
+#  -S オプションで状態を出力せず、ディレクトリ名だけを出力する。
+#    このオプションはシェルの backquote `` 機能から使うことを想定している。
+
+require "vasputils/calcinspector.rb"
+require "vasputils/vaspdir.rb"
+require "optparse"
+
+OPTIONS = {}
+OPTIONS[:target_status] = ["YET", "STARTED", "NEXT", "FINISHED"]
+tmp = []
+op = OptionParser.new
+op.on("-y", "--yet"      , "Show calcs of YET."     ){tmp << "YET"      }
+op.on("-s", "--started"  , "Show calcs of STARTED." ){tmp << "STARTED"  }
+op.on("-n", "--next"     , "Show calcs of NEXT."    ){tmp << "NEXT"     }
+op.on("-f", "--finished" , "Show calcs of FINISHED."){tmp << "FINISHED"}
+op.on("-S", "--no-status", "Not output status."     ){OPTIONS[:nostatus] = true}
+op.parse!(ARGV)
+OPTIONS[:target_status] = tmp unless tmp.empty?
+
+#pp OPTIONS
+
+dirs = ARGV
+dirs = Dir.glob("*").sort if ARGV.empty?
+dirs.each do |dir|
+	begin
+		vd = VaspDir.new(dir)
+		status = CalcInspector.inspect(vd)
+		#pp status
+		if OPTIONS[:target_status].include?(status)
+			printf("%10s ", status) unless OPTIONS[:nostatus]
+			printf("%s\n", dir)
+		end
+	rescue VaspDir::InitializeError
+		#pp "Something wrong."
+		next
+	end
+end

data/bin/lsvaspseries

@@ -0,0 +1,95 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+# vasp の計算シリーズの状態を表示する。
+#   - total ionic steps
+#   - total electronic steps
+#   - total elapsed time
+#   - final TOTEN
+#   - final cell parameter
+
+require "vasputils/calcseries.rb"
+require "vasputils/vaspdir.rb"
+require "optparse"
+
+## option analysis
+OPTIONS = Hash.new
+op = OptionParser.new
+op.on("-c", "--cutoff"          , "Show cutoff energy.    "){ OPTIONS[:cutoff ] = true }
+op.on("-e", "--electronic-steps", "Show electronic-steps." ){ OPTIONS[:e_steps] = true }
+op.on("-i", "--ionic-steps"     , "Show ionic-steps."      ){ OPTIONS[:i_steps] = true }
+op.on("-k", "--k-mesh"          , "Show k-mesh."           ){ OPTIONS[:k_mesh] = true }
+op.on("-K", "--irr-kpoints"     , "Show irreducible k-points."){ OPTIONS[:irr_k    ] = true }
+op.on("-l", "--lattice-const"   , "Show lattice constants."){ OPTIONS[:lattice] = true }
+op.on("-t", "--elapsed-time"    , "Show elapsed-time."     ){ OPTIONS[:time   ] = true }
+op.on("-T", "--time-estesp"     , "Show time per elec. steps."){ OPTIONS[:time_e] = true }
+op.on("-E", "--toten"           , "Show toten."            ){ OPTIONS[:toten  ] = true }
+op.on("-a", "--all"             , "Show all information."  ){
+	# -a で出力される順序はここで決まる。
+	OPTIONS[:cutoff    ] = true
+	OPTIONS[:irr_k     ] = true
+	OPTIONS[:k_mesh    ] = true
+	OPTIONS[:i_steps   ] = true
+	OPTIONS[:e_steps   ] = true
+	OPTIONS[:time      ] = true
+	OPTIONS[:time_e    ] = true
+	OPTIONS[:lattice   ] = true
+	OPTIONS[:toten     ] = true
+}
+
+op.parse!(ARGV)
+
+dirs = ARGV
+dirs = Dir.glob("*").sort if ARGV.empty?
+
+# generate headline
+items  = []
+OPTIONS.each do |key, val|
+	case key
+	when :cutoff    then items << sprintf("%6s", "ENCUT")
+	when :k_mesh    then items << sprintf("%2s, %2s, %2s, %8s, %8s, %8s",
+		"ka", "kb", "kc", "ka_shift", "kb_shift", "kc_shift")
+	when :irr_k     then items << sprintf("%4s", "ir_k")
+	when :lattice   then
+		items << sprintf("%9s, %9s, %9s, %9s, %9s, %9s",
+			"a", "b", "c", "alpha", "beta", "gamma")
+	when :e_steps   then items << sprintf("%4s", "el_s")
+	when :i_steps   then items << sprintf("%4s", "io_s")
+	when :time      then items << sprintf("%8s", "time")
+	when :time_e    then items << sprintf("%8s", "time/ele")
+	when :toten     then items << sprintf("%16s", "TOTEN")
+	end
+end
+items << "calc_name"
+puts "#" + items.join(", ")
+
+dirs.each do |dir|
+	begin
+		cs = CalcSeries.new(dir)
+		results = []
+		names = []
+		OPTIONS.each do |key, val|
+			case key
+			when :cutoff
+				results << sprintf("%6s", cs.finished_calc.incar["ENCUT"])
+			when :k_mesh    then results << sprintf("%2d, %2d, %2d, %8.5f, %8.5f, %8.5f",
+				* cs.finished_calc.kpoints[:mesh], * cs.finished_calc.kpoints[:shift])
+			when :irr_k     then results << sprintf("%4d", cs.finished_calc.outcar[:irreducible_kpoints])
+			when :lattice   then
+				results << sprintf("%9.5f, %9.5f, %9.5f, %9.5f, %9.5f, %9.5f",
+					* cs.finished_calc.contcar.axes.get_lattice_constants)
+			when :e_steps   then results << sprintf("%4d", cs.internal_steps)
+			when :i_steps   then results << sprintf("%4d", cs.external_steps)
+			when :time      then results << sprintf("%8d", cs.elapsed_time)
+			when :time_e    then results << sprintf("%8d", cs.elapsed_time/cs.internal_steps)
+			when :toten     then results << sprintf("%16.4f", cs.finished_calc.outcar[:totens][-1])
+			end
+		end
+		results << dir
+		puts " " + results.join(", ")
+	rescue
+		puts "Something wrong: #{dir}"
+		next
+	end
+
+end

data/bin/repeatvasp

@@ -0,0 +1,47 @@
+#! /usr/bin/ruby
+
+# Usage:
+#   repeatvasp calc_dir
+#   引数は計算ディレクトリで、1つのみ許可。
+#
+# 引数のディレクトリで計算を実行する。
+# 他所には投げず、localhost で実行する。
+#
+# 以下の条件を満たすとき、計算を繰り返す。
+# - cell shape, volume を含めた構造最適化で ionic step が 2以上。
+# - 構造最適化で ionic step が nsw と同じ。
+# - NEBM
+#
+# 計算の状態に応じたディレクトリ名をつける。
+# これらの名前は人間が分かり易くするためのものであって、
+# プログラムが計算の状態を把握するのはあくまで
+# 中のファイルを見て行うものとする。
+#
+# 状態を表す文字列をディレクトリ末尾に着ける機能はなし。
+# これをやると面倒。
+# この辺の確認は専用コマンドを用意することにする。
+#
+# 次の計算になるときに、
+# 元のディレクトリは末尾に 00 を付加する？
+# 次のディレクトリは末尾に 01 以降、インクリメント。
+# 01 からスタートしたときは、適当に 02 とか。
+#
+# 当面、計算ディレクトリは -try00 と末尾についていることを前提とする。
+
+require "vasputils/calcrepeater"
+require "vasputils/vaspdir"
+
+if ARGV.size != 1
+	puts "Number of arguments must be 1."
+	exit
+end
+
+vd = VaspDir.new(ARGV[0])
+cr = CalcRepeater.new
+begin
+	cr.repeat(vd)
+rescue VaspDir::PostfixMismatchError
+	puts "Postfix mismatches to #{VaspDir::POSTFIX}. Exit."
+rescue VaspDir::LockedError
+	puts "Lock file exist. Exit."
+end

data/bin/symposcar

@@ -0,0 +1,154 @@
+#!/usr/bin/env ruby
+#   Copyright (C) 2008 Atsushi Togo
+#   togo.atsushi@gmail.com
+#
+#   Little modified by Ippei Kishida. [2012-02-19]
+# 
+#   This program is free software; you can redistribute it and/or
+#   modify it under the terms of the GNU General Public License
+#   as published by the Free Software Foundation; either version 2
+#   of the License, or (at your option) any later version.
+#   
+#   This program is distributed in the hope that it will be useful,
+#   but WITHOUT ANY WARRANTY; without even the implied warranty of
+#   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+#   GNU General Public License for more details.
+#   
+#   You should have received a copy of the GNU General Public License
+#   along with this program; if not, write to
+#   the Free Software Foundation, Inc., 51 Franklin Street,
+#   Fifth Floor, Boston, MA 02110-1301, USA, or see
+#   http://www.gnu.org/copyleft/gpl.html
+#
+# Usage: symPoscar.rb [OPTION] [structure]
+# OPTION: -s, --symprec= : Symmetry check precision
+
+require 'optparse'
+require 'getspg.so'
+require 'poscar'
+require 'pp'
+include Getspg
+
+symprec = 1e-5
+nonewline = false
+pos_shift = [0,0,0]
+shift_string = false
+is_long_output = false
+is_operations = false
+is_dataset = false
+opt = OptionParser.new
+opt.on('-s', '--symprec=', 'Symmetry check precision') {|tmp| symprec = tmp.to_f}
+opt.on('-d', '--shift=', 'uniform shift of internal atomic positions') {|shift_string|}
+opt.on('-n', '--nonewline', 'Do not output the trailing newline') {nonewline = true}
+opt.on('-l', '--long', 'Long output') {is_long_output = true}
+opt.on('-o', '--operations', 'Symmetry operations') {is_operations = true}
+opt.on('-d', '--dataset', 'Dataset') {is_dataset = true}
+opt.parse!(ARGV)
+
+unless ARGV.size == 1
+	puts "Indicate only 1 POSCAR or CONTCAR. Exit."
+	exit
+end
+
+if shift_string
+  pos_shift = []
+  shift_string.split.each {|val| pos_shift << val.to_f }
+end
+cell = Vasp::Poscar.new(ARGV.shift).cell
+lattice = cell.axis.transpose
+names = (cell.atoms.collect {|atom| atom.name}).uniq
+position = []
+types = []
+names.each_with_index do |name, i|
+  cell.atoms.each do |atom|
+    if atom.name == name
+      apos = atom.position
+      position << [ apos[0]+pos_shift[0],
+                    apos[1]+pos_shift[1],
+                    apos[2]+pos_shift[2] ]
+      types << i+1
+    end
+  end
+end
+
+if is_long_output
+  spg, brv_lattice, brv_positions, brv_types, spgnum = 
+    refine_cell(types.size, lattice, position, types, symprec)
+end
+
+spgnum, spg, hall_symbol, t_mat, o_shift,
+rotations, translations, wyckoffs = get_dataset( lattice,
+                                                 position,
+                                                 types,
+                                                 symprec )
+
+if spgnum > 0
+  if nonewline
+    print "# #{spg.strip} (#{spgnum})"
+  else
+    puts "# #{spg.strip} (#{spgnum})"
+
+    if is_long_output
+      puts "----------- original -----------"
+      lattice.each do |vec|
+        printf("%10.5f %10.5f %10.5f\n", vec[0], vec[1], vec[2]);
+      end
+
+      puts "------------ final -------------"
+      brv_lattice.each do |vec|
+        printf("%10.5f %10.5f %10.5f\n", vec[0], vec[1], vec[2]);
+      end
+
+      brv_types.size.times do |i|
+        printf("%d: %d  %10.5f %10.5f %10.5f\n", i+1, brv_types[i], 
+               brv_positions[i][0], brv_positions[i][1], brv_positions[i][2]);
+      end
+    end
+
+    if is_dataset
+      # puts "------ transformation matrix -----"
+      # t_mat.each do |row|
+      #   printf("%10.5f %10.5f %10.5f\n", row[0], row[1], row[2]);
+      # end
+
+      # puts "---------- origin shift ----------"
+      # printf("%10.5f %10.5f %10.5f\n", o_shift[0], o_shift[1], o_shift[2]);
+      
+      puts "--------- Wyckoff position ----------"
+      wl = "abcdefghijklmnopqrstuvwxyz"
+      wyckoffs.each_with_index do |w, i|
+        pos = []
+        3.times do |j|
+          pos.push( position[i][j] - position[i][j].floor )
+        end
+        printf("%4d %2s  %s %8.5f %8.5f %8.5f\n",
+               i+1, cell.atoms[i].name, wl[w,1], pos[0], pos[1], pos[2])
+      end
+    end
+
+  end
+end
+
+#pp get_operations( types.size, lattice, position, types, symprec )
+
+if is_operations
+  rotations, translations = get_operations( types.size,
+                                            lattice,
+                                            position,
+                                            types,
+                                            symprec )
+end
+
+
+if is_operations
+  rotations.size.times do |i|
+    puts "- # ----#{i+1}----"
+    puts "  rotation:"
+    rotations[i].each do |row|
+      printf("    - [%2d, %2d, %2d]\n", row[0], row[1], row[2])
+    end
+    printf("  translation: [%f, %f, %f]\n",
+      translations[i][0], translations[i][1], translations[i][2])
+  end
+end
+

data/lib/vasputils/calcinspector.rb

@@ -0,0 +1,24 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+# 計算の状態を解析するクラス。
+class CalcInspector
+	#
+	#def initialize
+	#end
+
+	# 引数 calc で渡された計算の状態を返す。
+	# calc は以下のメソッドを持つ必要がある。
+	#   - started?
+	#   - normal_ended?
+	#   - to_be_continued?
+	def self.inspect(calc)
+		return "YET"       unless calc.started?
+		return "STARTED"   unless calc.normal_ended?
+		return "NEXT"      if     calc.to_be_continued?
+		return "FINISHED"
+		#raise "must not occur"
+	end
+
+end
+

data/lib/vasputils/calcrepeater.rb

@@ -0,0 +1,27 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+require "pp"
+#
+#
+#
+class CalcRepeater
+
+	#def initialize
+	#end
+
+	# 引数 calc は以下のメソッドを持っている必要がある。
+	def repeat(calc, io = STDOUT)
+		while true
+			io.puts "Calculation started: #{calc.name}."
+			calc.calculate
+			break unless calc.normal_ended?
+			break unless calc.to_be_continued?
+			calc.next
+			io.puts "Next to #{calc.name}."
+		end
+		calc.teardown
+		puts "Calculation got converged: #{calc.name}."
+	end
+end
+

data/lib/vasputils/calcseries.rb

@@ -0,0 +1,98 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+require "vasputils/vaspdir.rb"
+
+# 連続性のある一連の計算を表現するクラス。
+# 収束までに何度も繰り返すような計算を想定し、
+# これらが1つのディレクトリにまとめられていることを前提とする。
+# 計算は ASCII ソートできる順に
+# 名前が付けられていることを前提とする。
+class CalcSeries
+
+	class FinishedCalcError < Exception; end
+
+	#
+	def initialize(dir)
+		@dir = dir
+		@calculations = belonged_calculations
+	end
+
+	# Guess and return calculation type of an rgument 'file'.
+	# Argument 'file' can be a normal file or a directory.
+	# Return a calc class instance, e.g., VaspDir.
+	# If the file is unknown type, raise CalcSeriesUnknownCalcError 
+	#
+	# MEMO: Now, always return VaspDir.
+	def self.guess(file)
+		return VaspDir
+	end
+
+	# Return the final calculation that achieve to be convergence,
+	# whose state is FINISHED.
+	# If there are multiple FINISHED calcs,
+	# raise CalcSeriesMultipleFinishedError.
+	# If there is no FINISHED calcs, return nil.
+	def finished_calc
+		#tmp = @calculations.select{|calc| calc.normal_ended? && (! calc.to_be_continued?)}
+		tmp = @calculations.select{|calc| (! calc.to_be_continued?)}
+		#pp tmp.size
+		if tmp.size > 1
+			message = "Multiplicated finished calcs; #{tmp.map{|calc| calc.name}.join(', ')}"
+			raise FinishedCalcError, message
+		elsif tmp.size == 0
+			message = "No finished calc."
+			raise FinishedCalcError, message
+		else
+			return tmp[0]
+		end
+	end
+
+	# Return sum of number of 'internal steps',
+	# whish is returned by 'internal_steps' of each calculation.
+	# E.g., in case of VaspDir/try01, try02,...
+	#   internal_steps of VaspDir/try01 was 10, 
+	#   internal_steps of VaspDir/try02 was 3, 
+	def internal_steps
+		@calculations.inject(0.0) do |sum, item|
+			sum += item.internal_steps
+		end
+	end
+
+	# 外部ループ。
+	# vaspdir では ionic steps を返す筈。
+	def external_steps
+		@calculations.inject(0.0) do |sum, item|
+			sum += item.external_steps
+		end
+	end
+
+	# 全ての経過時間の合計を返す。
+	def elapsed_time
+		@calculations.inject(0.0) do |sum, item|
+			sum += item.elapsed_time
+		end
+	end
+
+	# 最後の計算が normal_ended ならば true。
+	# それ以外は false
+	def normal_ended?
+		begin
+			#pp finished_calc
+			return finished_calc.normal_ended?
+		rescue FinishedCalcError
+			return false
+		end
+	end
+
+	private
+
+	# 配下の calculation を VaspDir のようなインスタンスにし、
+	# それらをまとめた配列を返す。
+	def belonged_calculations
+		Dir.glob("#{@dir}/*").sort.map do |file|
+			self.class.guess(file).new(file)
+		end
+	end
+end
+