void segBlockFreeList(struct segBlock **pList)
/* Free up a list of segment blocks. */
{
struct segBlock *el, *next;
for (el = *pList; el != NULL; el = next)
{
next = el->next;
segBlockFree(&el);
}
*pList = NULL;
}
void segSplitPairwise(char *inSeg, char *outSeg)
/* segSplitPairwise - Split a segment file pairwise. */
{
struct segFile *sf = segOpen(inSeg);
FILE *of = mustOpen(outSeg, "w");
struct segBlock *sb;
struct segBlock newSb;
struct segComp *refComp, *sc, *refSc, *newSc;
char *splitName = NULL;
int sNum, bNum = 1;
segWriteStart(of);
while ((sb = segNext(sf)) != NULL)
{
/* Set ref if it wasn't set on the commandline. */
if (ref == NULL && ((ref = segFirstCompSpecies(sb, '.')) == NULL))
errAbort("ref is not set and the first segment block has no "
"components.");
/* Find and clone the reference species component. */
refComp = segFindCompSpecies(sb, ref, '.');
refSc = cloneSegComp(refComp);
sNum = 1;
for (sc = sb->components; sc != NULL; sc = sc->next)
{
if (sc == refComp)
continue;
/* Clone the current component. */
newSc = cloneSegComp(sc);
newSc->next = NULL;
/* Point the reference component to the newly cloned component. */
refSc->next = newSc;
/* Set up the segment block to hold the two components above. */
if (number)
{
if (sb->name == NULL)
{
AllocArray(sb->name, 6 + digitsBaseTen(bNum));
sprintf(sb->name, "Block%d", bNum);
}
AllocArray(splitName, strlen(sb->name) + digitsBaseTen(sNum) + 2);
sprintf(splitName, "%s.%d", sb->name, sNum++);
newSb.name = splitName;
}
else
newSb.name = sb->name;
newSb.val = sb->val;
newSb.components = refSc;
/* Output this new block. */
segWrite(of, &newSb);
if (number)
freeMem(splitName);
/* Free the cloned component. */
segCompFree(&newSc);
}
/* Free the cloned reference component. */
segCompFree(&refSc);
segBlockFree(&sb);
bNum++;
}
segWriteEnd(of);
carefulClose(&of);
segFileFree(&sf);
}
