4 extern void TomHydro(ProfPos *Prof, unsigned Length);
\r
6 // Apply hydrophobicity heuristic to a profile
\r
7 void Hydro(ProfPos *Prof, unsigned uLength)
\r
9 if (ALPHA_Amino != g_Alpha)
\r
14 TomHydro(Prof, uLength);
\r
18 if (0 == g_uHydrophobicRunLength)
\r
21 if (uLength <= g_uHydrophobicRunLength)
\r
24 unsigned uRunLength = 0;
\r
25 unsigned L2 = g_uHydrophobicRunLength/2;
\r
26 for (unsigned uColIndex = L2; uColIndex < uLength - L2; ++uColIndex)
\r
28 ProfPos &PP = Prof[uColIndex];
\r
29 bool bHydro = IsHydrophobic(PP.m_fcCounts);
\r
33 if (uRunLength >= g_uHydrophobicRunLength)
\r
35 Prof[uColIndex-L2].m_scoreGapOpen *= (SCORE) g_dHydroFactor;
\r
36 Prof[uColIndex-L2].m_scoreGapClose *= (SCORE) g_dHydroFactor;
\r