itkArray2D.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright NumFOCUS
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         https://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/
#ifndef itkArray2D_h
#define itkArray2D_h
 
#include "itkMacro.h"
#include "itkIntTypes.h"
#include "vnl/vnl_matrix.h"
 
namespace itk
{
template <typename TValue>
class ITK_TEMPLATE_EXPORT Array2D : public vnl_matrix<TValue>
{
public:
  using ValueType = TValue;
  using Self = Array2D;
  using VnlMatrixType = vnl_matrix<TValue>;

  Array2D() = default;

  Array2D(unsigned int numberOfRows, unsigned int numberOfCols);

  Array2D(unsigned int numberOfRows, unsigned int numberOfCols, const TValue & initialValue);

  Array2D(const Self &) = default;

  Array2D(Self && array) noexcept
    : vnl_matrix<TValue>(std::move(array))
  {
    // Note: GCC <= 9.5 does not yet support "defaulting" (`= default`) this `noexcept` move-constructor.
  }

  Array2D(const VnlMatrixType & matrix);

  Self &
  operator=(const Self &) = default;

  Self &
  operator=(Self && array) noexcept
  {
    // Note: GCC <= 9.5 does not yet support "defaulting" (`= default`) this `noexcept` move-assignment operator.
    this->VnlMatrixType::operator=(std::move(array));
    return *this;
  }

  Self &
  operator=(const VnlMatrixType & matrix);
 
  void
  Fill(const TValue & v)
  {
    this->fill(v);
  }

  [[nodiscard]] const TValue &
  GetElement(SizeValueType row, SizeValueType col) const
  {
    return this->operator()(static_cast<unsigned int>(row), static_cast<unsigned int>(col));
  }

  void
  SetElement(SizeValueType row, SizeValueType col, const TValue & value)
  {
    this->operator()(static_cast<unsigned int>(row), static_cast<unsigned int>(col)) = value;
  }

  void
  SetSize(unsigned int m, unsigned int n);

  ~Array2D() override = default;
};
 
template <typename TValue>
std::ostream &
operator<<(std::ostream & os, const Array2D<TValue> & arr)
{
  const unsigned int numberOfRows = arr.rows();
  const unsigned int numberOfColumns = arr.cols();
 
  for (unsigned int r = 0; r < numberOfRows; ++r)
  {
    os << '[';
    if (numberOfColumns >= 1)
    {
      const unsigned int lastColumn = numberOfColumns - 1;
      for (unsigned int c = 0; c < lastColumn; ++c)
      {
        os << arr(r, c) << ", ";
      }
      os << arr(r, lastColumn);
    }
    os << ']' << std::endl;
  }
 
  return os;
}
 
// declaration of specialization
template <>
ITKCommon_EXPORT std::ostream &
                 operator<<(std::ostream & os, const Array2D<float> & arr);
template <>
ITKCommon_EXPORT std::ostream &
                 operator<<(std::ostream & os, const Array2D<double> & arr);
 
} // namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkArray2D.hxx"
#endif
 
#endif