source code
# The Computer Language Benchmarks Game
# https://salsa.debian.org/benchmarksgame-team/benchmarksgame/
#
# regex-dna program contributed by jose fco. gonzalez
# optimized & parallelized by Rick Branson
# optimized for ruby2 by Aaron Tavistock
# converted from regex-dna program
# array not dictionary by Isaac Gouy
seq = STDIN.readlines.join
ilen = seq.size
seq.gsub!(/>.*\n|\n/,"")
clen = seq.length
MATCHERS = [
/agggtaaa|tttaccct/,
/[cgt]gggtaaa|tttaccc[acg]/,
/a[act]ggtaaa|tttacc[agt]t/,
/ag[act]gtaaa|tttac[agt]ct/,
/agg[act]taaa|ttta[agt]cct/,
/aggg[acg]aaa|ttt[cgt]ccct/,
/agggt[cgt]aa|tt[acg]accct/,
/agggta[cgt]a|t[acg]taccct/,
/agggtaa[cgt]|[acg]ttaccct/
]
threads = MATCHERS.map do |f|
Thread.new do
Thread.current[:result] = "#{f.source} #{seq.scan(f).size}"
end
end
threads.each do |t|
t.join
end
match_results = threads.map do |t|
t[:result]
end
threads.each do |t|
t.join
end
# ruby 1.8.7: to iterate in-order use array not dictionary
[
[/tHa[Nt]/, '<4>'], [/aND|caN|Ha[DS]|WaS/, '<3>'], [/a[NSt]|BY/, '<2>'],
[/<[^>]*>/, '|'], [/\|[^|][^|]*\|/, '-']
].each { |f,r| seq.gsub!(f,r) }
puts "#{match_results.join("\n")}\n\n#{ilen}\n#{clen}\n#{seq.length}"
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
Sun, 17 May 2020 02:45:24 GMT
COMMAND LINE:
/usr/bin/ruby regexredux.mri-2.mri 0 < regexredux-input5000000.txt
PROGRAM OUTPUT:
agggtaaa|tttaccct 356
[cgt]gggtaaa|tttaccc[acg] 1250
a[act]ggtaaa|tttacc[agt]t 4252
ag[act]gtaaa|tttac[agt]ct 2894
agg[act]taaa|ttta[agt]cct 5435
aggg[acg]aaa|ttt[cgt]ccct 1537
agggt[cgt]aa|tt[acg]accct 1431
agggta[cgt]a|t[acg]taccct 1608
agggtaa[cgt]|[acg]ttaccct 2178
50833411
50000000
27388361