ginac库矩阵运算

1.程序源码

#include <iostream>
#include <ginac/ginac.h>

using namespace std;
using namespace GiNaC;

int main()
{
        matrix A = matrix(2,2,lst(1,2,3,4));
        matrix B = matrix(2,2,lst(-1,0,2,1));
        matrix C = matrix(2,2,lst(8,4,2,1));

        //result = A * B - C * 2
        matrix D = A.mul(B).sub(C.mul_scalar(2));
        ex E = A * B - C * 2;

        cout << D << endl;
        cout << E.evalm() << endl;

        return 0;
}

2.编译源码

$ g++ -o example example.c -lginac -lcln -I/usr/local/include -L/usr/local/lib64

3.运行结果

./example
[[-13,-6],[1,2]]
[[-13,-6],[1,2]]