+++ /dev/null
-#include "muscle.h"\r
-#include "msa.h"\r
-#include "profile.h"\r
-#include "pwpath.h"\r
-\r
-SCORE GlobalAlign4(ProfPos *PA, unsigned uLengthA, ProfPos *PB,\r
- unsigned uLengthB, PWPath &Path);\r
-\r
-SCORE AlignTwoProfs(\r
- const ProfPos *PA, unsigned uLengthA, WEIGHT wA,\r
- const ProfPos *PB, unsigned uLengthB, WEIGHT wB,\r
- PWPath &Path, ProfPos **ptrPout, unsigned *ptruLengthOut)\r
- {\r
- assert(uLengthA < 100000);\r
- assert(uLengthB < 100000);\r
-\r
- float r = (float) uLengthA/ (float) (uLengthB + 1); // +1 to prevent div 0\r
- if (r < 1)\r
- r = 1/r;\r
-\r
- SCORE Score = GlobalAlign(PA, uLengthA, PB, uLengthB, Path);\r
-\r
- AlignTwoProfsGivenPath(Path, PA, uLengthB, wA/(wA + wB), PB, uLengthB, wB/(wA + wB),\r
- ptrPout, ptruLengthOut);\r
-\r
-#if HYDRO\r
- if (ALPHA_Amino == g_Alpha)\r
- Hydro(*ptrPout, *ptruLengthOut);\r
-#endif\r
- return Score;\r
- }\r