void inputoptions()
{
/* read the options information */
Char ch;
long i, extranum, cursp, curst;
if (!firstset) {
if (eoln(infile)) {
fscanf(infile, "%*[^\n]");
getc(infile);
}
fscanf(infile, "%ld%ld", &cursp, &curst);
if (cursp != spp) {
printf("\nERROR: INCONSISTENT NUMBER OF SPECIES IN DATA SET %4ld\n",
ith);
exxit(-1);
}
sites = curst;
}
for (i = 1; i <= sites; i++)
weight[i] = 1;
weightsum = sites;
extranum = 0;
fscanf(infile, "%*[ 0-9]");
readoptions(&extranum, "W");
for (i = 1; i <= extranum; i++) {
matchoptions(&ch, "W");
inputweights2(1, sites+1, &weightsum, weight, &weights, "RESTDIST");
}
} /* inputoptions */
void inputoptions()
{
/* input the information on the options */
long weightsum, maxfactsize, i, j, k, l, m;
if (data == genefreqs) {
k = 0;
l = 0;
for (i = 0; i < (loci); i++) {
if (all)
m = alleles[i];
else
m = alleles[i] - 1;
k++;
for (j = 1; j <= m; j++) {
l++;
factorr[l - 1] = k;
}
}
} else {
for (i = 1; i <= (sites); i++)
factorr[i - 1] = i;
}
if(factors){
seqboot_inputfactors();
}
for (i = 0; i < (sites); i++)
oldweight[i] = 1;
if (weights)
inputweights2(0, sites, &weightsum, oldweight, &weights, "seqboot");
if (factors && printdata) {
for(i = 0; i < sites; i++)
factor[i] = (char)('0' + (factorr[i]%10));
printfactors(outfile, sites, factor, " (least significant digit)");
}
if (weights && printdata)
printweights(outfile, 0, sites, oldweight, "Sites");
for (i = 0; i < (loci); i++)
how_many[i] = 0;
for (i = 0; i < (loci); i++)
where[i] = 0;
for (i = 1; i <= (sites); i++) {
how_many[factorr[i - 1] - 1]++;
if (where[factorr[i - 1] - 1] == 0)
where[factorr[i - 1] - 1] = i;
}
groups = factorr[sites - 1];
newgroups = 0;
newsites = 0;
maxfactsize = 0;
for(i = 0 ; i < loci ; i++){
if(how_many[i] > maxfactsize){
maxfactsize = how_many[i];
}
}
maxnewsites = groups * maxfactsize;
allocnew();
for (i = 0; i < groups; i++) {
if (oldweight[where[i] - 1] > 0) {
newgroups++;
newsites += how_many[i];
newwhere[newgroups - 1] = where[i];
newhowmany[newgroups - 1] = how_many[i];
}
}
} /* inputoptions */
