fasta benchmark ≈240MB N=25,000,000
Each chart bar shows how many times slower, one ↓ fasta program was, compared to the fastest program.
These are not the only programs that could be written. These are not the only compilers and interpreters. These are not the only programming languages.
Column × shows how many times more each program used compared to the benchmark program that used least.
| sort | sort | sort | ||||
| × | Program Source Code | CPU secs | Elapsed secs | Memory KB | Code B | ≈ CPU Load |
|---|---|---|---|---|---|---|
| 1.0 | Fortran Intel #4 | 2.70 | 2.70 | 244 | 1327 | 1% 0% 0% 100% |
| 1.5 | Haskell GHC #2 | 4.68 | 4.03 | 2,632 | 979 | 5% 86% 19% 5% |
| 1.5 | Ada 2005 GNAT #5 | 4.05 | 4.06 | 1,500 | 2186 | 1% 0% 100% 0% |
| 1.5 | C++ g++ #3 | 4.16 | 4.16 | 1,072 | 1286 | 0% 0% 100% 0% |
| 1.6 | C++ g++ #2 | 4.33 | 4.34 | 828 | 1105 | 0% 0% 100% 0% |
| 1.7 | Fortran Intel #3 | 4.46 | 4.46 | 508 | 1190 | 0% 0% 0% 100% |
| 1.7 | C++ g++ | 4.60 | 4.60 | 828 | 1033 | 0% 0% 100% 0% |
| 1.9 | Java 7 #4 | 5.21 | 5.14 | 18,284 | 1507 | 0% 1% 1% 100% |
| 1.9 | ATS | 5.20 | 5.20 | 364 | 2104 | 0% 100% 0% 0% |
| 2.2 | C gcc #4 | 5.96 | 5.97 | 372 | 1221 | 0% 0% 0% 100% |
| 2.2 | Fortran Intel | 5.97 | 5.97 | 508 | 1155 | 0% 0% 0% 100% |
| 2.3 | C++ g++ #4 | 6.20 | 6.21 | 372 | 1266 | 0% 0% 0% 100% |
| 2.3 | Scala #3 | 6.40 | 6.29 | 24,720 | 1053 | 0% 1% 1% 100% |
| 2.6 | Go | 6.94 | 6.94 | 1,272 | 1036 | 32% 67% 0% 0% |
| 2.7 | C gcc | 7.15 | 7.16 | 376 | 1185 | 0% 0% 100% 0% |
| 2.8 | Lisp SBCL #6 | 7.45 | 7.45 | 5,052 | 1751 | 0% 0% 100% 0% |
| 2.8 | C# Mono #2 | 7.47 | 7.48 | 14,748 | 1180 | 77% 0% 0% 23% |
| 3.0 | Java 7 #2 | 8.05 | 7.96 | 18,856 | 1240 | 1% 1% 0% 100% |
| 3.1 | F# Mono | 8.41 | 8.41 | 16,112 | 978 | 0% 0% 0% 100% |
| 3.2 | Ada 2005 GNAT | 8.73 | 8.74 | 1,496 | 1346 | 0% 0% 100% 0% |
| 3.7 | Lisp SBCL #3 | 9.91 | 9.92 | 4,788 | 1579 | 100% 0% 0% 0% |
| 3.7 | Dart | 9.96 | 9.97 | 291,008 | 1421 | 1% 0% 0% 100% |
| 3.8 | Scala | 10.43 | 10.31 | 27,832 | 1080 | 0% 0% 100% 1% |
| 3.9 | OCaml #6 | 10.50 | 10.52 | 198,060 | 1161 | 0% 0% 0% 100% |
| 3.9 | Clojure #5 | 11.91 | 10.59 | 63,312 | 1839 | 6% 68% 33% 6% |
| 4.1 | Pascal Free Pascal #4 | 11.17 | 11.18 | 8 | 1112 | 26% 0% 0% 74% |
| 4.9 | OCaml #3 | 13.31 | 13.32 | 2,848 | 1042 | 0% 0% 0% 100% |
| 6.4 | Lisp SBCL | 17.14 | 17.15 | 5,040 | 1419 | 0% 0% 100% 0% |
| 6.4 | Racket #3 | 17.32 | 17.32 | 24,100 | 1276 | 0% 0% 0% 100% |
| 6.5 | Lisp SBCL #2 | 17.48 | 17.49 | 5,044 | 1617 | 0% 0% 100% 0% |
| 11 | Erlang HiPE #2 | 29.43 | 29.43 | 14,164 | 1164 | 86% 0% 14% 0% |
| 18 | Racket | 49.56 | 49.57 | 18,708 | 1054 | 95% 5% 0% 0% |
| 22 | Erlang #2 | 60.56 | 60.58 | 12,284 | 1164 | 96% 0% 4% 0% |
| 24 | PHP #4 | 63.98 | 64.00 | 3,260 | 1110 | 0% 0% 100% 0% |
| 25 | Erlang HiPE | 66.89 | 66.90 | 13,536 | 1039 | 62% 38% 0% 0% |
| 35 | Erlang | 94.25 | 94.27 | 12,296 | 1039 | 100% 0% 0% 0% |
| 39 | Perl | 104.98 | 105.03 | 100,328 | 838 | 0% 100% 0% 0% |
| 40 | Ruby 2.0 #6 | 106.86 | 106.92 | 223,864 | 772 | 92% 0% 8% 0% |
| 53 | PHP #3 | 141.14 | 144.05 | 3,256 | 1030 | 91% 0% 1% 10% |
| 55 | Ruby 2.0 #4 | 145.07 | 149.58 | 223,872 | 904 | 30% 38% 14% 20% |
| 62 | Ruby 2.0 #5 | 162.95 | 166.65 | 6,160 | 987 | 38% 19% 2% 42% |
| 74 | Ruby 2.0 #2 | 198.71 | 198.80 | 250,028 | 732 | 97% 0% 3% 1% |
| 83 | Python 3 #2 | 219.23 | 222.96 | 6,208 | 788 | 31% 26% 3% 41% |
| 83 | Python 3 | 220.29 | 224.16 | 6,212 | 792 | 77% 1% 2% 21% |
| 84 | Ruby JRuby | 225.70 | 227.28 | 606,236 | 760 | 50% 48% 4% 2% |
| 107 | Perl #4 | 285.70 | 288.56 | 2,388 | 934 | 1% 1% 0% 99% |
| 142 | Perl #2 | 6 min | 6 min | 2,648 | 886 | 37% 63% 0% 0% |
| 146 | PHP #2 | 6 min | 6 min | 3,264 | 1006 | 41% 1% 7% 53% |
| Haskell GHC | Bad Output | 1421 | ||||
| Haskell GHC #4 | Bad Output | 1413 | ||||
| "wrong" (different) algorithm / less comparable programs | ||||||
| 0.8 | Perl #5 | 2.10 | 2.10 | 25,132 | 1113 | |
| 0.9 | C gcc #2 | 2.37 | 2.37 | 500 | 1169 | |
| 5.5 | Haskell GHC #3 | 16.69 | 14.83 | 2,112 | 1408 | |
fasta benchmark : Generate and write random DNA sequences
diff program output N = 1000 with this 10KB output file to check your program is correct before contributing.
We are trying to show the performance of various programming language implementations - so we ask that contributed programs not only give the correct result, but also use the same algorithm to calculate that result.
Each program should
- generate DNA sequences, by copying from a given sequence
- generate DNA sequences, by weighted random selection from 2 alphabets
- convert the expected probability of selecting each nucleotide into cumulative probabilities
- match a random number against those cumulative probabilities to select each nucleotide (use linear search or binary search)
- use this linear congruential generator to calculate a random number each time a nucleotide needs to be selected (don't cache the random number sequence)
IM = 139968 IA = 3877 IC = 29573 Seed = 42 Random (Max) Seed = (Seed * IA + IC) modulo IM = Max * Seed / IM
- write 3 sequences line-by-line in FASTA format
We'll use the generated FASTA file as input for other benchmarks (reverse-complement, k-nucleotide).
Random DNA sequences can be based on a variety of Random Models (554KB pdf). You can use Markov chains or independently distributed nucleotides to generate random DNA sequences online.