- GiNaC-list - ginac.de

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

Planning to use GiNaC
by Roberto Bagnara 07 Dec '01

07 Dec '01

Hi there, here at the University of Parma we are about to start a project on automatic program complexity analysis. Part of our work will consist in developing a solver for several classes of recurrence relations. From a superficial look at the GiNaC web pages we gathered the impression that your library would significantly simplify our job, since representation, manipulation, simplification and matching of symbolic expression would be readily available. Do you know of anyone using GiNaC to do something similar to what we have in mind? Any advice you would like to give to someone embarking in such an adventure? Thanks a lot, Roberto -- 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

SolverQ 2.0
by Eddy McHam 05 Dec '01

05 Dec '01

Hi, Just small link to the freetool that have some some part of symbolic comptuation: SolverQ 2.0. SolverQ is a program that solves equations. Equations are entered in natural form as text. It performs limited symbolic computation and presents graphic display of results. Together with installation comes about 40 lessons. As you move through these lessons you will learn some of the power of SolverQ. SolverQ represents a philosophy in which we feel the clarity of a problem formulation is paramount, and that solution details should be largely transparent to a user. SolverQ supports: solving nonlinear equations, repetitive solutions, graph plotting and exporting to EMF, complex and transcendental complex functions, symbolic simplification, symbolic differentiation, symbolic elimination, symbolic substitution, symbolic linearization (numeric), symbolic linearization (analytic), IF Statements, unconstrained and constrained optimization, regression, FFT's, differential equations, definite integral equations. Application is available as FREEWARE for Win 95/98/ME/NT/2000/XP and can be downloaded from: http://www.geocities.com/izudindz/SolverQ.zip Regards. _________________________________________________________________ Get your FREE download of MSN Explorer at http://explorer.msn.com/intl.asp

1 0