INNER CODE UNIT · C++
MinGlobalHSPScore
rcedgar/usearch12 · src/alnheuristics.cpp:40
MinGlobalHSPScore = MinGlobalHSPFractId*MinGlobalHSPLength*(float) oget_fltd(OPT_match, 1.0);
}
else
{
HSPFinderWordLength = 3;
// Avg BLOSUM62 score on the diagonal is 5.2, for comparison
const float * const *SubstMx = AP.SubstMx;
float MinDiagScore = 9e9f;
for (unsigned i = 0; i < 20; ++i)
{
byte c = g_LetterToCharAmino[i];
float Score = SubstMx[c][c];
if (Score < MinDiagScore)
MinDiagScore = Score;
}
MinGlobalHSPFractId = max((float) oget_fltd(OPT_id, 0.5f), 0.5f);