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 double T;
11 typedef DenseVector<Array<T> > DEVector;
12 typedef GeMatrix<FullStorage<T, ColMajor> > GEMatrix;
13
14 DEVector x(3);
15 x = 1, 2, 3;
16
17 GEMatrix A(3,3);
18 A = 1, 2, 3,
19 5, 6, 7,
20 5, 4, 3;
21
22 ///
23 /// We turn on logging of the closure evaluation, i.e. the evaluation
24 /// of linear algebra expressions that are coded through overloaded
25 /// operators.
26 ///
27 verbose::ClosureLog::start("mylogfile");
28
29 ///
30 /// Compute $x = A x$
31 ///
32 x = A*x;
33
34 cout << "x = " << x << endl;
35
36 ///
37 /// Stop logging.
38 ///
39 verbose::ClosureLog::stop();
40
41 return 0;
42 }
2 #include <iostream>
3
4 using namespace flens;
5 using namespace std;
6
7 int
8 main()
9 {
10 typedef double T;
11 typedef DenseVector<Array<T> > DEVector;
12 typedef GeMatrix<FullStorage<T, ColMajor> > GEMatrix;
13
14 DEVector x(3);
15 x = 1, 2, 3;
16
17 GEMatrix A(3,3);
18 A = 1, 2, 3,
19 5, 6, 7,
20 5, 4, 3;
21
22 ///
23 /// We turn on logging of the closure evaluation, i.e. the evaluation
24 /// of linear algebra expressions that are coded through overloaded
25 /// operators.
26 ///
27 verbose::ClosureLog::start("mylogfile");
28
29 ///
30 /// Compute $x = A x$
31 ///
32 x = A*x;
33
34 cout << "x = " << x << endl;
35
36 ///
37 /// Stop logging.
38 ///
39 verbose::ClosureLog::stop();
40
41 return 0;
42 }