From db92e1e50c2e0f23de028254f393f8deb84910f4 Mon Sep 17 00:00:00 2001
From: J-Donald Tournier <jdtournier@gmail.com>
Date: Thu, 2 May 2024 11:15:04 +0100
Subject: [PATCH] Adapter::Reslice: fix handling of complex types

---
 core/adapter/reslice.h | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/core/adapter/reslice.h b/core/adapter/reslice.h
index 792f75aa11..775b9ea5a8 100644
--- a/core/adapter/reslice.h
+++ b/core/adapter/reslice.h
@@ -43,8 +43,8 @@ typename std::enable_if<!std::is_same<value_type, bool>::value && std::is_integr
 }
 
 template <typename value_type>
-typename std::enable_if<std::is_floating_point<value_type>::value, value_type>::type inline normalise(
-    const default_type sum, const default_type norm) {
+typename std::enable_if<!std::is_same<value_type, bool>::value && !std::is_integral<value_type>::value,
+                        value_type>::type inline normalise(const default_type sum, const default_type norm) {
   return (sum * norm);
 }
 } // namespace
@@ -98,6 +98,7 @@ template <template <class ImageType> class Interpolator, class ImageType>
 class Reslice : public ImageBase<Reslice<Interpolator, ImageType>, typename ImageType::value_type> {
 public:
   using value_type = typename ImageType::value_type;
+  using internal_type = typename std::conditional<is_complex<value_type>::value, cdouble, double>::type;
 
   template <class HeaderType>
   Reslice(const ImageType &original,
@@ -175,7 +176,7 @@ class Reslice : public ImageBase<Reslice<Interpolator, ImageType>, typename Imag
     using namespace Eigen;
     if (oversampling) {
       Vector3d d(x[0] + from[0], x[1] + from[1], x[2] + from[2]);
-      default_type sum(0.0);
+      internal_type sum(0.0);
       Vector3d s;
       for (uint32_t z = 0; z < OS[2]; ++z) {
         s[2] = d[2] + z * inc[2];