CUDA CI test

.github/workflows/c-cpp.yml

@@ -67,6 +67,7 @@ jobs:
     env:
       CUDA_HOME: /usr/local/cuda/
       REQUIRE_CUDA: 1
+      FC: gfortran
     strategy:
       matrix:
         folder: [ epochX/cudacpp/ee_mumu.sa/SubProcesses/P1_Sigma_sm_epem_mupmum,
@@ -79,11 +80,8 @@ jobs:
     steps:
     - uses: actions/checkout@v2
     - name: make info
-      run: FC=`which gfortran`;
-           make FPTYPE=${{ matrix.precision }} -C ${{ matrix.folder }} info
+      run: make FPTYPE=${{ matrix.precision }} -C ${{ matrix.folder }} info
     - name: make
-      run: FC=`which gfortran`;
-           make FPTYPE=${{ matrix.precision }} -C ${{ matrix.folder }}
+      run: make FPTYPE=${{ matrix.precision }} -C ${{ matrix.folder }}
     - name: make check
-      run: FC=`which gfortran`;
-           make FPTYPE=${{ matrix.precision }} -C ${{ matrix.folder }} check
+      run: make FPTYPE=${{ matrix.precision }} -C ${{ matrix.folder }} check

epochX/cudacpp/CODEGEN/PLUGIN/CUDACPP_SA_OUTPUT/aloha/template_files/gpu/helas.cu

@@ -1,3 +1,4 @@
+! Copyright (C) 2010 The ALOHA Development team and Contributors.
 ! Copyright (C) 2010 The MadGraph5_aMC@NLO development team and contributors.
 ! Created by: J. Alwall (Sep 2010) for the MG5aMC CPP backend.
 !==========================================================================

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

@@ -15,6 +15,12 @@
 #include <iomanip>
 #include <iostream>
 
+#ifdef MGONGPUCPP_GPUIMPL
+using namespace mg5amcGpu;
+#else
+using namespace mg5amcCpu;
+#endif
+
 #ifndef MGONGPU_HARDCODE_PARAM
 
 // Initialize static instance

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

@@ -24,175 +24,194 @@
 
 #include "read_slha.h"
 
-class Parameters_%(model_name)s
+// NB: namespaces mg5amcGpu and mg5amcCpu includes types which are defined in different ways for CPU and GPU builds (see #318 and #725)
+#ifdef MGONGPUCPP_GPUIMPL
+namespace mg5amcGpu
+#else
+namespace mg5amcCpu
+#endif
 {
-public:
+  class Parameters_%(model_name)s
+  {
+  public:
+
+    static Parameters_%(model_name)s* getInstance();
 
-  static Parameters_%(model_name)s* getInstance();
+    // Define "zero"
+    double zero, ZERO;
 
-  // Define "zero"
-  double zero, ZERO;
+    %(independent_parameters)s
 
-  %(independent_parameters)s
+    %(independent_couplings)s
 
-  %(independent_couplings)s
+    %(dependent_parameters)s
 
-  %(dependent_parameters)s
+    %(dependent_couplings)s
 
-  %(dependent_couplings)s
+    // Set parameters that are unchanged during the run
+    void setIndependentParameters( SLHAReader& slha );
 
-  // Set parameters that are unchanged during the run
-  void setIndependentParameters( SLHAReader& slha );
+    // Set couplings that are unchanged during the run
+    void setIndependentCouplings();
 
-  // Set couplings that are unchanged during the run
-  void setIndependentCouplings();
+    // Set parameters that are changed event by event
+    //void setDependentParameters(); // now computed event-by-event (running alphas #373)
 
-  // Set parameters that are changed event by event
-  //void setDependentParameters(); // now computed event-by-event (running alphas #373)
+    // Set couplings that are changed event by event
+    //void setDependentCouplings(); // now computed event-by-event (running alphas #373)
 
-  // Set couplings that are changed event by event
-  //void setDependentCouplings(); // now computed event-by-event (running alphas #373)
+    // Print parameters that are unchanged during the run
+    void printIndependentParameters();
 
-  // Print parameters that are unchanged during the run
-  void printIndependentParameters();
+    // Print couplings that are unchanged during the run
+    void printIndependentCouplings();
 
-  // Print couplings that are unchanged during the run
-  void printIndependentCouplings();
+    // Print parameters that are changed event by event
+    //void printDependentParameters(); // now computed event-by-event (running alphas #373)
 
-  // Print parameters that are changed event by event
-  //void printDependentParameters(); // now computed event-by-event (running alphas #373)
+    // Print couplings that are changed event by event
+    //void printDependentCouplings(); // now computed event-by-event (running alphas #373)
 
-  // Print couplings that are changed event by event
-  //void printDependentCouplings(); // now computed event-by-event (running alphas #373)
+  private:
 
-private:
+    static Parameters_%(model_name)s* instance;
+  };
 
-  static Parameters_%(model_name)s* instance;
-};
+} // end namespace mg5amcGpu/mg5amcCpu
 
 #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
+// NB: namespaces mg5amcGpu and mg5amcCpu includes types which are defined in different ways for CPU and GPU builds (see #318 and #725)
+#ifdef MGONGPUCPP_GPUIMPL
+namespace mg5amcGpu
+#else
+namespace mg5amcCpu
+#endif
 {
-  // Constexpr implementation of sqrt (see https://stackoverflow.com/a/34134071)
-  double constexpr sqrtNewtonRaphson( double x, double curr, double prev )
-  {
-    return curr == prev ? curr : sqrtNewtonRaphson( x, 0.5 * ( curr + x / curr ), curr );
-  }
-  double constexpr constexpr_sqrt( double x )
+  // Hardcoded constexpr physics parameters
+  namespace Parameters_%(model_name)s // keep the same name rather than HardcodedParameters_%(model_name)s for simplicity
   {
-    return x >= 0 // && x < std::numeric_limits<double>::infinity() // avoid -Wtautological-constant-compare warning in fast math
-      ? sqrtNewtonRaphson( x, x, 0 )
-      : std::numeric_limits<double>::quiet_NaN();
-  }
+    // Constexpr implementation of sqrt (see https://stackoverflow.com/a/34134071)
+    double constexpr sqrtNewtonRaphson( double x, double curr, double prev )
+    {
+      return curr == prev ? curr : sqrtNewtonRaphson( x, 0.5 * ( curr + x / curr ), curr );
+    }
+    double constexpr constexpr_sqrt( double x )
+    {
+      return x >= 0 // && x < std::numeric_limits<double>::infinity() // avoid -Wtautological-constant-compare warning in fast math
+        ? 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 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 );
+    // 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;
+    %(hardcoded_independent_parameters)s
+    // Model couplings independent of aS
+    %(hardcoded_independent_couplings)s
+    // Model parameters dependent on aS
+    %(hardcoded_dependent_parameters)s
+    // Model couplings dependent on aS
+    %(hardcoded_dependent_couplings)s
+    // Print parameters that are unchanged during the run
+    void printIndependentParameters();
+
+    // Print couplings that are unchanged during the run
+    void printIndependentCouplings();
+
+    // Print parameters that are changed event by event
+    //void printDependentParameters(); // now computed event-by-event (running alphas #373)
+
+    // Print couplings that are changed event by event
+    //void printDependentCouplings(); // now computed event-by-event (running alphas #373)
   }
 
-  // Model parameters independent of aS
-  constexpr double zero = 0;
-  constexpr double ZERO = 0;
-  %(hardcoded_independent_parameters)s
-  // Model couplings independent of aS
-  %(hardcoded_independent_couplings)s
-  // Model parameters dependent on aS
-  %(hardcoded_dependent_parameters)s
-  // Model couplings dependent on aS
-  %(hardcoded_dependent_couplings)s
-  // Print parameters that are unchanged during the run
-  void printIndependentParameters();
-
-  // Print couplings that are unchanged during the run
-  void printIndependentCouplings();
-
-  // Print parameters that are changed event by event
-  //void printDependentParameters(); // now computed event-by-event (running alphas #373)
-
-  // Print couplings that are changed event by event
-  //void printDependentCouplings(); // now computed event-by-event (running alphas #373)
-}
+} // end namespace mg5amcGpu/mg5amcCpu
 
 #endif
 
 //==========================================================================
 
-namespace Parameters_%(model_name)s_dependentCouplings
+// NB: namespaces mg5amcGpu and mg5amcCpu includes types which are defined in different ways for CPU and GPU builds (see #318 and #725)
+#ifdef MGONGPUCPP_GPUIMPL
+namespace mg5amcGpu
+#else
+namespace mg5amcCpu
+#endif
 {
-  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
+  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 MGONGPUCPP_GPUIMPL
 #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 )
-  {
+    __host__ __device__ inline const DependentCouplings_sv computeDependentCouplings_fromG( const fptype_sv& G_sv )
+    {
 #ifdef MGONGPU_HARDCODE_PARAM
-    using namespace Parameters_%(model_name)s;
+      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;
+      // 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
+      {
+        const fptype_sv& G = G_sv;
+        // Model parameters dependent on aS
 %(dcoupsetdpar)s
-      // Model couplings dependent on aS
+        // Model couplings dependent on aS
 %(dcoupsetdcoup)s
-    }
+      }
 %(eftspecial2)s
-    return out;
-  }
+      return out;
+    }
 #ifdef MGONGPUCPP_GPUIMPL
 #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
+  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 MGONGPUCPP_GPUIMPL
-namespace mg5amcGpu
-#else
-namespace mg5amcCpu
-#endif
-{
 #pragma GCC diagnostic push
 #ifndef __clang__
 #pragma GCC diagnostic ignored "-Wunused-but-set-variable" // e.g. <<warning: variable ‘couplings_sv’ set but not used [-Wunused-but-set-variable]>>
@@ -212,7 +231,8 @@ namespace mg5amcCpu
     return;
   }
 #pragma GCC diagnostic pop
-}
+
+} // end namespace mg5amcGpu/mg5amcCpu
 
 //==========================================================================