- GiNaC-list - ginac.de

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

Why rearrange?
by gorb＠krasu.ru 04 Dec '01

04 Dec '01

symbol x("x"); symbol y("y"); symbol a("a"); symbol z("z"); ex e1=a-y+x+z; ex e2=a+x-y+z; cout <<"e1="<< e1 << endl; // e1=a-y+z+x cout <<"e2="<< e2 << endl; // e2=a-y+z+x What if I want to get the same order of the variables which was in expression on creation. How can I do this? And how can I even get ex e3=x+x; cout <<"e3="<< e3 << endl; // e3=x+x ? best regards, Dmitri Gorbenko

2 1

spam mails
by Hubert Spiesberger 26 Nov '01

26 Nov '01

hi ginac-fans, we apologize that the ginac-mailing list was abused to distribute spam mails with faked sender-address. please stay indulgent until we find a way to suppress this rubbish - we are working on it. on behalf of the ginac-authors: Hubert Spiesberger Institut fuer Physik WA ThEP Johannes-Gutenberg-Universitaet Mainz D-55099 Mainz Phone: +49-6131.39.23682 Fax: +49-6131.39.24611

2 1

OT: Annoying spam
by Roland Richter 26 Nov '01

26 Nov '01

Dear ginac moderators, it seems that the ginac list is abused by some idiots to distribute their sex spam mails (see the ginac-list archive for examples!) I guess nobody of the ginac community is responsible for this nonsense, but PLEASE do try to stop this abuse anyway, otherwise I'll have to consider to unsubscribe the list (and I don't want to do that). Regards, - Roland PS: Don't know whether this helps, here's the header of the last mail: Return-Path: <owner-ginac-list(a)zino.physik.uni-mainz.de> Received: from zino.physik.uni-mainz.de (root(a)zino.Physik.Uni-Mainz.DE [134.93.129.49]) by alijku04.edvz.uni-linz.ac.at (8.9.3/8.9.3) with ESMTP id MAA29132 for <roland.richter(a)students.uni-linz.ac.at>; Mon, 26 Nov 2001 12:46:23 +0100 Received: (from majordom@localhost) by zino.physik.uni-mainz.de (8.9.3/8.9.3/Debian 8.9.3-21) id MAA02064 for ginac-list-outgoing; Mon, 26 Nov 2001 12:45:17 +0100 X-Authentication-Warning: zino.physik.uni-mainz.de: majordom set sender to owner-ginac-list(a)ThEP.Physik.Uni-Mainz.DE using -f Received: from mailgate1.zdv.Uni-Mainz.DE (mailgate1.zdv.Uni-Mainz.DE [134.93.8.56]) by zino.physik.uni-mainz.de (8.9.3/8.9.3/Debian 8.9.3-21) with ESMTP id MAA02048; Mon, 26 Nov 2001 12:45:12 +0100 From: bkcurf(a)abacho.de Received: from femail37.sdc1.sfba.home.com (femail37.sdc1.sfba.home.com [24.254.60.31]) by mailgate1.zdv.Uni-Mainz.DE (8.11.0/8.10.2) with ESMTP id fAQBjBC11568; Mon, 26 Nov 2001 12:45:11 +0100 (MET) Received: from [65.4.82.71] by femail37.sdc1.sfba.home.com (InterMail vM.4.01.03.20 201-229-121-120-20010223) with SMTP id <20011126114419.LPFG9688.femail37.sdc1.sfba.home.com(a)[65.4.82.71]>; Mon, 26 Nov 2001 03:44:19 -0800 To: bkcurf(a)abacho.de Subject: DOWNLOADEN 20 neue Filme - IHR AKTUELLES KENNWORT = zDert456X Date: lun, 26 nov 2001 11:50:18 -1000 MIME-Version: 1.0 Content-Type: multipart/mixed; boundary="----=_NextPart_000_30C7_000000ED.00003772" X-Priority: 3 X-MSMail-Priority: Normal Message-Id: <20011126114419.LPFG9688.femail37.sdc1.sfba.home.com(a)[65.4.82.71]> Sender: owner-ginac-list(a)zino.physik.uni-mainz.de Precedence: bulk Reply-To: ginac-list(a)zino.physik.uni-mainz.de X-UIDL: 1=B!!`LR"!9G@!!Slj"!

2 1