1       2       3       4       5       6       7       8       9      10      11      12      13      14      15      16      17      18      19      20      21      22      23      24      25      26      27      28      29      30      31      32      33      34      35      36      37      38      39      40      41      42      43      44      45      46      47      48      49      50      51      52      53      54      55      56      57      58      59      60 #include  /// ///  With header __flens.cxx__ all of FLENS gets included. /// ///  :links:  __flens.cxx__ -> file:flens/flens.cxx #include  using namespace std; using namespace flens; typedef double   T; int main() {     ///     ///  Define some convenient typedefs for the matrix/vector types     ///  of our system of linear equations.     ///     typedef GeMatrix >           Matrix;     typedef DenseVector >              Vector;     ///     ///  We also need an extra vector type for the pivots.  The type of the     ///  pivots is taken for the system matrix.     ///     typedef Matrix::IndexType                   IndexType;     typedef DenseVector >      IndexVector;     ///     ///  Set up the baby problem ...     ///     const IndexType n = 4;     Matrix         A(n,n);     Vector         b(n);     IndexVector    piv(n);     A =  2,   3,  -1,   0,         -6,  -5,   0,   2,          2,  -5,   6,  -6,          4,   6,   2,  -3;     b = 20,        -33,        -43,         49;     cerr << "A = " << A << endl;     cerr << "b = " << b << endl;     ///     /// And solve it with __lapack::sv__     ///     /// :links: __lapack::sv__ -> doc:flens/lapack/ge/sv     lapack::sv(A, piv, b);     cerr << "x = " << b << endl; }