static int input(int y0, char*hlp, char*directive, char*text, int cur, int lim)
{
if(errTxt)
{
goto_xy(1,y0+1);
scrcolor(FGmain,BGmain);
clr_eol();
print("%s","Error: ");
scrcolor(Red,BGmain);
print("%s",errTxt);
scrcolor(FGmain,BGmain);
if(blind) { printf("%s\n",errTxt); sortie(110);}
be_be();
errTxt=NULL;
}
for(;;)
{ int rc;
goto_xy(1,y0);
scrcolor(FGmain,BGmain);
print("%s",directive);
rc=str_redact(text,cur, lim);
switch(rc)
{ case KB_UP:
case KB_DOWN:
case KB_PAGED:
case KB_PAGEU: prtcllist(rc); continue;
case KB_F1: show_help(hlp); continue;
}
return rc;
}
}
int informline(long curent, long total)
{ int xx;
int res=0;
static int Y;
static int X;
int xm=where_x(), ym=where_y(), fc=fColor,bc=bColor;
if (curent == 0)
{
Y=maxRow();
X=15;
/* goto_xy(15,Y); scrcolor(Black,Red);
memset(b,' ',50); b[50] = '\0';
*/
goto_xy(15,Y); scrcolor(White,Red);
print(" Calculation in progress. Calculation in progress.");
/* 12345678901234567890123456789012345678901234567890 */
/* print("%s",b); */
goto_xy(15,Y);
scrcolor(fc,bc);
escpressed();
goto exi;;
}
if(X>65 || X<15) X=15;
xx = 15 + (50 * curent) / total;
if (xx > 65) xx = 65;
scrcolor(Black,Red);
goto_xy(X,Y);
while (where_x() < xx) print("%c",'X');
if(xx !=X) { X=xx; res=escpressed();}
if (curent >= total)
{
scrcolor(White,Black);
goto_xy(1,Y); clr_eol();
goto exi;
}
exi:
scrcolor(fc,bc); goto_xy(xm,ym);
return res;
}
void menu0(int col,int row,char* label, char* menstr ,
void (**funcKey)(int) ,char** funcKeyMess, void ** hscr, int* kk)
{ int i, j, k, col1, npage,lastrow;
long lastpage;
int ink;
int ncol;
void * pscr;
int fkPos[11];
int height;
char fmt[20];
int lastLine;
/* colors */
int label_fg =Yellow;
int label_bg =Blue;
int help_fg1 =White;
int help_fg2 =Black;
int help_bg =DarkGray;
int box_fg =White;
int box_bg =DarkGray;
int star_fg =Red;
int page_fg =Black;
int actFunc_fg=Black;
int actFunc_bg=White;
int pasFunc_fg=White;
int pasFunc_bg=DarkGray;
/* save colors */
int fcolor_tmp=fColor;
int bcolor_tmp=bColor;
void *hscr_=NULL;
if(hscr==NULL) hscr=&hscr_;
lastLine=maxRow();
if (funcKey == NULL) for (i=0;i<11;i++) fkPos[i]=0; else
{ int xx;
scrcolor(FGmain,BGmain);
goto_xy(1,lastLine); clr_eol();
xx=0;
for (j=0;j<10;j++) { if(funcKey[j] && funcKeyMess[j])
xx=xx+4+strlen(funcKeyMess[j]);}
xx= (80 - xx )/2 ;
goto_xy(xx,lastLine);
for (i=0;i<10;i++)
{ fkPos[i]=where_x();
if (funcKey[i] && funcKeyMess[i])
{ scrcolor(help_fg1,help_bg); print(" F%i-",i+1);
scrcolor(help_fg2,help_bg); print(funcKeyMess[i]);
}
}
fkPos[10]=where_x();
}
clearTypeAhead();
if (*kk < 0) *kk = -(*kk);
ncol=menstr[0];
sprintf(fmt,"%%%d.%ds",ncol,ncol);
height=strlen(menstr)/ncol;
if(height==0) { *kk=0; return; }
if (row+height+1 >lastLine-2) height=lastLine-3-row;
lastpage = 1+ (strlen(menstr)/ncol -1)/height ;
if(label[0] ==0 || row == 1)
{ if (*hscr == NULL) get_text(col,row,col+ncol+1,row+2,hscr);}
else
{ char label_[STRSIZ];
int shft,sz;
if (*hscr == NULL) get_text(col,row-1,col+ncol+1,row+2,hscr);
for(i=0;i<ncol+2;i++) label_[i]=' ';
label_[ncol+2]=0;
sz=strlen(label);
if(sz >ncol+2) {shft=0;sz=ncol+2;} else shft=(ncol+2 -sz)/2;
memcpy(label_+shft,label,sz);
scrcolor(label_fg,label_bg);
goto_xy(col,row-1);
print(label_);
}
get_text(col,row + 3,col + ncol + 1,row + height + 1,&pscr);
if (*kk <= 0 || *kk > lastpage * height )
{ npage = 1;
k = 1;
}
else {
k = ((*kk) - 1) % height + 1;
npage = ((*kk) - 1) / height + 1;
}
col1 = col + 1;
label_1:
scrcolor(box_fg,box_bg);
chepbox(col,row,col + ncol + 1,row + height + 1);
scrcolor(star_fg,box_bg);
goto_xy(col+1,row); print("<");
/* goto_xy(col+1,row + height+1); print("?"); */
scrcolor(page_fg,box_bg);
if (npage > 1)
{
goto_xy(col + ncol - 2,row);
print("PgUp");
}
if (npage < lastpage)
{
goto_xy(col + ncol - 2,row + height + 1);
print("PgDn");
}
if(npage<lastpage) lastrow=height;
else lastrow = (strlen(menstr)/ncol)%height;
lastrow=0;
scrcolor(pasFunc_fg,pasFunc_bg);
for (j = 1; j <= height; j++)
{ int shift;
goto_xy(col + 1,row + j);
shift=1+(j-1 + (npage-1)*height)*ncol;
if(shift<strlen(menstr)) {print(fmt,menstr+shift );lastrow++;}
else print(fmt," ");
}
scrcolor(actFunc_fg,actFunc_bg);
if (k > lastrow) k = lastrow;
goto_xy(col + 1,row + k);
if (lastrow) print(fmt,menstr+1+(k-1+(npage-1)*height)*ncol);
while (1)
{ int jump=1,mousePos;
scrcolor(pasFunc_fg,pasFunc_bg);
ink = inkey();
/* mouse filter */
if ((ink==KB_MOUSE)&&(mouse_info.but1 == 2))
{
if (mouse_info.row == lastLine )
for(i=0; i<10;i++)
if ((mouse_info.col > fkPos[i]) && (mouse_info.col < fkPos[i+1]))
{ if (i==9)ink='0'; else ink='1'+i;}
if ( (mouse_info.col >= col ) && (mouse_info.col <=col+ncol+1) )
{ mousePos = mouse_info.row - row;
if (col+ncol+1-mouse_info.col <4)
{
if (mousePos==0) ink=KB_PAGEU;
if (mousePos==height+1) ink=KB_PAGED;
}
if ((mousePos == 0 ) && ( mouse_info.col-col <4)) ink=KB_ESC;
if ((mousePos < 0)&&(mousePos >= height))
{
if (mousePos > k) {ink=KB_DOWN; jump=mousePos - k;}
if (mousePos < k) {ink=KB_UP; jump=k - mousePos;}
if (mousePos==k ) ink=KB_ENTER;
}
}
}
/* end of filter */
if (lastrow == 0) goto label_3;
label_4:
switch (ink)
{
case KB_MOUSE:
if (mouse_info.but1 != 2) break;
if (mouse_info.row == lastLine )
for(i=0; i<10;i++)
if ((mouse_info.col > fkPos[i]) && (mouse_info.col < fkPos[i+1]))
{ if (i==9)ink='0'; else ink='1'+i;
goto label_4;
}
if ( (mouse_info.col < col ) || (mouse_info.col >col+ncol+1) ||
(mouse_info.row < row ) || (mouse_info.row >row+height+1) ) break;
mousePos = mouse_info.row - row;
if ((mousePos == 0 ) && ( mouse_info.col-col <4)) ink=KB_ESC;
if ((mousePos != 0)&&(mousePos != height+1))
{
if (mousePos > k) { ink=KB_DOWN; jump=mousePos - k;}
if (mousePos < k ) { ink=KB_UP; jump=k - mousePos;}
}
if (col+ncol+1-mouse_info.col <4)
{
if (mousePos==0) ink=KB_PAGEU;
if (mousePos==height+1) ink=KB_PAGED;
}
if (mousePos==k ) ink=KB_ENTER;
if (ink!=KB_MOUSE) goto label_4;
break;
case KB_DOWN:
if(k==lastrow)
{
if(npage < lastpage)
{ k=1;
npage++;
goto label_1;
} else { be_be(); break; }
}else {ink= KB_RIGHT; goto label_4;}
case KB_UP:
if(k==1)
{
if (npage > 1)
{
k=height;
npage--;
goto label_1;
}
else{ be_be(); break; }
}
else {ink= KB_LEFT; goto label_4;}
case KB_LEFT:
goto_xy(col1,row + k);
print(fmt,menstr+1+(k-1+(npage-1)*height)*ncol);
k = k - jump;
if (k == 0) k = lastrow;
scrcolor(actFunc_fg,actFunc_bg);
goto_xy(col1,row + k);
print(fmt,menstr+1+(k-1+(npage-1)*height)*ncol);
break;
case KB_RIGHT:
goto_xy(col1,row + k);
print(fmt,menstr+1+(k-1+(npage-1)*height)*ncol);
k = k + jump ;
if(k > lastrow) k = 1;
scrcolor(actFunc_fg,actFunc_bg);
goto_xy(col1,row + k);
print(fmt,menstr+1+(k-1+(npage-1)*height)*ncol);
break;
case KB_ENTER:
scrcolor(box_fg,box_bg);
chepbox(col,row,col+ncol+1,row+2);
put_text(&pscr);
goto_xy(col1,row + 1);
scrcolor(actFunc_fg,actFunc_bg);
print(fmt,menstr+1+(k-1+(npage-1)*height)*ncol);
*kk = (npage - 1) * height + k;
goto_xy(1,lastLine);scrcolor(FGmain,BGmain); clr_eol();
if(hscr==&hscr_) put_text(hscr);
refresh_scr();
escpressed();
return;
case KB_BACKSP:
case KB_ESC:
goto label_3;
case KB_PAGEU:
if (npage > 1)
{
npage--;
goto label_1;
}
else
be_be();
break;
case KB_PAGED:
if (npage < lastpage)
{
npage++;
goto label_1;
}
else
be_be();
break;
case '1': case'2': case '3': case '4': case '5':
case '6': case '7': case '8': case '9': case '0':
case KB_F1: case KB_F2: case KB_F3: case KB_F4: case KB_F5:
case KB_F6: case KB_F7: case KB_F8: case KB_F9: case KB_F10:
{ int fk;
if (funcKey==NULL) break;
if ( ink>='0' && ink <='9') { fk=ink-'0';if (fk==0) fk=10;}
else fk=ink-KB_F1+1;
if ((funcKey[fk-1]) != NULL)
{ void * saveHlp;
get_text(1,lastLine,maxCol(),lastLine,&saveHlp);
scrcolor(box_fg,box_bg);
goto_xy(col+1,row); print("%c",boxFrame[1]);
goto_xy(1,lastLine);scrcolor(FGmain,BGmain);clr_eol();
(*funcKey[fk-1])((npage-1)*height+k);
put_text(&saveHlp);
scrcolor(star_fg,box_bg);
goto_xy(col+1,row); print("<");
}
break;
}
case 6: /* ^F */
case 'f':
case 'F':
{ char name[32]="";
int key=correctStr(5,23,"Enter name(Esc for exit):",name,30,1);
if(key)
{ char * x=strstr(menstr+1,name);
if(x==NULL)
{
if(blind) sortie(121); else messanykey(10,10, "Not detected");
break;
}
k= ((x-menstr)-1)/menstr[0];
npage = k/ height + 1;
k = k% height + 1;
goto label_1;
}
} break;
}
}
label_3:
put_text(hscr);
put_text(&pscr);
scrcolor(fcolor_tmp,bcolor_tmp);
goto_xy(1,lastLine); clr_eol();
*kk = 0;
}
void mk_reduceprograms(void)
{
int ndel, ncalc, nrest, i;
long nrecord, naxu;
csdiagram csd;
unsigned ncalctot;
shortstr txt;
hlpcsptr gstlist, c;
vcsect vcs_copy;
s_listptr d_facts, df;
rmptr t_fact;
goto_xy(1,21); scrcolor(Yellow,Blue);
print(" REDUCE code generation \n");
scrcolor(Red,BGmain);
print(" Generated........\n");
print(" current diagram :\n");
scrcolor(Yellow,Blue);
print(" Press Esc to halt REDUCE codes generation ");
scrcolor(FGmain,BGmain);
diagrq=fopen(DIAGRQ_NAME,"rb");
ncalctot = 0;
menuq=fopen(MENUQ_NAME,"rb");
for(nsub=1;nsub<=subproc_sq;nsub++)
{
rd_menu(2,nsub,txt,&ndel,&ncalc,&nrest,&nrecord);
fseek(diagrq,nrecord*sizeof(csdiagram),SEEK_SET);
naxu = ndel + ncalc + nrest;
for (ndiagr = 1; ndiagr <= naxu; ndiagr++)
{
goto_xy(20,22); print("%u",ncalctot);
goto_xy(20,23); print("%u",ndiagr); clr_eol();
FREAD1(csd,diagrq);
if (csd.status != -1)
{
outFileOpen("%sresults%cp%d_%d.red",pathtouser,f_slash,nsub,ndiagr);
writeLabel('%');
writeF("%%\n");
transfdiagr(&csd,&vcs);
cwtarg(&vcs);
if (vcs.clrnum == 0)
{
writeF(
"%%------- Zero color factor --------\n");
writeF("totFactor_:=0$\n");
writeF("numerator_:=0$\n");
writeF("denominator_:=1$\n");
}
else
{
generateghosts(&vcs,&gstlist);
if (gstlist == NULL)
{
writeF( "%%------- non-existent diagram --------\n");
writeF("totFactor_:=0$\n");
writeF("numerator_:=0$\n");
writeF("denominator_:=1$\n");
}
else
{
goto_xy(40,23);
print("(%% %4d subdiagrams)",gstlist->maxnum);
writeF("%% The total number of diagrams %d\n",gstlist->maxnum);
preperdiagram();
head();
emitfactors();
diagramsrfactors(gstlist,&d_facts,&t_fact);
writeF("totFactor_:=%s$\n",rmonomtxt(*t_fact));
writeF("totFactor_:="
"totFactor_*SymmFact*AverFact*FermFact*ColorFact$\n");
clrvm(t_fact->n.v);
clrvm(t_fact->d.v);
free(t_fact);
writesubst();
writeF("numerator_:=0$\n");
c = gstlist;
df = d_facts;
vcs_copy = vcs;
while (c != NULL)
{
coloringvcs(c);
writeF("%% diagram number = %d\n", c->num);
DiagramToOutFile(&vcs,1,'%');
{int k; int sgn=c->sgn;
for(k=0;k<vcs.sizet;k++) sgn*=vertexes[k].lgrnptr->factor;
writeF(" GhostFact:=%d$\n",sgn);
}
findReversVert();
attachvertexes();
emitreducecode();
writeF(" numerator_:=numerator_ +(%s)*GhostFact*Vrt_1 $\n",
smonomtxt(df->monom));
writeF(" Clear Vrt_1,GhostFact$\n");
writeF("%%\n");
vcs = vcs_copy;
c = c->next;
df = df->next;
}
eraseslist(d_facts);
eraseghosts(gstlist);
vcs = vcs_copy;
emitdenoms();
writeF(" Clear p%d",nin + nout + 1);
for (i = nin + nout + 2; i <= 12; i++)
writeF(",p%d",i);
writeF("$\n");
writeF("%%\n");
}
}
writeF("End$\n");
outFileClose();
--(nrest);
++(ncalctot);
if (escpressed()) goto exi;
}
}
}
exi:
fclose(diagrq); fclose(menuq);
clrbox(1,21,70,24);
}
void paramdependence(r_func ff, char* procname, char* resultname)
{
double minprm, maxprm;
int npoints;
double memprm, stepprm;
unsigned count;
double f[201];
int Esc=0,mPos=1;
double prmval;
char txt[100];
char name[20];
double *vPos;
if(!selectParam(54,11,"Choose parameter",NULL,nin_int==2,1,0,&vPos,name,&mPos)) return;
memprm=*vPos;
minprm = memprm;
maxprm = minprm;
label1:
sprintf(txt,"'%s' min=",name);
if (!correctDouble(55,14,txt,&minprm,0)) return;
label2:
sprintf(txt,"'%s' max=",name);
if (!correctDouble(55,15,txt,&maxprm,0))
{
goto_xy(55,14);
clr_eol();
goto label1;
}
if (maxprm <= minprm)
{
messanykey(55,17,"Range check error");
goto_xy(55,15);
clr_eol();
goto label2;
}
label4: npoints = 101;
if (correctInt(55,16,"Number of points= ",&npoints,0))
{
if (npoints < 3)
{
messanykey(55,17,"Too few points!");
goto label4;
}
if (npoints > 201)
{
messanykey(55,17,"Too many points!");
goto label4;
}
}
else
{
goto_xy(55,15);
clr_eol();
goto label2;
}
goto_xy(55,14); clr_eol();
goto_xy(55,15); clr_eol();
goto_xy(55,16); clr_eol();
stepprm = (maxprm - minprm)/(npoints - 1);
informline(0,npoints);
stepprm = (maxprm - minprm) / (npoints - 1);
prmval=minprm;
err_code = 0;
for(count = 1; count <= npoints; count++)
{
*vPos=prmval;
err_code=checkParam();
if(err_code>1) break;
f[count-1] = ff();
if(err_code>1) break;
Esc=informline(count,npoints);
if(Esc) break;
prmval += stepprm;
}
if(err_code) errormessage();
strcpy(txt,name);
if (err_code <=1 && Esc==0) plot_1(minprm,maxprm ,npoints,f,NULL,procname,txt,resultname);
*vPos=memprm;
calcFunc_int();
}
static void clearstatistic(void)
{int i; for (i = 17; i < 24; i++) { goto_xy(1,i); clr_eol();} }
int viewresults(void)
{
int k,kmenu;
void * pscr = NULL;
shortstr newname;
int dirStat=checkDir("results");
if(dirStat==0){messanykey(10,15,"directory RESULTS is empty"); return 1;}
kmenu = 1;
label_1:
menu1(10,10,"","\010"
" View "
" Delete "
" Rename ","s_res",&pscr,&kmenu);
switch (kmenu)
{
case 0:
return 0;
case 1:
viewDir("results");
break;
case 2:
if(dirStat==2)
{ char mess[]="Can not clean dir 'results' because it contains the LOCK file";
if(blind) { printf("%s\n",mess); sortie(102);}
else { messanykey(3,13,mess); break;}
}
if ( mess_y_n( 6,13," Delete files ") )
{ struct dirent **namelist;
int n,i;
n = scandir("./results", &namelist, NULL, NULL);
for(i=0; i<n;i++)
{
char buff[100];
if(strcmp(namelist[i]->d_name,"aux") && strcmp(namelist[i]->d_name,"..") && strcmp(namelist[i]->d_name,".") )
{
sprintf(buff,"rm -rf results/%s",namelist[i]->d_name);
if(unlink(buff+7)) system(buff);
}
}
free(namelist);
}
put_text(&pscr);
return 1;
case 3:
strcpy(newname," ");
while(1)
{ void * pscr3;
get_text(1,maxRow(),maxCol(),maxRow(),&pscr3);
goto_xy(1,maxRow());
print("Enter new name: ");
k = str_redact(newname,1,30);
if (k == KB_ESC)
{ goto_xy(1,24);
clr_eol();
goto label_1;
}
if (k == KB_ENTER)
{
trim(newname);
if(rename("results",newname)==0)
{ char command[200];
mkdir("results",0755);
sprintf(command," cp -rp %s/aux results",newname);
system(command);
put_text(&pscr);
put_text(&pscr3);
return 1;
}
else messanykey(10,15," Can't rename the directory");
}
put_text(&pscr3);
}
}
goto label_1;
}
int qcdmen_(void)
{
void * pscr=NULL;
int mode;
int returnCode=0;
initStrFun(0);
L10:{ char strmen[]="\030"
" parton dist. alpha OFF "
" alpha(MZ)= ZZZZ "
" nf = NF "
" order= NNLO "
" mb(mb)= MbMb "
" Mtop(pole)= Mtp "
" Alpha(Q) plot "
" Qren = RRR "
" Qpdf1= FF1 "
" Qpdf2= FF2 "
" Qshow= FFS ";
if(alphaPDF)
{ int k=0;
//printf("alphaPDF=%d sf_alpha[0]=%p sf_alpha[1]=%p\n", alphaPDF,sf_alpha[0],sf_alpha[1]);
switch(alphaPDF)
{ case 1: if(sf_alpha[0]) k=1; else if(sf_alpha[1]) k=2; break;
case 2: if(sf_alpha[1]) k=2; else if(sf_alpha[0]) k=1; break;
}
if(k) improveStr(strmen,"OFF","pdf%d",k);
}
improveStr(strmen,"ZZZZ","%.4f", alphaMZ);
improveStr(strmen,"NF","%d",alphaNF);
if(alphaOrder==1) improveStr(strmen,"NNLO","%-.4s","LO");
else if(alphaOrder==2) improveStr(strmen,"NNLO","%-.4s","NLO");
else alphaOrder=3;
improveStr(strmen,"MbMb","%.3f", MbMb);
improveStr(strmen,"Mtp","%.2f", Mtp);
improveStr(strmen,"RRR","%-.16s", Rscale_str);
improveStr(strmen,"FF1","%-.16s", F1scale_str);
improveStr(strmen,"FF2","%-.16s", F2scale_str);
improveStr(strmen,"FFS","%-.16s", Sscale_str);
menu1(54,8,"QCD alpha",strmen,"n_alpha",&pscr,&mode);
}
switch (mode)
{ case 0: if(returnCode) init_alpha(); return returnCode;
case 1:
{ char alphaMen[100]="\006"
" OFF "
" pdf1 "
" pdf2 ";
int k=0;
menu1(54,12,"alpha",alphaMen,"",NULL,&k);
if(k)alphaPDF=k-1;
if(alphaPDF && !sf_alpha[alphaPDF-1]) messanykey(20,20,"WARNING! This pdf does not define alphaQCD ");
}
break;
case 2:
{ double alphaMZ_old=alphaMZ;
if(correctDouble(3,15,"Enter new value ",&alphaMZ,1)) returnCode=1;
if(alphaMZ>0 && alphaMZ<0.3) returnCode=1; else
{ alphaMZ=alphaMZ_old;
messanykey(5,15,"Your input is out of alphaMZ range");
}
}
break;
case 3:
{ int NF_old=alphaNF;
if(correctInt(3,15,"Enter new value ",&alphaNF,1))
{
if(alphaNF<=6 && alphaNF>=3) returnCode=1;
else { messanykey(5,15,"NF out of range"); alphaNF=NF_old;}
}
}
break;
case 4:
{ char lomen[]="\010"
" LO "
" NLO "
" NNLO ";
void *pscrlo=NULL;
int k=0;
menu1(52,12,"",lomen,"",&pscrlo,&k);
if(k) { alphaOrder=k; returnCode=1; put_text(&pscrlo);}
}
break;
case 5: correctDouble(3,15,"Enter new value ",&MbMb,1); break;
case 6: correctDouble(3,15,"Enter new value ",&Mtp,1); break;
case 7:
{ void * screen;
int i;
static double qMin=1, qMax=1000;
static int nPoints=100;
if(returnCode) init_alpha();
get_text(1,1,maxCol(),maxRow(),&screen);
if(correctDouble(40 ,15 ,"Q_min=",&qMin,0)&& qMin>=0.5
&& correctDouble(40 ,16 ,"Q_max=",&qMax,0)&& qMax>qMin
&& correctInt(33,17,"number of points=" ,&nPoints,0)
&& nPoints>3&& nPoints<=150)
{ double *f[3]={NULL,NULL,NULL};
double *ff[3]={NULL,NULL,NULL};
char buff[3][100];
char* Y[3]={buff[0],buff[1],buff[2]};
switch(alphaOrder)
{ case 1: sprintf(Y[0],"MSbar LO"); break;
case 2: sprintf(Y[0],"MSbar NLO"); break;
case 3: sprintf(Y[0],"MSbar NNLO"); break;
default:sprintf(Y[0],"MSbar");
}
int N,k;
f[0]=(double*) malloc(nPoints*sizeof(double));
int xLog= (qMin>0 && qMax/qMin >10)? 1 : 0;
for(i=0;i<nPoints;i++)
{ double z=(i+0.5)/(double)(nPoints),q;
if(xLog) q=pow(qMin,1-z)*pow(qMax,z); else q=qMin*(1-z)+qMax*z;
f[0][i]=alpha_0(q);
// printf("i=%d %E\n",i,f[0][i]);
}
N=1;
for(k=0;k<2;k++) if(sf_alpha[k])
{
char buff[300];
strFunName(k+1,Y[N]);
char *p=strstr(Y[N],"(proton");
if(p) p[0]=0; else
{ p=strstr(Y[N],"(anti-proton");
if(p) p[0]=0;
}
f[N]=(double*) malloc(nPoints*sizeof(double));
for(i=0;i<nPoints;i++)
{ double z=(i+0.5)/(double)(nPoints),q;
if(xLog) q=pow(qMin,1-z)*pow(qMax,z);
else q=qMin*(1-z)+qMax*z;
f[N][i]=(*sf_alpha[k])(q);
}
N++;
}
// printf("N=%d Y[0]=%s\n Y[1]=%s\n Y[2]=%s\n", N,Y[0],Y[1],Y[2]);
// plot_Nar(NULL, "Alpha(Q)", log10(qMin), log10(qMax),"log10(Q/GeV)", nPoints, N, f,ff,Y);
plot_Nar(NULL, "Alpha(Q)", qMin,qMax,"Q/GeV", nPoints, xLog, N, f,ff,Y);
for(k=0;k<N;k++) free(f[k]);
} else messanykey(40,18,
" Correct input is \n"
" 0.5<= Q_min <Q_max\n"
" number of points <=150 and >=4");
put_text(&screen);
} break;
case 8:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Renorm. scale: ");
if(str_redact(Rscale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str, mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
case 9:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Fct1.scale: ");
if(str_redact(F1scale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str, mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
case 10:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Fct1.scale: ");
if(str_redact(F2scale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str, mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
case 11:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Shworing scale: ");
if(str_redact(Sscale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str,mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
}
goto L10;
}
void outputter(SCREEN_Screen* scr, int terminator, GetCharFunction getf)
{
while (1) {
unsigned char c = getf();
if (c == terminator) {
return;
}
switch (c) {
case '\a': break; // TODO: alert?
case '\t': tab(scr); break;
case '\n': cursor_down(scr); break;
case '\b': cursor_left(scr); break;
case '\r': carriage_return(scr); break;
case 0x1b:
{
c = getf();
switch (c) {
case '[':
{
int x;
if (getnum(getf, &x, &c)) {
if (x == 1 && c == 'K') {
clr_bol(scr);
} else {
switch (c) {
case 'm': mode(scr, x); break;
case 'A': parm_up_cursor(scr, x); break;
case 'B': parm_down_cursor(scr, x); break;
case 'C': parm_right_cursor(scr, x); break;
case 'D': parm_left_cursor(scr, x); break;
case 'G': column_address(scr, (x-1)); break;
case 'L': parm_insert_line(scr, x); break;
case 'M': parm_delete_line(scr, x); break;
case 'P': parm_dch(scr, x); break;
case 'S': parm_index(scr, x); break;
case 'T': parm_rindex(scr, x); break;
case 'X': erase_chars(scr, x); break;
case '@': parm_ich(scr, x); break;
case 'd': row_address(scr, (x-1)); break;
case ';':
{
int y;
if (getnum(getf, &y, &c)) {
switch (c) {
case 'm':
mode(scr, x);
mode(scr, y);
break;
case 'H': cursor_address(scr, mkpos(y-1, x-1)); break;
default:
fprintf(stderr, "unhandled: ESC[%i;%i%c\n", x, y, c);
break;
}
} else {
fprintf(stderr, "unhandled: ESC[%i;%c\n", x, c);
}
} break;
default:
fprintf(stderr, "unahndled: ESC[%i%c\n", x, c);
break;
}
}
} else {
switch (c) {
case '@': insert_character(scr); break;
case 'A': cursor_up(scr); break;
case 'B': cursor_down(scr); break;
case 'C': cursor_right(scr); break;
case 'D': cursor_left(scr); break;
case 'H': cursor_home(scr); break;
case 'I': tab(scr); break;
case 'J': clr_eos(scr); break;
case 'K': clr_eol(scr); break;
case 'L': insert_line(scr); break;
case 'M': delete_line(scr); break;
case 'P': delete_character(scr); break;
case 'm': exit_attribute_mode(scr); break;
default:
fprintf(stderr, "unhandled: ESC[%c\n", c);
break;
}
}
} break;
case 'M':
scroll_reverse(scr);
break;
default:
fprintf(stderr, "unhandled: ESC%c\n", c);
break;
}
} break;
default:
{
wchar_t wc;
if ((c & 0x80) == 0x00) {
wc = c;
} else if ((c & 0xE0) == 0xC0) {
char c1 = 0x1F & c;
char c2 = 0x3F & getf();
wc = (c1 << 6) | c2;
} else if ((c & 0xF0) == 0xE0) {
char c1 = 0x0F & c;
char c2 = 0x3F & getf();
char c3 = 0x3F & getf();
wc = (c1 << 12) | (c2 << 6) | c3;
} else if ((c & 0xF8) == 0xF0) {
char c1 = 0x07 & c;
char c2 = 0x3F & getf();
char c3 = 0x3F & getf();
char c4 = 0x3F & getf();
wc = (c1 << 18) | (c2 << 12) | (c3 << 6) | c4;
} else if ((c & 0xFC) == 0xF8) {
char c1 = 0x03 & c;
char c2 = 0x3F & getf();
char c3 = 0x3F & getf();
char c4 = 0x3F & getf();
char c5 = 0x3F & getf();
wc = (c1 << 24) | (c2 << 18) | (c3 << 12) | (c4 << 6) | c5;
} else if ((c & 0xFE) == 0xFC) {
char c1 = 0x01 & c;
char c2 = 0x3F & getf();
char c3 = 0x3F & getf();
char c4 = 0x3F & getf();
char c5 = 0x3F & getf();
char c6 = 0x3F & getf();
wc = (c1 << 30) | (c2 << 24) | (c3 << 18) | (c4 << 12) | (c5 < 6) | c6;
} else {
fprintf(stderr, "bad utf-8 sequence: c=0x%02x\n", c);
wc = '\0';
}
put_char(scr, wc);
} break;
}
}
}
int main(int argc,char** argv)
{
/*===================================
0 - First start.
1 - Model menu.
2 - Enter process menu.
3 - Feynman diagrams menu; squaring.
4 - Squared diagram menu; symbolic calculation.
5 - Write results; new process.
10 -Restart symbolic calculations
==========================================================*/
/* 0-Start; 1-Restart; 2-Heap Error,3-Edit Model,4-UserBreak */
void *pscr1=NULL,*pscr2=NULL,*pscr3=NULL,*pscr4=NULL,*pscr5=NULL;
int k1=1,k2=1,k3=1,k4=1,k5=1;
int n;
int pid=0;
char LOCKtxt[]="Directory 'results/' contains the .lock file created by n_calchep.\n"
"To continue you may a)Close the n_calchep session, or b)Rename 'results/',\n";
blind=0;
strcpy(pathtocomphep,argv[0]);
for(n=strlen(pathtocomphep)-1; n && pathtocomphep[n]!=f_slash; n--);
pathtocomphep[n-3]=0;
for ( n=1;n<argc;n++)
{ if (strcmp(argv[n],"-blind")==0 && n<argc-1 )
{ blind=1;
inkeyString=argv[++n];
}
if (strcmp(argv[n],"+blind")==0 ) blind=2;
}
if(!writeLockFile(".lock"))
{ fprintf(stderr,"locked by other s_calchep. See .lock file\n");
exit(100);
}
strcpy(pathtouser,"");
sprintf(pathtohelp,"%shelp%c",pathtocomphep,f_slash);
outputDir="results/";
{ char * icon=(char *) malloc(strlen(pathtocomphep)+10);
int err;
sprintf(icon,"%sicon",pathtocomphep);
err=start1("CalcHEP/symb",icon,"calchep.ini",&xw_error);
if(err && blind==0)
{ printf("Error:You have launched interactive session for version compiled without X11 library.\n");
printf("Options: a) Use blind session; b) Recompile CalcHEP with X11\n");
exit(66);
}
free(icon);
}
fillModelMenu();
f3_key[0]=f3_key_prog; f3_mess[0]="Model";
f3_key[1]=f4_key_prog; f3_mess[1]="Diagrams";
f3_key[2]=f5_key_prog; f3_mess[2]="Switches";
f3_key[3]=f6_key_prog; f3_mess[3]="Results";
f3_mess[4]="Del";
f3_mess[5]="UnDel";
f3_key[6]=f9_key_prog; f3_mess[6]="Ref";
f3_key[7]=f10_key_prog; f3_mess[7]="Quit";
restoreent(&menulevel);
if(!blind && menulevel<2) cheplabel();
switch (menulevel)
{
case 10:
case 6:
case 5:
case 4:
case 3: readModelFiles("./models",n_model);
modelinfo();
loadModel(0,forceUG);
processinfo();
diagramsinfo();
case 2: k1=n_model;
case 1: break;
}
switch (menulevel)
{
case 2: menuhelp();
goto label_20;
case 3: goto label_31;
case 4: goto label_40;
case 5: goto label_50;
case 6:
case 10: goto restart2;
}
label_10: /* Menu2(ModelMenu): */
f3_key[0]=NULL; /*models*/
f3_key[1]=NULL; /*diagrams*/
menulevel = 1;
forceUG=0;
menuhelp();
for(;;)
{
showheap();
k1=n_model;
menu1(56,4,"",modelmenu,"s_1",&pscr1,&k1);
n_model=k1;
if(n_model == 0)
{
if( mess_y_n(56,4,"Quit session")) { saveent(menulevel); goto exi; }
}
else if(n_model == maxmodel+1)
{
clrbox(1,4,55,18);
makenewmodel();
menuhelp();
}
else if (n_model > 0)
{
put_text(&pscr1);
goto label_20;
}
}
label_20: /* Menu3:Enter Process */
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel = 2;
saveent(menulevel);
readModelFiles("./models",n_model);
modelinfo();
k2 = 1;
do
{ char strmen[]="\026"
" Enter Process "
" Force Unit.Gauge OFF "
" Edit model "
" Delete model ";
if(forceUG)improveStr(strmen,"OFF","ON");
menu1(56,4,"",strmen,"s_2_*",&pscr2,&k2);
switch (k2)
{
case 0: goto_xy(1,1); clr_eol(); goto label_10;
case 2: forceUG=!forceUG; modelinfo(); break;
case 3: editModel(1); break;
case 4:
if ( deletemodel(n_model) )
{
goto_xy(1,1);
clr_eol();
n_model=1;
fillModelMenu();
goto label_10;
}
else readModelFiles("./models",n_model);
}
} while (k2 != 1);
loadModel(0,forceUG);
label_21:
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel=2;
errorcode=enter(); /* Enter a process */
newCodes=0;
showheap();
if (errorcode) /* 'Esc' pressed */ { menuhelp(); goto label_20;}
errorcode=construct(); /* unSquared diagrams */
if (errorcode)
{ if(blind)
{ printf("Processes of this type are absent\n"); sortie(111); }
else
{ messanykey(5,22,"Processes of this type are absent");
clrbox(1,19,80,24);
goto label_21;
}
}
else if(!blind)
{ int dirStat=checkDir("results");
if(dirStat!=0)
{ messanykey( 10,10,"There are files in directory 'results/'.\n"
"To continue you has to clean or rename this directory.");
viewresults();
if(checkDir("results")!=0) goto label_21;
}
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
goto label_31;
}
label_30: /* Menu4: Squaring,...*/
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
label_31:
f3_key[0]=f3_key_prog;
f3_key[1]=NULL;
menulevel=3;
k3 = 1;
do
{
menu1(56,4,"","\026"
" View diagrams "
/* " Amplitude calculation" */
" Squaring technique "
" Write down processes ","s_squa_*",&pscr3,&k3);
switch (k3)
{
case 0: clrbox(1,2,55,11); menuhelp(); goto label_20;
case 1: viewfeyndiag(1); break;
case 3: { FILE*f=fopen("results/list_prc.txt","w");
int k,ndel,ncalc,nrest;
char process[100];
long recpos;
menup=fopen(MENUP_NAME,"r");
for(k=1;;k++)
{ int err=rd_menu(1,k,process,&ndel,&ncalc,&nrest,&recpos);
if(!err) break;
trim(process);
fprintf(f,"%s\n",process);
}
fclose(f);
fclose(menup);
messanykey(20,14,"See file 'results/list_prc.txt'");
} break;
/* case 2: messanykey(10,10,"Not implemented yet"); Amplitudes(); */
}
} while (k3 != 2);
if (!squaring()) goto label_30; /* process is absent */
clear_tmp();
saveent(menulevel);
restoreent(&menulevel);
label_40: /* Menu5: View squared diagrams..... */
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
sq_diagramsinfo(); /* ???????? */
k4=1;
saveent(menulevel);
pscr4=NULL;
for(;;)
{ int res;
menu1(56,4,"","\026"
" View squared diagrams"
" Symbolic calculations"
" Make&Launch n_calchep"
" Make n_calchep "
" REDUCE program "
,"s_calc_*",&pscr4,&k4);
res=checkDir("results");
if(res==1)
{
int n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id) unLockFile(n_calchep_id); else res=2;
}
switch (k4)
{ case 0:
if (mess_y_n(50,3,"Return to previous menu?"))goto label_30;
break;
case 1:
viewsqdiagr(); break;
case 2: /* Compute all diagrams */
restart2:
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
calcallproc();
sq_diagramsinfo();
if(!continuetest()) break;
showheap();
put_text(&pscr4);
break;
case 3:
{ static char keystr[12]="]{{[{";
if(res==2) messanykey(3,10,LOCKtxt); else inkeyString=keystr;
}
break;
case 4:
if(res==2) messanykey(3,10,LOCKtxt); else
{ char command[100];
FILE *f=fopen("results/EXTLIB","w");
fprintf(f,"EXTLIB=\"%s\"\nexport EXTLIB\n",EXTLIB);
fclose(f);
saveent(menulevel);
{ char command[100];
sprintf(command,"cd results; $CALCHEP/bin/make_VandP ../models %d",n_model);
system(command);
}
finish();
sortie(22);
}
case 5: mk_reduceprograms(); break;
}
if(k4==2 && continuetest()) { put_text(&pscr4); break; }
}
label_50:
k5=1;
pscr5=NULL;
menulevel=5;
saveent(menulevel);
for(;;)
{ int n_calchep_id;
menu1(56,4,"","\026"
" C code "
" C-compiler "
" Edit Linker "
" REDUCE code "
" MATHEMATICA code "
" FORM code "
" Enter new process " ,"s_out_*",&pscr5,&k5);
if((k5==1||k5==2)&&pid)
{ int epid=waitpid(pid,NULL,WNOHANG);
if(epid) pid=0; else
{
messanykey(3,10,"This option is frozen while n_calchep runs or is compiled");
continue;
}
}
switch (k5)
{ case 0: goto label_40; break;
case 1:
c_prog(); newCodes=0; saveent(menulevel);
break;
case 2:
if(newCodes) { c_prog(); newCodes=0; saveent(menulevel);}
pid=fork();
if(pid==0)
{
char * command;
n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id)
{
command=(char*)malloc(strlen(pathtocomphep)+strlen(EXTLIB)+100);
sprintf(command,"EXTLIB=\"%s\"\n"
"export EXTLIB\n"
"cd results\n"
"xterm -e %s/make__n_calchep %d",
EXTLIB,pathtocomphep,n_model);
system(command);
sprintf(command,"cd results; ./n_calchep");
unLockFile(n_calchep_id);
system(command);
free(command);
}
exit(0);
}
break;
case 3: if(edittable(1,1,&modelTab[4],1," ",0))
{ char fName[STRSIZ];
readEXTLIB();
sprintf(fName,"%smodels%c%s%d.mdl",pathtouser,f_slash,mdFls[4],n_model);
writetable( &modelTab[4],fName);
}
break;
case 4: makeReduceOutput(); break;
case 5: makeMathOutput(); break;
case 6: makeFormOutput(); break;
case 7:
put_text(&pscr5);
clrbox(1,2,55,11);
menuhelp();
goto label_20;
}
}
exi:
finish();
sortie(0);
return 0;
}
int main(int argc,char** argv)
{
/*===================================
0 - First start.
1 - Model menu.
2 - Enter process menu.
3 - Feynman diagrams menu; squaring.
4 - Squared diagram menu; symbolic calculation.
5 - Write results; new process.
10 -Restart symbolic calculations
==========================================================*/
/* 0-Start; 1-Restart; 2-Heap Error,3-Edit Model,4-UserBreak */
void *pscr1=NULL,*pscr2=NULL,*pscr3=NULL,*pscr4=NULL,*pscr5=NULL;
int k1=1,k2=1,k3=1,k4=1,k5=1;
int n;
int pid=0;
char LOCKtxt[]="Directory 'results/' contains the .lock file created by n_calchep.\n"
"To continue you may a)Close the n_calchep session, or b)Rename 'results/',\n";
blind=0;
strcpy(pathtocalchep,argv[0]);
for(n=strlen(pathtocalchep)-1; n && pathtocalchep[n]!=f_slash; n--);
pathtocalchep[n-3]=0;
for ( n=1;n<argc;n++)
{ if (strcmp(argv[n],"-blind")==0 && n<argc-1 )
{ blind=1;
inkeyString=argv[++n];
}
if (strcmp(argv[n],"+blind")==0 ) blind=2;
if (strcmp(argv[n],"--version")==0 ) { printf("%s\n", VERSION_); exit(0);}
}
if(!writeLockFile(".lock"))
{ fprintf(stderr,"locked by other s_calchep. See .lock file\n");
exit(100);
}
strcpy(pathtouser,"");
sprintf(pathtohelp,"%shelp%c",pathtocalchep,f_slash);
outputDir="results/";
{ char * icon=(char *) malloc(strlen(pathtocalchep)+20);
char title[30];
int err;
sprintf(icon,"%s/include/icon",pathtocalchep);
sprintf(title,"CalcHEP_%s/symb", VERSION);
err=start1(title,icon,"calchep.ini",&xw_error);
if(err && blind==0)
{ printf("Error:You have launched interactive session for version compiled without X11 library.\n");
printf(" Presumably X11 development package is not installed in your computer.\n");
printf(" In this case directory /usr/include/X11/ is empty.\n");
printf("Options: a) Use blind session; b) Update Linux and recompile CalcHEP \n");
printf("Name of needed package\n");
printf(" libX11-devel Fedora/Scientific Linux/CYGWIN/Darwin(MAC)\n");
printf(" libX11-dev Ubuntu/Debian\n");
printf(" xorg-x11-devel SUSE\n");
exit(66);
}
free(icon);
}
fillModelMenu();
f3_key[0]=f3_key_prog; f3_mess[0]="Model";
f3_key[1]=f4_key_prog; f3_mess[1]="Diagrams";
f3_key[2]=f5_key_prog; f3_mess[2]="Switches";
f3_key[3]=f6_key_prog; f3_mess[3]="Results";
f3_mess[4]="Del";
f3_mess[5]="UnDel";
f3_key[6]=f9_key_prog; f3_mess[6]="Ref";
f3_key[7]=f10_key_prog; f3_mess[7]="Quit";
restoreent(&menulevel);
if(!blind && menulevel<2) cheplabel();
switch (menulevel)
{
case 10:
case 6:
case 5:
case 4:
case 3: readModelFiles("./models",n_model);
modelinfo();
loadModel(0,forceUG);
processinfo();
diagramsinfo();
k1=n_model;
break;
case 2: k1=n_model;
readModelFiles("./models",n_model);
break;
}
switch (menulevel)
{
case 2: menuhelp();
goto label_20;
case 3: goto label_31;
case 4: goto label_40;
case 5: goto label_50;
case 6:
case 10: goto restart2;
}
label_10: /* Menu2(ModelMenu): */
f3_key[0]=NULL; /*models*/
f3_key[1]=NULL; /*diagrams*/
menulevel = 1;
forceUG=0;
menuhelp();
for(;;)
{
showheap();
k1=n_model;
menu1(56,4,"",modelmenu,"s_1",&pscr1,&k1);
if(k1 == 0)
{
if( mess_y_n(56,4,"Quit session")) {n_model=0; saveent(menulevel); goto exi; }
}
else if(k1 == maxmodel+1)
{
clrbox(1,4,55,18);
makenewmodel();
menuhelp();
}
else if (k1 > 0)
{ int err=0;
put_text(&pscr1);
if(k1!=n_model || ldModelStatus==0) { err=readModelFiles("./models",k1);}
n_model=k1;
if(err){ if(blind) sortie(133); else goto label_10;} else goto label_20;
}
}
label_20: /* Menu3:Enter Process */
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel = 2;
saveent(menulevel);
modelinfo();
k2 = 1;
do
{ char strmen[]="\026"
" Enter Process "
" Force Unit.Gauge= OFF"
" Edit model "
" Numerical Evaluation "
"======================"
" Delete model ";
if(forceUG)improveStr(strmen,"OFF","ON");
menu1(56,4,"",strmen,"s_2_*",&pscr2,&k2);
switch (k2)
{
case 0: goto_xy(1,1); clr_eol(); goto label_10;
case 2: forceUG=!forceUG; modelinfo(); break;
case 3: editModel(1); break;
case 4: numcheck();
case 5: break;
case 6:
if(deletemodel(n_model))
{
goto_xy(1,1);
clr_eol();
n_model=1;
ldModelStatus=0;
fillModelMenu();
goto label_10;
}
}
} while (k2 != 1);
if(!loadModel(0,forceUG)) goto label_20;
label_21:
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel=2;
checkAuxDir(n_model);
errorcode=enter(); /* Enter a process */
newCodes=0;
showheap();
if (errorcode) /* 'Esc' pressed */ { menuhelp(); goto label_20;}
errorcode=construct(); /* unSquared diagrams */
if (errorcode)
{ if(blind)
{ printf("Processes of this type are absent\n"); sortie(111); }
else
{ messanykey(5,22,"Processes of this type are absent");
clrbox(1,19,80,24);
goto label_21;
}
}
else if(!blind)
{ int dirStat=checkDir("results");
if(dirStat!=0)
{ messanykey( 10,10,"There are files in directory 'results/'.\n"
"To continue you has to clean or rename this directory.");
viewresults();
if(checkDir("results")!=0) goto label_21;
}
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
goto label_31;
}
label_30: /* Menu4: Squaring,...*/
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
label_31:
f3_key[0]=f3_key_prog;
f3_key[1]=NULL;
menulevel=3;
k3 = 1;
do
{
menu1(56,4,"","\026"
" View diagrams "
/* " Amplitude calculation" */
" Square diagrams "
" Write down processes ","s_squa_*",&pscr3,&k3);
switch (k3)
{
case 0: clrbox(1,2,55,11); menuhelp(); goto label_20;
case 1: viewfeyndiag(1); break;
case 3: { FILE*f=fopen("results/list_prc.txt","w");
int k,ndel,ncalc,nrest;
char process[100];
long recpos;
menup=fopen(MENUP_NAME,"r");
for(k=1;;k++)
{ int err=rd_menu(1,k,process,&ndel,&ncalc,&nrest,&recpos);
if(!err) break;
trim(process);
fprintf(f,"%s\n",process);
}
fclose(f);
fclose(menup);
messanykey(20,14,"See file 'results/list_prc.txt'");
} break;
/* case 2: messanykey(10,10,"Not implemented yet"); Amplitudes(); */
}
} while (k3 != 2);
if (!squaring()) goto label_30; /* process is absent */
clear_tmp();
saveent(menulevel);
restoreent(&menulevel);
label_40: /* Menu5: View squared diagrams..... */
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
sq_diagramsinfo(); /* ???????? */
k4=1;
saveent(menulevel);
pscr4=NULL;
for(;;)
{ int res;
menu1(56,4,"","\026"
" View squared diagrams"
" Symbolic calculations"
" Make&Launch n_calchep"
" Make n_calchep "
" REDUCE program "
,"s_calc_*",&pscr4,&k4);
res=checkDir("results");
if(res==1)
{
int n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id) unLockFile(n_calchep_id); else res=2;
}
switch (k4)
{ case 0:
if (mess_y_n(50,3,"Return to previous menu?"))goto label_30;
break;
case 1:
viewsqdiagr(); break;
case 2: /* Compute all diagrams */
restart2:
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
if(!nPROCSS ) calcallproc(); else
{ int *pids=malloc(sizeof(int)*nPROCSS);
int *pipes=malloc(2*sizeof(int)*nPROCSS);
int **qd=malloc(sizeof(int*)*nPROCSS);
int totD=sqDiagList(qd, nPROCSS);
int totC,nProc;
int k;
fflush(NULL);
for(k=0;k<nPROCSS;k++)
{ int* kpipe=pipes+2*k;
pipe(kpipe);
pids[k]=fork();
if(pids[k]==0)
{
close(kpipe[0]);
calcWithFork(k,qd[k],kpipe[1]);
exit(0);
}
}
for(k=0;k<nPROCSS;k++) close(pipes[2*k+1]);
infoLine(0.);
for(nProc=nPROCSS,totC=0;nProc;)
{ int one;
int err;
nProc=0;
for(k=0;k<nPROCSS;k++) if(pids[k])
{ if(waitpid(pids[k],NULL,WNOHANG)==0)
{ nProc++;
if(read(pipes[2*k],&one,sizeof(int)))totC+=one;
} else
{ for(;read(pipes[2*k],&one,sizeof(int))>0;) totC+=one;
pids[k]=0;
}
}
if(infoLine((double)totC/(double)totD)) for(k=0;k<nPROCSS;k++) if(pids[k]) kill(pids[k],SIGUSR1);
}
infoLine(2);
for(k=0;k<nPROCSS;k++)
{ char ctlgName[100];
char command[200];
sprintf(ctlgName,"%s_%d",CATALOG_NAME,k);
if(access(ctlgName,R_OK)==0)
{ sprintf(command," cat %s >> %s", ctlgName,CATALOG_NAME);
system(command);
unlink(ctlgName);
}
}
if(totC) newCodes=1;
updateMenuQ();
for(k=0;k<nPROCSS;k++) close(pipes[2*k]);
free(pids);
free(pipes);
for(k=0;k<nPROCSS;k++) free(qd[k]);
free(qd);
}
sq_diagramsinfo();
if(!continuetest()) break;
showheap();
put_text(&pscr4);
break;
case 3:
{ static char keystr[12]="]{{[{";
if(res==2) messanykey(3,10,LOCKtxt); else inkeyString=keystr;
}
break;
case 4:
if(res==2) messanykey(3,10,LOCKtxt); else
{ char command[100];
saveent(menulevel);
chdir("results");
makeVandP(0,"../models",n_model,2,pathtocalchep);
finish();
sortie(22);
}
case 5: mk_reduceprograms(); break;
}
if(k4==2 && continuetest()) { put_text(&pscr4); break; }
}
label_50:
k5=1;
pscr5=NULL;
menulevel=5;
saveent(menulevel);
sq_diagramsinfo();
for(;;)
{ int n_calchep_id;
menu1(56,4,"","\026"
" C code "
" C-compiler "
" Edit Linker "
" REDUCE code "
" MATHEMATICA code "
" FORM code "
" Enter new process " ,"s_out_*",&pscr5,&k5);
if((k5==1||k5==2)&&pid)
{ int epid=waitpid(pid,NULL,WNOHANG);
if(epid) pid=0; else
{
messanykey(3,10,"This option is frozen while n_calchep runs or is compiled");
continue;
}
}
switch (k5)
{ case 0: goto label_40; break;
case 1:
c_prog(); newCodes=0; saveent(menulevel);
break;
case 2:
if(newCodes) { c_prog(); newCodes=0; saveent(menulevel);}
n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id)
{
if(nPROCSS)
{ chdir("results");
makeVandP(0,"../models",n_model,2,pathtocalchep);
pCompile();
if(access("./n_calchep",X_OK)==0)
{
fflush(NULL);
pid=fork();
if(pid==0)
{
system("./n_calchep");
exit(0);
}
} else messanykey(15,15,"n_calchep is not generated");
chdir("../");
} else
{ fflush(NULL);
pid=fork();
if(pid==0)
{ if(chdir("results")==0)
{ char*command=malloc(100+strlen(pathtocalchep));
makeVandP(0,"../models",n_model,2,pathtocalchep);
sprintf(command,"xterm -e %s/sbin/make__n_calchep %d", pathtocalchep,n_model);
system(command);
free(command);
system("./n_calchep");
}
exit(0);
}
}
unLockFile(n_calchep_id);
}
break;
case 3: if(edittable(1,1,&modelTab[4],1," ",0))
{ char fName[STRSIZ];
sprintf(fName,"%smodels%c%s%d.mdl",pathtouser,f_slash,mdFls[4],n_model);
writetable( &modelTab[4],fName);
}
break;
case 4: makeReduceOutput(); break;
case 5: makeMathOutput(); break;
case 6: makeFormOutput(); break;
case 7:
put_text(&pscr5);
clrbox(1,2,55,11);
menuhelp();
goto label_20;
}
}
exi:
finish();
sortie(0);
return 0;
}
static void prtcllist(int key)
{
char fullname[STRSIZ];
char hlp[60];
char p1[60], p2[60];
int i, j, pnum;
linelist ln;
int tabMax,tabSz;
static int nTot,nFirst;
tabMax=ycons -7;
if (key==0)
{
scrcolor(FGmain,BGmain);
for (i = 2; i <= 24; i++)
{ goto_xy(1,i);
clr_eol();
}
goto_xy(14,3);
scrcolor(Blue,BGmain);
print("List of particles (antiparticles)");
nTot=0;
nFirst=1;
}
else
{
if (nTot <= 3 *tabMax ) return;
switch (key)
{
case KB_DOWN : nFirst+=3; break;
case KB_UP : nFirst-=3; break;
case KB_PAGED: nFirst +=3*tabMax; break;
case KB_PAGEU: nFirst -=3*tabMax; break;
}
if (nFirst <1) nFirst=1;
if (nTot-nFirst+3<3*tabMax ) nFirst=1+3*((nTot+2)/3) -3*tabMax;
clrbox(1,4,79,5+tabMax);
}
goto_xy(3,5); scrcolor(FGmain,BGmain);
for(i=0,ln=prtcls_tab.strings; ln; ln=ln->next )
{ sscanf(ln->line,"%[^|]%*c%[^|]%*c%[^|]%*c%*[^|]%*c%*[^|]%*c%*[^|]%*c%*[^|]%*c%*[^|]%*c%[^|]",
fullname,p1,p2,hlp);
trim(p1);
locateinbase(p1,&pnum);
trim(hlp);
if (prtclbase[pnum-1].top != NULL && strcmp(hlp,"*") != 0)
{
i++;
if (i>=nFirst && (i-nFirst)/3 <tabMax )
{
print("%s",p1);
if (strcmp(p1,p2) == 0) print(" "); else print("(%s)",p2);
trim(fullname);
print("- %s",fullname);
j = i % 3;
if (j == 0) goto_xy(3,where_y() + 1);
else goto_xy(3 + 26 * j,where_y());
}
}
}
nTot=i;
tabSz=MIN((nTot+2)/3,tabMax);
chepbox(1,4,79,5+tabSz);
if (nFirst >1 ) { goto_xy(72,4); print("PgUp"); }
if (nFirst+3*tabSz <= nTot ) { goto_xy(72,5+tabMax); print("PgDn"); }
scrcolor(FGmain,BGmain);
}
int enter(void)
{
int i, y0,scat;
int redres;
char ** items=NULL;
char * errpos=NULL;
char * arrpos=NULL;
int curh=0;
prtcllist(0);
scrcolor(Red,BGmain);
y0 = ycons;
errTxt=NULL;
label_1:
if(y0<maxRow()) y0++;
for(;y0>ycons;y0--){goto_xy(1,y0); clr_eol();}
if(arrpos) strncpy(arrpos,"->",2);
redres=input(y0,"s_ent_1", "Enter process: ",processch,errpos?errpos-processch+1:1,SSTRLEN);
switch (redres)
{
case KB_F1: /* Help */
show_help("s_ent_1");
goto label_1;
case KB_ESC:
clrbox(1,2,maxCol(),24);
return 1;
} /* Case */
arrpos=strstr(processch,"->");
if(arrpos) strncpy(arrpos,", ",2);
else { errTxt="'->' is absent "; goto label_1; }
if( strchr(processch,',')<arrpos ) scat=1; else scat=0;
if(items) free(items);
items=stritems(" ,",processch);
for(i=0,nin=0,nout=0,n_x=0,curh=0; items[i]; i++)
{ char name[100];
sscanf(items[i],"%[^, ]", name);
if(strlen(name)>7) {errTxt="Too long name"; break;}
if(strlen(name)==0){errTxt="Empty item"; break;}
if(items[i]>arrpos && strlen(name) == 3 && name[1] == '*' &&
(name[2] == 'X' || name[2] == 'x') &&
isdigit(name[0])) {n_x += name[0]-'0'; nout +=name[0]-'0';}
else
{
if(curh==MAXINOUT-1){errTxt="Too many particles"; break;}
if(items[i]>arrpos) nout++; else nin++;
if(nin>2){errTxt="Too many incoming particles"; break;}
if(enter_h(&y0,name,curh,scat)) break;
curh++;
}
}
errpos=items[i];
free(items); items=NULL;
if(errpos && !errTxt)
{ sprintf(err_Txt,"wrong definition of %d-th particle",i+1);
errTxt=err_Txt;
}
if(errTxt) goto label_1;
if(nout<2)
{ errTxt="The number of outgoing particles should exceed 1";
goto label_1;
}
if(nin<1)
{ errTxt="Incoming particle(s) is not defined";
goto label_1;
}
if(errpos||nin+nout>MAXINOUT)
{if(!errTxt) errTxt=errtxt; goto label_1; }
y0++;
if(y0>=maxRow()-1) {y0=maxRow()-1; goto_xy(1,y0); clr_eol();}
if(restrict_p(y0,1, "Exclude diagrams with ","s_ent_4",limpch,liminsp)
) goto label_1;
for(nilprtcl(LimQ), i=0;i<liminsp[i].who;i++) if(liminsp[i].how <0)
{
liminsp[i].how =-liminsp[i].how;
addlim(LimQ,liminsp[i].who,liminsp[i].how,0);
}
y0++;
if(y0>=maxRow()-1) {y0=maxRow()-1; goto_xy(1,y0); clr_eol();}
if(n_x && restrict_p(y0,0,"Exclude X-particles ","s_ent_5",deloutch,limout)
) goto label_1;
for(i=nin+nout;i<MAXINOUT;i++) hadrons[i].name[0]=0;
strncpy(arrpos,"->",2);
return 0;
}
static int enter_h(int * y,char* name,int num,int scat)
{ int i,m,j=0;
int redres;
shortstr hadrch;
char ** items;
char * errpos=NULL;
locateinbase(name,&j);
if(j)
{
if(pseudop(j)) return -1;
strcpy(hadrons[num].name,name);
strcpy(hadrons[num].contents,name);
hadrons[num].parton[0] = j;
hadrons[num].pow = 1;
if(!scat && num==0 && strcmp(prtclbase1[j].massidnt,"0")==0)
{ errTxt="Decay of massless particle.";
return -1;
}
if(name[strlen(name)-1]=='%')
{
if( !scat || nout||strcmp(prtclbase1[j].massidnt,"0") ||
(prtclbase1[j].spin!=1 && prtclbase1[j].spin!=2) ||
strchr("LR",prtclbase1[j].hlp)
)
{ errTxt="This particle can not be polarized";
return -1;
} else hadrons[num].polarized[0]=1;
} else hadrons[num].polarized[0]=0;
return 0;
}
for(i=0;i<num;i++)
if(!strcmp(hadrons[i].name,name))
{ strcpy(hadrons[num].name,name);
strcpy(hadrons[num].contents,hadrons[i].contents);
hadrons[num].pow=hadrons[i].pow;
for(j=0;j<hadrons[num].pow;j++)
{ hadrons[num].parton[j]=hadrons[i].parton[j];
if(nout) hadrons[num].polarized[j]=0;
else hadrons[num].polarized[j]=hadrons[i].polarized[j];
}
return 0;
}
hadrch[0]=0;
for(i=num;i<MAXINOUT;i++) if(!strcmp(hadrons[i].name,name))
{
strcpy(hadrch,hadrons[i].contents);
break;
}
if(*y>=maxRow()-1) { goto_xy(1,*y); clr_eol();} else (*y)++;
do
{
m=errpos? errpos-hadrch+1: 0;
do
{ char direction[100];
sprintf(direction,"composit '%s' consists of: ",name);
redres=input(*y, "s_ent_2", direction, hadrch, m , SSTRLEN);
if(redres==KB_ESC) return 1;
} while (redres!=KB_ENTER && redres!=KB_ESC);
if (redres == KB_ESC || strcmp(hadrch,"") == 0) return 1;
items=stritems(" ,",hadrch);
for(m=0,hadrons[num].pow=0; items[m]; m++)
{ char name[100];
if(hadrons[num].pow>=100) {errTxt="too many partons";break;}
sscanf(items[m],"%[^ ,]",name);
locateinbase(name,&j);
if (j==0 || pseudop(j))
{ errTxt= "This particle is absent in the model"; break;}
if(!scat && num==0 && strcmp(prtclbase1[j].massidnt,"0")==0)
{ errTxt="Decay of massless particle.";
break;
}
if(name[strlen(name)-1]=='%')
{
if( !scat || nout||strcmp(prtclbase1[j].massidnt,"0") ||
(prtclbase1[j].spin!=1 && prtclbase1[j].spin!=2) ||
strchr("LR",prtclbase1[j].hlp)
)
{ errTxt="This particle can not be polarized";
break;
} else hadrons[num].polarized[hadrons[num].pow]=1;
} else hadrons[num].polarized[hadrons[num].pow]=0;
hadrons[num].parton[hadrons[num].pow++]=j;
}
errpos=items[m];
if(!errpos)
{ for(i=0;i<hadrons[num].pow;i++)
for(j=i+1;j<hadrons[num].pow;j++)
if(hadrons[num].parton[i]==hadrons[num].parton[j])
{ errpos=items[j];
errTxt="duplicate parton";
}
}
strcpy(hadrons[num].name,name);
strcpy(hadrons[num].contents,hadrch);
free(items);
}
while(errpos);
return 0;
}
int viewresults(void)
{
int k,kmenu;
void * pscr = NULL;
shortstr newname;
int dirStat=checkDir("results");
if(dirStat==0){messanykey(10,15,"directory RESULTS is empty"); return 1;}
kmenu = 1;
label_1:
menu1(10,10,"","\010"
" View "
" Delete "
" Rename ","s_res",&pscr,&kmenu);
switch (kmenu)
{
case 0:
return 0;
case 1:
viewDir("results");
break;
case 2:
if(dirStat==2)
{ char mess[]="Can not clean dir 'results' because it contains the LOCK file";
if(blind) { printf("%s\n",mess); sortie(102);}
else { messanykey(3,13,mess); break;}
}
if ( mess_y_n( 6,13," Delete files ") ) system("rm -r results; mkdir results");
put_text(&pscr);
return 1;
case 3:
strcpy(newname," ");
while(1)
{ void * pscr3;
get_text(1,maxRow(),maxCol(),maxRow(),&pscr3);
goto_xy(1,maxRow());
print("Enter new name: ");
k = str_redact(newname,1,30);
if (k == KB_ESC)
{ goto_xy(1,24);
clr_eol();
goto label_1;
}
if (k == KB_ENTER)
{
trim(newname);
if(rename("results",newname)==0)
{
mkdir("results",-1);
put_text(&pscr);
put_text(&pscr3);
return 1;
}
else messanykey(10,15," Can't rename the directory");
}
put_text(&pscr3);
}
}
goto label_1;
}
