Subversion Repositories wimsdev

Rev

Rev 3879 | Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
10 reyssat 1
/*    Copyright (C) 1998-2003 XIAO, Gang of Universite de Nice - Sophia Antipolis
2
 *
3
 *  This program is free software; you can redistribute it and/or modify
4
 *  it under the terms of the GNU General Public License as published by
5
 *  the Free Software Foundation; either version 2 of the License, or
6
 *  (at your option) any later version.
7
 *
8
 *  This program is distributed in the hope that it will be useful,
9
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
10
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
11
 *  GNU General Public License for more details.
12
 *
13
 *  You should have received a copy of the GNU General Public License
14
 *  along with this program; if not, write to the Free Software
15
 *  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
16
 */
17
 
18
        /* Interface gp to wims */
19
 
20
/* This program is now limited to pari version 2.??. */
21
 
22
/*************** Customization: change values hereafter ****************/
23
 
24
        /* gp prompt string */
25
#define gpprompt "\n? "
26
        /* This is the good bye string of gp, signaling end of output. */
27
#define goodbyestring "Good bye!"
28
        /* This is PARI home page. To be kept up to date. */
29
#define homepage "http://pari.math.u-bordeaux.fr/"
30
        /* String to tell the program to quit. */
31
#define quitstring "\nquit\n"
32
        /* The way to print a string in the program. */
33
#define stringprinter "print(\"%s\")\n"
34
        /* limit of input/output file sizes */
35
int fsizelim=131072;
36
int precision=28;       /* default */
37
 
38
char *nameofcmd="gp";
39
char *gprcenv="GPRC";
40
char *gprc="../.gprc";
41
char header[]="default(output,0)\n\
42
alias(ln,log)\n\
43
alias(pi,Pi)\n\
44
alias(euler,Euler)\n\
45
alias(PI,Pi)\n\
46
i=I\n\
47
e=exp(1)\n\
48
E=exp(1)\n\
49
{ sec(x)=return(1/cos(x)) }\n\
50
{ csc(x)=return(1/sin(x)) }\n\
51
{ lg(x) =log(x)/log(10) }\n\
52
{ log2(x) =log(x)/log(10) }\n\
53
alias(log10,lg)\n\
54
alias(sgn,sign)\n\
55
alias(ch,cosh)\n\
56
alias(sh,sinh)\n\
57
alias(th,tanh)\n\
58
alias(arccos,acos)\n\
59
alias(arcsin,asin)\n\
60
alias(tg,tan)\n\
61
alias(arctan,atan)\n\
62
alias(arctg,atan)\n\
63
alias(Argch,acosh)\n\
64
alias(Argsh,asinh)\n\
65
alias(Argth,atanh)\n\
66
alias(argch,acosh)\n\
67
alias(argsh,asinh)\n\
68
alias(argth,atanh)\n\
69
alias(cot,cotan)\n\
70
alias(ctg,cotan)\n\
71
alias(rint,round)\n\
72
alias(RANDOM,random)\n\
73
alias(SOLVE,solve)\n\
74
alias(parirandom,random)\n\
75
alias(parisolve,solve)\n\
76
";
77
 
78
struct {
79
    char *wname;    char *defaultval;    char *gpset;
80
} setups[]={
81
      {"w_pari_precision",      "20",   "\\p "},
82
      {"w_pari_serieslength",   "8",    "\\ps "}
83
};
84
 
85
        /* names which are not allowed */
86
char *illegal[]={
87
        "alias","allocatemem","default","extern",
88
      "getrand","getstack","input","psdraw",
89
      "read","system","install"
90
};
91
 
92
        /* name parts which are not allowed */
93
char *illpart[]={
94
    "plot", "write", "help"
95
};
96
 
97
/***************** Nothing should need change hereafter *****************/
98
 
99
#define progname "pari"
100
#include "common.h"
101
#include "common.c"
102
 
103
int pariray=0;
104
 
105
        /* check for security violations in command string */
106
void check_parm(char *pm)
107
{
108
    char *p;
109
    for(p=pm;*p!=0;p++) {
110
          /* Underscore replacement */
111
        if(*p=='_') {*p='K'; continue;}
112
          /* no escape commands. */
113
        if(*p!='\n') continue;
114
        while(*p!=0 && isspace(*p)) p++;
115
        if(*p=='\\' && *(p+1)!='\\' && *(p+1)!='v' && *(p+1)!='p') *p='.';
116
        if(*p=='\\') p++;
117
    }
118
    find_illegal(pm);
119
}
120
 
121
        /* process and print gp output */
122
void output(char *p)
123
{
124
    int i,n;
125
    char *pp, *pe, *pt;
126
 
127
    if(pariray) {puts(p); return;}
128
    pp=strstr(p,gpprompt); if(pp==NULL) return;
129
    pe=strstr(pp,goodbyestring);
130
    if(pe>=pp) *pe=0;
131
    while(pp!=NULL) {
132
        pp++;
133
        pe=strstr(pp,gpprompt);
134
        if(pe>=pp) *pe=0;
135
        pp=strchr(pp,'\n');
136
        if(pp==NULL) {
137
            emptyline:
138
            puts(""); pp=pe; continue;
139
        }
140
        pp++; n=strlen(pp);
141
        if(n==0) goto emptyline;
142
                /* make every output one-line */
143
        for(i=0;i<n;i++) {
144
            if(*(pp+i)=='\n') {
145
                if(*(pp+i+1)!='%') *(pp+i)=' ';
146
                else {*(pp+i)=0; break;}
147
            }
148
        }
149
          /* strip leading and trailing spaces */
150
        while(isspace(*pp) && pp<pe) pp++;
151
        pt=pp+strlen(pp)-1;
152
        while(isspace(*pt) && pt>pp) *pt--=0;
153
          /* remove parentheses of matrix output */
154
        if(memcmp(pp,"Mat(",4)==0 && *pt==')' && find_matching(pp+4,')')==pt) {
155
            *(pt--)=0; pp+=4;
156
        }
157
        if(*pp=='[' && *pt==']' && find_matching(pp+1,']')==pt) {
158
            *(pt--)=0; pp++;
159
        }
160
        strip_zeros(pp);
161
        puts(pp); pp=pe;
162
    }
163
}
164
 
165
void about(void)
166
{
167
    char *p;
168
 
169
    prepabout("\\v\nquit\n",outputfname,NULL);
170
    if(readabout()>0) {
171
        p=strchr(aboutbuf,'\n'); if(p!=NULL) *p=0;
172
        strip_trailing_spaces(aboutbuf);
173
        printf("<A HREF=\"%s\">%s</A>",homepage,aboutbuf);
174
    }
175
}
176
 
177
char *dynsetup(char *ptr, char *end)
178
{
179
    int i;
180
    char *p, *pp;
181
 
182
    snprintf(ptr,end-ptr,"\nsetrand(%u)\n",seed&(0x7FFFFFFF));
183
    ptr+=strlen(ptr);
184
    for(i=0;i<SETUP_NO;i++) {
185
        p=getenv(setups[i].wname);
186
        if(p!=NULL) for(pp=p;*pp;pp++) if(!isspace(*pp) && !isalnum(*pp)) p="";
187
        if(p==NULL || *p==0) p=setups[i].defaultval;
188
        snprintf(ptr,end-ptr,"%s%s\n",setups[i].gpset,p);
189
        ptr+=strlen(ptr);
190
        if(strstr(setups[i].wname,"pari_precision")!=NULL)
191
          precision=atoi(p);
192
        if(precision<0) precision=-precision;
193
    }
194
    if(pariray) snprintf(ptr,end-ptr,"\\e0\n");
195
    else snprintf(ptr,end-ptr,"\\e1\n");
196
    ptr+=strlen(ptr);
197
    return ptr;
198
}
199
 
200
int main(int argc,char *argv[])
201
{
202
    char *p;
203
 
204
    p=getenv("pari_ray");
205
    if(p!=NULL && *p!=0) {pariray=1; fsizelim=4*1024*1024;}
206
    setenv(gprcenv,gprc,1);
207
    prepare1();
208
    setenv("GPTMPDIR",tmp_dir,1);
209
    run();
210
    return 0;    
211
}
212