diff --git a/Modules/Numerics/FEM/include/itkImageToRectilinearFEMObjectFilter.hxx b/Modules/Numerics/FEM/include/itkImageToRectilinearFEMObjectFilter.hxx
index d1d528dfccd..a08e3878b6a 100644
--- a/Modules/Numerics/FEM/include/itkImageToRectilinearFEMObjectFilter.hxx
+++ b/Modules/Numerics/FEM/include/itkImageToRectilinearFEMObjectFilter.hxx
@@ -299,7 +299,7 @@ ImageToRectilinearFEMObjectFilter<TInputImage>::PrintSelf(std::ostream & os, Ind
   os << indent << "NumberOfElements: " << m_NumberOfElements << std::endl;
   os << indent << "PixelsPerElement: " << m_PixelsPerElement << std::endl;
 
-  itkPrintSelfObject(Material);
+  itkPrintSelfObjectMacro(Material);
   itkPrintSelfObjectMacro(Element);
 }
 
diff --git a/Modules/Registration/FEM/include/itkFEMFiniteDifferenceFunctionLoad.hxx b/Modules/Registration/FEM/include/itkFEMFiniteDifferenceFunctionLoad.hxx
index 3b29821d7b3..bf20560678e 100644
--- a/Modules/Registration/FEM/include/itkFEMFiniteDifferenceFunctionLoad.hxx
+++ b/Modules/Registration/FEM/include/itkFEMFiniteDifferenceFunctionLoad.hxx
@@ -355,7 +355,15 @@ FiniteDifferenceFunctionLoad<TMoving, TFixed>::PrintSelf(std::ostream & os, Inde
   os << indent << "SolutionIndex2: " << m_SolutionIndex2 << std::endl;
   os << indent << "Gamma: " << m_Gamma << std::endl;
 
-  itkPrintSelfObjectMacro(Solution);
+  os << indent << "Solution: ";
+  if (m_Solution != nullptr)
+  {
+    os << m_Solution << std::endl;
+  }
+  else
+  {
+    os << "(null)" << std::endl;
+  }
 
   os << indent << "GradSigma: " << static_cast<typename NumericTraits<Float>::PrintType>(m_GradSigma) << std::endl;
   os << indent << "Sign: " << m_Sign << std::endl;
diff --git a/Modules/Registration/FEM/include/itkFEMRegistrationFilter.hxx b/Modules/Registration/FEM/include/itkFEMRegistrationFilter.hxx
index c23898a5e2b..a6c1d88c221 100644
--- a/Modules/Registration/FEM/include/itkFEMRegistrationFilter.hxx
+++ b/Modules/Registration/FEM/include/itkFEMRegistrationFilter.hxx
@@ -1468,7 +1468,8 @@ FEMRegistrationFilter<TMovingImage, TFixedImage, TFemObject>::PrintSelf(std::ost
   os << indent << "FileCount: " << m_FileCount << std::endl;
   os << indent << "CurrentLevel: " << m_CurrentLevel << std::endl;
   os << indent << "CurrentLevelImageSize: "
-     << static_cast<typename NumericTraits<typename FixedImageType::SizeType>::PrintType>(m_CurrentLevel) << std::endl;
+     << static_cast<typename NumericTraits<typename FixedImageType::SizeType>::PrintType>(m_CurrentLevelImageSize)
+     << std::endl;
   os << indent << "WhichMetric: " << m_WhichMetric << std::endl;
 
   os << indent << "MeshPixelsPerElementAtEachResolution: " << m_MeshPixelsPerElementAtEachResolution << std::endl;
@@ -1526,7 +1527,11 @@ FEMRegistrationFilter<TMovingImage, TFixedImage, TFemObject>::PrintSelf(std::ost
   itkPrintSelfObjectMacro(Metric);
   itkPrintSelfObjectMacro(FEMObject);
 
-  os << indent << "LandmarkArray: " << m_LandmarkArray << std::endl;
+  os << indent << "LandmarkArray: ";
+  for (const auto & elem : m_LandmarkArray)
+  {
+    os << indent.GetNextIndent() << "[" << &elem - &*(m_LandmarkArray.begin()) << "]: " << *elem << std::endl;
+  }
 
   itkPrintSelfObjectMacro(Interpolator);