reverse-complement Matz's Ruby program
source code
# The Computer Language Benchmarks Game
# https://salsa.debian.org/benchmarksgame-team/benchmarksgame/
#
# Contributed by Aaron Tavistock
def process_segment(segment)
(header, *body, tail) = segment.split("\n")
# This is an ugly way to handle the first row, but returning nil causes the fork to implode
return '' if !tail
body << tail if tail != '>'
sequence = body.join
sequence.reverse!
sequence.tr!(
'wsatugcyrkmbdhvnATUGCYRKMBDHVN',
'WSTAACGRYMKVHDBNTAACGRYMKVHDBN'
)
results = [">#{header}"]
idx = 0
length = sequence.length
while idx < length
results << sequence[idx,60]
idx += 60
end
results.join("\n")
end
def forking_worker(segment)
reader, writer = IO.pipe
pid = Process.fork do
begin
reader.close
results = process_segment(segment)
Marshal.dump(results, writer) if results
ensure
writer.close
end
end
writer.close
begin
results = Marshal.load(reader)
ensure
reader.close
end
Process.waitpid(pid)
results
end
threads = []
$stdin.each_line('>') do |segment|
threads << Thread.new do
if (RUBY_PLATFORM == 'java')
Thread.current[:output] = process_segment(segment)
else
Thread.current[:output] = forking_worker(segment)
end
end
end
threads.each(&:join)
threads.each do |thread|
if !thread[:output].empty? # Artifact of non-empty results from forking worker
puts thread[:output]
end
end
notes, command-line, and program output
NOTES:
64-bit Ubuntu quad core
ruby 1.8.7 (2008-08-11 patchlevel 72) [x86_64-linux]
So old that I haven't been able to get rubygems to work
no backport, no gmp
Sat, 04 Jan 2020 18:52:31 GMT
COMMAND LINE:
/usr/bin/ruby revcomp.mri 0 < revcomp-input250000.txt
PROGRAM FAILED
PROGRAM OUTPUT:
revcomp.mri:7: syntax error, unexpected ',', expecting '='
(header, *body, tail) = segment.split("\n")
^