[GiNaC-list] solving linear systems with ginac

[Sheplyakov Alexei]

Hello,

On Fri, Dec 16, 2005 at 01:22:50PM +0100, Kent-Andre Mardal wrote:
> > 
> > Could you please post the actual code which has this problem?
> > 
 
> I am not able to give you the code

Well, let's try the other way around.

Attached program (which calculates the inverse of your matrix) seems
to Work for me (TM). 

ginac-config --version
1.3.3
g++ --version
g++ (GCC) 4.0.3 20051201 (prerelease) (Debian 4.0.2-5)
Copyright (C) 2005 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

N.B. Passing solve_algo::gauss is essential, so GiNaC heuristics _does_
need to be improved.

So, I'd like to see an example (minimal) which demonstrates the problem
you have mentioned.

> (or that would involve downloading
> and kompiling the library SyFi which I guess you do not want.) 

Why not? 

-- 
All science is either physics or stamp collecting.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_matr.cpp.gz
Type: application/octet-stream
Size: 671 bytes
Desc: not available
Url : http://www.cebix.net/pipermail/ginac-list/attachments/20051216/11344b01/test_matr.cpp.obj
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 189 bytes
Desc: Digital signature
Url : http://www.cebix.net/pipermail/ginac-list/attachments/20051216/11344b01/attachment.pgp