void
SplitQuery_SetEffectiveSearchSpace(CRef<CBlastOptions> options,
CRef<IQueryFactory> full_query_fact,
CRef<SInternalData> full_data)
{
_ASSERT(full_data);
_ASSERT(full_data->m_SeqSrc);
// If the effective search options have been set, we don't need to
// recompute those...
if (options->GetEffectiveSearchSpace() != 0) {
return;
}
const BlastSeqSrc* seqsrc = full_data->m_SeqSrc->GetPointer();
Int8 total_length = BlastSeqSrcGetTotLenStats(seqsrc);
if (total_length <= 0)
total_length = BlastSeqSrcGetTotLen(seqsrc);
Int4 num_seqs = BlastSeqSrcGetNumSeqsStats(seqsrc);
if (num_seqs <= 0)
num_seqs = BlastSeqSrcGetNumSeqs(seqsrc);
CEffectiveSearchSpaceCalculator calc(full_query_fact, *options,
num_seqs, total_length,
full_data->m_ScoreBlk->GetPointer());
BlastQueryInfo* qinfo = full_data->m_QueryInfo;
_ASSERT(qinfo);
vector<Int8> eff_searchsp;
for (size_t index = 0; index <= (size_t)qinfo->last_context; index++) {
eff_searchsp.push_back(calc.GetEffSearchSpaceForContext(index));
}
options->SetEffectiveSearchSpace(eff_searchsp);
}
