BUG: Fix warp filter exception and apply clang-format to src/

The ContinuousBorderWarpImageFilter incorrectly threw an
exception when TransformPhysicalPointToContinuousIndex returned
false (point outside image domain). Out-of-bounds points are
expected in warp operations — the subsequent clamping to boundary
indices provides the continuous-border behavior. This caused all
14 registration tests to fail with "Failed to warp moving image."

Also apply clang-format 19 to VariationalRegistrationMain.cxx
and win32_compatibility/getopt.h (lint CI failures).

Author: hjmjohnson

include/itkContinuousBorderWarpImageFilter.hxx

@@ -70,11 +70,9 @@ ContinuousBorderWarpImageFilter<TInputImage, TOutputImage, TDisplacementField>::
         point[j] += displacement[j];
       }
 
-      // project point into image region
-      if (!inputPtr->TransformPhysicalPointToContinuousIndex(point, contIndex))
-      {
-        itkExceptionMacro("Physical point " << point << " is outside the image domain");
-      }
+      // project point into image region (out-of-bounds is expected;
+      // the clamping below provides continuous-border behavior)
+      inputPtr->TransformPhysicalPointToContinuousIndex(point, contIndex);
 
       for (unsigned int j = 0; j < ImageDimension; j++)
       {

src/VariationalRegistrationMain.cxx

@@ -640,7 +640,8 @@ main(int argc, char * argv[])
   FunctionType::Pointer function;
   switch (forceType)
   {
-    case 0: {
+    case 0:
+    {
       DemonsFunctionType::Pointer demonsFunction = DemonsFunctionType::New();
       switch (forceDomain)
       {
@@ -658,7 +659,8 @@ main(int argc, char * argv[])
       function = demonsFunction;
     }
     break;
-    case 1: {
+    case 1:
+    {
       SSDFunctionType::Pointer ssdFunction = SSDFunctionType::New();
       switch (forceDomain)
       {
@@ -676,7 +678,8 @@ main(int argc, char * argv[])
       function = ssdFunction;
     }
     break;
-    case 2: {
+    case 2:
+    {
       NCCFunctionType::Pointer    nccFunction = NCCFunctionType::New();
       NCCFunctionType::RadiusType r;
       for (unsigned int dim = 0; dim < NCCFunctionType::ImageDimension; dim++)
@@ -719,19 +722,22 @@ main(int argc, char * argv[])
   RegularizerType::Pointer regularizer;
   switch (regularizerType)
   {
-    case 0: {
+    case 0:
+    {
       GaussianRegularizerType::Pointer gaussRegularizer = GaussianRegularizerType::New();
       gaussRegularizer->SetStandardDeviations(std::sqrt(regulVar));
       regularizer = gaussRegularizer;
     }
     break;
-    case 1: {
+    case 1:
+    {
       DiffusionRegularizerType::Pointer diffRegularizer = DiffusionRegularizerType::New();
       diffRegularizer->SetAlpha(regulAlpha);
       regularizer = diffRegularizer;
     }
     break;
-    case 2: {
+    case 2:
+    {
 #if defined(ITK_USE_FFTWD) || defined(ITK_USE_FFTWF)
       ElasticRegularizerType::Pointer elasticRegularizer = ElasticRegularizerType::New();
       elasticRegularizer->SetMu(regulMu);
@@ -742,7 +748,8 @@ main(int argc, char * argv[])
 #endif
     }
     break;
-    case 3: {
+    case 3:
+    {
 #if defined(ITK_USE_FFTWD) || defined(ITK_USE_FFTWF)
       CurvatureRegularizerType::Pointer curvatureRegularizer = CurvatureRegularizerType::New();
       curvatureRegularizer->SetAlpha(regulAlpha);
@@ -768,17 +775,20 @@ main(int argc, char * argv[])
   RegistrationFilterType::Pointer regFilter;
   switch (searchSpace)
   {
-    case 0: {
+    case 0:
+    {
       regFilter = RegistrationFilterType::New();
       break;
     }
-    case 1: {
+    case 1:
+    {
       DiffeomorphicRegistrationFilterType::Pointer diffeoRegFilter = DiffeomorphicRegistrationFilterType::New();
       diffeoRegFilter->SetNumberOfExponentiatorIterations(numberOfExponentiatorIterations);
       regFilter = diffeoRegFilter;
       break;
     }
-    case 2: {
+    case 2:
+    {
       SymmetricDiffeomorphicRegistrationFilterType::Pointer symmDiffeoRegFilter =
         SymmetricDiffeomorphicRegistrationFilterType::New();
       symmDiffeoRegFilter->SetNumberOfExponentiatorIterations(numberOfExponentiatorIterations);

src/win32_compatibility/getopt.h

@@ -43,8 +43,8 @@
 #ifdef _WIN32
 /* from <sys/cdefs.h> */
 #  ifdef __cplusplus
-#    define __BEGIN_DECLS                                                                                              \
-      extern "C"                                                                                                       \
+#    define __BEGIN_DECLS \
+      extern "C"          \
       {
 #    define __END_DECLS }
 #  else