+++ /dev/null
-#include "muscle.h"\r
-#include "profile.h"\r
-\r
-extern void TomHydro(ProfPos *Prof, unsigned Length);\r
-\r
-// Apply hydrophobicity heuristic to a profile\r
-void Hydro(ProfPos *Prof, unsigned uLength)\r
- {\r
- if (ALPHA_Amino != g_Alpha)\r
- return;\r
-\r
- if (g_bTomHydro)\r
- {\r
- TomHydro(Prof, uLength);\r
- return;\r
- }\r
-\r
- if (0 == g_uHydrophobicRunLength)\r
- return;\r
-\r
- if (uLength <= g_uHydrophobicRunLength)\r
- return;\r
-\r
- unsigned uRunLength = 0;\r
- unsigned L2 = g_uHydrophobicRunLength/2;\r
- for (unsigned uColIndex = L2; uColIndex < uLength - L2; ++uColIndex)\r
- {\r
- ProfPos &PP = Prof[uColIndex];\r
- bool bHydro = IsHydrophobic(PP.m_fcCounts);\r
- if (bHydro)\r
- {\r
- ++uRunLength;\r
- if (uRunLength >= g_uHydrophobicRunLength)\r
- {\r
- Prof[uColIndex-L2].m_scoreGapOpen *= (SCORE) g_dHydroFactor;\r
- Prof[uColIndex-L2].m_scoreGapClose *= (SCORE) g_dHydroFactor;\r
- }\r
- }\r
- else\r
- uRunLength = 0;\r
- }\r
- }\r