FilesExpand file tree

master

/

threaded_consume.py

Copy path

Blame

More file actions

Blame

More file actions

Latest commit

 

History

History

History

30 lines (24 loc) · 622 Bytes

master

/

threaded_consume.py

Copy path

Top

File metadata and controls

• Code
• Blame

30 lines (24 loc) · 622 Bytes

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

from threading import Thread

import time

import khmer

class ConsumeThread(Thread):

def __init__(self, wordsize, genome):

self.wordsize = wordsize

self.genome = genome

Thread.__init__(self)

def run(self):

self.kt = khmer.consume_genome(self.wordsize, self.genome)

genome = open('/tmp/all2.dna').read()

nthreads = 2

length = len(genome)

genome1 = genome[:length/2]

genome2 = genome[length/2:]

t1 = ConsumeThread(5, genome1)

t2 = ConsumeThread(5, genome2)

t1.start()

t2.start()

t1.join()

t2.join()

master_kt = khmer.new_ktable(5)

master_kt.update(t1.kt)

master_kt.update(t2.kt)