func benchmarkAnno(b *testing.B) {
var configs shared.Config
if _, err := toml.DecodeFile("example/conf.toml", &configs); err != nil {
panic(err)
}
out := bufio.NewWriter(ioutil.Discard)
Js, _ := xopen.Ropen("example/custom.js")
jbytes, _ := ioutil.ReadAll(Js)
js_string := string(jbytes)
srcs, err := configs.Sources()
if err != nil {
log.Fatal(err)
}
a := api.NewAnnotator(srcs, js_string, false, true)
files, err := a.SetupStreams()
if err != nil {
log.Fatal(err)
}
for n := 0; n < b.N; n++ {
q, err := bix.New("example/query.vcf.gz", 1)
if err != nil {
log.Fatal(err)
}
bx, err := q.Query(nil)
if err != nil {
log.Fatal(err)
}
a.UpdateHeader(q)
//files := []string{"example/cadd.sub.txt.gz", "example/query.vcf.gz", "example/fitcons.bed.gz"}
//files := []string{"example/query.vcf.gz", "example/fitcons.bed.gz"}
queryables := make([]interfaces.Queryable, len(files))
for i, f := range files {
q, err := bix.New(f, 1)
if err != nil {
log.Fatal(err)
}
queryables[i] = q
}
fn := func(v interfaces.Relatable) {
a.AnnotateEnds(v, api.INTERVAL)
}
stream := irelate.PIRelate(4000, 20000, bx, false, fn, queryables...)
for interval := range stream {
fmt.Fprintf(out, "%s\n", interval)
}
out.Flush()
}
}
func benchmarkAnno(b *testing.B) {
var configs shared.Config
if _, err := toml.DecodeFile("example/conf.toml", &configs); err != nil {
panic(err)
}
out := bufio.NewWriter(ioutil.Discard)
Lua, _ := xopen.Ropen("example/custom.lua")
lbytes, _ := ioutil.ReadAll(Lua)
l_string := string(lbytes)
srcs, err := configs.Sources()
if err != nil {
log.Fatal(err)
}
empty := make([]api.PostAnnotation, 0)
for n := 0; n < b.N; n++ {
if err != nil {
log.Fatal(err)
}
a := api.NewAnnotator(srcs, l_string, false, true, empty)
qrdr, err := xopen.Ropen("example/query.vcf.gz")
if err != nil {
log.Fatal(err)
}
qstream, query, err := parsers.VCFIterator(qrdr)
queryables, err := a.Setup(query)
if err != nil {
log.Fatal(err)
}
//files := []string{"example/cadd.sub.txt.gz", "example/query.vcf.gz", "example/fitcons.bed.gz"}
//files := []string{"example/query.vcf.gz", "example/fitcons.bed.gz"}
fn := func(v interfaces.Relatable) {
a.AnnotateEnds(v, api.INTERVAL)
}
stream := irelate.PIRelate(6000, 20000, qstream, false, fn, queryables...)
for interval := range stream {
fmt.Fprintf(out, "%s\n", interval)
}
out.Flush()
}
}