ITK/Examples/Math/Matrix: Difference between revisions

From KitwarePublic
< ITK‎ | Examples
Jump to navigationJump to search
(Created page with "How do you resize a matrix at runtime? How do you compute the inverse and get back an itkMatrix instead of a vnl_matrix? ==Matrix.cxx== <source lang="cpp"> #include <itkMatrix.h...")
 
(Deprecated content that is moved to sphinx)
 
(3 intermediate revisions by 2 users not shown)
Line 1: Line 1:
How do you resize a matrix at runtime? How do you compute the inverse and get back an itkMatrix instead of a vnl_matrix?
{{warning|1=The media wiki content on this page is no longer maintained.  The examples presented on the https://itk.org/Wiki/*  pages likely require ITK version 4.13 or earlier releases.  In many cases, the examples on this page no longer conform to the best practices for modern ITK versions.
}}


==Matrix.cxx==
[https://itk.org/ITKExamples[ITK Sphinx Examples]]
<source lang="cpp">
#include <itkMatrix.h>
#include <itkVector.h>
 
#include <iostream>
 
void Construct();
void ConstructRunTimeDims();
void Multiply();
void Inverse();
 
int main()
{
  Construct();
  Multiply();
  return 0;
}
 
void Construct()
{
  typedef itk::Matrix<double, 3, 3> MatrixType;
  MatrixType M;
  M(0,0) = 1.0;
  M(0,1) = 2.0;
  M(0,2) = 3.0;
  M(1,0) = 4.0;
  M(1,1) = 5.0;
  M(1,2) = 6.0;
  M(2,0) = 7.0;
  M(2,1) = 8.0;
  M(2,2) = 9.0;
 
  std::cout << "M: " << M << std::endl;
}
 
void ConstructRunTimeDims()
{
  /*
  int matrixSize = 3;
  typedef itk::Matrix<double, matrixSize, matrixSize> MatrixType;
  MatrixType M;
  M(0,0) = 1.0;
  M(0,1) = 2.0;
  M(0,2) = 3.0;
  M(1,0) = 4.0;
  M(1,1) = 5.0;
  M(1,2) = 6.0;
  M(2,0) = 7.0;
  M(2,1) = 8.0;
  M(2,2) = 9.0;
 
  std::cout << "M: " << M << std::endl;
  */
}
 
void Multiply()
{
  typedef itk::Matrix<double, 3, 3> MatrixType;
  MatrixType M;
  M(0,0) = 1.0;
  M(0,1) = 2.0;
  M(0,2) = 3.0;
  M(1,0) = 4.0;
  M(1,1) = 5.0;
  M(1,2) = 6.0;
  M(2,0) = 7.0;
  M(2,1) = 8.0;
  M(2,2) = 9.0;
 
  std::cout << "M: " << M << std::endl;
 
  typedef itk::Vector<double, 3> VectorType;
  VectorType V;
  V[0] = 1.0;
  V[1] = 2.0;
  V[2] = 3.0;
 
  std::cout << "V: " << V << std::endl;
 
  std::cout << "MV: " << M*V << std::endl;
}
 
void Inverse()
{
 
}
</source>
 
 
{{ITKCMakeLists|Matrix}}

Latest revision as of 19:43, 3 June 2019

Warning: The media wiki content on this page is no longer maintained. The examples presented on the https://itk.org/Wiki/* pages likely require ITK version 4.13 or earlier releases. In many cases, the examples on this page no longer conform to the best practices for modern ITK versions.

[ITK Sphinx Examples]