The Computer Language
Benchmarks Game

regex-redux C# .NET Core #8 program

source code

/* The Computer Language Benchmarks Game
   http://benchmarksgame.alioth.debian.org/
 *
 * contributed by Jimmy Tang
 * modified by Josh Goldfoot (2016)
 * modified by Jan de Vaan (compile regex, small stuff)
 */
using System;
using System.Text;
using System.Text.RegularExpressions;
using System.Threading.Tasks;

class regexredux
{
    static string readStdIn(out int seqLength, out int inputLength)
    {
        StringBuilder sb = new StringBuilder();
        int commentLength = 0;
        String line;
        
        while ((line = Console.ReadLine()) != null)
        {
            if (line[0] == '>')
                commentLength += line.Length + 1;
            else
            {
                sb.Append(line);
                commentLength += 1;
            }
        }
        seqLength = sb.Length;
        inputLength = seqLength + commentLength; 
        return sb.ToString();
    }

    static void Main()
    {

        string[] variants = {
           "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"
        };
         
        int seqLength, initialLength;
        var sequence = readStdIn(out seqLength, out initialLength);
        var newSequenceLength = Task.Run(() =>
            {
                var table = new int['Z'];
                table['D'] = "(a|g|t)".Length - 1;
                table['H'] = "(a|c|t)".Length - 1;
                table['K'] = "(g|t)".Length - 1;
                table['M'] = "(a|c)".Length - 1;
                table['N'] = "(a|c|g|t)".Length - 1;
                table['R'] = "(a|g)".Length - 1;
                table['S'] = "(c|g)".Length - 1;
                table['V'] = "(a|c|g)".Length - 1;
                table['W'] = "(a|t)".Length - 1;
                table['Y'] = "(c|t)".Length - 1;
                table['B'] = "(c|g|t)".Length - 1;

                var r = new Regex("[WYKMSRBDVHN]", RegexOptions.Compiled);

                int length = sequence.Length;

                for (Match m = r.Match(sequence); m.Success; m = m.NextMatch())
                {
                    length += table[m.Value[0]];
                }
                
                return length;
            });

        var output = new string[variants.Length];
        Parallel.For(0, variants.Length, i =>
        {
            Regex r = new Regex(variants[i], RegexOptions.Compiled);            
            output[i] = r.ToString() + " " + r.Matches(sequence).Count;
        });

        foreach (var s in output)
            Console.WriteLine(s);
        
        Console.WriteLine("\n{0}\n{1}\n{2}", initialLength, seqLength, newSequenceLength.Result);        
    }
}
    

notes, command-line, and program output

NOTES:
64-bit Ubuntu quad core
2.0.0 cdcd1928c9
"System.GC.Server": true


Mon, 14 Aug 2017 20:23:06 GMT

MAKE:
cp regexredux.csharpcore-8.csharpcore Program.cs
cp Include/csharpcore/tmp.csproj .
cp Include/csharpcore/runtimeconfig.template.json .
mkdir obj
cp Include/csharpcore/tmp.csproj.nuget.g.props ./obj
cp Include/csharpcore/tmp.csproj.nuget.g.targets ./obj
/usr/bin/dotnet build -c Release
Microsoft (R) Build Engine version 15.3.409.57025 for .NET Core
Copyright (C) Microsoft Corporation. All rights reserved.

  tmp -> /home/dunham/benchmarksgame_quadcore/regexredux/tmp/bin/Release/netcoreapp2.0/tmp.dll

Build succeeded.
    0 Warning(s)
    0 Error(s)

Time Elapsed 00:00:03.61

6.32s to complete and log all make actions

COMMAND LINE:
/usr/bin/dotnet ./bin/Release/netcoreapp2.0/tmp.dll 0 < regexredux-input50000.txt

UNEXPECTED OUTPUT 

13c13
< 668262
---
> 273927

PROGRAM OUTPUT:
agggtaaa|tttaccct 3
[cgt]gggtaaa|tttaccc[acg] 12
a[act]ggtaaa|tttacc[agt]t 43
ag[act]gtaaa|tttac[agt]ct 27
agg[act]taaa|ttta[agt]cct 58
aggg[acg]aaa|ttt[cgt]ccct 16
agggt[cgt]aa|tt[acg]accct 15
agggta[cgt]a|t[acg]taccct 18
agggtaa[cgt]|[acg]ttaccct 20

508411
500000
668262