1 #include <flens/flens.cxx>
2 #include <iostream>
3
4 using namespace flens;
5 using namespace std;
6
7 int
8 main()
9 {
10 typedef GeMatrix<FullStorage<double, ColMajor> > GEMatrix;
11
12 GEMatrix A(4,4);
13 A = 1, 2, 3, 4,
14 5, 6, 7, 8,
15 9, 8, 7, 6,
16 5, 4, 3, 20;
17
18 cout << "A = " << A << endl;
19
20 ///
21 /// ...
22 ///
23 auto U = A.upper();
24 cout << "U = " << U << endl;
25
26 ///
27 /// ...
28 ///
29 auto SU = U.symmetric();
30 cout << "SU = " << SU << endl;
31
32 ///
33 /// ...
34 ///
35 auto SL = A.lower().symmetric();
36 cout << "SL = " << SL << endl;
37
38 ///
39 /// ...
40 ///
41 auto L = A.lowerUnit();
42 cout << "L = " << L << endl;
43
44 return 0;
45 }
2 #include <iostream>
3
4 using namespace flens;
5 using namespace std;
6
7 int
8 main()
9 {
10 typedef GeMatrix<FullStorage<double, ColMajor> > GEMatrix;
11
12 GEMatrix A(4,4);
13 A = 1, 2, 3, 4,
14 5, 6, 7, 8,
15 9, 8, 7, 6,
16 5, 4, 3, 20;
17
18 cout << "A = " << A << endl;
19
20 ///
21 /// ...
22 ///
23 auto U = A.upper();
24 cout << "U = " << U << endl;
25
26 ///
27 /// ...
28 ///
29 auto SU = U.symmetric();
30 cout << "SU = " << SU << endl;
31
32 ///
33 /// ...
34 ///
35 auto SL = A.lower().symmetric();
36 cout << "SL = " << SL << endl;
37
38 ///
39 /// ...
40 ///
41 auto L = A.lowerUnit();
42 cout << "L = " << L << endl;
43
44 return 0;
45 }