tools/clitkMatrixMultiply.cxx

   1 /*=========================================================================
   2   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
   3
   4   Authors belong to:
   5   - University of LYON              http://www.universite-lyon.fr/
   6   - Léon Bérard cancer center       http://www.centreleonberard.fr
   7   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
   8
   9   This software is distributed WITHOUT ANY WARRANTY; without even
  10   the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  11   PURPOSE.  See the copyright notices for more information.
  12
  13   It is distributed under dual licence
  14
  15   - BSD        See included LICENSE.txt file
  16   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  17 ===========================================================================**/
  18
  19 // clitk
  20 #include "clitkMatrixMultiply_ggo.h"
  21 #include "clitkTransformUtilities.h"
  22 #include "clitkIO.h"
  23 #include "clitkMatrix.h"
  24
  25 //--------------------------------------------------------------------
  26 int main(int argc, char * argv[])
  27 {
  28   // Init command line
  29   GGO(clitkMatrixMultiply, args_info);
  30   CLITK_INIT;
  31
  32   // Read matrix
  33   itk::Matrix<double, 4, 4> matrix1, matrix2;
  34   try {
  35     matrix1 = clitk::ReadMatrix3D(args_info.input1_arg);
  36     matrix2 = clitk::ReadMatrix3D(args_info.input2_arg);
  37   }
  38   catch (itk::ExceptionObject & err) {
  39     std::cerr << "Error reading "
  40               << args_info.input1_arg
  41               << " and "
  42               << args_info.input2_arg
  43               << std::endl;
  44     std::cerr << err.GetDescription() << std::endl;
  45     exit(-1);
  46   }
  47
  48   itk::Matrix<double, 4, 4> matrix(matrix2.GetVnlMatrix() * matrix1.GetVnlMatrix());
  49
  50   // Print matrix
  51   std::ofstream os;
  52   clitk::openFileForWriting(os, args_info.output_arg);
  53   os << clitk::Get4x4MatrixDoubleAsString(matrix, 16);
  54   os.close();
  55
  56   return EXIT_SUCCESS;
  57 }// end main
  58
  59 //--------------------------------------------------------------------