1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
|
#include "libMUSCLE/muscle.h"
#include "libMUSCLE/msa.h"
#include "libMUSCLE/tree.h"
#include "libMUSCLE/seqvect.h"
namespace muscle {
#if DEBUG
static TLS<SeqVect *> g_ptrMuscleSeqVect(0);
static TLS<MSA> MuscleInputMSA;
void SetMuscleInputMSA(MSA &msa)
{
MuscleInputMSA.get().Copy(msa);
}
void SetMuscleSeqVect(SeqVect &v)
{
g_ptrMuscleSeqVect.get() = &v;
}
void ValidateMuscleIdsSeqVect(const MSA &msa)
{
const unsigned uSeqCount = msa.GetSeqCount();
for (unsigned uSeqIndex = 0; uSeqIndex < uSeqCount; ++uSeqIndex)
{
const unsigned uId = msa.GetSeqId(uSeqIndex);
const char *ptrNameMSA = msa.GetSeqName(uSeqIndex);
const char *ptrName = g_ptrMuscleSeqVect.get()->GetSeqName(uId);
if (0 != strcmp(ptrNameMSA, ptrName))
Quit("ValidateMuscleIdsSeqVect, names don't match");
}
}
void ValidateMuscleIdsMSA(const MSA &msa)
{
const unsigned uSeqCount = msa.GetSeqCount();
for (unsigned uSeqIndex = 0; uSeqIndex < uSeqCount; ++uSeqIndex)
{
const unsigned uId = msa.GetSeqId(uSeqIndex);
const char *ptrNameMSA = msa.GetSeqName(uSeqIndex);
const char *ptrName = MuscleInputMSA.get().GetSeqName(uId);
if (0 != strcmp(ptrNameMSA, ptrName))
{
Log("Input MSA:\n");
MuscleInputMSA.get().LogMe();
Log("MSA being tested:\n");
msa.LogMe();
Log("Id=%u\n", uId);
Log("Input name=%s\n", ptrName);
Log("Test name=%s\n", ptrNameMSA);
Quit("ValidateMuscleIdsMSA, names don't match");
}
}
}
void ValidateMuscleIds(const MSA &msa)
{
if (0 != g_ptrMuscleSeqVect.get())
ValidateMuscleIdsSeqVect(msa);
else if (0 != MuscleInputMSA.get().GetSeqCount())
ValidateMuscleIdsMSA(msa);
else
Quit("ValidateMuscleIds, ptrMuscleSeqVect=0 && 0 == MuscleInputMSA.get().SeqCount()");
}
void ValidateMuscleIdsSeqVect(const Tree &tree)
{
const unsigned uSeqCount = g_ptrMuscleSeqVect->GetSeqCount();
const unsigned uNodeCount = tree.GetNodeCount();
for (unsigned uNodeIndex = 0; uNodeIndex < uNodeCount; ++uNodeIndex)
{
if (!tree.IsLeaf(uNodeIndex))
continue;
const unsigned uId = tree.GetLeafId(uNodeIndex);
if (uId >= uSeqCount)
{
tree.LogMe();
Quit("Leaf with node index %u has id=%u, there are %u seqs",
uNodeIndex, uId, uSeqCount);
}
const char *ptrNameTree = tree.GetLeafName(uNodeIndex);
const char *ptrName = g_ptrMuscleSeqVect.get()->GetSeqName(uId);
if (0 != strcmp(ptrNameTree, ptrName))
Quit("ValidateMuscleIds: names don't match");
}
}
void ValidateMuscleIdsMSA(const Tree &tree)
{
const unsigned uNodeCount = tree.GetNodeCount();
for (unsigned uNodeIndex = 0; uNodeIndex < uNodeCount; ++uNodeIndex)
{
if (!tree.IsLeaf(uNodeIndex))
continue;
const unsigned uId = tree.GetLeafId(uNodeIndex);
const char *ptrNameTree = tree.GetLeafName(uNodeIndex);
const char *ptrName = MuscleInputMSA.get().GetSeqName(uId);
if (0 != strcmp(ptrNameTree, ptrName))
Quit("ValidateMuscleIds: names don't match");
}
}
void ValidateMuscleIds(const Tree &tree)
{
if (0 != g_ptrMuscleSeqVect.get())
ValidateMuscleIdsSeqVect(tree);
else if (0 != MuscleInputMSA.get().GetSeqCount())
ValidateMuscleIdsMSA(tree);
else
Quit("ValidateMuscleIds, ptrMuscleSeqVect=0 && 0 == MuscleInputMSA.get().SeqCount");
}
#endif
}
|
