INNER CODE UNIT · C++

min_no_symbols

refresh-bio/FAMSA · src/core/io_service.cpp:116

		uint32_t min_no_symbols = static_cast<uint32_t>(no_seqs * rare_column_threshold);

		for (size_t i = 0; i < global_column_freq.size(); ++i)
			accepted_columns[i] = global_column_freq[i] >= min_no_symbols;
	}

	// *** Save sequences
	seq_id = 0;
	CLimitedPriorityQueue<vector<uint8_t>> lpq(no_threads, 5 * no_threads);

	for (int i = 0; i < no_threads; ++i)
		v_threads.emplace_back([pack_size, gzip_level, rare_column_threshold, &seq_id, &sequences, &lpq, &accepted_columns] {
		int no_seqs = static_cast<int>(sequences.size());
		string s_tmp;

		libdeflate_compressor* compressor = gzip_level >= 0 ? libdeflate_alloc_compressor(gzip_level) : nullptr;
		vector<uint8_t> gz_vec, raw_vec;

View source record →

📰 Research Paper
Loading…
⏳ Fetching content…