g++ [flags ...] file ... -l /isip/tools/lib/$ISIP_BINARY/lib_math_matrix.a

#include <MatrixFloat.h>

MatrixFloat(long nrows = DEF_SIZE, long ncols = DEF_SIZE,
             long type = DEF_TYPE);
MatrixFloat(const MatrixFloat& matrix);

float data0[6] = {1.2, 2.5, 3.6, 4.8, 5.1, 6.7};
float data1[6] = {4.2, 6.1, 7.0, 2.5, 3.3, 9.5};

MatrixFloat val0;
MatrixFloat val1(3, 3, Integral::LOWER_TRIANGULAR);
val0.assign(data0, 2, 3)
val1.assign(data1, 3, 3);

• none

• This example demonstrates basic linear algebra operations using the matrix class
  by solving the equation: w = x^-1y + z.

         
  // declare the operands: note that these matrices are of different types.
  // the class will handle type conversion automatically.
  //
  MatrixFloat x;
  x.assign(3, 3, L"1.0, 2.0, 3.0, 6.0, 5.0, 4.0, 7.0, 8.0, 9.0", Integral::FULL);
   
  MatrixFloat y;
  y.assign(3, 3, L"1.0, 1.0, 1.0", Integral::DIAGONAL);
    
  MatrixFloat z;
  z.setDimensions(3, 3);
  z.setValue(1, 2, 27.0);
         
  // declare an output matrix: not that we don't need to declare the size
  // of this matrix.
  //
  MatrixFloat w;
    
  // compute a matrix inverse in place: this destroys the input.
  //
  x.inverse();
    
  // compute the product of the inverse of x and y, and add z
  //
  x.mult(y);       
  w.add(x, z);
  
  // print the output
  //
  w.debug(L"the answer is");
  

  The output that results from this example is:
  
  

      [ 0.190476 0.52381 -0.333333  ]
  w = [ -0.952381 -0.619048 27.6667 ]
      [ 0.904762 0.238095 -0.333333 ]