void
work(void){
int c;
struct keyw *lptr;
char *pt;
int cc;
int ct;
while ((c = getch()) != Beof){
switch (c){
case '{':
if ((lptr = lookup(lastlook,p)) != 0){
if (lptr->type == ELSE)gotelse();
else if(lptr->type == DO)gotdo();
else if(lptr->type == STRUCT)structlev++;
}
if(++clev >= &ind[CLEVEL-1]){
fprint(2,"too many levels of curly brackets\n");
clev = &ind[CLEVEL-1];
}
clev->pdepth = 0;
clev->tabs = (clev-1)->tabs;
clearif(clev);
if(strict && clev->tabs > 0)
putspace(' ',NO);
putch(c,NO);
getnl();
if(keyflag == DATADEF){
OUT;
}
else {
OUTK;
}
clev->tabs++;
pt = getnext(0); /* to handle initialized structures */
if(*pt == '{'){ /* hide one level of {} */
while((c=getch()) != '{')
if(c == Beof)error("{");
putch(c,NO);
if(strict){
putch(' ',NO);
eatspace();
}
keyflag = SINIT;
}
continue;
case '}':
pt = getnext(0); /* to handle initialized structures */
if(*pt == ','){
if(strict){
putspace(' ',NO);
eatspace();
}
putch(c,NO);
putch(*pt,NO);
*pt = '\0';
ct = getnl();
pt = getnext(0);
if(*pt == '{'){
OUT;
while((cc = getch()) != '{')
if(cc == Beof)error("}");
putch(cc,NO);
if(strict){
putch(' ',NO);
eatspace();
}
getnext(0);
continue;
}
else if(strict || ct){
OUT;
}
continue;
}
else if(keyflag == SINIT && *pt == '}'){
if(strict)
putspace(' ',NO);
putch(c,NO);
getnl();
OUT;
keyflag = DATADEF;
*pt = '\0';
pt = getnext(0);
}
outs(clev->tabs);
if(--clev < ind)clev = ind;
ptabs(clev->tabs);
putch(c,NO);
lbegin = 0;
lptr=lookup(pt,lastplace+1);
c = *pt;
if(*pt == ';' || *pt == ','){
putch(*pt,NO);
*pt = '\0';
lastplace=pt;
}
ct = getnl();
if((dolevel && clev->tabs <= dotabs[dolevel]) || (structlev )
|| (lptr != 0 &&lptr->type == ELSE&& clev->pdepth == 0)){
if(c == ';'){
OUTK;
}
else if(strict || (lptr != 0 && lptr->type == ELSE && ct == 0)){
putspace(' ',NO);
eatspace();
}
else if(lptr != 0 && lptr->type == ELSE){
OUTK;
}
if(structlev){
structlev--;
keyflag = DATADEF;
}
}
else {
OUTK;
if(strict && clev->tabs == 0){
if((c=getch()) != '\n'){
Bputc(output, '\n');
Bputc(output, '\n');
unget(c);
}
else {
lineno++;
Bputc(output, '\n');
if((c=getch()) != '\n')unget(c);
else lineno++;
Bputc(output, '\n');
}
}
}
if(lptr != 0 && lptr->type == ELSE && clev->pdepth != 0){
UNBUMP;
}
if(lptr == 0 || lptr->type != ELSE){
clev->iflev = 0;
if(dolevel && docurly[dolevel] == NO && clev->tabs == dotabs[dolevel]+1)
clev->tabs--;
else if(clev->pdepth != 0){
UNBUMP;
}
}
continue;
case '(':
paren++;
if ((lptr = lookup(lastlook,p)) != 0){
if(!(lptr->type == TYPE || lptr->type == STRUCT))keyflag=KEYWORD;
if (strict){
putspace(lptr->punc,NO);
opflag = 1;
}
putch(c,NO);
if (lptr->type == IF)gotif();
}
else {
putch(c,NO);
lastlook = p;
opflag = 1;
}
continue;
case ')':
if(--paren < 0)paren = 0;
putch(c,NO);
if((lptr = lookup(lastlook,p)) != 0){
if(lptr->type == TYPE || lptr->type == STRUCT)
opflag = 1;
}
else if(keyflag == DATADEF)opflag = 1;
else opflag = 0;
outs(clev->tabs);
pt = getnext(1);
if ((ct = getnl()) == 1 && !strict){
if(dolevel && clev->tabs <= dotabs[dolevel])
resetdo();
if(clev->tabs > 0 && (paren != 0 || keyflag == 0)){
if(join){
eatspace();
putch(' ',YES);
continue;
} else {
OUT;
split = 1;
continue;
}
}
else if(clev->tabs > 0 && *pt != '{'){
BUMP;
}
OUTK;
}
else if(strict){
if(clev->tabs == 0){
if(*pt != ';' && *pt != ',' && *pt != '(' && *pt != '['){
OUTK;
}
}
else {
if(keyflag == KEYWORD && paren == 0){
if(dolevel && clev->tabs <= dotabs[dolevel]){
resetdo();
eatspace();
continue;
}
if(*pt != '{'){
BUMP;
OUTK;
}
else {
*pt='\0';
eatspace();
unget('{');
}
}
else if(ct){
if(paren){
if(join){
eatspace();
} else {
split = 1;
OUT;
}
}
else {
OUTK;
}
}
}
}
else if(dolevel && clev->tabs <= dotabs[dolevel])
resetdo();
continue;
case ' ':
case '\t':
if ((lptr = lookup(lastlook,p)) != 0){
if(!(lptr->type==TYPE||lptr->type==STRUCT))
keyflag = KEYWORD;
else if(paren == 0)keyflag = DATADEF;
if(strict){
if(lptr->type != ELSE){
if(lptr->type == TYPE){
if(paren != 0)putch(' ',YES);
}
else
putch(lptr->punc,NO);
eatspace();
}
}
else putch(c,YES);
switch(lptr->type){
case CASE:
outs(clev->tabs-1);
continue;
case ELSE:
pt = getnext(1);
eatspace();
if((cc = getch()) == '\n' && !strict){
unget(cc);
}
else {
unget(cc);
if(checkif(pt))continue;
}
gotelse();
if(strict) unget(c);
if(getnl() == 1 && !strict){
OUTK;
if(*pt != '{'){
BUMP;
}
}
else if(strict){
if(*pt != '{'){
OUTK;
BUMP;
}
}
continue;
case IF:
gotif();
continue;
case DO:
gotdo();
pt = getnext(1);
if(*pt != '{'){
eatallsp();
OUTK;
docurly[dolevel] = NO;
dopdepth[dolevel] = clev->pdepth;
clev->pdepth = 0;
clev->tabs++;
}
continue;
case TYPE:
if(paren)continue;
if(!strict)continue;
gottype(lptr);
continue;
case STRUCT:
gotstruct();
continue;
}
}
else if (lbegin == 0 || p > string)
if(strict)
putch(c,NO);
else putch(c,YES);
continue;
case ';':
putch(c,NO);
if(paren != 0){
if(strict){
putch(' ',YES);
eatspace();
}
opflag = 1;
continue;
}
outs(clev->tabs);
pt = getnext(0);
lptr=lookup(pt,lastplace+1);
if(lptr == 0 || lptr->type != ELSE){
clev->iflev = 0;
if(clev->pdepth != 0){
UNBUMP;
}
if(dolevel && docurly[dolevel] == NO && clev->tabs <= dotabs[dolevel]+1)
clev->tabs--;
/*
else if(clev->pdepth != 0){
UNBUMP;
}
*/
}
getnl();
OUTK;
continue;
case '\n':
if ((lptr = lookup(lastlook,p)) != 0){
pt = getnext(1);
if (lptr->type == ELSE){
if(strict)
if(checkif(pt))continue;
gotelse();
OUTK;
if(*pt != '{'){
BUMP;
}
}
else if(lptr->type == DO){
OUTK;
gotdo();
if(*pt != '{'){
docurly[dolevel] = NO;
dopdepth[dolevel] = clev->pdepth;
clev->pdepth = 0;
clev->tabs++;
}
}
else {
OUTK;
if(lptr->type == STRUCT)gotstruct();
}
}
else if(p == string)Bputc(output, '\n');
else {
if(clev->tabs > 0 &&(paren != 0 || keyflag == 0)){
if(join){
putch(' ',YES);
eatspace();
continue;
} else {
OUT;
split = 1;
continue;
}
}
else if(keyflag == KEYWORD){
OUTK;
continue;
}
OUT;
}
continue;
case '"':
case '\'':
putch(c,NO);
while ((cc = getch()) != c){
if(cc == Beof)
error("\" or '");
putch(cc,NO);
if (cc == '\\'){
putch(getch(),NO);
}
if (cc == '\n'){
outs(clev->tabs);
lbegin = 1;
count = 0;
}
}
putch(cc,NO);
opflag=0;
if (getnl() == 1){
unget('\n');
}
continue;
case '\\':
putch(c,NO);
putch(getch(),NO);
continue;
case '?':
question = 1;
gotop(c);
continue;
case ':':
if (question == 1){
question = 0;
gotop(c);
continue;
}
putch(c,NO);
if(structlev)continue;
if ((lptr = lookup(lastlook,p)) != 0){
if (lptr->type == CASE)outs(clev->tabs - 1);
}
else {
lbegin = 0;
outs(clev->tabs);
}
getnl();
OUTK;
continue;
case '/':
if ((cc = getch()) == '/') {
putch(c,NO);
putch(cc,NO);
cpp_comment(YES);
OUT;
lastlook = 0;
continue;
}
else if (cc != '*') {
unget(cc);
gotop(c);
continue;
}
putch(c,NO);
putch(cc,NO);
cc = comment(YES);
if(getnl() == 1){
if(cc == 0){
OUT;
}
else {
outs(0);
Bputc(output, '\n');
lbegin = 1;
count = 0;
}
lastlook = 0;
}
continue;
case '[':
putch(c,NO);
ct = 0;
while((c = getch()) != ']' || ct > 0){
if(c == Beof)error("]");
putch(c,NO);
if(c == '[')ct++;
if(c == ']')ct--;
}
putch(c,NO);
continue;
case '#':
putch(c,NO);
while ((cc = getch()) != '\n'){
if(cc == Beof)error("newline");
if (cc == '\\'){
putch(cc,NO);
cc = getch();
}
putch(cc,NO);
}
putch(cc,NO);
lbegin = 0;
outs(clev->tabs);
lbegin = 1;
count = 0;
continue;
default:
if (c == ','){
opflag = 1;
putch(c,YES);
if (strict){
if ((cc = getch()) != ' ')unget(cc);
if(cc != '\n')putch(' ',YES);
}
}
else if(isop(c))gotop(c);
else {
if(isalnum(c) && lastlook == 0)lastlook = p;
if(isdigit(c)){
putch(c,NO);
while(isdigit(c=Bgetc(input))||c == '.')putch(c,NO);
if(c == 'e'){
putch(c,NO);
c = Bgetc(input);
putch(c, NO);
while(isdigit(c=Bgetc(input)))putch(c,NO);
}
Bungetc(input);
}
else putch(c,NO);
if(keyflag != DATADEF)opflag = 0;
}
}
}
}
/* \fcnfh
Reads partition info
*/
static _Bool
db_part(char **name,
unsigned short *nT,
PREC_TEMP **Tp,
unsigned short *niso,
char ***isonamesp,
PREC_MASS **massp,
PREC_Z ***Zp,
PREC_CS ***CSp)
{
char line[maxline], *lp, *lp2, rc;
//skip comments and blank lines
while((rc=fgetupto(line,maxline,fp)) == '#' ||
rc == '\n')
currline++;
if(!rc) earlyend(dbfilename, currline);
//read name, number of temps, and number of isotopes
if((*name = readstr_sp_alloc(line, &lp, '_')) == NULL)
mpallocerror(0);
checkprepost(lp, 0, *lp==' '||*lp=='\t',*lp=='\0');
int rn;
rn = getnl(2, ' ', lp, niso, nT);
checkprepost(lp, rn!=2, 0, 0);
//allocate temperature, mass, cross section and partition info
PREC_CS **CS;
PREC_Z **Z;
char **inames;
PREC_MASS *mass;
PREC_TEMP *T;
inames = *isonamesp = (char **)calloc(*niso, sizeof(char *));
mass = *massp = (PREC_MASS *)calloc(*niso, sizeof(PREC_ZREC));
T = *Tp = (PREC_TEMP *)calloc(*nT, sizeof(PREC_ZREC));
CS = *CSp = (PREC_CS **)calloc(*niso, sizeof(PREC_CS *));
Z = *Zp = (PREC_Z **)calloc(*niso, sizeof(PREC_ZREC *));
*Z = (PREC_Z *)calloc(*niso**nT,sizeof(PREC_ZREC));
*CS = (PREC_CS *)calloc(*niso**nT,sizeof(PREC_CS));
//get isotope name and mass
while((rc=fgetupto(lp2=line,maxline,fp)) == '#' || rc == '\n')
currline++;
if(!rc) earlyend(dbfilename, currline);
for(int i=0 ; i<*niso ; i++){
Z[i] = Z[0] + *nT*i;
CS[i] = CS[0] + *nT*i;
if((inames[i]=readstr_sp_alloc(lp2,&lp,'_')) == NULL)
mpallocerror(0);
//get mass
mass[i] = strtod(lp,&lp2);
if(i == *niso-1)
break;
checkprepost(lp2, lp==lp2, *lp2==' '||*lp2=='\t',*lp2=='\0');
}
//Last isotope has to be followed by an end of string
checkprepost(lp2, 0, *lp2==' '||*lp2=='\t', *lp2!='\0');
//get for each temperature a line with temperature, partfcn and
//cross-sect info
for(long t=0 ; t<*nT ; t++){
while((rc=fgetupto(lp=line, maxline, fp)) == '#' ||
rc == '\n')
currline++;
if (!rc) earlyend(dbfilename, currline);
while (*lp==' ') lp++;
T[t] = strtod(lp, &lp);
checkprepost(lp, *lp=='\0', *lp==' '||*lp=='\t', *lp=='\0');
for(unsigned short i=0 ; i<*niso ; i++){
Z[i][t] = strtod(lp, &lp);
checkprepost(lp, *lp=='\0', *lp==' '||*lp=='\t', *lp=='\0');
}
unsigned short i;
for(i=0 ; i<*niso-1 ; i++){
CS[i][t] = strtod(lp, &lp);
checkprepost(lp, *lp=='\0', *lp==' '||*lp=='\t', *lp=='\0');
}
CS[i][t] = strtod(lp, &lp);
checkprepost(lp, 0, *lp==' '||*lp=='\t', *lp!='\0');
}
partitionread = 1;
if (--ndb) return 1;
else return 0;
}
