- GiNaC-list - ginac.de

symmetrize function
by Ruben Zilibowitz 22 Jul '14

22 Jul '14

The C++ sample code here http://www.ginac.de/tutorial/Symmetrization.html#Symmetrization is not working for me with the latest GiNaC (version 1.6.2). I get the following error main.cpp:58:27: error: no member named 'symmetrize' in 'GiNaC::indexed' cout << indexed(A, i, j).symmetrize(i,j) << endl; Is the tutorial outdated? Regards, Ruben

2 1

bug in series()
by Mario Prausa 18 Jun '14

18 Jun '14

Hi, there seems to be a bug in the series() implementation. The attached code causes a segfault. As log(-x) doesn't depend on the expansion parameter ep, series() should return log(-x). The bug only occurs, if x is defined as possymbol, so it might have something to do with the branch cut of the logarithm. I reproduced the bug with GiNaC 1.6.2 and the current git version. Regards, Mario

1 0

Re: [GiNaC-list] [PATCH] add configure option to enable deterministic sorting
by Valery Yundin 29 May '14

29 May '14

Hi, On 29 May 2014 07:31, Alexei Sheplyakov <alexei.sheplyakov(a)gmail.com> wrote: > I'm terribly sorry, but everything which makes the hash function worse > (slower or produces more collisions) is not acceptable. If you insist on that regardless of how it affects overall performance, it is not a technical problem and cannot be fixed with a patch. In any case, if my previous patches are not acceptable for whatever reasons GiNaC developers might have, I suggest to close this discussion by adding a configure option. So users can decide for themselves if (non-)determinism is worth few per-cent speed difference. Below is the simplest patch one can imagine. By default it keeps things as they are now. If user specifies --enable-determinism then make_hash_seed simply returns 0. Objects with the same calchash() implementation *and* the same contents *but* different types will have a hash collision, which be resolved in base::compare by using typeid_1.before(typeid_2) (on Linux it is a simple strcmp). Since it happens quite rarely (depends on a particular code of course), it is better than always using crc32 as in windows version of your code. In fact it can sometimes be faster than current "fast" hash. Best, Valery

1 0

[PATCH] fix bug/feature causing non-deterministic output
by Valery Yundin 29 May '14

29 May '14

Hi, Please see attached patch which fixes rather peculiar "feature" of GiNaC - non-determinism. One can think of many arguments why non-determinism of CAS is a bad idea. Shortly speaking, it is simply unacceptable. Especially when there is no good reason for it (and I can hardly think of a good reason anyway). At the moment GiNaC relies on OS-assigned memory addresses of type-names to define its term-ordering, which results in different output from the same program. When this bug was brought up in the mailing list (several times) the developers called it a "design decision". Well, it is time to stop the insanity. The attached patch fixes the bug at zero cost to the performance. Here are the results of "/usr/bin/time make check" (second run with pre-warmed caches) unpatched non-deterministic: 132.57 s fixed deterministic: 132.62 s Keep in mind that it wouldn't make GiNaC canonical ordering entirely predictable, since it still depends on the order of symbol declaration and such (btw, fixing that wouldn't be too hard either). But at least for the same program and compiler the result must be identical. PS I added FNV1a hash, one could just keep using crc32 which is already there. It makes no difference since the hash function is now called only once for each type (before it was called for every object). PPS I'm attaching patches against master and release_1-6-2 since I hope that developers recognize it as a bug and update stable releases as well. With best regards, Valery Yundin.

3 5

Re: [GiNaC-list] [PATCH] fix bug/feature causing non-deterministic output
by Valery Yundin 25 May '14

25 May '14

Hi, Updated patch, rev 3. Overloaded gethashseed() in a few classes which I missed before. Also modified gethashseed() to use "fast" hash for derived classes which do not overload it and made it independent of mangling scheme for non-template classes. Removed some debugging output. IMO "fast" hash can be removed entirely and if classes do not overload gethashseed() it will be taken care of by typeid. But since it is relevant only for custom user-defined classes, I don't care much. I kept it because GiNaC developers seem to have sentimental attachment to it. Best, Valery PS patch chunks have only 1 context line, because this list has post limit of 40kB With best regards, Mr. Valery Yundin. On 25 May 2014 18:44, Valery Yundin <yuvalery(a)gmail.com> wrote: > Hi, > > Updated patch, rev 3. > Overloaded gethashseed() in a few classes which I missed before. Also > modified gethashseed() to use "fast" hash for derived classes which do > not overload it and made it independent of mangling scheme for > non-template classes. Removed some debugging output. > > IMO "fast" hash can be removed entirely and if classes do not overload > gethashseed() it will be taken care of by typeid. But since it is > relevant only for custom user-defined classes, I don't care much. I > kept it since GiNaC developers seem to have sentimental attachment to > it. > > Best, > Valery

1 0

Possible bug in factor()
by tluthe＠physik.uni-bielefeld.de 18 May '14

18 May '14

Hi, there appears to be a bug in the multivariate factor() routine. For some expressions it seems to get stuck in an infinite loop, without any result or error messages appearing. The simplest expression I could find to reproduce this behaviour consistently is the expanded form of (3*x+2*y)*(x+y+x*y). For some other expressions (e.g. (x+2*y)*(x+y+x*y)) the problem occurs randomly, but not all the time. I'm guessing this is due to the randomness of the internal ordering of symbols. A minimal example code exhibiting the problem is attached. I've reproduced this bevahiour on openSuse and Ubuntu machines using GiNaC 1.6.2 and CLN 1.3.2. Regards, Thomas

2 1

Example Input Program Hanging Up
by Casey Conger 03 Feb '14

03 Feb '14

Hello all, I’m having trouble with running the example interactive program from the tutorial. I simply cut and pasted it, but after it asks for “an expression involving ‘x’ “ it just hangs and never does anything. The problem seems to be in the line ‘ex e = reader(cin)’, I got the other examples to run just fine, but the reader(cin) seems to be problematic. My setup: Macbook Air1.5GHZ i7, OSX 10.9.1, ginac 1.6.2 My program: #include <iostream> #include <string> #include <stdexcept> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { cout << "Enter an expression containing ’x’: " << flush; parser reader; try { ex e = reader(cin); symtab table = reader.get_syms(); symbol x = table.find("x") != table.end() ? ex_to<symbol>(table["x"]) : symbol("x"); cout << "The derivative of " << e << " with respect to x is "; cout << e.diff(x) << "." << endl; } catch (exception &p) { cerr << p.what() << endl; } } Commandline: Caseys-MacBook-Air:ginac casey$ g++ ginacPlay.cpp -o ginacPlay.out -lginac -lcln Caseys-MacBook-Air:ginac casey$ ./ginacPlay.out Enter an expression containing ’x’: x ^C …It just hangs there for several minutes.

2 1

Stay in touch with me through LinkedIn
by Jose Antonio Garcia Peiro 29 Jan '14

29 Jan '14

LinkedIn ------------ Me gustaría añadirte a mi red profesional en LinkedIn. -Jose Antonio Jose Antonio Garcia Peiro Estudiante de Universidad Miguel Hernández de Elche Murcia y alrededores, España Confirm that you know Jose Antonio Garcia Peiro: https://www.linkedin.com/e/-qmwzwc-hr0w8crr-5d/isd/4281101300/r-1XX1c3/?hs=… -- You are receiving Invitation to Connect emails. Click to unsubscribe: http://www.linkedin.com/e/-qmwzwc-hr0w8crr-5d/kGy7Yg-6HHAoBPMg-Y4n5eWUQMwWs… (c) 2012 LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Problems running ginac on OSX 10.9
by Casey Conger 19 Jan '14

19 Jan '14

Hello, I apologize in advance because I am extremely bad with anything beyond simple computing. My Setup: OSX 10.9.1, gcc4.8 My Situation: After some trouble compiling, I finally was able make/ make install with no errors. Now when I try to run one of the example programs, it doesn’t compile. I see all of the header files located in /usr/local/include/ginac. I try to compile with g++ derivative.cpp -I /usr/local/include or g++ derivative.cpp -L /usr/local/include but I get the following: Caseys-MacBook-Air:examples casey$ g++ derivative.cpp -I/usr/local/include/ Undefined symbols for architecture x86_64: "GiNaC::library_init::library_init()", referenced from: ___cxx_global_var_init in derivative-sy6zT0.o "GiNaC::library_init::~library_init()", referenced from: ___cxx_global_var_init in derivative-sy6zT0.o "GiNaC::add_unarchiver::add_unarchiver()", referenced from: ___cxx_global_var_init18 in derivative-sy6zT0.o "GiNaC::add_unarchiver::~add_unarchiver()", referenced from: ___cxx_global_var_init18 in derivative-sy6zT0.o "GiNaC::idx_unarchiver::idx_unarchiver()", referenced from: ___cxx_global_var_init22 in derivative-sy6zT0.o "GiNaC::idx_unarchiver::~idx_unarchiver()", referenced from: ___cxx_global_var_init22 in derivative-sy6zT0.o "GiNaC::lst_unarchiver::lst_unarchiver()", referenced from: ___cxx_global_var_init2 in derivative-sy6zT0.o "GiNaC::lst_unarchiver::~lst_unarchiver()", referenced from: ___cxx_global_var_init2 in derivative-sy6zT0.o "GiNaC::mul_unarchiver::mul_unarchiver()", referenced from: ___cxx_global_var_init19 in derivative-sy6zT0.o "GiNaC::mul_unarchiver::~mul_unarchiver()", referenced from: ___cxx_global_var_init19 in derivative-sy6zT0.o "GiNaC::fail_unarchiver::fail_unarchiver()", referenced from: ___cxx_global_var_init4 in derivative-sy6zT0.o "GiNaC::fail_unarchiver::~fail_unarchiver()", referenced from: ___cxx_global_var_init4 in derivative-sy6zT0.o "GiNaC::su3d_unarchiver::su3d_unarchiver()", referenced from: ___cxx_global_var_init34 in derivative-sy6zT0.o "GiNaC::su3d_unarchiver::~su3d_unarchiver()", referenced from: ___cxx_global_var_init34 in derivative-sy6zT0.o "GiNaC::su3f_unarchiver::su3f_unarchiver()", referenced from: ___cxx_global_var_init33 in derivative-sy6zT0.o "GiNaC::su3f_unarchiver::~su3f_unarchiver()", referenced from: ___cxx_global_var_init33 in derivative-sy6zT0.o "GiNaC::su3t_unarchiver::su3t_unarchiver()", referenced from: ___cxx_global_var_init32 in derivative-sy6zT0.o "GiNaC::su3t_unarchiver::~su3t_unarchiver()", referenced from: ___cxx_global_var_init32 in derivative-sy6zT0.o "GiNaC::color_unarchiver::color_unarchiver()", referenced from: ___cxx_global_var_init30 in derivative-sy6zT0.o "GiNaC::color_unarchiver::~color_unarchiver()", referenced from: ___cxx_global_var_init30 in derivative-sy6zT0.o "GiNaC::ncmul_unarchiver::ncmul_unarchiver()", referenced from: ___cxx_global_var_init10 in derivative-sy6zT0.o "GiNaC::ncmul_unarchiver::~ncmul_unarchiver()", referenced from: ___cxx_global_var_init10 in derivative-sy6zT0.o "GiNaC::power_unarchiver::power_unarchiver()", referenced from: ___cxx_global_var_init8 in derivative-sy6zT0.o "GiNaC::power_unarchiver::~power_unarchiver()", referenced from: ___cxx_global_var_init8 in derivative-sy6zT0.o "GiNaC::matrix_unarchiver::matrix_unarchiver()", referenced from: ___cxx_global_var_init6 in derivative-sy6zT0.o "GiNaC::matrix_unarchiver::~matrix_unarchiver()", referenced from: ___cxx_global_var_init6 in derivative-sy6zT0.o "GiNaC::su3one_unarchiver::su3one_unarchiver()", referenced from: ___cxx_global_var_init31 in derivative-sy6zT0.o "GiNaC::su3one_unarchiver::~su3one_unarchiver()", referenced from: ___cxx_global_var_init31 in derivative-sy6zT0.o "GiNaC::symbol_unarchiver::symbol_unarchiver()", referenced from: ___cxx_global_var_init11 in derivative-sy6zT0.o "GiNaC::symbol_unarchiver::~symbol_unarchiver()", referenced from: ___cxx_global_var_init11 in derivative-sy6zT0.o "GiNaC::unarchive_table_t::unarchive_table_t()", referenced from: ___cxx_global_var_init1 in derivative-sy6zT0.o "GiNaC::unarchive_table_t::~unarchive_table_t()", referenced from: ___cxx_global_var_init1 in derivative-sy6zT0.o "GiNaC::varidx_unarchiver::varidx_unarchiver()", referenced from: ___cxx_global_var_init23 in derivative-sy6zT0.o "GiNaC::varidx_unarchiver::~varidx_unarchiver()", referenced from: ___cxx_global_var_init23 in derivative-sy6zT0.o "GiNaC::get_default_reader()", referenced from: _main in derivative-sy6zT0.o "GiNaC::indexed_unarchiver::indexed_unarchiver()", referenced from: ___cxx_global_var_init15 in derivative-sy6zT0.o "GiNaC::indexed_unarchiver::~indexed_unarchiver()", referenced from: ___cxx_global_var_init15 in derivative-sy6zT0.o "GiNaC::numeric_unarchiver::numeric_unarchiver()", referenced from: ___cxx_global_var_init7 in derivative-sy6zT0.o "GiNaC::numeric_unarchiver::~numeric_unarchiver()", referenced from: ___cxx_global_var_init7 in derivative-sy6zT0.o "GiNaC::pseries_unarchiver::pseries_unarchiver()", referenced from: ___cxx_global_var_init16 in derivative-sy6zT0.o "GiNaC::pseries_unarchiver::~pseries_unarchiver()", referenced from: ___cxx_global_var_init16 in derivative-sy6zT0.o "GiNaC::spinidx_unarchiver::spinidx_unarchiver()", referenced from: ___cxx_global_var_init24 in derivative-sy6zT0.o "GiNaC::spinidx_unarchiver::~spinidx_unarchiver()", referenced from: ___cxx_global_var_init24 in derivative-sy6zT0.o "GiNaC::clifford_unarchiver::clifford_unarchiver()", referenced from: ___cxx_global_var_init35 in derivative-sy6zT0.o "GiNaC::clifford_unarchiver::~clifford_unarchiver()", referenced from: ___cxx_global_var_init35 in derivative-sy6zT0.o "GiNaC::constant_unarchiver::constant_unarchiver()", referenced from: ___cxx_global_var_init3 in derivative-sy6zT0.o "GiNaC::constant_unarchiver::~constant_unarchiver()", referenced from: ___cxx_global_var_init3 in derivative-sy6zT0.o "GiNaC::diracone_unarchiver::diracone_unarchiver()", referenced from: ___cxx_global_var_init36 in derivative-sy6zT0.o "GiNaC::diracone_unarchiver::~diracone_unarchiver()", referenced from: ___cxx_global_var_init36 in derivative-sy6zT0.o "GiNaC::function_unarchiver::function_unarchiver()", referenced from: ___cxx_global_var_init20 in derivative-sy6zT0.o "GiNaC::function_unarchiver::~function_unarchiver()", referenced from: ___cxx_global_var_init20 in derivative-sy6zT0.o "GiNaC::integral_unarchiver::integral_unarchiver()", referenced from: ___cxx_global_var_init5 in derivative-sy6zT0.o "GiNaC::integral_unarchiver::~integral_unarchiver()", referenced from: ___cxx_global_var_init5 in derivative-sy6zT0.o "GiNaC::symmetry_unarchiver::symmetry_unarchiver()", referenced from: ___cxx_global_var_init17 in derivative-sy6zT0.o "GiNaC::symmetry_unarchiver::~symmetry_unarchiver()", referenced from: ___cxx_global_var_init17 in derivative-sy6zT0.o "GiNaC::wildcard_unarchiver::wildcard_unarchiver()", referenced from: ___cxx_global_var_init14 in derivative-sy6zT0.o "GiNaC::wildcard_unarchiver::~wildcard_unarchiver()", referenced from: ___cxx_global_var_init14 in derivative-sy6zT0.o "GiNaC::possymbol_unarchiver::possymbol_unarchiver()", referenced from: ___cxx_global_var_init13 in derivative-sy6zT0.o "GiNaC::possymbol_unarchiver::~possymbol_unarchiver()", referenced from: ___cxx_global_var_init13 in derivative-sy6zT0.o "GiNaC::tensdelta_unarchiver::tensdelta_unarchiver()", referenced from: ___cxx_global_var_init25 in derivative-sy6zT0.o "GiNaC::tensdelta_unarchiver::~tensdelta_unarchiver()", referenced from: ___cxx_global_var_init25 in derivative-sy6zT0.o "GiNaC::diracgamma_unarchiver::diracgamma_unarchiver()", referenced from: ___cxx_global_var_init37 in derivative-sy6zT0.o "GiNaC::diracgamma_unarchiver::~diracgamma_unarchiver()", referenced from: ___cxx_global_var_init37 in derivative-sy6zT0.o "GiNaC::minkmetric_unarchiver::minkmetric_unarchiver()", referenced from: ___cxx_global_var_init27 in derivative-sy6zT0.o "GiNaC::minkmetric_unarchiver::~minkmetric_unarchiver()", referenced from: ___cxx_global_var_init27 in derivative-sy6zT0.o "GiNaC::realsymbol_unarchiver::realsymbol_unarchiver()", referenced from: ___cxx_global_var_init12 in derivative-sy6zT0.o "GiNaC::realsymbol_unarchiver::~realsymbol_unarchiver()", referenced from: ___cxx_global_var_init12 in derivative-sy6zT0.o "GiNaC::relational_unarchiver::relational_unarchiver()", referenced from: ___cxx_global_var_init9 in derivative-sy6zT0.o "GiNaC::relational_unarchiver::~relational_unarchiver()", referenced from: ___cxx_global_var_init9 in derivative-sy6zT0.o "GiNaC::spinmetric_unarchiver::spinmetric_unarchiver()", referenced from: ___cxx_global_var_init28 in derivative-sy6zT0.o "GiNaC::spinmetric_unarchiver::~spinmetric_unarchiver()", referenced from: ___cxx_global_var_init28 in derivative-sy6zT0.o "GiNaC::tensmetric_unarchiver::tensmetric_unarchiver()", referenced from: ___cxx_global_var_init26 in derivative-sy6zT0.o "GiNaC::tensmetric_unarchiver::~tensmetric_unarchiver()", referenced from: ___cxx_global_var_init26 in derivative-sy6zT0.o "GiNaC::diracgamma5_unarchiver::diracgamma5_unarchiver()", referenced from: ___cxx_global_var_init38 in derivative-sy6zT0.o "GiNaC::diracgamma5_unarchiver::~diracgamma5_unarchiver()", referenced from: ___cxx_global_var_init38 in derivative-sy6zT0.o "GiNaC::diracgammaL_unarchiver::diracgammaL_unarchiver()", referenced from: ___cxx_global_var_init39 in derivative-sy6zT0.o "GiNaC::diracgammaL_unarchiver::~diracgammaL_unarchiver()", referenced from: ___cxx_global_var_init39 in derivative-sy6zT0.o "GiNaC::diracgammaR_unarchiver::diracgammaR_unarchiver()", referenced from: ___cxx_global_var_init40 in derivative-sy6zT0.o "GiNaC::diracgammaR_unarchiver::~diracgammaR_unarchiver()", referenced from: ___cxx_global_var_init40 in derivative-sy6zT0.o "GiNaC::fderivative_unarchiver::fderivative_unarchiver()", referenced from: ___cxx_global_var_init21 in derivative-sy6zT0.o "GiNaC::fderivative_unarchiver::~fderivative_unarchiver()", referenced from: ___cxx_global_var_init21 in derivative-sy6zT0.o "GiNaC::tensepsilon_unarchiver::tensepsilon_unarchiver()", referenced from: ___cxx_global_var_init29 in derivative-sy6zT0.o "GiNaC::tensepsilon_unarchiver::~tensepsilon_unarchiver()", referenced from: ___cxx_global_var_init29 in derivative-sy6zT0.o "GiNaC::basic::basic(GiNaC::basic const&)", referenced from: GiNaC::symbol::symbol(GiNaC::symbol const&) in derivative-sy6zT0.o "GiNaC::parser::parser(std::__1::map<std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> >, GiNaC::ex, std::__1::less<std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> > >, std::__1::allocator<std::__1::pair<std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> > const, GiNaC::ex> > > const&, bool, std::__1::map<std::__1::pair<std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> >, unsigned long>, GiNaC::ex (*)(std::__1::vector<GiNaC::ex, std::__1::allocator<GiNaC::ex> > const&), GiNaC::PrototypeLess, std::__1::allocator<std::__1::pair<std::__1::pair<std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> >, unsigned long> const, GiNaC::ex (*)(std::__1::vector<GiNaC::ex, std::__1::allocator<GiNaC::ex> > const&)> > > const&)", referenced from: _main in derivative-sy6zT0.o "GiNaC::parser::~parser()", referenced from: _main in derivative-sy6zT0.o "GiNaC::parser::operator()(std::__1::basic_istream<char, std::__1::char_traits<char> >&)", referenced from: _main in derivative-sy6zT0.o "GiNaC::symbol::symbol(std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> > const&)", referenced from: _main in derivative-sy6zT0.o "GiNaC::_num0_bp", referenced from: GiNaC::ex::ex() in derivative-sy6zT0.o "GiNaC::operator<<(std::__1::basic_ostream<char, std::__1::char_traits<char> >&, GiNaC::ex const&)", referenced from: _main in derivative-sy6zT0.o "GiNaC::ex::diff(GiNaC::symbol const&, unsigned int) const", referenced from: _main in derivative-sy6zT0.o "vtable for GiNaC::basic", referenced from: GiNaC::basic::~basic() in derivative-sy6zT0.o NOTE: a missing vtable usually means the first non-inline virtual member function has no definition. "vtable for GiNaC::symbol", referenced from: GiNaC::symbol::~symbol() in derivative-sy6zT0.o GiNaC::symbol::symbol(GiNaC::symbol const&) in derivative-sy6zT0.o NOTE: a missing vtable usually means the first non-inline virtual member function has no definition. ld: symbol(s) not found for architecture x86_64 clang: error: linker command failed with exit code 1 (use -v to see invocation) It seems like it just can’t find these functions, but I have no idea why. Thank you in advance for your help. -Casey

1 0

Randomness when using lsolve
by Federico Bergero 16 Jan '14

16 Jan '14

Hello everybody I'm new to GiNaC. I'm trying to use it to solve linear algebraic system of equations and running the example from the tutorial https://www.ginac.de/tutorial/Solving-linear-systems-of-equations.html#Solv… symbol a("a"), b("b"), x("x"), y("y"); lst eqns, vars; eqns = a*x+b*y==3, x-y==b; vars = x, y; cout << lsolve(eqns, vars) << endl; While it works perfectly when I run it several times I get different solutions (they are all equivalent though). For example: {x==(3+b^2)*(a+b)^(-1),y==-(a+b)^(-1)*(-3+a*b)} or {x==(b+a)^(-1)*(3+b^2),y==-(-3+b*a)*(b+a)^(-1)} what is going on there? Is there a way to have a deterministic solution? Thanks in advance. Fede -- I find television very educational. Every time someone switches it on, I go into another room and read a good book.

2 2