- GiNaC-list - ginac.de

Use of STL iterators in GiNaC
by Michael Goffioul 24 Apr '09

24 Apr '09

Hi, I'm trying to compile ginac-1.5.1 with VC++ under Windows. I know "I should get myself a better compiler", but that being said, I've encountered runtime problems due to usage of STL iterators that might be wrong, essentially on std::vector. Problems I've noticed are: dereferencing end(), incrementing end(), decrementing begin()... I'm currently tracking down these issues by running the test suite in the source package, but this might not catch all of them. If you are interested, I can send a patch once I succeed to run the test suite completely. Bye, Michael.

8 25

Re: [GiNaC-list] [CLN-list] msvc support patch
by Richard B. Kreckel 20 Apr '09

20 Apr '09

Hi! Bruno Haible wrote: > Michael Goffioul wrote: >> - struct/class, const/non-const: at some places, CLN foward declares >> objects as being "class", while there are actually "struct"; VC++ generates >> different mangled names for struct and class, so you end up with unresolved >> references. > > What does the C++ standard say about this? Is CLN right, or is MSVC++ right? I don't think it is worth checking what the standard says about this. Even if it permits such sloppiness, they should be avoided. I'll prepare a patch for CLN. -richy. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

1 0

Re: [GiNaC-list] msvc support patch
by Michael Goffioul 20 Apr '09

20 Apr '09

On Mon, Apr 20, 2009 at 12:33 PM, Bruno Haible <bruno(a)clisp.org> wrote: > Michael Goffioul wrote: >> This link will hopefully clarify this. >> >> http://www.gnu.org/licenses/gpl-faq.html#WindowsRuntimeAndGPL > > This FAQ does not make the distinction between whether the installer > contains the run-time libraries or not. Based on the verbiage in the GPL > (explained in the other mail) I therefore think that this FAQ assumes > that the installer does not contain the run-time libraries. (Just for the record, I didn't find back this verbiage in GPLv3). Grand. I'll strip symbolic package from future octave binary releases. Michael.

1 0

Re: [GiNaC-list] msvc support patch
by Michael Goffioul 20 Apr '09

20 Apr '09

On Mon, Apr 20, 2009 at 3:22 AM, Bruno Haible <bruno(a)clisp.org> wrote: > Michael Goffioul wrote: >> CLN/GiNaC compilation with VC++ is part of the making of Windows >> package for octave. > > There is a problem with the Windows installer of Octave: It violates the > GNU GPL, under which Octave is distributed. > > The download octave-3.0.3-setup.exe [1] contains an installable part called > "Microsoft C/C++ runtime libraries required by Octave." According to the > packaging scripts to which you pointed us [2], these runtime libraries are > the files msvcr80.dll, msvcp80.dll, msvcm80.dll. Because these files are > not normally distributed with the Windows OS and because they are not > distributed with source code, it's a breach of the GPL. (Of section 3 of > GPLv2, to be precise.) This link will hopefully clarify this. If you still don't agree, there's no point in discussing the rest. http://www.gnu.org/licenses/gpl-faq.html#WindowsRuntimeAndGPL Michael.

1 0

Release GiNaC 1.5.1
by Jens Vollinga 25 Feb '09

25 Feb '09

Hi everybody, GiNaC 1.5.1 is out and available. The changes are: * Fixed compiling problems with GNU g++ 3.4. * Work around weird C++ RTTI behaviour on woe32 so GiNaC works on MinGW again. As always, this release can be downloaded from http://www.ginac.de/Download.html Best wishes, Jens

1 0

Substitution troubles
by Jeremy Jay 25 Feb '09

25 Feb '09

Hello all, I'm new to GiNaC, but so far I've been very impressed by its power. But, I seem to be having a problem with the substitutions, can anyone tell me what I'm doing wrong? A short example is the following code: ---- #include <iostream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; ////////////////////////////////////////////////////////////// // delta( c_i, c_j ) = 1 if c_i==c_j otherwise 0 DECLARE_FUNCTION_2P(delta); REGISTER_FUNCTION(delta, dummy()); ////////////////////////////////////////////////////////////// int main(int argc, char **argv) { symbol q("q"), c0("c0"), c1("c1"), c2("c2"); ex poly( (1-delta(c0,c1))*(1-delta(c0,c2)) ); poly=poly.expand(); cout << "START : " << poly << endl; poly = poly.subs( 1 == q ); cout << "sub 1 => q : " << poly << endl; return 0; } ---- which gives the output: ---- START : 1-delta(c0,c2)-delta(c0,c1)+delta(c0,c2)*delta(c0,c1) sub 1 => q : 1-delta(c0,c2)-delta(c0,c1)+delta(c0,c2)*delta(c0,c1) ---- which is not what I would expect. I've also tried using the 'numeric' type instead of the literal 1's, but get the same effect. Am I missing something? BTW, for anyone interesed in the point of this, I'm trying to implement the algorithm in this paper: http://www.iop.org/EJ/abstract/1367-2630/11/2/023001/ Thanks in advance! Jeremy

5 10

Polynomial companion matrix
by James Jackson 18 Feb '09

18 Feb '09

Hi, As far as I can see, there is no method available to compute the companion matrix of a polynomial. May I suggest the addition of a method along these lines?: // Creates the Companion matrix for a given polynomial GiNaC::matrix CreateCompanion(const GiNaC::ex &poly, const GiNaC::ex &symbol) { if(poly.is_polynomial(symbol)) { int rank = poly.degree(symbol); GiNaC::matrix mat(rank, rank); GiNaC::ex majorCoeff = poly.coeff(symbol, rank); for(int i = 0; i < rank; ++i) { for(int j = 0; j < rank; ++j) { if(i - j == 1) { mat(i, j) = 1; } else if(i == 0) { mat(i, j) = -1 * poly.coeff(symbol, rank - j - 1) / majorCoeff; } } } return mat; } else { // Throw exception or return empty matrix } } Regards, James.

2 3

Release GiNaC 1.5.0
by Jens Vollinga 17 Feb '09

17 Feb '09

Hi everybody, GiNaC 1.5.0 is out and available. The changes are: * Added polynomial factorization. * New, faster (recursive descent) expression parser. * Faster GCD computation. * gcd() now allows the user to override (some of the) heuristics. * Writing extensions classes is simpler now: - Replaced custom RTTI by standard C++ RTTI. - Definition of archiving methods now optional when declaring a new GiNaC class. * Fixed recursion in polynomial divide that caused a significant slowdown in sqrfree(). * Improved lsolve() of systems containing non-numeric coefficients. * Removed implicit conversions from cl_N to numeric. * Improved configuration and compatibility. This is a new branch that introduces some incompatibilities on the source code level! If you have defined your own GiNaC classes, you will have to change their definitions slightly (cf. section 6.5.6 in http://www.ginac.de/tutorial/Adding-classes.html#Adding-classes) As always, this release can be downloaded from http://www.ginac.de/Download.html Best wishes, Jens

2 2

Calculating Matrix eigenvalues?
by James Jackson 17 Feb '09

17 Feb '09

Hi, I'm wondering if the Matrix class permits the calculation of eigenvectors / values? I'm guessing, reading the docs, that this isn't implemented, but thought I'd check! Regards, James.

2 1

Not factoring expression
by James Jackson 17 Feb '09

17 Feb '09

Hi, I've just started using GiNaC (and very nice it is too!), however I have some strange behavior with an expression not being simplified. One output of a linear equation solver is the expression: (C_2 *s ^ 2 *L_2 *V_1 + C_2 *L_1 *C_1 *s ^ 4 *L_2 *V_1 + L_1 *C_1 *s ^ 2 *V_1 + V_1 )*(1 +L_1*C_1*s^2+C_2*L_1*s^2+C_2*L_1*C_1*s^4*L_2+C_2*s^2*L_2)^(-1)*V_1^(-1) What is happening is that the V_1 terms in the first bracket are not being factorised, and therefore no cancellation with the final (V_1)^(-1) term is occuring. All my symbols are created with a SymbolFactory pattern, so I can be sure each V_1 is the same object. Can anyone suggest why GiNaC isn't simplifying this to cancel the V_1 terms? Regards, James.

3 3