Откопал интересную библиотеку Blitz.

Никто не пробовал решить с ее помощью решить СЛАУ или хотя бы привести матрицу к треугольному виду (в арсенале видел срезы, но упростят ли они жизнь или нет пока непонял)?

Пока догадался что нужно юзать Array

    blitz::Array<double,2> A(2,3);
    blitz::firstIndex i;
    blitz::secondIndex j;
    A = where(i!=j,10+i+j,1);
    std::cout << A << std::endl;

А есть ли подобные библиотеки?