- GiNaC-list - ginac.de

Problems
by Jonathan Heinen 31 Jan '07

31 Jan '07

Hi, i want to use ginac as a lib under Mac OS X 10.4..8. I compiled cln 1.1.13 without problems. But compiling ginac 1.3.6 gave me the well known ginsh_parser.yy error: if g++ -DHAVE_CONFIG_H -I. -I. -I.. -I./../ginac -I../ginac - DIN_GINAC -I/usr/local/include -g -O2 -MT ginsh_parser.o -MD -MP - MF ".deps/ginsh_parser.Tpo" -c -o ginsh_parser.o ginsh_parser.cc; \ then mv -f ".deps/ginsh_parser.Tpo" ".deps/ginsh_parser.Po"; else rm - f ".deps/ginsh_parser.Tpo"; exit 1; fi ginsh_parser.yy:50:6: error: missing binary operator before token "wrapper" ginsh_parser.yy:56:6: error: missing binary operator before token "wrapper" ginsh_parser.yy:860:6: error: missing binary operator before token "wrapper" ginsh_parser.yy:870:6: error: missing binary operator before token "wrapper" ginsh_parser.yy:922:6: error: missing binary operator before token "wrapper" ginsh_parser.yy: In function 'char** fcn_completion(const char*, int, int)': ginsh_parser.yy:859: error: invalid conversion from 'const char*' to 'char*' ginsh_parser.yy:869: error: invalid conversion from 'const char*' to 'char*' make[2]: *** [ginsh_parser.o] Error 1 make[1]: *** [all-recursive] Error 1 make: *** [all] Error 2 I'm using gcc --version i686-apple-darwin8-gcc-4.0.1 (GCC) 4.0.1 (Apple Computer, Inc. build 5363) Copyright (C) 2005 Free Software Foundation, Inc. and yacc --version bison (GNU Bison) 2.3 Written by Robert Corbett and Richard Stallman. Is there another thing i have to instal? O shold i use a different version of gcc/Yacc? Thanks Jonathan

6 19

Bug report: Can't construct an expression from a string
by Mitch Watrous 20 Jan '07

20 Jan '07

Hello. After the latest version of GiNaC was installed on my machine, my programs are no longer able to construct an expression from a string. I am using version 1.3.6-1.fc5 of GiNaC. My operating system is Fedora Core 5, and I am using version 4.1.1 of the gcc compiler. Here is the code that causes the error, which is based on the example code from the GiNaC tutorial section on expression input: #include <ginac/ginac.h> int main() { using namespace GiNaC; symbol x("x"), y("y"); ex e("2*x+sin(y)", lst(x, y)); return 0; } Here is the error message I get when I execute the program main: terminate called after throwing an instance of 'std::runtime_error' what(): syntax error, unexpected ')' at ) Aborted Here is my Makefile: LDFLAGS = -L/usr/lib CPPFLAGS = -g -I/usr/include CC = g++ LIBS = -lcln -lginac main_OBJ = main.o .PHONY: all clean all: $(MAKE) main clean: -rm $(main_OBJ) -rm main -rm -rf ./tmp main: $(main_OBJ) $(CC) $(LDFLAGS) $^ -o $@ $(LIBS) %.o: %.cpp $(CC) -c $(CPPFLAGS) $(CFLAGS) $< Thanks for your help, Mitch Watrous

2 1

Re: GiNaC-list Digest, Vol 22, Issue 5
by Eric Collins 06 Jan '07

06 Jan '07

Thank you for your assistance. I was able to properly serialize the ex class using char arrays, facilitated by the I/O capabilities of the archive class. Eric M. Collins On 12/24/06, ginac-list-request(a)ginac.de <ginac-list-request(a)ginac.de> wrote: > > Send GiNaC-list mailing list submissions to > ginac-list(a)ginac.de > > To subscribe or unsubscribe via the World Wide Web, visit > https://www.cebix.net/mailman/listinfo/ginac-list > or, via email, send a message with subject or body 'help' to > ginac-list-request(a)ginac.de > > You can reach the person managing the list at > ginac-list-owner(a)ginac.de > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of GiNaC-list digest..." > > > Today's Topics: > > 1. Serialization of datatypes (Eric Collins) > 2. Re: Serialization of datatypes (Chris.Dams(a)mi.infn.it) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 23 Dec 2006 10:14:51 -0600 > From: "Eric Collins" <emc2(a)ra.msstate.edu> > Subject: [GiNaC-list] Serialization of datatypes > To: ginac-list(a)ginac.de > Message-ID: > <805f78150612230814k3689154epd1d3f1bb5905992(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Would anyone be interested in helping me figure out how to serialize the > datatypes in GiNaC/CLN/GMP? I've looked at the source headers for some of > the datatypes and the use of reference counted pointers everywhere makes > this seem like a non-trivial task. > > I need to be able to pack all data relevant to a particular instantiation > of > a datatype (e.g. numeric or pseries), into a contiguous array. This data > would then be passed to another process where it would be unpacked into a > new instantiation of the same datatype. > > Any help, or suggestions of where to begin looking for relevant > information > would be much appreciated. > > Eric M. Collins > -- > Ad Astra Per Aspera >

1 0

(const) ostream&
by Harald Devos 27 Dec '06

27 Dec '06

Hi, the print functions declared in ginac/print.h all take "std::ostream &" as an argument. I wonder if this could be changed to "const std::ostream &". Else functions receiving a "const std::ostream &" (in any program using the GiNaC library) cause troubles when calling the print functions of the ginac library. Kind regards and best wishes for 2007. Harald Devos. -------------------------------------------------------------------------- Harald Devos Electronics and Information Systems Department, Faculty of Engineering, Ghent University, St. Pietersnieuwstraat 41, O O B-9000 Ghent, Belgium O O work phone: +32-9-264.34.51 fax: +32-9-264.35.94 O email: hmdevos(a)elis.UGent.be

2 1

Serialization of datatypes
by Eric Collins 23 Dec '06

23 Dec '06

Would anyone be interested in helping me figure out how to serialize the datatypes in GiNaC/CLN/GMP? I've looked at the source headers for some of the datatypes and the use of reference counted pointers everywhere makes this seem like a non-trivial task. I need to be able to pack all data relevant to a particular instantiation of a datatype (e.g. numeric or pseries), into a contiguous array. This data would then be passed to another process where it would be unpacked into a new instantiation of the same datatype. Any help, or suggestions of where to begin looking for relevant information would be much appreciated. Eric M. Collins -- Ad Astra Per Aspera

2 1

Release of GiNaC 1.3.6
by Jens Vollinga 13 Dec '06

13 Dec '06

Hello everybody, GiNaC 1.3.6 is out and available. The changes are: - Better worst case behavior in heuristic gcd. - Improved collecting for sparse multivariate polynomials. - Improved collect_common_factors. - Fixed bug in numerical evaluation of multiple polylogs. - Fixed numerical integration of complex functions. - Fixed bug in tensor::replace_contr_index(). As always, this release can be downloaded from http://www.ginac.de/Download.html Best wishes, Jens

1 0

Internal error: statement in file ./real/conv/cl_R_to_double.cc, line 55 has been reached!!
by Pedro Cruz 08 Dec '06

08 Dec '06

Hi, -------------------- this output -------------------- Pi=3.1415926535897932385 Internal error: statement in file ./real/conv/cl_R_to_double.cc, line 55 has been reached!! Please send the authors of the program a description how you produced this error! -------------------- has been produced by -------------------- #include <iostream> #include <sstream> #include <stdexcept> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { cout << "Pi=" << Pi.evalf() << endl; if (is_a<constant>(Pi)) { double f = ex_to<numeric>(Pi).to_double(); cout << f << endl; } return 0; } --------------------------------------------- on debian (etch testing) with libginac1.3c2a. --------------------------------------------- It was a first example i tried to do in order to understand ginac/cln concepts: <ex> <constant> and <numeric>. How to convert from Pi to C++ double type in one C++ line ? Thanks, Pedro Cruz

2 1

GiNaC as a polynomials manipulator?
by Francesco Biscani 03 Dec '06

03 Dec '06

Hello, I'd like to know if GiNaC can be used as a polynomials manipulator for celestial mechanics. It's needed to multiplicate and add polynomials by polynomials. An important requisite is to be able to extract the degree of the monomials in order to decide whether that particular monomial should be eliminated from the polynomial or not. Example: x^2*y+x^2*y^2+x^3*z If the truncation threshold is set to order 4 this polynomial gets truncated to x^2*y I imagine that a function to place polynomials in some kind of "canonical" form is needed, and then it is needed to inspect each monomial to see if it will be kept or not. Is this possible with GiNaC? Thanks and best regards, Francesco

2 1

Fwd: compute groebner basis?
by fabian＠fseoane.net 03 Dec '06

03 Dec '06

Hello! I'm new to GiNaC and my question is: Has GiNaC got a routine to compute Groebner basis? If the answer is no, then if i tried to implement one, would it be welcome ? Greetings ------------- Fabian

2 1

Problem with scalar_products : Scalar product of two vectors in index notation
by aravind b 09 Nov '06

09 Nov '06

Hi, Hope everyone's doing fine ! Please go through this code : #include <iostream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { symbol i_sym("i"), j_sym("j"), e("e"), a("a"), b("b") ; idx i(i_sym, 3), j( j_sym, 3) ; //defining indices i and j scalar_products sp; sp.add( indexed(e, i), indexed(e, j), delta_tensor(i,j) ); // e.i* e.j = delta_tensor.i.j cout << (indexed(e,i)*indexed(e,j)).simplify_indexed(sp) <<endl ; return 0; } does NOT give delta.j.i instead just gives e.i * e.j - why? Basically, i want to define an orthonormal basis vector set {e.i}. If i have vector_a = indexed(a, i) * indexed(e, i) and vector_b = indexed(b, j) * indexed(e, j) , then i want adotb to reduce to a.i * b.i ( ofcoure after i apply simplify_indexed(sp) to adotb ) Thanks a lot for going through this mail.

2 3