- GiNaC-list - ginac.de

Request for a sample source file to create 8x8 matrix of real numbers.
by Ananda Murthy R S 19 Jan '02

19 Jan '02

Hello: I hereby request anybody to give me sample code to create an 8x8 matrix of real numbers using GiNaC and to print it. It would help me a lot. Thanks in advance for taking the trouble of helping me. Sincerely, Anand

2 1

Invalid conversions and error messages.
by Bob McElrath 19 Jan '02

19 Jan '02

I wrote a program that GiNaC doesn't like, and it told me so. (I'm becoming somewhat of an expert at this) ;) #include <ginac/ginac.h> #include <iostream.h> using namespace GiNaC; using namespace std; int main() { symbol a("a"); ex blah = a; numeric blahn = ex_to<numeric>(blah); // The line below generates the error: // Internal error: statement in file ./float/misc/cl_F_digits.cc, line 25 has been reached!! // Please send the authors of the program a description how you produced this error! cout << "blahn = " << blahn << endl; // If the cout is commented out, the line below generates the error: // Internal error: statement in file ./real/elem/cl_R_minusp.cc, line 30 has been reached!! // Please send the authors of the program a description how you produced this error! if(blahn.is_negative()) { cout << "blahn is negative." << endl; } } Looks to me like a cln error. But GiNaC should probably throw an exception for the invalid conversion of "a" to a numeric. Cheers, -- Bob Bob McElrath (rsmcelrath(a)students.wisc.edu) Univ. of Wisconsin at Madison, Department of Physics

2 1

Feature request: clearing lsts
by Roberto Bagnara 19 Jan '02

19 Jan '02

I did not find in the documentation a way to erase all the elements of an lst. More formally, what is needed is the functionality provided by, e.g., void clear(GiNaC::lst& l) { for (unsigned n = l.nops(); n-- > 0; ) l.remove_first(); assert(l.nops() == 0); } Something like that turns out to be very useful in several situations. For an example, consider the following: if (match(e, pattern_1, substitution)) ... else if (clear(substitution), match(e, pattern_2, substitution)) ... else if (clear(substitution), match(e, pattern_3, substitution)) ... Please, consider the addition of such a method to GiNaC::lst. You may not want to call it clear() though, as this is perhaps too STL-ish (someone may be induced to think that lst is an STL container and we do not want that ;-) All the best Roberto -- Prof. Roberto Bagnara Computer Science Group Department of Mathematics, University of Parma, Italy http://www.cs.unipr.it/~bagnara/ mailto:bagnara@cs.unipr.it

2 1

Unable to create 3x8 matrix having real numbers.
by Ananda Murthy R S 13 Jan '02

13 Jan '02

Hello: The program given below is not compiling. //Creating a 3x8 Matrix #include <iostream> #include <ginac/ginac.h> using namespace GiNaC; using namespace std; int main() { matrix A(3,8, lst(-2.902280128, 4.902280128, 3.698135576 , -1.091710896, -1.339597088, -3.833588112, 0.1354525359, 0.1354525359, 14.57598539, 8.663631092, 4.471485536, -6.163631092, -1.031929416, 1.692145555, -11.85191042, 0.5280631913, -13.01717766, 5.242845433, 6.808655713, 0.3658471056, -1.250335699, -9.401666250, 5.593010537, -0.1838540200)); cout << "A = " << A << endl; return 0; } But if I change it to a 2 x 8 matrix as shown below it is compiling. //Creating a 2x8 Matrix #include <iostream> #include <ginac/ginac.h> using namespace GiNaC; using namespace std; int main() { matrix A(2,8, lst(-2.902280128, 4.902280128, 3.698135576 , -1.091710896, -1.339597088, -3.833588112, 0.1354525359, 0.1354525359, 14.57598539, 8.663631092, 4.471485536, -6.163631092, -1.031929416, 1.692145555, -11.85191042, 0.5280631913)); cout << "A = " << A << endl; return 0; } Can somebody tell me what is the problem? Is there any limit on the dimensions of the matrix that can be created in GiNaC? I am runnilng GiNaC on Linux-Mandrake 8.1. Thanks in advance, Sincerely, Ananda Murthy R S

1 0

Unable to create 3x8 matrix having real numbers.
by Ananda Murthy R S 13 Jan '02

13 Jan '02

Hello: The program given below is not compiling. //Creating a 3x8 Matrix #include <iostream> #include <ginac/ginac.h> using namespace GiNaC; using namespace std; int main() { matrix A(3,8, lst(-2.902280128, 4.902280128, 3.698135576 , -1.091710896, -1.339597088, -3.833588112, 0.1354525359, 0.1354525359, 14.57598539, 8.663631092, 4.471485536, -6.163631092, -1.031929416, 1.692145555, -11.85191042, 0.5280631913, -13.01717766, 5.242845433, 6.808655713, 0.3658471056, -1.250335699, -9.401666250, 5.593010537, -0.1838540200)); cout << "A = " << A << endl; return 0; } But if I change it to a 2 x 8 matrix as shown below it is compiling. //Creating a 2x8 Matrix #include <iostream> #include <ginac/ginac.h> using namespace GiNaC; using namespace std; int main() { matrix A(2,8, lst(-2.902280128, 4.902280128, 3.698135576 , -1.091710896, -1.339597088, -3.833588112, 0.1354525359, 0.1354525359, 14.57598539, 8.663631092, 4.471485536, -6.163631092, -1.031929416, 1.692145555, -11.85191042, 0.5280631913)); cout << "A = " << A << endl; return 0; } Can somebody tell me what is the problem? Is there any limit on the dimensions of the matrix that can be created in GiNaC? I am runnilng GiNaC on Linux-Mandrake 8.1. Thanks in advance, Sincerely, Ananda Murthy R S

1 0

Unable to create 3x8 matrix having real numbers.
by Ananda Murthy R S 13 Jan '02

13 Jan '02

Hello: The program given below is not compiling. //Creating a 3x8 Matrix #include <iostream> #include <ginac/ginac.h> using namespace GiNaC; using namespace std; int main() { matrix A(3,8, lst(-2.902280128, 4.902280128, 3.698135576 , -1.091710896, -1.339597088, -3.833588112, 0.1354525359, 0.1354525359, 14.57598539, 8.663631092, 4.471485536, -6.163631092, -1.031929416, 1.692145555, -11.85191042, 0.5280631913, -13.01717766, 5.242845433, 6.808655713, 0.3658471056, -1.250335699, -9.401666250, 5.593010537, -0.1838540200)); cout << "A = " << A << endl; return 0; } But if I change it to a 2 x 8 matrix as shown below it is compiling. //Creating a 2x8 Matrix #include <iostream> #include <ginac/ginac.h> using namespace GiNaC; using namespace std; int main() { matrix A(2,8, lst(-2.902280128, 4.902280128, 3.698135576 , -1.091710896, -1.339597088, -3.833588112, 0.1354525359, 0.1354525359, 14.57598539, 8.663631092, 4.471485536, -6.163631092, -1.031929416, 1.692145555, -11.85191042, 0.5280631913)); cout << "A = " << A << endl; return 0; } Can somebody tell me what is the problem? Is there any limit on the dimensions of the matrix that can be created in GiNaC? I am runnilng GiNaC on Linux-Mandrake 8.1. Thanks in advance, Sincerely, Ananda Murthy R S

2 1

evalf()/Digits bugs
by Bob McElrath 08 Jan '02

08 Jan '02

Some interesting (and wrong) calculations, probably rounding somewhere: (1)<mcelrath@hawk:/u/mcelrath> ginsh ginsh - GiNaC Interactive Shell (GiNaC V1.0.3) __, _______ Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, (__) * | Germany. This is free software with ABSOLUTELY NO WARRANTY. ._) i N a C | You are welcome to redistribute it under certain conditions. <-------------' For details type `warranty;'. Type ?? for a list of help topics. > evalf(1/3); 0.33333333333333333334 > Digits=30; 30 > evalf(1/3); 0.333333333333333333333333333333333333334 > Digits=40; 40 > evalf(1/3); 0.3333333333333333333333333333333333333333333333334 > Digits=50; 50 > evalf(1/3); 0.33333333333333333333333333333333333333333333333333333333336 > Note also that this is on page 14 of the tutorial, and the tutorial does not have the unusual final digit... Cheers, -- Bob Bob McElrath (rsmcelrath(a)students.wisc.edu) Univ. of Wisconsin at Madison, Department of Physics

3 6

ROOT 3.02/06 and GiNaC 1.0.3
by Hans Peter W�rmli 08 Jan '02

08 Jan '02

I am an unexperienced user of root, cint or ginac, which latest available versions I installed on my Linux 2.4.2 (ROOT 3.02/06 and GiNaC 1.0.3). I tried to interprete the simple GiNaC hello.cc example (from the Tutorial) with cint and root/cint, but the interpreter complains that library_init() is declared, but not defined. It is defined in utils.cpp. I tried to "gSystem->Load("libginac.so");" in root, but if I' lucky, the new cint error message in root is "~library_init() declared, but not defined" (~library_init() does in fact no seem to be defined). Could anybody give me a hint, what I should do? Thanks Hans Peter

2 3

Just finished compiling version 1.0.2?
by Richard B. Kreckel 21 Dec '01

21 Dec '01

Well, here comes the next load to keep your machine busy. Version 1.0.3 contains two bugfixes without compromising binary compatiblilty: * Fixed a bug where quo() would call vector::reserve() with a negative argument (reported by Pearu Peterson). * Fixed several bugs in code generation (reported by Do Hoang Son). We apologize for the inconvenience... -richy. -- Richard B. Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

1 0

After lifting the "no-bugs" ordinance...
by Richard B. Kreckel 19 Dec '01

19 Dec '01

...a couple of bugs were indeed discovered and fixed wich leeds us to the release of GiNaC-1.0.2, available at the usual places. It does not break compatibility with prior 1.0.n-versions. Enjoy -richy. -- Richard B. Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

1 0