partial fix for EFT processes after merging alphas (only workd with HRDCOD=1 for now)

epochX/cudacpp/CODEGEN/PLUGIN/CUDACPP_SA_OUTPUT/madgraph/iolibs/template_files/cpp_model_parameters_h.inc

@@ -13,74 +13,7 @@
 
 //==========================================================================
 
-namespace Parameters_%(model_name)s_dependentCouplings
-{
-  constexpr size_t ndcoup = %(ndcoup)d; // #couplings that vary event by event because they depend on the running alphas QCD
-%(idcoup)s
-  struct DependentCouplings_sv
-  {
-%(dcoupdecl)s
-  };
-#pragma GCC diagnostic push
-#pragma GCC diagnostic ignored "-Wunused-variable"  // e.g. <<warning: unused variable ‘mdl_G__exp__2’ [-Wunused-variable]>>
-#pragma GCC diagnostic ignored "-Wunused-parameter" // e.g. <<warning: unused parameter ‘G’ [-Wunused-parameter]>>
-#ifdef __CUDACC__
-#pragma nv_diagnostic push
-#pragma nv_diag_suppress 177 // e.g. <<warning #177-D: variable "mdl_G__exp__2" was declared but never referenced>>
-#endif
-  __host__ __device__ inline const DependentCouplings_sv computeDependentCouplings_fromG( const fptype_sv& G )
-  {
-    // Model parameters dependent on aS
-  %(dcoupsetdpar)s
-    // Model couplings dependent on aS
-    DependentCouplings_sv out;
-  %(dcoupsetdcoup)s
-    return out;
-  }
-#ifdef __CUDACC__
-#pragma GCC diagnostic pop
-#pragma nv_diagnostic pop
-#endif
-}
-
-//==========================================================================
-
-namespace Parameters_%(model_name)s_independentCouplings
-{
-  constexpr size_t nicoup = %(nicoup)d; // #couplings that are fixed for all events because they do not depend on the running alphas QCD
-%(iicoup)s
-}
-
-//==========================================================================
-
-#ifdef __CUDACC__
-namespace mg5amcGpu
-#else
-namespace mg5amcCpu
-#endif
-{
-#pragma GCC diagnostic push
-#pragma GCC diagnostic ignored "-Wunused-but-set-variable" // e.g. <<warning: variable ‘couplings_sv’ set but not used [-Wunused-but-set-variable]>>
-  // Compute the output couplings (e.g. gc10 and gc11) from the input gs
-  template<class G_ACCESS, class C_ACCESS>
-  __device__ inline void
-  G2COUP( const fptype gs[],
-          fptype couplings[] )
-  {
-    mgDebug( 0, __FUNCTION__ );
-    using namespace Parameters_%(model_name)s_dependentCouplings;
-    const fptype_sv& gs_sv = G_ACCESS::kernelAccessConst( gs );
-    DependentCouplings_sv couplings_sv = computeDependentCouplings_fromG( gs_sv );
-%(dcoupaccessbuffer)s%(dcoupkernelaccess)s%(dcoupcompute)s
-    mgDebug( 1, __FUNCTION__ );
-    return;
-  }
-#pragma GCC diagnostic pop
-}
-
-//==========================================================================
-
-#ifndef MGONGPU_HARDCODE_PARAM
+#ifndef MGONGPU_HARDCODE_PARAM // this is only supported in SM processes (e.g. not in EFT models) for the moment (#439)%(efterror)s
 
 #include "read_slha.h"
 
@@ -132,23 +65,43 @@ private:
 
 #else
 
+#include <cassert>
 #include <limits>
 
 // Hardcoded constexpr physics parameters
 namespace Parameters_%(model_name)s // keep the same name rather than HardcodedParameters_%(model_name)s for simplicity
 {
   // Constexpr implementation of sqrt (see https://stackoverflow.com/a/34134071)
-  double constexpr detailSqrtNewtonRaphson( double x, double curr, double prev )
+  double constexpr sqrtNewtonRaphson( double x, double curr, double prev )
   {
-    return curr == prev ? curr : detailSqrtNewtonRaphson( x, 0.5 * ( curr + x / curr ), curr );
+    return curr == prev ? curr : sqrtNewtonRaphson( x, 0.5 * ( curr + x / curr ), curr );
   }
-  double constexpr sqrtNR( double x )
+  double constexpr constexpr_sqrt( double x )
   {
     return x >= 0 // && x < std::numeric_limits<double>::infinity() // avoid -Wtautological-constant-compare warning in fast math
-      ? detailSqrtNewtonRaphson( x, x, 0 )
+      ? sqrtNewtonRaphson( x, x, 0 )
       : std::numeric_limits<double>::quiet_NaN();
   }
 
+  // Constexpr implementation of floor (see https://stackoverflow.com/a/66146159)
+  constexpr int constexpr_floor( double d )
+  {
+    const int i = static_cast<int>( d );
+    return d < i ? i - 1 : i;
+  }
+
+  // Constexpr implementation of pow
+  constexpr double constexpr_pow( double base, double exp )
+  {
+    // NB(1): this implementation of constexpr_pow requires exponent >= 0
+    assert( exp >= 0 ); // NB would fail at compile time with "error: call to non-‘constexpr’ function ‘void __assert_fail'"
+    // NB(2): this implementation of constexpr_pow requires an integer exponent
+    const int iexp = constexpr_floor( exp );
+    assert( static_cast<double>( iexp ) == exp ); // NB would fail at compile time with "error: call to non-‘constexpr’ function ‘void __assert_fail'"
+    // Iterative implementation of pow if exp is a non negative integer
+    return iexp == 0 ? 1 : base * constexpr_pow( base, iexp - 1 );
+  }
+
   // Model parameters independent of aS
   constexpr double zero = 0;
   constexpr double ZERO = 0;
@@ -176,4 +129,82 @@ namespace Parameters_%(model_name)s // keep the same name rather than HardcodedP
 
 //==========================================================================
 
+namespace Parameters_%(model_name)s_dependentCouplings
+{
+  constexpr size_t ndcoup = %(ndcoup)d; // #couplings that vary event by event because they depend on the running alphas QCD
+%(idcoup)s
+  struct DependentCouplings_sv
+  {
+%(dcoupdecl)s
+  };
+#pragma GCC diagnostic push
+#pragma GCC diagnostic ignored "-Wunused-variable"  // e.g. <<warning: unused variable ‘mdl_G__exp__2’ [-Wunused-variable]>>
+#pragma GCC diagnostic ignored "-Wunused-parameter" // e.g. <<warning: unused parameter ‘G’ [-Wunused-parameter]>>
+#ifdef __CUDACC__
+#pragma nv_diagnostic push
+#pragma nv_diag_suppress 177 // e.g. <<warning #177-D: variable "mdl_G__exp__2" was declared but never referenced>>
+#endif
+  __host__ __device__ inline const DependentCouplings_sv computeDependentCouplings_fromG( const fptype_sv& G_sv )
+  {
+#ifdef MGONGPU_HARDCODE_PARAM
+    using namespace Parameters_%(model_name)s;
+#endif
+    // NB: hardcode cxtype cI(0,1) instead of cxtype (or hardcoded cxsmpl) mdl_complexi (which exists in Parameters_%(model_name)s) because:
+    // (1) mdl_complexi is always (0,1); (2) mdl_complexi is undefined in device code; (3) need cxsmpl conversion to cxtype in code below
+    const cxtype cI( 0., 1. );
+    DependentCouplings_sv out;
+%(eftspecial1)s
+    {
+      const fptype_sv& G = G_sv;
+      // Model parameters dependent on aS
+%(dcoupsetdpar)s
+      // Model couplings dependent on aS
+%(dcoupsetdcoup)s
+    }
+%(eftspecial2)s
+    return out;
+  }
+#ifdef __CUDACC__
+#pragma GCC diagnostic pop
+#pragma nv_diagnostic pop
+#endif
+}
+
+//==========================================================================
+
+namespace Parameters_%(model_name)s_independentCouplings
+{
+  constexpr size_t nicoup = %(nicoup)d; // #couplings that are fixed for all events because they do not depend on the running alphas QCD
+%(iicoup)s
+}
+
+//==========================================================================
+
+#ifdef __CUDACC__
+namespace mg5amcGpu
+#else
+namespace mg5amcCpu
+#endif
+{
+#pragma GCC diagnostic push
+#pragma GCC diagnostic ignored "-Wunused-but-set-variable" // e.g. <<warning: variable ‘couplings_sv’ set but not used [-Wunused-but-set-variable]>>
+  // Compute the output couplings (e.g. gc10 and gc11) from the input gs
+  template<class G_ACCESS, class C_ACCESS>
+  __device__ inline void
+  G2COUP( const fptype gs[],
+          fptype couplings[] )
+  {
+    mgDebug( 0, __FUNCTION__ );
+    using namespace Parameters_%(model_name)s_dependentCouplings;
+    const fptype_sv& gs_sv = G_ACCESS::kernelAccessConst( gs );
+    DependentCouplings_sv couplings_sv = computeDependentCouplings_fromG( gs_sv );
+%(dcoupaccessbuffer)s%(dcoupkernelaccess)s%(dcoupcompute)s
+    mgDebug( 1, __FUNCTION__ );
+    return;
+  }
+#pragma GCC diagnostic pop
+}
+
+//==========================================================================
+
 #endif // Parameters_%(model_name)s_H

epochX/cudacpp/CODEGEN/PLUGIN/CUDACPP_SA_OUTPUT/madgraph/iolibs/template_files/gpu/mgOnGpuCxtypes.h

@@ -138,6 +138,12 @@ operator*( const FP& a, const cxsmpl<FP>& b )
   return cxsmpl<FP>( a, 0 ) * b;
 }
 
+inline __host__ __device__ constexpr cxsmpl<float>
+operator*( const double& a, const cxsmpl<float>& b )
+{
+  return cxsmpl<float>( a, 0 ) * b;
+}
+
 template<typename FP>
 inline __host__ __device__ constexpr cxsmpl<FP>
 operator/( const cxsmpl<FP>& a, const cxsmpl<FP>& b )