1 #include "G__ci.h" /* G__atpause is defined in G__ci.h */
3 #if (!defined(G__CINTVERSION)) || (G__CINTVERSION < 501438)
4 #error You need at least cint 5.14.38 to compile GiNaC-cint. Download it via http from root.cern.ch/root/Cint.html or via ftp from ftpthep.physik.uni-mainz.de/pub/cint
5 #endif // (!defined(G__CINTVERSION)) || (G__CINTVERSION < 501438)
12 #include "ginac/ginac.h"
16 #ifndef NO_NAMESPACE_GINAC
17 using namespace GiNaC;
18 #endif // ndef NO_NAMESPACE_GINAC
20 extern "C" G__value G__exec_tempfile G__P((char *file));
21 extern "C" void G__store_undo_position(void);
23 #define PROMPT1 "GiNaC> "
26 #ifdef OBSCURE_CINT_HACK
31 string ToString(const T & t)
34 ostrstream(buf,sizeof(buf)) << t << ends;
38 basic * ex::last_created_or_assigned_bp = 0;
39 basic * ex::dummy_bp = 0;
40 long ex::last_created_or_assigned_exp = 0;
42 #endif // def OBSCURE_CINT_HACK
44 G__value exec_tempfile(string const & command);
45 char * process_permanentfile(string const & command);
46 void process_tempfile(string const & command);
48 void helpmessage(void);
49 string preprocess(char const * const line, bool & comment, bool & single_quote,
50 bool & double_quote, unsigned & open_braces);
52 void sigterm_handler(int n);
53 void initialize(void);
54 void initialize_cint(void);
56 bool is_command(string const & command, string & preprocessed,
57 string const & comparevalue, bool substr=false);
58 bool readlines(istream * is, string & allcommands);
59 bool readfile(string const & filename, string & allcommands);
60 void savefile(string const & filename, string const & allcommands);
62 typedef list<char *> cplist;
64 bool redirect_output=false;
66 G__value exec_tempfile(string const & command)
69 char *tmpfilename = tempnam(NULL,"ginac");
71 fout.open(tmpfilename);
72 fout << "{" << endl << command << endl << "}" << endl;
74 G__store_undo_position();
75 retval = G__exec_tempfile(tmpfilename);
76 G__security_recover(stdout);
82 char * process_permanentfile(string const & command)
84 char *tmpfilename = tempnam(NULL,"ginac");
85 cout << "creating file " << tmpfilename << endl;
87 fout.open(tmpfilename);
88 fout << command << endl;
90 G__store_undo_position();
91 G__loadfile(tmpfilename);
92 G__security_recover(stdout);
96 void process_tempfile(string const & command)
98 #ifdef OBSCURE_CINT_HACK
99 static G__value ref_symbol = exec_tempfile("symbol ginac_cint_internal_symbol; ginac_cint_internal_symbol;");
100 static G__value ref_constant = exec_tempfile("constant ginac_cint_internal_constant; ginac_cint_internal_constant;");
101 static G__value ref_function = exec_tempfile("sin(ginac_cint_internal_symbol);");
102 static G__value ref_power = exec_tempfile("power(ex(ginac_cint_internal_symbol),ex(ginac_cint_internal_symbol));");
103 static G__value ref_numeric = exec_tempfile("numeric ginac_cint_internal_numeric; ginac_cint_internal_numeric;");
104 static G__value ref_ex = exec_tempfile("ex ginac_cint_internal_ex; ginac_cint_internal_ex;");
105 static bool basic_type_warning_already_displayed = false;
106 #endif // def OBSCURE_CINT_HACK
108 G__value retval = exec_tempfile(command);
110 #ifdef OBSCURE_CINT_HACK
112 #define TYPES_EQUAL(A,B) (((A).type==(B).type) && ((A).tagnum==(B).tagnum))
114 static unsigned out_count = 0;
115 if (TYPES_EQUAL(retval,ref_ex)) {
116 string varname = "Out"+ToString(++out_count);
117 if (retval.obj.i!=ex::last_created_or_assigned_exp) {
118 // an ex was returned, but this is not the ex which was created last
119 // => this is not a temporary ex, but one that resides safely in memory
121 // cout << "warning: using ex from retval (experimental)" << endl;
122 ex::dummy_bp=((ex *)(void *)(retval.obj.i))->bp;
123 exec_tempfile("ex "+varname+"(*ex::dummy_bp);");
124 } else if (ex::last_created_or_assigned_bp_can_be_converted_to_ex()) {
126 //for (int i=4-int(log10(out_count)); i>0; --i)
128 exec_tempfile("ex "+varname+"(*ex::last_created_or_assigned_bp);");
130 cout << "warning: last_created_or_assigned_bp modified 0 or not evaluated or not dynallocated" << endl;
132 exec_tempfile(string()+"LLLAST=LLAST;\n"
134 +"LAST="+varname+";\n"
135 +"if (ginac_cint_internal_redirect_output&&"
136 +" ginac_cint_internal_fout.good()) {"
137 +" ginac_cint_internal_fout << \""+varname+" = \" << "+varname+" << endl << endl;"
139 +" cout << \""+varname+" = \" << "+varname+" << endl << endl;"
141 } else if (TYPES_EQUAL(retval,ref_symbol)||
142 TYPES_EQUAL(retval,ref_constant)||
143 TYPES_EQUAL(retval,ref_function)||
144 TYPES_EQUAL(retval,ref_power)||
145 TYPES_EQUAL(retval,ref_numeric)) {
146 if (!basic_type_warning_already_displayed) {
148 <<"WARNING: The return value of the last expression you entered was a symbol," << endl
149 << "constant, function, power or numeric, which cannot be safely displayed." << endl
150 << "To force the output, cast it explicitly to type 'ex' or use 'cout'," << endl
151 << "for example (assume 'x' is a symbol):" << endl
152 << PROMPT1 "ex(x);" << endl
153 << "OutX = x" << endl << endl
154 << PROMPT1 "cout << x << endl;" << endl
155 << "x" << endl << endl
156 << "This warning will not be shown again." << endl;
157 basic_type_warning_already_displayed = true;
160 #endif // def OBSCURE_CINT_HACK
163 //----+----#----+----#----+----#----+----#----+----#----+----#----+----#----+----#
164 // __, _______ Cint: Copyright 1995-2000 Masaharu Goto and Agilent Technologies, JP." << endl
167 cout << "Welcome to GiNaC-cint (" << PACKAGE << " V" << VERSION << ")" << endl;
168 cout << " __, _______ GiNaC: (C) 1999-2000 Johannes Gutenberg University Mainz," << endl
169 << " (__) * | Germany. Cint C/C++ interpreter: (C) 1995-2000 Masaharu" << endl
170 << " ._) i N a C | Goto and Agilent Technologies, Japan. This is free software" << endl
171 << "<-------------' with ABSOLUTELY NO WARRANTY. For details, type `.warranty'" << endl
172 << "Type .help for help." << endl
177 void helpmessage(void)
179 cout << "GiNaC-cint recognizes some special commands which start with a dot:" << endl << endl
180 << " .cint switch to cint interactive mode (see cint" << endl
181 << " documentation for further details)" << endl
182 << " .function define the body of a function (necessary due to a" << endl
183 << " cint limitation)" << endl
184 << " .help the text you are currently reading" << endl
185 << " .q, .quit, .exit, .bye quit GiNaC-cint" << endl
186 << " .read filename read a file from disk and execute it in GiNaC-cint" << endl
187 << " (recursive call is possible)" << endl
188 << " .redirect [filename] redirect 'OutXY = ...' output to a file" << endl
189 << " (.redirect alone redirects output back to console)" << endl
190 << " .restart restart GiNaC-cint (does not re-read command line" << endl
192 << " .save filename save the commands you have entered so far in a file" << endl
193 << " .silent suppress 'OutXY = ...' output (variables are still" << endl
194 << " accessible)" << endl
195 << " .warranty information on redistribution and warranty" << endl
196 << " .> [filename] same as .redirect [filename]" << endl << endl
197 << "Instead of '.cmd' you can also write '//GiNaC-cint.cmd' to be compatible with" << endl
198 << "programs that will be compiled later." << endl
199 << "Additionally you can exit GiNaC-cint with quit; exit; or bye;" << endl
204 void warrantymessage(void)
206 cout << "GiNaC is free software; you can redistribute it and/or modify it under the" << endl
207 << "the terms of the GNU General Public License as published by the Free Software" << endl
208 << "Foundation; either version 2 of the License, or (at your option) any later" << endl
209 << "version." << endl
210 << "This program is distributed in the hope that it will be useful, but WITHOUT" << endl
211 << "ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS" << endl
212 << "FOR A PARTICULAR PURPOSE. See the GNU General Public License for more" << endl
213 << "details." << endl
214 << "You should have received a copy of the GNU General Public License along with" << endl
215 << "this program. If not, write to the Free Software Foundation, 675 Mass Ave," << endl
216 << "Cambridge, MA 02139, USA." << endl << endl;
217 cout << "Cint and associated tools are copyright by Agilent Technologies Japan Company" << endl
218 << "and Masaharu Goto <MXJ02154@niftyserve.or.jp>." << endl
219 << "Source code, binary executable or library of Cint and associated tools can be" << endl
220 << "used, modified and distributed with no royalty for any purpose provided that" << endl
221 << "the copyright notice appear in all copies and that both that copyright notice" << endl
222 << "and this permission notice appear in supporting documentation." << endl
223 << "Agilent Technologies Japan and the author make no representations about the" << endl
224 << "suitability of this software for any purpose. It is provided \"AS IS\"" << endl
225 << "without express or implied warranty." << endl;
229 string preprocess(char const * const line, bool & comment, bool & single_quote,
230 bool & double_quote, unsigned & open_braces)
232 // "preprocess" the line entered to be able to decide if the command shall be
233 // executed directly or more input is needed or this is a special command
235 // ALL whitespace will be removed
236 // all comments (/* */ and //) will be removed
237 // open and close braces ( { and } ) outside strings will be counted
240 cout << "line=" << line << endl;
241 cout << "comment=" << comment << ", single_quote=" << single_quote
242 << ",double_quote=" << double_quote << ", open_braces=" << open_braces
251 bool asterisk = false;
252 while ((line[pos]!='\0')&&!end) {
254 // last character was a \, ignore this one
257 // last character was a /, test if * or /
259 if (line[pos]=='/') {
261 } else if (line[pos]=='*') {
265 preprocessed += line[pos];
267 } else if (asterisk) {
268 // last character was a *, test if /
270 if (line[pos]=='/') {
272 } else if (line[pos]=='*') {
282 // whitespace: ignore
285 // escape character, ignore next
289 if ((!single_quote)&&(!comment)) {
290 double_quote = !double_quote;
294 if ((!double_quote)&&(!comment)) {
295 single_quote = !single_quote;
299 if ((!single_quote)&&(!double_quote)&&(!comment)) {
304 if ((!single_quote)&&(!double_quote)&&(!comment)&&(open_braces>0)) {
315 preprocessed += line[pos];
322 cout << "preprocessed=" << preprocessed << endl;
323 cout << "comment=" << comment << ", single_quote=" << single_quote
324 << ",double_quote=" << double_quote << ", open_braces=" << open_braces
333 for (cplist::iterator it=filenames.begin(); it!=filenames.end(); ++it) {
334 cout << "removing file " << *it << endl;
340 void sigterm_handler(int n)
345 void initialize(void)
351 signal(SIGTERM,sigterm_handler);
355 void initialize_cint(void)
357 G__init_cint("cint"); /* initialize cint */
359 // no longer needed as of cint 5.14.31:
360 // exec_tempfile("#include <string>\n");
362 #ifndef NO_NAMESPACE_GINAC
363 exec_tempfile("using namespace GiNaC;");
364 #endif // ndef NO_NAMESPACE_GINAC
366 exec_tempfile("ex LAST,LLAST,LLLAST;\n");
367 exec_tempfile("bool ginac_cint_internal_redirect_output=false;\n");
368 exec_tempfile("ofstream ginac_cint_internal_fout;\n");
373 cout << "Restarting GiNaC-cint." << endl;
378 void redirect(string const & filename)
381 cout << "Redirecting output back to console..." << endl;
382 exec_tempfile( string()
383 +"ginac_cint_internal_redirect_output=false;\n"
384 +"ginac_cint_internal_fout.close();");
386 cout << "Redirecting output to " << filename << "..." << endl;
387 exec_tempfile( string()
388 +"ginac_cint_internal_redirect_output=true;\n"
389 +"ginac_cint_internal_fout.open(\""+filename+"\");\n");
393 bool is_command(string const & command, string & preprocessed,
394 string const & comparevalue, bool substr)
396 bool single_quote = false;
397 bool double_quote = false;
398 bool comment = false;
399 unsigned open_braces = 0;
400 if ((preprocessed=="."+comparevalue)||
401 substr&&(preprocessed.substr(0,comparevalue.length()+1)==
405 if ((command=="//GiNaC-cint."+comparevalue+"\n")||
406 substr&&(command.substr(0,comparevalue.length()+13)==
407 "//GiNaC-cint."+comparevalue)) {
408 preprocessed = preprocess(command.substr(12).c_str(),comment,
409 single_quote,double_quote,open_braces);
415 bool readlines(istream * is, string & allcommands)
418 char prompt[G__ONELINE];
423 bool next_command_is_function = false;
424 bool single_quote = false;
425 bool double_quote = false;
426 bool comment = false;
427 unsigned open_braces = 0;
429 while ((!quit)&&(!eof)) {
430 strcpy(prompt,PROMPT1);
431 bool end_of_command = false;
434 while (!end_of_command) {
436 line = G__input(prompt);
438 getline(*is,linebuffer);
439 line = linebuffer.c_str();
443 preprocessed += preprocess(line,comment,single_quote,double_quote,open_braces);
444 if ((open_braces==0)&&(!single_quote)&&(!double_quote)&&(!comment)) {
445 unsigned l = preprocessed.length();
447 (preprocessed[0]=='#')||
448 (preprocessed[0]=='.')||
449 (preprocessed[l-1]==';')||
450 (preprocessed[l-1]=='}')) {
451 end_of_command = true;
454 strcpy(prompt,PROMPT2);
456 if ((preprocessed=="quit;")||
457 (preprocessed=="exit;")||
458 (preprocessed=="bye;")||
459 (is_command(command,preprocessed,"quit"))||
460 (is_command(command,preprocessed,"exit"))||
461 (is_command(command,preprocessed,"bye"))||
462 (is_command(command,preprocessed,"q"))) {
464 } else if (is_command(command,preprocessed,"function")) {
465 cout << "next expression can be a function definition" << endl;
466 next_command_is_function=true;
467 } else if (is_command(command,preprocessed,"cint")) {
468 cout << endl << "switching to cint interactive mode" << endl;
469 cout << "'h' for help, 'q' to quit, '{ statements }' or 'p [expression]' to evaluate" << endl;
471 cout << "back from cint" << endl;
472 } else if (is_command(command,preprocessed,"help")) {
474 } else if (is_command(command,preprocessed,"read",true)) {
475 quit=readfile(preprocessed.substr(5),allcommands);
476 } else if (is_command(command,preprocessed,"save",true)) {
477 command = "/* "+command+" */"; // we do not want the .save command itself in saved files
478 savefile(preprocessed.substr(5),allcommands);
479 } else if (is_command(command,preprocessed,"restart")) {
481 } else if (is_command(command,preprocessed,"redirect",true)) {
482 redirect(preprocessed.substr(9));
483 } else if (is_command(command,preprocessed,">",true)) {
484 redirect(preprocessed.substr(2));
485 } else if (is_command(command,preprocessed,"silent")) {
486 redirect("/dev/null");
487 } else if (is_command(command,preprocessed,"warranty")) {
489 /* test for more special commands
490 } else if (preprocessed==".xyz") {
491 cout << "special command (TBD): " << command << endl;
493 } else if (command.substr(0,2)=="#!") {
494 // ignore lines which indicate that this file is executed as a script
496 // cout << "now processing: " << command << endl;
497 if (next_command_is_function) {
498 next_command_is_function = false;
499 filenames.push_back(process_permanentfile(command));
501 process_tempfile(command);
505 // test for end of file if reading from a stream
508 // save commands only when reading from keyboard
509 allcommands += command;
516 bool readfile(string const & filename, string & allcommands)
518 cout << "Reading commands from file " << filename << "." << endl;
521 fin.open(filename.c_str());
523 quit = readlines(&fin,allcommands);
525 cout << "Cannot open " << filename << " for reading." << endl;
531 void savefile(string const & filename, string const & allcommands)
533 cout << "Saving commands to file " << filename << "." << endl;
535 fout.open(filename.c_str());
539 cout << "Cannot save commands to " << filename << "." << endl;
542 cout << "Cannot open " << filename << " for writing." << endl;
547 int main(int argc, char ** argv)
553 bool argsexist = argc>1;
556 allcommands = "/* Files given as command line arguments:\n";
560 while (argc && !quit) {
561 allcommands += *argv;
563 quit=readfile(*argv,allcommands);
568 allcommands += "*/\n";
572 readlines(NULL,allcommands);