INNER CODE UNIT · Python
n_read_pairs
HadrienG/InSilicoSeq · iss/app.py:81
n_read_pairs = n_reads // 2
chunk_size = -((n_read_pairs) // -args.cpus) # this is ceildiv, see https://stackoverflow.com/a/17511341
logger.debug("Chunk size: %s" % chunk_size)
# Divide the work of generating n_reads for each record into chunks
work_chunks = generate_work_divider(
fasta_file,
readcount_dic,
abundance_dic,
n_reads,
args.coverage,
args.coverage_file,
error_model,
args.output,
chunk_size,
)
# Generate reads for each chunk in parallel