ba1e.c3 (690B) - raw

      1 // clump finding
      2 
      3 module ba1e;
      4 import std::io;
      5 import std::io::file;
      6 import std::collections;
      7 
      8 import util;
      9 
     10 fn int main(String[] args) {
     11 	if (args.len != 2) {
     12 		io::eprintn("Please supply path to data file");
     13 		return 1;
     14 	}
     15 	File f = file::open(args[1], "rb")!!;
     16 	String genome = io::treadline(&f)!!;
     17 	String ints = io::treadline(&f)!!;
     18 	String[] ints_split = ints.tsplit(" ");
     19 	int k = ints_split[0].to_integer(int)!!;
     20 	int region_len = ints_split[1].to_integer(int)!!; // L
     21 	int clump_threshold = ints_split[2].to_integer(int)!!; // t
     22 	String[] clumps = util::clump_finding(genome, k, region_len, clump_threshold);
     23 	foreach (clump : clumps) {
     24 		io::printf("%s ", clump);
     25 	}
     26 	return 0;
     27 }