C++23 compatibility fixes

Author: dancasarin

include/kfr/base/expression.hpp

@@ -147,7 +147,7 @@ template <typename E>
 concept expression_argument = input_expression<E> && expression_traits<E>::explicit_operand;
 
 template <typename... E>
-concept expression_arguments = ((input_expression<E> && expression_traits<E>::explicit_operand) || ...);
+concept expression_arguments = (input_expression<E> && ...) && (expression_argument<E> || ...);
 
 template <typename T>
 concept input_output_expression = input_expression<T> && output_expression<T>;

include/kfr/base/simd_expressions.hpp

@@ -506,64 +506,64 @@ KFR_INTRINSIC expression_make_function<fn::satsub, E1, E2> satsub(E1&& x, E2&& y
 }
 
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator+=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator+=(E1&& e1, E2&& e2)
 {
     process(e1, operator+(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator-=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator-=(E1&& e1, E2&& e2)
 {
     process(e1, operator-(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator*=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator*=(E1&& e1, E2&& e2)
 {
     process(e1, operator*(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator/=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator/=(E1&& e1, E2&& e2)
 {
     process(e1, operator/(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator%=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator%=(E1&& e1, E2&& e2)
 {
     process(e1, operator%(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator|=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator|=(E1&& e1, E2&& e2)
 {
     process(e1, operator|(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator&=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator&=(E1&& e1, E2&& e2)
 {
     process(e1, operator&(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator^=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator^=(E1&& e1, E2&& e2)
 {
     process(e1, operator^(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator<<=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator<<=(E1&& e1, E2&& e2)
 {
     process(e1, operator<<(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 template <input_output_expression E1, input_expression E2>
-KFR_INTRINSIC E1& operator>>=(E1&& e1, E2&& e2)
+KFR_INTRINSIC std::remove_reference_t<E1> operator>>=(E1&& e1, E2&& e2)
 {
     process(e1, operator>>(e1, e2));
-    return e1;
+    return std::forward<E1>(e1);
 }
 
 } // namespace KFR_ARCH_NAME

src/dft/ft.hpp

@@ -1517,12 +1517,18 @@ KFR_INTRINSIC void cwrite_transposed(cbool_t<transposed>, complex<T>* ptr, vec<T
 template <size_t I, size_t radix, typename T, size_t N, size_t width = N / 2>
 KFR_INTRINSIC vec<T, N> mul_tw(cbool_t<false>, const vec<T, N>& x, const complex<T>* twiddle)
 {
-    return I == 0 ? x : cmul(x, cread<width>(twiddle + width * (I - 1)));
+    if constexpr (I == 0)
+        return x;
+    else
+        return cmul(x, cread<width>(twiddle + width * (I - 1)));
 }
 template <size_t I, size_t radix, typename T, size_t N, size_t width = N / 2>
 KFR_INTRINSIC vec<T, N> mul_tw(cbool_t<true>, const vec<T, N>& x, const complex<T>* twiddle)
 {
-    return I == 0 ? x : cmul_conj(x, cread<width>(twiddle + width * (I - 1)));
+    if constexpr (I == 0)
+        return x;
+    else
+        return cmul_conj(x, cread<width>(twiddle + width * (I - 1)));
 }
 
 // Non-final