February 2013 - GiNaC-list

cmake msvc build
by Philipp Kraus 25 Feb '13

25 Feb '13

Hello, I love GiNaC and I can build it with MSYS very well. My question is for a feature request is: Can a cmake build added, so that I can build GiNaC native under MSVC? I think this can be a nice solution. Thanks Phil

3 3

subs does not work
by Zhao Li 21 Feb '13

21 Feb '13

Here I am trying to make a substitution. For testing, the function "has" works as shown below. But the "subs" just does not implement the substitution. cout << has(temp, JF1(wild(1), GA(wild(5)), wild(2))*JF1(wild(3), GA(wild(5)), wild(4)) , has_options::algebraic ) << endl; cout << temp.subs( JF1(wild(1),GA(wild(5)),wild(2))*JF1(wild(3),GA(wild(5)),wild(4)) == 2, subs_options::algebraic ) << endl; The functions "JF1" and "GA" are both declared using DECLARE_FUNCTION and commutative. DECLARE_FUNCTION_1P( GA ); REGISTER_FUNCTION( GA, dummy() ); Best regards, Zhao Theory Division Institute of High Energy Physics Chinese Academy of Sciences

2 1

Creating a number of FUNCTION_1P
by Dale Lukas Peterson 20 Feb '13

20 Feb '13

I am trying to declare N GiNaC 1 parameter functions (DECLARE_FUNCTION_1P), each of which represents a different function of the same single parameter (in my case, this parameter represents time). My primary motivation for doing this is to be able to easily create a bunch of symbol-like objects that all depend on time so that I can trust that the chain rule will work properly when differentiating with respect to this single parameter. I think this has come up before, but I haven't seen a way to declare a vector or list of these functions without N explicit statements to DECLARE_FUNCTION_1P. Is there a simple way to do this in GiNaC? Thanks, Luke

1 1

Compilation under macos + question on matrices
by Laurent Bartholdi 16 Feb '13

16 Feb '13

Dear all, I'm new to ginac, and very eager to learn more about it. I installed the latest (1.6.2) version on a recent (mountain lion 12.2.1) system, using the automatic packaging system 'homebrew'. 1. There must be a bug somewhere along the line, because the most naive, simple program from the manual (and all other examples I tried) cause a malloc error: % cat x.C #include <iostream> #include <ginac/ginac.h> int main() { GiNaC::symbol v("v"); std::cout << v << std::endl; return 0; } % g++ x.C -lginac && ./a.out v a.out(31763) malloc: *** error for object 0x7fff7b419570: pointer being freed was not allocated *** set a breakpoint in malloc_error_break to debug Abort trap: 6 2. Setting matrix entries causes crash: ... GiNaC::matrix m(10,10); m[1][2] = 1; ... causes "libc++abi.dylib: terminate called throwing an exception" 3. I can't find the tensor product of matrices (aka Kronecker product) in the documentation. If M in an axb matrix and N is a cxd matrix, this would be the (ac)x(bd) matrix with entries the products of entries of M and N. While it's easy to implement myself, I'm surprised not to find it in the library. Did I miss something? I'm not sure I know well how this mailing list works, so I'd appreciate a reply directly to my email address. Thanks in advance, Laurent -- Prof. Dr. Laurent Bartholdi \ laurent.bartholdi<at>gmail<dot>com G.-A. Universität zu Göttingen \ Phone: +49 551 39 7826 Bunsenstraße 3-5 \ Secr: +49 551 39 7752 D-37073 Göttingen, Germany \ Fax: +49 551 39 22674

1 0