/* Copyright (C) 1998-2003 XIAO, Gang of Universite de Nice - Sophia Antipolis
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
*/
/* Interface gp to wims */
/* This program is now limited to pari version 2.??. */
/*************** Customization: change values hereafter ****************/
#include "common.h"
/* gp prompt string */
#define gpprompt "\n? "
/* This is the good bye string of gp, signaling end of output. */
#define goodbyestring "Good bye!"
/* This is PARI home page. To be kept up to date. */
#define homepage "http://pari.math.u-bordeaux.fr/"
/* String to tell the program to quit. */
char *quitstring="\nquit\n";
/* The way to print a string in the program. */
char *stringprinter="print(\"%s\")\n";
/* limit of input/output file sizes */
int fsizelim=131072;
int precision=28; /* default */
char *nameofcmd="gp";
char *gprcenv="GPRC";
char *gprc="../.gprc";
char header[]="default(output,0)\n\
default(secure,1)\n\
alias(ln,log)\n\
alias(pi,Pi)\n\
alias(euler,Euler)\n\
alias(PI,Pi)\n\
i=I\n\
e=exp(1)\n\
E=exp(1)\n\
sec(x)=1/cos(x);\n\
csc(x)=1/sin(x);\n\
lg(x)=log(x)/log(10);\n\
log2(x)=log(x)/log(10);\n\
alias(log10,lg)\n\
alias(sgn,sign)\n\
alias(ch,cosh)\n\
alias(sh,sinh)\n\
alias(th,tanh)\n\
alias(arccos,acos)\n\
alias(arcsin,asin)\n\
alias(tg,tan)\n\
alias(arctan,atan)\n\
alias(arctg,atan)\n\
alias(Argch,acosh)\n\
alias(Argsh,asinh)\n\
alias(Argth,atanh)\n\
alias(argch,acosh)\n\
alias(argsh,asinh)\n\
alias(argth,atanh)\n\
alias(cot,cotan)\n\
alias(ctg,cotan)\n\
alias(rint,round)\n\
alias(RANDOM,random)\n\
alias(SOLVE,solve)\n\
alias(parirandom,random)\n\
alias(parisolve,solve)\n\
";
struct {
char *wname; char *defaultval; char *gpset;
} setups[]={
{"w_pari_precision", "20", "\\p "},
{"w_pari_serieslength", "8", "\\ps "}
};
/* names which are not allowed */
char *illegal[]={
};
int illegal_no=(sizeof(illegal)/sizeof(illegal[0]));
/* name parts which are not allowed */
char *illpart[]={
"plot", "write", "help"
};
int illpart_no=(sizeof(illpart)/sizeof(illpart[0]));
/***************** Nothing should need change hereafter *****************/
char *progname="pari";
int pariray=0;
/* check for security violations in command string */
void check_parm(char *pm)
{
char *p;
for(p=pm;*p!=0;p++) {
/* Underscore replacement */
if(*p=='_') {*p='K'; continue;}
/* no escape commands. */
if(*p!='\n') continue;
if(*p=='\\' && *(p+1)!='\\' && *(p+1)!='v' && *(p+1)!='p') *p='.';
if(*p=='\\') p++;
}
find_illegal(pm);
}
/* process and print gp output */
void output(char *p)
{
int i,n;
char *pp, *pe, *pt;
if(pariray
) {puts(p
); return;}
pp
=strstr(p
,gpprompt
); if(pp
==NULL
) return;
if(pe>=pp) *pe=0;
while(pp!=NULL) {
pp++;
if(pe>=pp) *pe=0;
if(pp==NULL) {
emptyline:
puts(""); pp
=pe
; continue;
}
if(n==0) goto emptyline;
/* make every output one-line */
for(i=0;i<n;i++) {
if(*(pp+i)=='\n') {
if(*(pp+i+1)!='%') *(pp+i)=' ';
else {*(pp+i)=0; break;}
}
}
/* strip leading and trailing spaces */
while(isspace(*pt
) && pt
>pp
) *pt
--=0;
/* remove parentheses of matrix output */
if(memcmp(pp
,"Mat(",4)==0 && *pt
==')' && find_matching2
(pp
+4,')')==pt
) {
*(pt--)=0; pp+=4;
}
if(memcmp(pp
,"Vecsmall(",9)==0 && *pt
==')' && find_matching2
(pp
+9,')')==pt
) {
*(pt--)=0; pp+=9;
}
if(*pp=='[' && *pt==']' && find_matching2(pp+1,']')==pt) {
*(pt--)=0; pp++;
}
strip_zeros(pp);
}
}
void about(void)
{
char *p;
prepabout("\\v\nquit\n",outputfname,NULL);
if(readabout()>0) {
p
=strchr(aboutbuf
,'\n'); if(p
!=NULL
) *p
=0;
strip_trailing_spaces2(aboutbuf);
printf("<a target=\"wims_external\" href=\"%s\">%s</a>",homepage
,aboutbuf
);
}
}
char *dynsetup(char *ptr, char *end)
{
int i;
char *p, *pp;
if(wseed
!= NULL
) snprintf(ptr
,end
-ptr
,"\nsetrand(%u)\n",atoi(wseed
)&(0x7FFFFFFF));
else snprintf(ptr
,end
-ptr
,"\nsetrand(%u);\n",seed
&(0x7FFFFFFF));
for(i=0;i<SETUP_NO;i++) {
if(p
!=NULL
) for(pp
=p
;*pp
;pp
++) if(!isspace(*pp
) && !isalnum(*pp
)) p
="";
if(p==NULL || *p==0) p=setups[i].defaultval;
snprintf(ptr
,end
-ptr
,"%s%s\n",setups
[i
].
gpset,p
);
if(strstr(setups
[i
].
wname,"pari_precision")!=NULL
)
if(precision<0) precision=-precision;
}
if(pariray
) snprintf(ptr
,end
-ptr
,"\\e0\n");
return ptr;
}
int main(int argc,char *argv[])
{
char *p;
if(p!=NULL && *p!=0) {pariray=1; fsizelim=4*1024*1024;}
setenv(gprcenv,gprc,1);
prepare1();
setenv("GPTMPDIR",tmp_dir,1);
run();
return 0;
}