]> Creatis software - clitk.git/blob - tools/clitkMatrixMultiply.cxx
Merge branch 'master' of git.creatis.insa-lyon.fr:clitk
[clitk.git] / 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 //--------------------------------------------------------------------