1 // This file is part of Eigen, a lightweight C++ template library 2 // for linear algebra. 3 // 4 // Copyright (C) 2009 Benoit Jacob <jacob.benoit.1 (at) gmail.com> 5 // 6 // This Source Code Form is subject to the terms of the Mozilla 7 // Public License v. 2.0. If a copy of the MPL was not distributed 8 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 10 #include "binary_library.h" 11 #include "stdio.h" 12 13 void demo_MatrixXd() 14 { 15 struct C_MatrixXd *matrix1, *matrix2, *result; 16 printf("*** demo_MatrixXd ***\n"); 17 18 matrix1 = MatrixXd_new(3, 3); 19 MatrixXd_set_zero(matrix1); 20 MatrixXd_set_coeff(matrix1, 0, 1, 2.5); 21 MatrixXd_set_coeff(matrix1, 1, 0, 1.4); 22 printf("Here is matrix1:\n"); 23 MatrixXd_print(matrix1); 24 25 matrix2 = MatrixXd_new(3, 3); 26 MatrixXd_multiply(matrix1, matrix1, matrix2); 27 printf("Here is matrix1*matrix1:\n"); 28 MatrixXd_print(matrix2); 29 30 MatrixXd_delete(matrix1); 31 MatrixXd_delete(matrix2); 32 } 33 34 // this helper function takes a plain C array and prints it in one line 35 void print_array(double *array, int n) 36 { 37 struct C_Map_MatrixXd *m = Map_MatrixXd_new(array, 1, n); 38 Map_MatrixXd_print(m); 39 Map_MatrixXd_delete(m); 40 } 41 42 void demo_Map_MatrixXd() 43 { 44 struct C_Map_MatrixXd *map; 45 double array[5]; 46 int i; 47 printf("*** demo_Map_MatrixXd ***\n"); 48 49 for(i = 0; i < 5; ++i) array[i] = i; 50 printf("Initially, the array is:\n"); 51 print_array(array, 5); 52 53 map = Map_MatrixXd_new(array, 5, 1); 54 Map_MatrixXd_add(map, map, map); 55 Map_MatrixXd_delete(map); 56 57 printf("Now the array is:\n"); 58 print_array(array, 5); 59 } 60 61 int main() 62 { 63 demo_MatrixXd(); 64 demo_Map_MatrixXd(); 65 } 66
