/* outerqn2.c
 * Copyright (c)  MIT/JGOFS                                     *
   From outerqn, renamed to outergn2.
   Modified 8/8/95 WHOI/GLOBEC Georges Bank.  (rcg) Increase token size
	from 880 to 160
   This is a version of outer which allows
   -  quoted strings "This is a string" - will not generate an anchor
       for these in html mode.
   -  does not generate links to flat listings or other options
  grf   02/16/95
*/
#include <stdio.h>
#include <string.h>
#include <malloc.h>
#include "../libsrc/htmlescape.c"

double atof();
double strtod();

#define TOKEN 160
#define NPREFIX 2
#define NINFIX 6
#define NTEST 12
#define NCOMPS 10

int ntotal;
int numsel;
int maxlevel;
int reqlevel;
int cl;
int minlevel;
int htmlflag=0;
int flatflag=0;
int *level=NULL;

static int anchor;
static char *outnames=NULL;
static char objsav[1024],parsav[1024];

char *getenv();

/* 
===========================================
Subroutines required:
===========================================
*/
int ioopen_();
  /*
  int ioopen_(s,nparams,ntotal)
  char s[][80];
  int *nparams,*ntotal;
    s[0..nparams-1]: parameter strings. Inner sets s[j][0]=0 for any 
    strings which it processes; others will be processed by outer. Thus
    selection/projections would normally be ignored by inner.
    nparams: number of parameter strings
    ntotal (returned): total number of variable names
  */
int ioreadrec_();
  /*
  int ioreadrec_(level)
  int *level;
    Read record at appropriate level. Return 0 if end at that level. Return
    1 if ok.
  */
ioclose_();
  /*
  ioclose_()
    Close files
  */
iovalreal_();
  /*
  iovalreal_(vn,f)
  int *vn;
  float *f;
    Return real value (f) for variable indexed by vn. -9999 for strings
  */
iovalstr_();
  /*
  iovalstr_(vn,tmp)
  int *vn;
  char *tmp;
   Return string value (tmp) for variable indexed by vn.
  */
int iovarlevel_();
  /*
  int iovarlevel_(vn)
  int *vn;
    Return level corresponding to variable indexed by vn.   
  */
int ioattrout_();
  /*
  ioattrout_(vn, str)
  int *vn;
  char *str;
    Output next attribute for variable indexed by vn. 0=none left.
  */
ioname_();
  /*
  ioname_(vn,s)
  int *vn;
  char *s;
    Return name (s) corresponding to variable number vn.
  */
int iocommout_();
  /*
  int iocommout_(s)
  char *s;
    Return next comment string. 0=none left.
  */


quit(s1,s2)
char *s1,*s2;
{
printf("&x %s %s\n",s1,s2);
if(htmlflag)printf("</pre>\n");
if(level)free(level);
if(outnames)free(outnames);
ioclose_();
exit(1);
}
  
/*
P-code form:
128..255 exec test(n-128) [negative #'s on Sun]
1=or,2=and,3=not
0: end
*/

static char *prefix[]={"!","not"};
static char *infix[]={"||","or","|","&&","and","&"};
static char *testop[]={"<","","=","==",">","","<=","","<>","!=",">=",""};
static struct {
  int testlev;
  int testvar;
  int testcode;
  float testval;
  char teststr[TOKEN];
  int testres;
  } tst[NCOMPS];
static int tstcnt=0;
static int tstproc[NCOMPS];
static int tstproccnt=0;

parse(s)
char *s;
{
int opstkp;
float f;
int opstack[20];
int opprior[20];
char *t,*t1,tok[TOKEN],tmp[TOKEN];
int i,state,typ,finaland;

opstkp= 0;
opstack[0]=0;
opprior[0]= -1;
state=0;

if (tstproccnt==0) finaland=0; else finaland=1;

strcat(s,"$");
t=s;
while(t != NULL){
  if(isalpha(*t) || *t=='#'){
    i=strcspn(t," ()$=<>|&!");
    strncpy(tok,t,i);
    tok[i]='\0';
    t += i;
    typ=0;
    }
  else if(isdigit(*t) || *t=='.' || *t=='+' || *t=='-'){
    f=strtod(t,&t1);
    typ=1;
    t=t1;
    }
  else if(*t == '('){
    typ=2;
    t++;
    }
  else if(*t == ')'){
    typ=3;
    t++;
    }
  else if(*t == '$'){
    typ=4;
    t=NULL;
    }
  else {
    tok[0]= *t;
    if((int)strspn(t+1,"&|=>")>0){
      tok[1]= *(++t);
      tok[2]='\0';
      }
    else
      tok[1]='\0';
    t++;
    typ=5;
  };
  switch(state*10+typ){
    case 0: /* variable,not */
      for(i=0;i<ntotal;i++){
	ioname_(&i,tmp);
	if(!strcmp(tok,tmp))break;
        };
      if(i== ntotal){
        if(strcmp(tok,prefix[1])) quit("Error on conditions",tok);
        opstack[++opstkp]=3;
        opprior[opstkp]=10;
        break;
        };
      tst[tstcnt].testvar=i;
      tst[tstcnt].testlev=iovarlevel_(&i);
      state=1;
      tstproc[tstproccnt++]=tstcnt+128;
      break;
    case 2: /* ( */
      opstack[++opstkp]=4;
      opprior[opstkp]=1;
      break;
    case 5: /* ! */
      for (i=0;i<NPREFIX;i++)if(strcmp(tok,prefix[i])==0)break;
      if(i==NPREFIX) quit("Bad prefix operator",tok);
      opstack[++opstkp]=3;
      opprior[opstkp]=10;
      break;
    case 15: /* ==, etc. */
      for(i=0;i<NTEST;i++)if(strcmp(tok,testop[i])==0)break;
      if(i==NTEST)quit("Bad condition",tok);
      tst[tstcnt].testcode=(i/2)+1;
      state=2;
      break;
    case 20: /* string comparison */
      strcpy(tst[tstcnt].teststr,tok);
      tst[tstcnt].testcode *= -1;
      tstcnt++;
      state=3;
      break;
    case 21: /* float comparison */
      tst[tstcnt++].testval=f;
      state=3;
      break;
    case 30:
    case 35: /* and or */
      for(i=0;i<NINFIX;i++)if(strcmp(tok,infix[i])==0)break;
      if(i==NINFIX) quit("Bad infix operator",tok);
      i=(i/3)+1;
      while(i+1<=opprior[opstkp]){
	tstproc[tstproccnt++]=opstack[opstkp--];
	};
      opstack[++opstkp]=i;
      opprior[opstkp]=i+1;
      state=0;
      break;
    case 33: /* ) */
      while(1<opprior[opstkp]){
        tstproc[tstproccnt++]=opstack[opstkp--];
	};
      opstkp--;
      break;
    case 34: /* $ */
      while(0<=opprior[opstkp]){
        tstproc[tstproccnt++]=opstack[opstkp--];
	};
      if(finaland)tstproc[tstproccnt++]=2;
#ifdef DEBUG
      printf("finaland = %d, tstproccnt = %d\n",finaland,tstproccnt);
      for(i=0;i<tstproccnt;i++){
	printf("tstproc %d ",tstproc[i]);
	if(tstproc[i]>=128)
	  printf("%d %d %f %s\n",tst[tstproc[i]-128].testvar,
		 tst[tstproc[i]-128].testcode,tst[tstproc[i]-128].testval,
		 tst[tstproc[i]-128].teststr);
	else printf("\n");
      };
#endif
      return;
    default:
      quit("Parsing error","");
    };
  if(*t==' ')t++;
  };
}

int test(cl)
int cl;
{
float f;
char tmp[TOKEN];
int i,j,m;
int stk[NCOMPS],stkp;

#ifdef DEBUG
  printf("entered test at level %d\n",cl);
#endif
stkp = -1;
stk[0]=1;
for(i=0;i<tstproccnt;i++){
  if(tstproc[i] >= 128){
    j=tstproc[i]-128;
    if(tst[j].testlev < cl)stk[++stkp]=tst[j].testres;
    else if(tst[j].testlev>cl)stk[++stkp]=1;
    else {
      m=0;
      if(tst[j].testcode<0){
	iovalstr_(&(tst[j].testvar),tmp);
	switch (tst[j].testcode){
	case -2:if(strcmp(tmp,tst[j].teststr)==0)m=1;break;
	case -5:if(strcmp(tmp,tst[j].teststr)!=0)m=1;break;
	};
      } else {
	iovalreal_(&(tst[j].testvar),&f);
#ifdef DEBUG
 	printf("-- %f %d %f -- ",f,tst[j].testcode,tst[j].testval);
#endif
	switch (tst[j].testcode){
	case 1:if(f<tst[j].testval)m=1;break;
	case 2:if(f==tst[j].testval)m=1;break;
	case 3:if(f>tst[j].testval)m=1;break;
	case 4:if(f<=tst[j].testval)m=1;break;
	case 5:if(f!=tst[j].testval)m=1;break;
	case 6:if(f>=tst[j].testval)m=1;break;
	};
      };
    tst[j].testres=m;
    stk[++stkp]=m;    
    };
  } else {
    switch(tstproc[i]){
    case 1: stk[stkp-1]=stk[stkp-1] || stk[stkp];stkp--;break;
    case 2: stk[stkp-1]=stk[stkp-1] && stk[stkp];stkp--;break;
    case 3: stk[stkp]=1-stk[stkp];
    };
  };
};
#ifdef DEBUG
  printf("test result = %d\n",stk[0]);
#endif
return stk[0];
}

outvar()
{
  int i,j;
  char tmp[80],attr[1024];
  printf("Content-type: text/plain\n\n"); 
  i=1;
  while(iocommout_(tmp)){
    if(i){printf("&c***********************\n");i=0;};
    printf("%s\n",tmp);
  };
  for(i=0;i<=maxlevel;i++){
    printf("&v%d======================\n",i);
    for(j=0;j<ntotal;j++)if(*(level+j)==i){
      strcpy(attr,"[");
      while(ioattrout_(&j,tmp)){strcat(attr,tmp);strcat(attr,";");};
      if(strcmp(attr,"[")) attr[strlen(attr)-1]=']'; else attr[0]=0;
      strcpy(tmp,outnames+j*TOKEN);
      strcat(tmp,attr);
      printf("%s\t",tmp);
    };
    printf("\n");
  };
  printf("&r=======================\n",i);
}

output(firsttime)
int firsttime;
{
  int i,j;
  char tmp[80];
  if(minlevel<maxlevel || firsttime){
    i=1;
    while(iocommout_(tmp)){
      if(i){printf("&c***********************\n");i=0;};
      printf("%s\n",tmp);
    };
    for(i=minlevel;i<=maxlevel;i++){
      printf("&d%d----------------------\n",i);
      for(j=0;j<ntotal;j++)if(*(level+j)==i){
	iovalstr_(&j,tmp);
	printf("%s\t",tmp);
      };
      printf("\n");
    };
    minlevel=maxlevel;
  } else {
    for(j=0;j<ntotal;j++)if(*(level+j)==maxlevel){
      iovalstr_(&j,tmp);
      printf("%s\t",tmp);
    };
    printf("\n");
  };
}

outvarflat()
{
  int j;
  char tmp[80],attr[1024];
  printf("Content-type: text/plain\n\n"); 
  for(j=0;j<ntotal;j++)if(*(level+j)>=0){
    strcpy(attr,"[");
    while(ioattrout_(&j,tmp)){strcat(attr,tmp);strcat(attr,";");};
    if(strcmp(attr,"[")) attr[strlen(attr)-1]=']'; else attr[0]=0;
    strcpy(tmp,outnames+j*TOKEN);
    strcat(tmp,attr);
    printf("%s\t",tmp);
  };
  printf("\n");
}

outputflat(firsttime)
int firsttime;
{
  int j;
  char tmp[80];
  for(j=0;j<ntotal;j++)if(*(level+j)>=0){
    iovalstr_(&j,tmp);
    printf("%s\t",tmp);
  };
  printf("\n");
  minlevel=maxlevel;
}

outvarhtml()
{
  char tmp[80],*sp;
  printf("Content-type: text/html\n\n");
  htmlescape(parsav);
  printf("<title>%s -- Level %d</title>\n",objsav,reqlevel);
  printf("<h1>%s --%s-- Level %d</h1>\n",objsav,parsav,reqlevel);
/*
  printf("<a href=\"/jg/serv%s.html%d?%s\">List all at next level</a> | ",
	 objsav,reqlevel+1,parsav);
  printf("<a href=\"/jg/serv%s.flat%d?%s\">Flat listing</a> | ",
	 objsav,reqlevel,parsav);
  printf("<a href=\"%s/jg/otheropt%s%s.html%d?%s\">Other options</a><p>\n",
         OPTIONHOST,MYADDR,objsav,reqlevel,parsav);
*/
  printf("<pre>\n");
  if (*parsav) strcat(parsav,",");
  while(iocommout_(tmp)) {
    htmlescape(tmp);
    printf("# %s\n",tmp);
  };
}

outputhtml(firsttime)
int firsttime;
{
  int i,j;
  char tmp[80];
  if(minlevel<maxlevel || firsttime){
    while(iocommout_(tmp)) {
      htmlescape(tmp);
      printf("# %s\n",tmp);
    };
    for(i=minlevel;i<=maxlevel;i++){
      printf("=========================\n");
      for(j=0;j<ntotal;j++)if(*(level+j)==i){
	printf("%s\t",outnames+j*TOKEN);
      };
      printf("\n");
      printf("-------------------------\n");
      for(j=0;j<ntotal;j++)if(*(level+j)==i){
	iovalstr_(&j,tmp);
	if(*tmp == '\"'){
	  tmp[strlen(tmp)-1]=0;
	  printf("%s\t",tmp+1);
	} else if(i==reqlevel && anchor)
	  printf("<a href=\"/jg/serv%s.html%d?%s%s=%s\">%s</a>\t",
            objsav,reqlevel+1,parsav,outnames+j*TOKEN,tmp,tmp);
	else printf("%s\t",tmp);
      };
      printf("\n");
    };
    minlevel=maxlevel;
  } else {
    for(j=0;j<ntotal;j++)if(*(level+j)==maxlevel){
      iovalstr_(&j,tmp);
      if(*tmp == '\"'){
	tmp[strlen(tmp)-1]=0;
	printf("%s\t",tmp+1);
      } else if(anchor)
	printf("<a href=\"/jg/serv%s.html%d?%s%s=%s\">%s</a>\t",
	       objsav,reqlevel+1,parsav,outnames+j*TOKEN,tmp,tmp);
      else printf("%s\t",tmp);
    };
    printf("\n");
  };
}

main(argc,argv)
int argc;
char *argv[];
{
  int i,j,k,nparams,noprojflag;
  char tmp[80],stmp[80],*pp,*sp;

  reqlevel=999;
  pp=getenv("PATH_INFO");
  if(pp) {
    if(sp=strstr(pp,".html")) htmlflag=1;
    else if(sp=strstr(pp,".flat")) flatflag=1;
    if(htmlflag || flatflag) {
      reqlevel=atoi(sp+5);
      *sp=0;
      strcpy(objsav,pp);
      pp=getenv("QUERY_STRING");
      if(pp) strcpy(parsav,pp); else *parsav=0;
    };
  };
  level=NULL;
/*
  system("printenv > /d3/glenn/jg/temp");
*/


  nparams=argc-1;
  k=ioopen_(argv+1,&nparams,&ntotal);
  if(k>reqlevel) anchor=1; else anchor=0;
  noprojflag=1;numsel= -1;
  outnames=(char *)malloc(TOKEN*ntotal);
  for(j=0;j<ntotal;j++)ioname_(&j,outnames+j*TOKEN);

  level= (int *)malloc(ntotal*sizeof(int));
  for(i=0;i<ntotal;i++)*(level+i)= -iovarlevel_(&i)-1;

  for (i=1;i<argc;i++) {
    strcpy(stmp,argv[i]);
    if(*stmp){
      if(strpbrk(stmp,"<>=")) parse(stmp);
      else {
	for(k=0;k<ntotal;k++){if(strcmp(stmp,outnames+k*TOKEN)==0)break;};
	if(k==ntotal)quit("Bad Name",stmp);
	if(*(level+k)<0 && *(level+k) >= -reqlevel-1) 
	  *(level+k)= - *(level+k)-1;
	noprojflag=0;
      };
    };
  };

  if(noprojflag)for(i=0;i<ntotal;i++)
    if(*(level+i) >= -reqlevel-1) *(level+i)= - *(level+i)-1;

  maxlevel=0;
  for(j=0;j<ntotal;j++)
    if(*(level+j)>maxlevel)maxlevel= *(level+j);
/*
  for(j=0;j<tstproccnt;j++)
    if(tstproc[j] >= 128)
      if(tst[j].testlev > maxlevel) maxlevel=tst[j].testlev;
*/
  if(reqlevel>maxlevel) reqlevel=maxlevel;

#ifdef DEBUG
  printf("maxlevel %d\n",maxlevel);
#endif

  minlevel=maxlevel;
  cl=0;
  k=1;
  if(htmlflag) outvarhtml();
  else if (flatflag) outvarflat();
  else outvar();

  while(1){
    if(ioreadrec_(&cl)){
#ifdef DEBUG
      printf("succ. read at lev %d\n",cl);
#endif
      if(test(cl))
	if(cl<maxlevel){
	  if(cl<minlevel)minlevel=cl;
	  cl++;
	} else {
	  if(htmlflag)outputhtml(k); 
	  else if(flatflag)outputflat(k); 
	  else output(k);
	  k=0;
	};
    } else {
#ifdef DEBUG
      printf("failed read at lev %d\n",cl);
#endif
      if(cl == 0)break;
      cl--;
    };
  };
  ioclose_();
  free(level);
  if(htmlflag)printf("</pre>\n"); 
  else if(!flatflag)printf("&e**** End of object ****\n");
}