- GiNaC-list - ginac.de

type casting
by Marko Riedel 12 Feb '23

12 Feb '23

Hi all, I have a set of custom objects that is declared like this: struct cindsym_is_less { bool operator() (const cindsymbol &lh, const cindsymbol &rh) const { return lh.compare(rh) < 0; } }; typedef std::set<cindsymbol, cindsym_is_less> cindsymset; Now I want to collect all custom objects that appear in an expression tree into this kind of set. I tried this: void indets_recurse(ex src, cindsymset & coll) { if(is_a<add>(src) || is_a<mul>(src) || is_a<power>(src)){ for(size_t idx = 0; idx < src.nops(); idx++) indets_recurse(src.op(idx), coll); return; } if(is_a<cindsymbol>(src)){ coll.insert(src); } } However this will produce the error "no matching function call" because src is an expression and not a cindsymbol. I know that I have my kind of symbol because is_a<cindsymbol> returned true. How can I get at the symbol wrapped in the expression? Best regards, Marko

2 4

symbol equality
by Marko Riedel 09 Feb '23

09 Feb '23

Greetings to all. I have been subclassing "symbol" for my project, where I need to attach an integer value to a symbol. I am not sure this is the best way to go about it but I would like to try. Everything works so far but now I have moved on to substituting variables in polynomials. Here is my sample program. #include <iostream> #include <ginac/ginac.h> #include "cindsymbol.h" #include "util.h" using namespace std; using namespace GiNaC; int main() { cindsymbol R = cindsymbol("R",0), G = cindsymbol("G",0), B = cindsymbol("B",0), X = cindsymbol("X",0); ex pl = X + pow(X, 3) + R + G + B; cout << R << G << B << " " << pl << endl; cout << pl.subs({X}, {pow(X,4)}) << endl; return 0; } The output of this program is as follows (the print functions from the subclass are being invoked correctly and do not print the zero value): RGB X+R+G+B+X^3 4*X^4+X^12 Obviously GiNaC thinks that R=X, B=X and G=X. How do I get this comparison to work? Do I need to override some function in my subclass? The function "cindsymbol" goes as follows (copied from the tutorial): static map<string, cindsymbol> directory; cindsymbol cindsym(const std::string & initname, unsigned cl) { std::string key = initname + "|" + std::to_string(cl); map<std::string, cindsymbol>::iterator it = directory.find(key); if (it != directory.end()) return it->second; else{ std::pair<std::string, cindsymbol> ent = make_pair(key, cindsymbol(initname, cl)); return directory.insert(ent).first->second; } } Can you advise me here? Thanks! Best regards, Marko Riedel

2 1

GiNaC 1.8.6 released
by Richard B. Kreckel 09 Feb '23

09 Feb '23

Hi, GiNaC 1.8.6 has been tagged and released. The News: * Fix wrong numeric info on transcendental functions. * Fix crash of evaluation of binomial(n, k) with negative integer n, k. Happy hacking! -richy. PS: If anyone knows how to push a commit together with the tags using git, do let me know. Yes, I've googled it up and tried several things but it never worked. -- Richard B. Kreckel <https://in.terlu.de/~kreckel/>

2 1

binomial coefficients
by Marko Riedel 07 Feb '23

07 Feb '23

Greetings. I would like to point out that binomial(-5, -7) will produce a segmentation fault. Compare to Maple which says that when binomial(n,k) has k<0 but n-k>=0 we should use binomial(n, n-k) which gives 15 in this case. Best regards, Marko

2 1

Segmentation Fault when calculate cos(cos(1/2))
by 翁湛阳 07 Feb '23

07 Feb '23

Dear GiNaC Team, When I tried to use GiNaC to calculate cos(cos(1/2)),I got a Segmentation Fault. My environment is as following: CPU: Intel i7-10750H Windows10 21H2 gcc version 12.2.0 (MinGW-W64 x86_64-ucrt-posix-seh) GiNaC 1.8.5 CLN 1.3.6 GMP 6.2.1 ArchLinux(in VMwave) Kernel: 6.1.9-arch1-1 GiNaC 1.8.4 installed by pacman CLN 1.3.6 installed by pacman gcc version 12.2.1 20230111(GCC) Relative code file is attached The program crash at line 8.I tried to compile and run the same program at ArchLinux(GCC 12.2.1 GiNaC 1.8.4 CLN 1.3.6),and also got a segmentation fault. Besides , I tried to use ginsh to execute “cos(cos(1/2));” in ArchLinux ,it gave a error: Internal error: statement in file float/misc/cl_F_extendsqrt.cc, line 30 has been reached!! The same error also occurred at Windows,but it just exit without the information above. After that, I built some lower version of GiNaC in Windows10 and re-run my code.I found the bug disappeared at GiNaC 1.8.3 . After reading the commit history I rollback ginac\inifcns_trans.cpp to version 1.8.3 and rebuilt ,the bug disappeared. All my best, -Weng Zhanyang.

2 1

subscripted variables
by Marko Riedel 06 Feb '23

06 Feb '23

Greetings to all. I am at this time considering to use GiNaC for some computations from group theory. I am not a C++ programmer. I would like to ask for assistance with the following task. I have a polynomial in several variables like a[1], a[2], a[5] etc. I understand I can initialize these by creating a symbol object with the variable name. The problem is, I need an easy and conformal means of accessing the integer in the variable name. This could be done by matching the name to a regular expression and extracting the value but I imagine there must be a better way. I am also looking for a means of extracting a set of variables that are present in a polynomial, or the degree these variables have. Furthermore, another question, does the C++ standard library have a means of memoizing functions? Thank you for your assistance. Marko Riedel

2 2

GiNaC 1.8.5 released
by Richard B. Kreckel 01 Jan '23

01 Jan '23

Happy New Year! As in 2022, New Year's Day brings a new GiNaC release. The only change is in factor(): * Speed up multivariate polynomial factorization; fix it in some rare corner cases where it didn't previously terminate. ...and happy hacking! -richy. -- Richard B. Kreckel <https://in.terlu.de/~kreckel/>

1 0

Failed Tests
by Robin Marzucca 21 Sep '22

21 Sep '22

Dear GiNaC Team, I just pulled the latest version of GiNaC from git and followed the installation steps as described. When running ‘make check’ I found two failed tests, the logs to which I attach. I just got a new laptop with an Apple Silicon processor, so the issue might still be related to that somehow, but most of the checks went through. Cheers, Robin

2 1

GiNaC 1.8.4 released
by Richard B. Kreckel 19 Sep '22

19 Sep '22

Hi, GiNaC 1.8.4 has been tagged and released with some unplanned git noise. The News: * Complete the support for sqrfree_parfrac(). * Add info methods for transcendental functions. -richy. -- Richard B. Kreckel <https://in.terlu.de/~kreckel/>

1 0

Potential issue in comparator in test_antipode
by Yuri Gribov 06 Jun '22

06 Jun '22

Hi, The comparison in `node::operator<` in time_antipode.cpp seems to be missing a dereference operation in return statement: if (!(*vert==*n.vert)) return (vert<n.vert); This causes the comparator to be non-transitive: (gdb) p *keys[5] < *keys[4] $5 = true (gdb) p *keys[4] < *keys[8] $6 = true (gdb) p *keys[5] < *keys[8] $7 = false which violates the strict weak ordering requirements of std::map (https://en.cppreference.com/w/cpp/named_req/Compare) Such violations may in practice result in non-deterministic behavior or crashes. The issue was found with [SortChecker++ tool](https://github.com/yugr/sortcheckxx). Best regards, Yuri Gribov

2 1