INNER CODE UNIT · C++

size

refresh-bio/FAMSA · src/msa_refinement.cpp:23

	size_t size = profile_to_refine->data.front()->gapped_size;
	size_t card = profile_to_refine->data.size();

	dest_prof_id.clear();

	if (!valid_gap_stats)
		profile_to_refine->GetGapStats(gap_stats);

	vector<pair<size_t, size_t>> tmp;

	for (size_t i = 1; i <= size; ++i)
	{
		int x = (int)min(gap_stats[i], card - gap_stats[i]);
		if (x > 0)
			tmp.emplace_back(i, x);
	}

	stable_sort(tmp.begin(), tmp.end(), [](const pair<size_t, size_t>& x, const pair<size_t, size_t>& y) {

View source record →