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

#include <MMatrix.h>

MMatrix();
MMatrix(const MMatrix& arg);
MMatrix(long nrows, long ncols = DEF_SIZE,
        Integral::MTYPE type = Integral::DEF_MTYPE);
boolean assign(TIntegral value);
boolean assign(long num_rows, long num_cols, long* arg, 
               Integral::MTYPE type = Integral::DEF_MTYPE);
boolean setCapacity(const MMatrix& prototype_matrix,
                    boolean preserve_values = true,
                    Integral::MTYPE type = Integral::UNCHANGED);
boolean setDimensions(const MMatrix& prototype_matrix,
                      boolean preserve_values = true,
                      Integral::MTYPE type = Integral::UNCHANGED);
boolean mult(MMatrix& m1, MMatrix& m2);
boolean multv(TVector& output_vector, const TVector& input_vector);
boolean vmult(TVector& output_vector, const TVector& input_vector);
double determinant() const;
boolean inverse();
boolean transpose();

MMatrix<Double, double> mat0(3, 3, Integral::FULL);
MMatrix<Double, double> mat2;

double det;

mat0.assign(3, 3, L"2.1, 8.3, 4.5, 4.6, 0, 8.1, 5.5, 3.2, 1.7");
MMatrix<Double, double> mat1(4, 4, Integral::SYMMETRIC);

MVector<Double, double> vec_in;
MVector<Double, double> vec_out1;
MVector<Double, double> vec_out2;

mat1.assign(4, 4, L"1.1, 4.5, 6.5, 3.6, 0.3, 6.7, 8.1, 2.3, 5.5, 9.2",
            Integral::SYMMETRIC);

vec_in.assign(L"2.2, 4.3, 6.5");
mat1.setDimensions(3, 3); 

mat0.multv(vec_out1, vec_in);
mat1.vmult(vec_out2, vec_in);

mat0.transpose();
mat2.mult(mat0, mat1);

det = mat2.determinant();
mat2.inverse();
mat0.setCapacity(0);

• MMatrix is never used directly. See the other matrix types for examples of computational methods. Here, we present an example that demonstrations the basic I/O capabilities:

  long data[9] = {1, 2, 3, 4, 5, 6, 7, 8, 9};
         
  MMatrix<Long, long> mat1;
  mat1.assign(3, 3, data);
  MMatrix<Long, long> mat2;
  
  String file_name(L"foo");
  Sof tmp_file;
  tmp_file.open(file_name, File::WRITE_ONLY, File::TEXT);
         
  // write the values
  //
  mat1.write(tmp_file, (long)0);
  
  // close the files
  //
  tmp_file.close();
         
  // open the files in read mode
  //
  tmp_file.open(file_name);
         
  // read the values back in
  //
  if (!mat2.read(tmp_file, (long)0) || (!mat2.almostEqual(mat1))) {
         mat1.debug(L"mat g1");
         mat2.debug(L"mat2");
         return Error::handle(mat2.name(), L"read-text", Error::TEST,
                              __FILE__,__LINE__);
  }
  tmp_file.close();