July 2024 - GiNaC-list

Patch for ginac-1.8.7 (MacOS 14.6,Apple Silicon)
by 足立聡 02 Aug '24

02 Aug '24

Dear Developers of GiNaC, [0] Thank you very much for developing and maintaining GiNaC, which is a very nice system to write C++ programs for symbolic computation. [1] During recent several months, I had a trouble to build ginac-1.8.7 on my MacBookPro (Apple Silicon). The configuration and compilation steps succeeded, but the check step failed at two tests as follows: ./configure --> checking for a BSD-compatible install... /Ext/sw/bin/ginstall -c ... Configuration of GiNaC 1.8.7 done. Now type "make". make --> Making all in ginac ... /Ext/sw/bin/makeinfo --no-split --html ./ginac-examples.texi make check --> Making check in ginac ... PASS: exam_paranoia ... PASS: exam_color FAIL: exam_clifford FAIL: exam_archive PASS: exam_structure ... PASS: time_parser ============================================================================ Testsuite summary for GiNaC 1.8.7 ============================================================================ # TOTAL: 61 # PASS: 59 # SKIP: 0 # XFAIL: 0 # FAIL: 2 # XPASS: 0 # ERROR: 0 ============================================================================ See check/test-suite.log Please report to ginac-list(a)ginac.de ============================================================================ My environment is as follows: Hardware: MacBookPro (Apple Silicon,14inch,2023) OS: MacOS 14.6 Compiler: Xcode 15.4 g++ --version --> Apple clang version 15.0.0 (clang-1500.3.9.4) Target: arm64-apple-darwin23.6.0 Thread model: posix InstalledDir: /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin [2] I have investigated what caused this trouble and found that the phenomenon described in "ginac/hash_seed.h" occurs not only in Win32 but also in my MacOS 14.6. typeid(ex.bp).name() returns always correct C string as value but its address varies from a call to another call. Accordingly, the problem can be fixed as follows: *** ginac/hash_seed.h.org Sun Jan 1 03:47:20 2023 --- ginac/hash_seed.h Tue Jul 30 23:02:36 2024 *************** *** 8,14 **** #include <typeinfo> #include <cstring> #include "utils.h" ! #ifdef _WIN32 #define GINAC_HASH_USE_MANGLED_NAME 1 #endif --- 8,14 ---- #include <typeinfo> #include <cstring> #include "utils.h" ! #if defined(_WIN32) || defined(__APPLE__) #define GINAC_HASH_USE_MANGLED_NAME 1 #endif [3] Now, everything is fine. make distclean ./configure ... make ... make check --> Making check in ginac ... PASS: exam_paranoia PASS: exam_heur_gcd PASS: exam_match PASS: exam_parser PASS: exam_numeric PASS: exam_relational PASS: exam_powerlaws PASS: exam_collect PASS: exam_inifcns PASS: exam_inifcns_nstdsums PASS: exam_inifcns_elliptic PASS: exam_differentiation PASS: exam_polygcd PASS: exam_collect_common_factors PASS: exam_normalization PASS: exam_sqrfree PASS: exam_factor PASS: exam_pseries PASS: exam_matrices PASS: exam_lsolve PASS: exam_indexed PASS: exam_color PASS: exam_clifford PASS: exam_archive PASS: exam_structure PASS: exam_misc PASS: exam_pgcd PASS: exam_mod_gcd PASS: exam_chinrem_gcd PASS: exam_function_exvector PASS: exam_real_imag PASS: check_numeric PASS: check_inifcns PASS: check_matrices PASS: check_lsolve PASS: check_cra PASS: time_dennyfliegner PASS: time_gammaseries PASS: time_vandermonde PASS: time_toeplitz PASS: time_lw_A PASS: time_lw_B PASS: time_lw_C PASS: time_lw_D PASS: time_lw_E PASS: time_lw_F PASS: time_lw_G PASS: time_lw_H PASS: time_lw_IJKL PASS: time_lw_M1 PASS: time_lw_M2 PASS: time_lw_N PASS: time_lw_O PASS: time_lw_P PASS: time_lw_Pprime PASS: time_lw_Q PASS: time_lw_Qprime PASS: time_antipode PASS: time_fateman_expand PASS: time_uvar_gcd PASS: time_parser ============================================================================ Testsuite summary for GiNaC 1.8.7 ============================================================================ # TOTAL: 61 # PASS: 61 # SKIP: 0 # XFAIL: 0 # FAIL: 0 # XPASS: 0 # ERROR: 0 sudo make install ... [4] There may be other people who run into the same problem as I. In that case, please utilize this information. Finally, GiNaC is a really nice system to write programs for symbolic computation. Thank you very much for publishing GiNaC to help us. Sincerely yours, Satoshi Adachi

2 1

I want to Add Integrating by Parts in GiNaC Symbolic Integration
by DS Glanzsche 01 Aug '24

01 Aug '24

Dear all, I just learned about GiNaC, and I want to add integration by parts for a bit of complex integration like integral of *cos nx * x*, for definite and indefinite integral. I think polynomial integration in GiNaC works amazing, but if we can add all other symbolic integration it will be better. Should I look and modify the source code *integral.cpp*? I never really modified open source code before in my life. Thanks.

2 3

Memory leaks
by tstk3268 30 Jul '24

30 Jul '24

HI, GiNaC seems to have memory leaks that show up in a no-op program under Windows/MSVC using the _CrtSetDbgFlag mechanism. This seems to be due to static initializers. Is this a known and intended feature? It makes developing a proper leak free application pretty difficult.

2 1

Access denied when cloning git repository
by Matthias Volk 08 Jul '24

08 Jul '24

Hello, I am trying the clone the GiNaC and CLN git repositories. Both commands git clone git://www.ginac.de/cln.git and git clone git://www.ginac.de/ginac.git fail with the following error message: fatal: remote error: access denied or repository not exported: /ginac.git Am I using a wrong command or are the git repositories currently not available? Best regards, Matthias

2 2

degree function
by Santos Jha 04 Jul '24

04 Jul '24

Greetings, I am new to GiNac. I was trying to use the degree function, but I got a weird result. My initial expression is -4*(-3+x)*(-1+x)+9/2*(-2+x)*(-1+x)+1/2*(-2+x)*(-3+x) After the "expand" function call, it returns 9/2*x^2-5/2*x-4*x^2+16*x-27/2*x+1/2*x^2 My goal is to add/subtract the coefficient of equal degree terms. I could not find any function for that. ( if it is already there, please point me) I wanted to write a function to achieve it. To do so. I am taking each term and try to see if for (size_t i = 0; i != poly.nops(); ++i) { // Here poly is polynomial as above ex subExpres=poly.op(i); // I get individual terms cout << "is polynomial="<< is_polynomial(subExpres,symb) <<endl; // GiNaC::ex pow2=pow(symb,2); cout << subExpres<< " degree=" << subExpres.degree(symb) << endl; } It determines the subexpression as a polynomial but can not determine the degree. E.g is polynomial=1 9/2*x^2 degree=0 // Here degree returned is wrong. Your thoughts will be appreciated. Regards, Santos

2 2