Merge lp:~maddevelopers/mg5amcnlo/lo_small_width into lp:~mg5core1/mg5amcnlo/2.6.4

[Olivier Mattelaer]

[lp:~maddevelopers/mg5amcnlo/lo_small_width]

311. By olivier-mattelaer on 2018-10-26

better writing of matrix.f to avoid to recompute the fake width systematically

312. By olivier-mattelaer on 2018-10-26

change small_width default value + Updates notes

313. By olivier-mattelaer on 2018-10-31

merge with 2.6.4

314. By olivier-mattelaer on 2018-10-31

fixing an issue with multiple particle with 0 width

315. By olivier-mattelaer on 2018-11-08

fixing tests

316. By olivier-mattelaer on 2018-11-08

merging with latest 2.6.4

317. By olivier-mattelaer on 2018-11-09

secure check in auto/scan case

=== modified file 'MadSpin/interface_madspin.py'
--- MadSpin/interface_madspin.py	2018-10-31 10:07:53 +0000
+++ MadSpin/interface_madspin.py	2018-11-09 08:24:41 +0000
@@ -1067,8 +1067,13 @@
                     if len(hepmc_output) == 0:
                         hepmc_output.append(lhe_parser.Particle(event=hepmc_output))
                         hepmc_output[0].color2 = 0
+                        hepmc_output[0].status = -1
+                        hepmc_output.nexternal+=1
                     decayed_particle = lhe_parser.Particle(particle, hepmc_output)
+                    decayed_particle.mother1 = hepmc_output[0]
+                    decayed_particle.mother2 = hepmc_output[0]
                     hepmc_output.append(decayed_particle)
+                    hepmc_output.nexternal+=1
                     decayed_particle.add_decay(decay)
             # change the weight associate to the event
             if self.options['new_wgt'] == 'cross-section':
@@ -1086,7 +1091,6 @@
             else:
                 hepmc_output.wgt = event.wgt
                 hepmc_output.nexternal = len(hepmc_output) # the append does not update nexternal
-                hepmc_output.assign_mother()
                 output_lhe.write(str(hepmc_output))
         else:
             if counter==0:
@@ -1436,7 +1440,7 @@
             if self.options['fixed_order']:
                 production, counterevt= production[0], production[1:]
             if curr_event and curr_event % 1000 == 0 and float(str(curr_event)[1:]) ==0:
-                print "decaying event number %s. Efficiency: %s [%s s]" % (curr_event, 1/self.efficiency, time.time()-start)
+                logger.info("decaying event number %s. Efficiency: %s [%s s]" % (curr_event, 1/self.efficiency, time.time()-start))
             while 1:
                 nb_try +=1
                 decays = self.get_decay_from_file(production, evt_decayfile, nb_event-curr_event)
 
=== modified file 'Template/LO/Source/run.inc'
--- Template/LO/Source/run.inc	2018-02-01 22:30:56 +0000
+++ Template/LO/Source/run.inc	2018-11-09 08:24:41 +0000
@@ -91,4 +91,7 @@
       double precision mxxmin4pdg(0:25) 
       logical mxxpart_antipart(1:25)
       common/TO_PDG_SPECIFIC_CUT/pdg_cut, ptmin4pdg,ptmax4pdg, Emin4pdg, Emax4pdg, etamin4pdg,
-     &etamax4pdg, mxxmin4pdg,mxxpart_antipart
 \ No newline at end of file
+     &etamax4pdg, mxxmin4pdg,mxxpart_antipart
+
+       double precision small_width_treatment
+       common/narrow_width/small_width_treatment
 \ No newline at end of file
 
=== modified file 'Template/LO/Source/transpole.f'
--- Template/LO/Source/transpole.f	2016-06-15 22:46:42 +0000
+++ Template/LO/Source/transpole.f	2018-11-09 08:24:41 +0000
@@ -28,13 +28,24 @@
       double precision z,zmin,zmax,xmin,xmax,ez
       double precision pole1,width1,x,xc
       double precision a,b
+c
+c     small width treatment
+c
+      double precision small_width_treatment
+      common/narrow_width/small_width_treatment
 c-----
 c  Begin Code
 c-----
       pole=pole1
       width=width1
+
       x = x1
       if (pole .gt. 0d0) then
+         if (width.lt.pole*small_width_treatment)then
+            width = pole * small_width_treatment
+            jac = jac * width/width1
+         endif
+
          zmin = atan((-pole)/width)/width
          zmax = atan((1d0-pole)/width)/width
          if (x .gt. del .and. x .lt. 1d0-del) then
@@ -175,7 +186,11 @@
 c
       double precision pole1,width1,y1,jac
       real*8 x
-
+c
+c     small width treatment
+c
+      double precision small_width_treatment
+      common/narrow_width/small_width_treatment
 c
 c     Local
 c
@@ -191,6 +206,10 @@
       width=width1
       y = y1
       if (pole .gt. 0d0) then                   !BW 
+         if (width.lt.pole*small_width_treatment)then
+            width = pole * small_width_treatment
+            jac = jac * width/width1
+         endif
          zmin = atan((-pole)/width)/width
          zmax = atan((1d0-pole)/width)/width
          z = atan((y-pole)/width)/width
 
=== modified file 'Template/LO/SubProcesses/myamp.f'
--- Template/LO/SubProcesses/myamp.f	2016-12-17 16:26:56 +0000
+++ Template/LO/SubProcesses/myamp.f	2018-11-09 08:24:41 +0000
@@ -1,91 +1,3 @@
-      double precision function testamp(p)
-c*****************************************************************************
-c     Approximates matrix element by propagators
-c*****************************************************************************
-      implicit none
-c
-c     Constants
-c     
-      include 'genps.inc'
-      include 'maxconfigs.inc'
-      include 'nexternal.inc'
-      double precision   zero
-      parameter (zero = 0d0)
-c
-c     Arguments
-c
-      double precision p(0:3,nexternal)
-c      integer iconfig
-c
-c     Local
-c
-      double precision xp(0:3,-nexternal:nexternal)
-      double precision mpole(-nexternal:0),shat,tsgn
-      integer i,j,iconfig
-
-      double precision prmass(-nexternal:0,lmaxconfigs)
-      double precision prwidth(-nexternal:0,lmaxconfigs)
-      integer pow(-nexternal:0,lmaxconfigs)
-      logical first_time
-c
-c     Global
-c
-      integer iforest(2,-max_branch:-1,lmaxconfigs)
-      common/to_forest/ iforest
-      integer            mapconfig(0:lmaxconfigs), this_config
-      common/to_mconfigs/mapconfig, this_config
-      
-      include 'coupl.inc'
-c
-c     External
-c
-      double precision dot
-
-      save prmass,prwidth,pow
-      data first_time /.true./
-c-----
-c  Begin Code
-c-----      
-      iconfig = this_config
-      if (first_time) then
-c         include 'props.inc'
-         first_time=.false.
-      endif
-
-      do i=1,nexternal
-         mpole(-i)=0d0
-         do j=0,3
-            xp(j,i)=p(j,i)
-         enddo
-      enddo
-c      mpole(-3) = 174**2
-c      shat = dot(p(0,1),p(0,2))/(1800)**2
-      shat = dot(p(0,1),p(0,2))/(10)**2
-c      shat = 1d0
-      testamp = 1d0
-      tsgn    = +1d0
-      do i=-1,-(nexternal-3),-1              !Find all the propagotors
-         if (iforest(1,i,iconfig) .eq. 1) tsgn=-1d0
-         do j=0,3
-            xp(j,i) = xp(j,iforest(1,i,iconfig))
-     $           +tsgn*xp(j,iforest(2,i,iconfig))
-         enddo
-         if (prwidth(i,iconfig) .ne. 0d0 .and. .false.) then
-            testamp=testamp/((dot(xp(0,i),xp(0,i))
-     $                        -prmass(i,iconfig)**2)**2
-     $         -(prmass(i,iconfig)*prwidth(i,iconfig))**2)
-         else
-            testamp = testamp/((dot(xp(0,i),xp(0,i)) -
-     $                          prmass(i,iconfig)**2)
-     $                          **(pow(i,iconfig)))
-         endif
-        testamp=testamp*shat**(pow(i,iconfig))
-c        write(*,*) i,iconfig,pow(i,iconfig),prmass(i,iconfig)
-      enddo
-c      testamp = 1d0/dot(xp(0,-1),xp(0,-1))
-      testamp=abs(testamp)
-c      testamp = 1d0
-      end
 
       logical function cut_bw(p)
 c*****************************************************************************
@@ -114,6 +26,7 @@
 
       double precision prmass(-nexternal:0,lmaxconfigs)
       double precision prwidth(-nexternal:0,lmaxconfigs)
+      double precision prwidth_tmp(-nexternal:0,lmaxconfigs)
       integer pow(-nexternal:0,lmaxconfigs)
       logical first_time, onshell
       double precision xmass
@@ -152,13 +65,14 @@
 c
       double precision dot
 
-      save prmass,prwidth,pow
+      save prmass,prwidth,pow,prwidth_tmp
       data first_time /.true./
 c-----
 c  Begin Code
 c-----      
       cut_bw = .false.    !Default is we passed the cut
       iconfig = this_config
+
       if (first_time) then
          include 'props.inc'
          nbw = 0
@@ -209,10 +123,12 @@
 c            write(*,*) 'xmass',xmass,prmass(i,iconfig)
 c
 c           Here we set if the BW is "on-shell" for LesHouches
-c
+c            
+            prwidth_tmp(i,iconfig) = max(prwidth(i,iconfig), prmass(i,iconfig)*small_width_treatment)
+ 
             onshell = (abs(xmass - prmass(i,iconfig)) .lt.
-     $           bwcutoff*prwidth(i,iconfig).and.
-     $           (prwidth(i,iconfig)/prmass(i,iconfig).lt.0.1d0.or.
+     $           bwcutoff*prwidth_tmp(i,iconfig).and.
+     $           (prwidth_tmp(i,iconfig)/prmass(i,iconfig).lt.0.1d0.or.
      $            gForceBW(i,iconfig).eq.1))
             if(onshell)then
 c     Remove on-shell forbidden s-channels (gForceBW=2) (JA 2/10/11)
@@ -264,10 +180,10 @@
 c     For decay-chain syntax use BWcutoff here too (22/12/14)
             if (gForceBW(i, iconfig).eq.1) then
                onshell = (abs(xmass - prmass(i,iconfig)) .lt.
-     $           bwcutoff*prwidth(i,iconfig))
+     $           bwcutoff*prwidth_tmp(i,iconfig))
             else
                onshell = (abs(xmass - prmass(i,iconfig)) .lt.
-     $           5d0*prwidth(i,iconfig))
+     $           5d0*prwidth_tmp(i,iconfig))
             endif
 
             if (onshell .and. (lbw(nbw).eq. 2) .or.
@@ -313,6 +229,7 @@
 
       double precision prmass(-nexternal:0,lmaxconfigs)
       double precision prwidth(-nexternal:0,lmaxconfigs)
+      double precision prwidth_tmp(-nexternal:0,lmaxconfigs)
       integer pow(-nexternal:0,lmaxconfigs)
 
       integer idup(nexternal,maxproc,maxsproc)
@@ -382,11 +299,21 @@
 
 c-----
 c  Begin Code
-c-----      
+c-----     
+      iconfig = this_config
+c     needs to be initialise to avoid segfault
+      do i = -nexternal,-1
+         prwidth(i,iconfig) = 0
+         prmass(i,iconfig) =0
+      enddo 
       include 'props.inc'
 c      etmin = 10
       nt = 0
-      iconfig = this_config
+      do i = -nexternal,-1
+         prwidth_tmp(i,iconfig) = max(prwidth(i,iconfig), prmass(i,iconfig)*small_width_treatment)
+      enddo
+
+
       mtot = 0d0
       etot = 0d0   !Total energy needed
       spmass = 0d0 !Keep track of BW masses for shat
@@ -423,7 +350,7 @@
 c     Look for identical particles to map radiation processes
       call idenparts(iden_part, iforest(1,-max_branch,iconfig),
      $     sprop(1,-max_branch,iconfig), gForceBW(-max_branch,iconfig),
-     $     prwidth(-nexternal,iconfig))
+     $     prwidth_tmp(-nexternal,iconfig))
 
 c     Start loop over propagators
       do i=-1,-(nexternal-3),-1
@@ -446,15 +373,15 @@
               xm(i)=max(xm(i),max(xqcutij(l1,l2),0d0))
             endif
 c            write(*,*) 'iconfig,i',iconfig,i
-c            write(*,*) prwidth(i,iconfig),prmass(i,iconfig)
-            if (prwidth(i,iconfig) .gt. 0 ) then
+c            write(*,*) prwidth_tmp(i,iconfig),prmass(i,iconfig)
+            if (prwidth_tmp(i,iconfig) .gt. 0 ) then
                nbw=nbw+1
 c              JA 6/8/2011 Set xe(i) for resonances
                if (gforcebw(i,iconfig).eq.1) then
-                  xm(i) = max(xm(i), prmass(i,iconfig)-bwcutoff*prwidth(i,iconfig))
+                  xm(i) = max(xm(i), prmass(i,iconfig)-bwcutoff*prwidth_tmp(i,iconfig))
                   bwcut_for_PS(i) = bwcutoff
                else if (lbw(nbw).eq.1) then
-                  xm(i) = max(xm(i), prmass(i,iconfig)-5d0*prwidth(i,iconfig))
+                  xm(i) = max(xm(i), prmass(i,iconfig)-5d0*prwidth_tmp(i,iconfig))
                   bwcut_for_PS(i) = 5d0
                else
                   bwcut_for_PS(i) = 5d0
@@ -465,19 +392,19 @@
 c     Either: required onshell and daughter masses too large
 c     Or: forced and daughter masses too large
 c     Or: required offshell and forced
-            if(prwidth(i,iconfig) .gt. 0.and.
+            if(prwidth_tmp(i,iconfig) .gt. 0.and.
      $         (lbw(nbw).eq.1.and.
-     $          (prmass(i,iconfig)+bwcut_for_PS(i)*prwidth(i,iconfig).lt.xm(i)
-     $           .or.prmass(i,iconfig)-bwcut_for_PS(i)*prwidth(i,iconfig).gt.dsqrt(stot))
+     $          (prmass(i,iconfig)+bwcut_for_PS(i)*prwidth_tmp(i,iconfig).lt.xm(i)
+     $           .or.prmass(i,iconfig)-bwcut_for_PS(i)*prwidth_tmp(i,iconfig).gt.dsqrt(stot))
      $          .or.gforcebw(i,iconfig).eq.1.and.
-     $              prmass(i,iconfig)+bwcutoff*prwidth(i,iconfig).lt.xm(i)
+     $              prmass(i,iconfig)+bwcutoff*prwidth_tmp(i,iconfig).lt.xm(i)
      $          .or.lbw(nbw).eq.2.and.gforcebw(i,iconfig).eq.1))
      $        then
 c     Write results.dat and quit
                call write_null_results()
                stop
             endif
-            if (prwidth(i,iconfig) .gt. 0 .and. lbw(nbw) .le. 1) then         !B.W.
+            if (prwidth_tmp(i,iconfig) .gt. 0 .and. lbw(nbw) .le. 1) then         !B.W.
                if (i .eq. -(nexternal-(nincoming+1))) then  !This is s-hat
                   j = 3*(nexternal-2)-4+1    !set i to ndim+1
 c-----
@@ -489,14 +416,14 @@
      $                 .or. lbw(nbw).eq.1) then
                      write(*,*) 'Setting PDF BW',j,nbw,prmass(i,iconfig)
                      spole(j)=prmass(i,iconfig)*prmass(i,iconfig)/stot
-                     swidth(j) = prwidth(i,iconfig)*prmass(i,iconfig)/stot
+                     swidth(j) = prwidth(i,iconfig)*prmass(i,iconfig)/stot ! keep the real width here (important for the jacobian)
                   endif
-               else if((prmass(i,iconfig)+bwcut_for_PS(i)*prwidth(i,iconfig)).ge.xm(i)
+               else if((prmass(i,iconfig)+bwcut_for_PS(i)*prwidth_tmp(i,iconfig)).ge.xm(i)
      $                  .and. iden_part(i).eq.0 .or. lbw(nbw).eq.1) then
 c              JA 02/13 Only allow BW if xm below M+5*Gamma
                   write(*,*) 'Setting BW',i,nbw,prmass(i,iconfig)
                   spole(-i)=prmass(i,iconfig)*prmass(i,iconfig)/stot
-                  swidth(-i) = prwidth(i,iconfig)*prmass(i,iconfig)/stot
+                  swidth(-i) = prwidth(i,iconfig)*prmass(i,iconfig)/stot ! keep the real width here (important for the jacobian)
                endif
 c     JA 4/1/2011 Set grid in case there is no BW (radiation process)
                if (swidth(-i) .eq. 0d0 .and.
@@ -509,7 +436,7 @@
 c     Set spmass for BWs
                if (swidth(-i) .ne. 0d0)
      $              spmass=spmass-xm(i) +
-     $              max(xm(i),prmass(i,iconfig)-bwcut_for_PS(i)*prwidth(i,iconfig))
+     $              max(xm(i),prmass(i,iconfig)-bwcut_for_PS(i)*prwidth_tmp(i,iconfig))
             else                                  !1/x^pow
               a=prmass(i,iconfig)**2/stot
 c     JA 4/1/2011 always set grid
@@ -530,7 +457,7 @@
                  endif
               endif
               if (xo.eq.0d0) xo=1d0/stot
-c              if (prwidth(i, iconfig) .eq. 0d0.or.iden_part(i).gt.0) then 
+c              if (prwidth_tmp(i, iconfig) .eq. 0d0.or.iden_part(i).gt.0) then 
               call setgrid(-i,xo,a,1)
 c              else 
 c                 write(*,*) 'Using flat grid for BW',i,nbw,
@@ -610,7 +537,7 @@
             spole(i)= -2.0d0    ! 1/s pole
             write(*,*) "Transforming s_hat 1/s ",i,xo, smin, stot
          else
-            write(*,*) "Transforming s_hat BW ",spole(i),swidth(i)
+            write(*,*) "Transforming s_hat BW ",spole(i), max(swidth(i), spole(i)*small_width_treatment)
          endif
       endif
 
 
=== modified file 'UpdateNotes.txt'
--- UpdateNotes.txt	2018-10-17 13:25:07 +0000
+++ UpdateNotes.txt	2018-11-09 08:24:41 +0000
@@ -1,12 +1,18 @@
 Update notes for MadGraph5_aMC@NLO (in reverse time order)
 
 2.6.4 (//)
+      OM: add specific treatement for small width (at LO only and not for loop-induced)
+          if the width is smaller than 1e-6 times the mass, a fake width (at that value) is used for the
+          numerical evaluation of the matrix-element. S-channel resonances are re-scaled according to 
+          narrow-width approximation to return the correct total cross-section (the distribution of events 
+          will on the other hand follow the new width). 
+          The parameter '1e-6' can be changed by adding to (LO) run_card the parameter: "small_width_treatment"
+      OM: add a new command "install looptools" to trigger the question that is automatically trigger 
+          the first time a loop computation is needed.
       RF: Fixed a bug when using TopDrawer plots for f(N)LO runs, where the combination of the plots could lead
           to completely wrong histograms/distributions in case of high-precision runs.
       OM: Fix some MLM crash for some processes (in particular BSM processes with W'). 
       OM: Fix a bug in the reweighting due to the new lhe format (the one avoiding some issue with py8)
-      OM: add a new command "install looptools" to trigger the question that is automatically trigger 
-          the first time a loop computation is needed.
       OM: Fix a behavior for negative mass, the width was set to negative in the param_card automatically
           making the Parton-shower (and other code) to crash since this does not follow the convention.
       OM: Change compiler flag to support Mojave.
 
=== modified file 'madgraph/core/helas_objects.py'
--- madgraph/core/helas_objects.py	2017-08-10 13:09:28 +0000
+++ madgraph/core/helas_objects.py	2018-11-09 08:24:41 +0000
@@ -4569,6 +4569,13 @@
         return sum([d.get('wavefunctions') for d in \
                        self.get('diagrams')], [])
 
+
+    def get_all_mass_widths(self):
+        """Gives a list of all widths used by this ME (from propagator)"""
+
+        return set([(d.get('mass'),d.get('width')) for d in self.get_all_wavefunctions()])
+
+
     def get_all_amplitudes(self):
         """Gives a list of all amplitudes for this ME"""
 
 
=== modified file 'madgraph/interface/common_run_interface.py'
--- madgraph/interface/common_run_interface.py	2018-11-08 22:22:58 +0000
+++ madgraph/interface/common_run_interface.py	2018-11-09 08:24:41 +0000
@@ -3242,7 +3242,19 @@
         - Check that no width are too small (raise a warning if this is the case)
         3) if dependent is on True check for dependent parameter (automatic for scan)"""
         
-        return self.static_check_param_card(path, self, run=run, dependent=dependent)
+        self.static_check_param_card(path, self, run=run, dependent=dependent)
+        
+        card = param_card_mod.ParamCard(path)
+        for param in card['decay']:
+            width = param.value
+            if width == 0:
+                continue
+            try:
+                mass = card['mass'].get(param.lhacode).value
+            except Exception:
+                continue
+        
+        
         
     @staticmethod
     def static_check_param_card(path, interface, run=True, dependent=False, 
@@ -5680,6 +5692,30 @@
                  self.run_card['mass_ion1'] != self.run_card['mass_ion2']):
                 raise Exception, "Heavy ion profile for both beam are different but the symmetry used forbids it. \n Please generate your process with \"set group_subprocesses False\"."
             
+            # check the status of small width status from LO
+            for param in self.param_card['decay']:
+                width = param.value
+                if width == 0 or isinstance(width,str):
+                    continue
+                try:
+                    mass = self.param_card['mass'].get(param.lhacode).value
+                except Exception:
+                    continue
+                if isinstance(mass,str):
+                    continue
+                
+                if mass:
+                    if abs(width/mass) < self.run_card['small_width_treatment']:
+                        logger.warning("Particle %s will use a fake width  ( %s instead of %s ).\n" +
+                          "Cross-section will be rescaled according to NWA if needed."  +
+                          "To force exact treatment reduce the value of 'small_width_treatment' parameter of the run_card",
+                          param.lhacode[0], mass*self.run_card['small_width_treatment'], width)
+                    elif abs(width/mass) < 1e-12:
+                        logger.error('The width of particle %s is too small for an s-channel resonance (%s). If you have this particle in an s-channel, this is likely to create numerical instabilities .', param.lhacode[0], width)
+                    if CommonRunCmd.sleep_for_error:
+                        time.sleep(5)
+                        CommonRunCmd.sleep_for_error = False
+
 
         ########################################################################
         #       NLO specific check
 
=== modified file 'madgraph/iolibs/export_v4.py'
--- madgraph/iolibs/export_v4.py	2018-11-08 22:22:58 +0000
+++ madgraph/iolibs/export_v4.py	2018-11-09 08:24:41 +0000
@@ -12,6 +12,7 @@
 # For more information, visit madgraph.phys.ucl.ac.be and amcatnlo.web.cern.ch
 #
 ################################################################################
+from madgraph.iolibs.helas_call_writers import HelasCallWriter
 """Methods and classes to export matrix elements to v4 format."""
 
 import copy
@@ -93,12 +94,28 @@
     #    - None, madgraph do nothing for initialisation
     exporter = 'v4'
     # language of the output 'v4' for Fortran output
-    #                        'cpp' for C++ output 
+    #                        'cpp' for C++ output
     
     
     def __init__(self, dir_path = "", opt=None):
-        # cmd_options is a dictionary with all the optional argurment passed at output time 
-        return
+        # cmd_options is a dictionary with all the optional argurment passed at output time
+        
+        # Activate some monkey patching for the helas call writer.
+        helas_call_writers.HelasCallWriter.customize_argument_for_all_other_helas_object = \
+                self.helas_call_writer_custom
+        
+
+    # helper function for customise helas writter
+    @staticmethod
+    def custom_helas_call(call, arg):
+        """static method to customise the way aloha function call are written
+        call is the default template for the call
+        arg are the dictionary used for the call
+        """
+        return call, arg
+    
+    helas_call_writer_custom = lambda x,y,z: x.custom_helas_call(y,z)
+
 
     def copy_template(self, model):
         return
@@ -160,6 +177,9 @@
         #place holder to pass information to the run_interface
         self.proc_characteristic = banner_mod.ProcCharacteristic()
         
+        # call mother class
+        super(ProcessExporterFortran,self).__init__(dir_path, opt)
+        
         
     #===========================================================================
     # process exporter fortran switch between group and not grouped
@@ -3423,6 +3443,7 @@
         return s_and_t_channels
 
 
+
 #===============================================================================
 # ProcessExporterFortranME
 #===============================================================================
@@ -3431,7 +3452,16 @@
     MadEvent format."""
 
     matrix_file = "matrix_madevent_v4.inc"
-
+    
+    # helper function for customise helas writter
+    @staticmethod
+    def custom_helas_call(call, arg):
+        if arg['mass'] == '%(M)s,%(W)s,':
+            arg['mass'] = '%(M)s, fk_%(W)s,'
+        elif '%(W)s' in arg['mass']:
+            raise Exception
+        return call, arg
+    
     def copy_template(self, model):
         """Additional actions needed for setup of Template
         """
@@ -3455,6 +3485,7 @@
         
         
 
+    
 
 
     #===========================================================================
@@ -3945,10 +3976,29 @@
         # Extract helas calls
         helas_calls = fortran_model.get_matrix_element_calls(\
                     matrix_element)
+        
 
         replace_dict['helas_calls'] = "\n".join(helas_calls)
 
 
+        #adding the support for the fake width (forbidding too small width)
+        mass_width = matrix_element.get_all_mass_widths()
+        width_list = set([e[1] for e in mass_width])
+        
+        replace_dict['fake_width_declaration'] = \
+            ('  double precision fk_%s \n' * len(width_list)) % tuple(width_list)
+        replace_dict['fake_width_declaration'] += \
+            ('  save fk_%s \n' * len(width_list)) % tuple(width_list)
+        fk_w_defs = []
+        one_def = ' fk_%(w)s = SIGN(MAX(ABS(%(w)s), ABS(%(m)s*small_width_treatment)), %(w)s)'     
+        for m, w in mass_width:
+            if w == 'zero':
+                if ' fk_zero = 0d0' not in fk_w_defs: 
+                    fk_w_defs.append(' fk_zero = 0d0')
+                continue    
+            fk_w_defs.append(one_def %{'m':m, 'w':w})
+        replace_dict['fake_width_definitions'] = '\n'.join(fk_w_defs)
+
         # Extract version number and date from VERSION file
         info_lines = self.get_mg5_info_lines()
         replace_dict['info_lines'] = info_lines
 
=== modified file 'madgraph/iolibs/helas_call_writers.py'
--- madgraph/iolibs/helas_call_writers.py	2016-03-03 23:17:06 +0000
+++ madgraph/iolibs/helas_call_writers.py	2018-11-09 08:24:41 +0000
@@ -21,11 +21,13 @@
 import aloha.aloha_writers as aloha_writers
 import aloha
 from madgraph import MadGraph5Error
+import madgraph.various.misc as misc
 
 class HelasWriterError(Exception):
     """Class for the error of this module """
     pass
 
+
 #===============================================================================
 # HelasCallWriter
 #===============================================================================
@@ -41,6 +43,19 @@
     # Dictionaries from spin states to letters in Helas call
     mother_dict = {1: 'S', 2: 'O', -2: 'I', 3: 'V', 5: 'T', 4:'OR', -4:'IR',
                    99:'P'}
+    
+    @staticmethod
+    def customize_argument_for_all_other_helas_object(call, arg):
+        """ Place holder for PLUGIN/...
+            (used by madevent output for small width handling) 
+        """
+        return call, arg
+        
+    @staticmethod
+    def default_customize_argument_for_all_other_helas_object(call,arg):
+        return call,arg  
+    #customize_argument_for_all_other_helas_object = fct_customize_argument_for_all_other_helas_object
+    #default_customize_argument_for_all_other_helas_object = fct_customize_argument_for_all_other_helas_object
 
     def default_setup(self):
 
@@ -85,6 +100,11 @@
 
         return ['model', 'wavefunctions', 'amplitudes']
 
+
+
+    
+     
+
     def get_loop_amp_helas_calls(self, matrix_element):
         """Return a list of strings, corresponding to the Helas calls
         for building loop amplitudes (AMPL) only."""
@@ -1197,6 +1217,7 @@
                  arg['second_line'] = ampl+"="+ampl+"*(%(uvct)s)"           
 
         # ALL ARGUMENT FORMATTED ###############################################
+        call, arg = HelasCallWriter.customize_argument_for_all_other_helas_object(call, arg)
         # Store the result.
         call = call % arg
         # Now we have a line correctly formatted
@@ -1209,6 +1230,8 @@
         else:
             self.add_amplitude(argument.get_call_key(), call_function)
 
+
+
     def get_loop_amplitude_helas_calls(self, loop_matrix_element):
         """ Returns a list of strings corresponding to the Helas calls for each 
         loop amplitude of this loop matrix element. This function is placed in
 
=== modified file 'madgraph/iolibs/template_files/matrix_loop_induced_madevent_group.inc'
--- madgraph/iolibs/template_files/matrix_loop_induced_madevent_group.inc	2017-11-01 10:08:52 +0000
+++ madgraph/iolibs/template_files/matrix_loop_induced_madevent_group.inc	2018-11-09 08:24:41 +0000
@@ -39,6 +39,7 @@
 C   When set negative, the security above is removed
 	DOUBLE PRECISION MULTICHANNEL_THRES
 	PARAMETER (MULTICHANNEL_THRES=1.0d-5)
+
 c
 c global (due to reading writting) 
 c
 
=== modified file 'madgraph/iolibs/template_files/matrix_madevent_group_v4.inc'
--- madgraph/iolibs/template_files/matrix_madevent_group_v4.inc	2017-11-01 10:08:52 +0000
+++ madgraph/iolibs/template_files/matrix_madevent_group_v4.inc	2018-11-09 08:24:41 +0000
@@ -68,6 +68,9 @@
  
     REAL*8 POL(2)
     COMMON/TO_POLARIZATION/ POL
+
+	double precision small_width_treatment
+	common/narrow_width/small_width_treatment
     
     INTEGER          ISUM_HEL
     LOGICAL                    MULTI_CHANNEL
@@ -87,6 +90,7 @@
 C ----------
 C BEGIN CODE
 C ----------
+
     NTRY(IMIRROR)=NTRY(IMIRROR)+1
     THIS_NTRY(IMIRROR) = THIS_NTRY(IMIRROR)+1 
     DO I=1,NEXTERNAL
@@ -255,6 +259,11 @@
 C   Needed for v4 models
     COMPLEX*16 DUM0,DUM1
     DATA DUM0, DUM1/(0d0, 0d0), (1d0, 0d0)/
+
+	%(fake_width_declaration)s
+	logical first
+	data first /.true./
+	save first
 C
 C FUNCTION
 C
@@ -265,6 +274,9 @@
     Double Precision amp2(maxamps), jamp2(0:maxflow)
     common/to_amps/  amp2,       jamp2
     include 'coupl.inc'
+
+	double precision small_width_treatment
+	common/narrow_width/small_width_treatment
 C  
 C COLOR DATA
 C  
@@ -272,6 +284,12 @@
 C ----------
 C BEGIN CODE
 C ----------
+if (first) then
+	first=.false.
+	%(fake_width_definitions)s
+endif
+
+
 %(helas_calls)s
 %(jamp_lines)s
 
 
=== modified file 'madgraph/iolibs/template_files/matrix_madevent_v4.inc'
--- madgraph/iolibs/template_files/matrix_madevent_v4.inc	2017-11-01 10:08:52 +0000
+++ madgraph/iolibs/template_files/matrix_madevent_v4.inc	2018-11-09 08:24:41 +0000
@@ -49,6 +49,7 @@
     INTEGER IDUM, NGOOD, IGOOD(NCOMB), JHEL, J, JJ
     REAL     XRAN1
     EXTERNAL XRAN1
+
 C  
 C GLOBAL VARIABLES
 C  
@@ -58,7 +59,7 @@
     CHARACTER*101        HEL_BUFF
     COMMON/TO_HELICITY/  HEL_BUFF
     
-    REAL*8 POL(2)
+	REAL*8 POL(2)
     COMMON/TO_POLARIZATION/ POL
     
     INTEGER          ISUM_HEL
@@ -68,6 +69,7 @@
     DATA IDUM /-1/
     DATA XTRY, XREJ, NGOOD /0,0,0/
     SAVE YFRAC, IGOOD, JHEL
+
 %(helicity_lines)s
 %(den_factor_line)s
 C ----------
@@ -235,6 +237,11 @@
 C   Needed for v4 models
     COMPLEX*16 DUM0,DUM1
     DATA DUM0, DUM1/(0d0, 0d0), (1d0, 0d0)/
+
+	%(fake_width_declaration)s
+	logical first
+	data first /.true./
+	save first
 C
 C FUNCTION
 C
@@ -245,6 +252,9 @@
     Double Precision amp2(maxamps), jamp2(0:maxflow)
     common/to_amps/  amp2,       jamp2
     include 'coupl.inc'
+
+	double precision small_width_treatment
+	common/narrow_width/small_width_treatment
 C  
 C COLOR DATA
 C  
@@ -252,6 +262,11 @@
 C ----------
 C BEGIN CODE
 C ----------
+if (first) then
+	first=.false.
+	%(fake_width_definitions)s
+endif
+
 %(helas_calls)s
 %(jamp_lines)s
 
 
=== modified file 'madgraph/iolibs/template_files/matrix_madweight_group_v4.inc'
--- madgraph/iolibs/template_files/matrix_madweight_group_v4.inc	2012-02-29 12:53:46 +0000
+++ madgraph/iolibs/template_files/matrix_madweight_group_v4.inc	2018-11-09 08:24:41 +0000
@@ -22,6 +22,7 @@
 C ARGUMENTS 
 C  
       REAL*8 P(0:3,NEXTERNAL),ANS
+
 C  
 C LOCAL VARIABLES 
 C  
 
=== modified file 'madgraph/various/banner.py'
--- madgraph/various/banner.py	2018-10-19 08:43:52 +0000
+++ madgraph/various/banner.py	2018-11-09 08:24:41 +0000
@@ -2952,6 +2952,7 @@
         self.add_param('survey_splitting', -1, hidden=True, include=False, comment="for loop-induced control how many core are used at survey for the computation of a single iteration.")
         self.add_param('survey_nchannel_per_job', 2, hidden=True, include=False, comment="control how many Channel are integrated inside a single job on cluster/multicore")
         self.add_param('refine_evt_by_job', -1, hidden=True, include=False, comment="control the maximal number of events for the first iteration of the refine (larger means less jobs)")
+        self.add_param('small_width_treatment', 1e-6, hidden=True, comment="generation where the width is below VALUE times mass will be replace by VALUE times mass for the computation. The cross-section will be corrected assuming NWA. Not used for loop-induced process")
         
         # parameter allowing to define simple cut via the pdg
         # Special syntax are related to those. (can not be edit directly)
 
=== modified file 'madgraph/various/lhe_parser.py'
--- madgraph/various/lhe_parser.py	2018-06-21 08:57:21 +0000
+++ madgraph/various/lhe_parser.py	2018-11-09 08:24:41 +0000
@@ -167,6 +167,7 @@
         """associate to this particle the decay in the associate event"""
         
         return self.event.add_decay_to_particle(self.event_id, decay_event)
+
             
     def __repr__(self):
         return 'Particle("%s", event=%s)' % (str(self), self.event)
@@ -1263,26 +1264,28 @@
                 
             if 'part' == status:
                 part = Particle(line, event=self)
-                if part.E != 0:
+                if part.E != 0 or part.status==-1:
                     self.append(part)
                 elif self.nexternal:
-                        self.nexternal-=1
+                    self.nexternal-=1
             else:
                 if '</event>' in line:
                     line = line.replace('</event>','',1)
                 self.tag += '%s\n' % line
                 
         self.assign_mother()
-        
+    
+    
     def assign_mother(self):
         """convert the number in actual particle"""
-        #Security if not incoming particle. Define a fake particle and set all particle as 
-        # decaying from that fake particle
+        #Security if not incoming particle. Define a fake particle 
         if all(p.status != -1 for p in self):
+            if not self.nexternal:
+                return
             if self.warning_order:
                 Event.warning_order = False
                 logger.warning("Weird format for lhe format: no incoming particle... adding a fake one")
-            
+                raise Exception
             mother = Particle(event=self)
             mother.status = -1
             mother.pid = 0
@@ -1291,8 +1294,10 @@
             mother.event_id = 0
             self.nexternal += 1
             for p in self[1:]:
-                p.mother1 = 1
-                p.mother2 = 1
+                if isinstance(p.mother1, int) and p.mother1 > 1:
+                    p.mother1 += 1
+                if isinstance(p.mother2, int) and p.mother2 > 1:
+                    p.mother2 += 1
                 p.event_id += 1
             
         
 
=== modified file 'tests/acceptance_tests/test_cmd_madloop.py'
--- tests/acceptance_tests/test_cmd_madloop.py	2018-04-24 22:39:50 +0000
+++ tests/acceptance_tests/test_cmd_madloop.py	2018-11-09 08:24:41 +0000
@@ -492,7 +492,7 @@
         # Select the Tensor Integral to include in the test
         misc.deactivate_dependence('pjfry', cmd = self.interface, log='stdout')
         misc.deactivate_dependence('samurai', cmd = self.interface, log='stdout')        
-        misc.activate_dependence('golem', cmd = self.interface, log='stdout')
+        misc.deactivate_dependence('golem', cmd = self.interface, log='stdout')
         misc.activate_dependence('ninja', cmd = self.interface, log='stdout',MG5dir=MG5DIR)
         misc.activate_dependence('collier', cmd = self.interface, log='stdout',MG5dir=MG5DIR)
 
@@ -533,7 +533,7 @@
         # Select the Tensor Integral to include in the test
         misc.deactivate_dependence('pjfry', cmd = interface, log='stdout')
         misc.deactivate_dependence('samurai', cmd = interface, log='stdout')        
-        misc.activate_dependence('golem', cmd = interface, log='stdout')
+        misc.deactivate_dependence('golem', cmd = interface, log='stdout')
         misc.activate_dependence('ninja', cmd = interface, log='stdout',MG5dir=MG5DIR)
 
         run_cmd('generate g g > t t~ [virt=QCD]')
 
=== modified file 'tests/acceptance_tests/test_export_fks.py'
--- tests/acceptance_tests/test_export_fks.py	2017-05-15 10:39:31 +0000
+++ tests/acceptance_tests/test_export_fks.py	2018-11-09 08:24:41 +0000
@@ -75,7 +75,7 @@
         # Select the Tensor Integral to include in the test
         misc.deactivate_dependence('pjfry', cmd = interface, log='stdout')
         misc.deactivate_dependence('samurai', cmd = interface, log='stdout')        
-        misc.activate_dependence('golem', cmd = interface, log='stdout')
+        misc.deactivate_dependence('golem', cmd = interface, log='stdout')
         misc.activate_dependence('ninja', cmd = interface, log='stdout',MG5dir=MG5DIR)        
         
         run_cmd('import model %s' % model)
 
=== removed file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%GOLEM_interface.f'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%GOLEM_interface.f	2016-08-07 07:16:07 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%GOLEM_interface.f	1970-01-01 00:00:00 +0000
@@ -1,748 +0,0 @@
-      SUBROUTINE GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-C     
-C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
-C     By the MadGraph5_aMC@NLO Development Team
-C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
-C     
-C     Interface between MG5 and Golem95.
-C     The Golem95 version should be higher than 1.3.0.
-C     It supports RANK = NLOOPLINE + 1 tensor integrals when 1 <
-C      NLOOPLINE < 6.
-C     
-C     Process: d~ u > w+ QED<=1 WEIGHTED<=2 [ all = QCD ]
-C     Process: s~ c > w+ QED<=1 WEIGHTED<=2 [ all = QCD ]
-C     
-C     
-C     MODULES
-C     
-      USE MATRICE_S
-      USE FORM_FACTOR_TYPE, ONLY: FORM_FACTOR
-      USE PRECISION_GOLEM, ONLY: KI
-      USE TENS_COMB
-      USE TENS_REC
-      USE FORM_FACTOR_1P, ONLY: A10
-      USE FORM_FACTOR_2P, ONLY: A20
-      USE FORM_FACTOR_3P, ONLY: A30
-      USE FORM_FACTOR_4P, ONLY: A40
-      USE FORM_FACTOR_5P, ONLY: A50
-      USE FORM_FACTOR_6P, ONLY: A60
-      IMPLICIT NONE
-C     
-C     CONSTANTS 
-C     
-      INTEGER    NEXTERNAL
-      PARAMETER (NEXTERNAL=3)
-      LOGICAL CHECKPCONSERVATION
-      PARAMETER (CHECKPCONSERVATION=.TRUE.)
-      REAL*8 NORMALIZATION
-      PARAMETER (NORMALIZATION = 1.D0/(16.D0*3.14159265358979323846D0*
-     $ *2))
-      REAL(KI),DIMENSION(0:3),PARAMETER::NULL_VEC = (/0.0_KI,0.0_KI
-     $ ,0.0_KI,0.0_KI/)
-C     GOLEM_RUN_MODE = 1: Use directly MadLoop tensorial coefficients
-C     GOLEM_RUN_MODE = 2: Reconstruct the tensorial coefficeints
-C      directly from 
-C     numerator using golem internal reconstruction routine
-C     GOLEM_RUN_MODE = 3: Cross-checked reconstructed coefficients
-C      against
-C     MadLoop internal ones.
-      INTEGER GOLEM_RUN_MODE
-      PARAMETER (GOLEM_RUN_MODE=1)
-C     The following is the acceptance threshold used for
-C      GOLEM_RUN_MODE = 3
-      REAL*8 COEF_CHECK_THRS
-      DATA COEF_CHECK_THRS/1.0D-13/
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-
-      LOGICAL PASS_COEF_CHECK
-C     
-C     ARGUMENTS 
-C     
-      INTEGER NLOOPLINE, RANK
-      REAL*8 PL(0:3,NLOOPLINE)
-      REAL*8 PCT(0:3,0:NLOOPLINE-1), ABSPCT(0:3)
-      REAL*8 REF_P
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX*16 M2L(NLOOPLINE)
-      COMPLEX(KI) M2LGOLEM(NLOOPLINE)
-      COMPLEX*16 RES(3)
-      LOGICAL STABLE
-C     
-C     LOCAL VARIABLES 
-C     
-      INTEGER I, J, K
-      TYPE(FORM_FACTOR)::RES_GOLEM
-
-      COMPLEX(KI)::COEFFS0,COEFFS0_REC
-      TYPE(COEFF_TYPE_1)::COEFFS1,COEFFS1_REC
-      TYPE(COEFF_TYPE_2)::COEFFS2,COEFFS2_REC
-      TYPE(COEFF_TYPE_3)::COEFFS3,COEFFS3_REC
-      TYPE(COEFF_TYPE_4)::COEFFS4,COEFFS4_REC
-      TYPE(COEFF_TYPE_5)::COEFFS5,COEFFS5_REC
-      TYPE(COEFF_TYPE_6)::COEFFS6,COEFFS6_REC
-
-C     The pinch propagator optimization is not used, so for now it is
-C     always 0.
-      INTEGER PINCH
-C     
-C     EXTERNAL FUNCTIONS
-C     
-      COMPLEX(KI) GOLEM_LOOPNUM
-      EXTERNAL GOLEM_LOOPNUM
-      LOGICAL COMPARE_COEFS_0
-      LOGICAL COMPARE_COEFS_1
-      LOGICAL COMPARE_COEFS_2
-      LOGICAL COMPARE_COEFS_3
-      LOGICAL COMPARE_COEFS_4
-      LOGICAL COMPARE_COEFS_5
-      LOGICAL COMPARE_COEFS_6
-C     
-C     GLOBAL VARIABLES
-C     
-      INCLUDE 'coupl.inc'
-      INTEGER CTMODE
-      REAL*8 LSCALE
-      COMMON/CT/LSCALE,CTMODE
-
-      INTEGER ID,SQSOINDEX,R
-      COMMON/LOOP/ID,SQSOINDEX,R
-
-      LOGICAL CTINIT, TIRINIT, GOLEMINIT, SAMURAIINIT, NINJAINIT
-     $ ,COLLIERINIT
-      COMMON/REDUCTIONCODEINIT/CTINIT, TIRINIT,GOLEMINIT,SAMURAIINIT
-     $ ,NINJAINIT,COLLIERINIT
-
-      INTEGER NLOOPGROUPS
-      PARAMETER (NLOOPGROUPS=1)
-      INTEGER NSQUAREDSO
-      PARAMETER (NSQUAREDSO=1)
-      INCLUDE 'loop_max_coefs.inc'
-
-      COMPLEX*16 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
-      COMMON/LCOEFS/LOOPCOEFS
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     The CT initialization is also performed here if not done already
-C      because it calls MPINIT of OneLOop which is necessary on some
-C      system
-      IF (CTINIT) THEN
-        CTINIT=.FALSE.
-        CALL INITCT()
-      ENDIF
-
-C     INITIALIZE GOLEM IF NEEDED
-      IF (GOLEMINIT) THEN
-        GOLEMINIT=.FALSE.
-        CALL INITGOLEM()
-      ENDIF
-
-C     No stability test intrisic to Golem95 now
-      STABLE=.TRUE.
-
-C     This initialization must be done for each reduction because we
-C     have not setup anyoptimization using pinched propagators yet.
-      CALL INITGOLEM95(NLOOPLINE)
-      PINCH = 0
-
-C     YOU CAN FIND THE DETAILS ABOUT THE DIFFERENT CTMODE AT THE
-C      BEGINNING OF THE FILE CTS_CUTS.F90 IN THE CUTTOOLS DISTRIBUTION
-
-C     CONVERT THE MASSES TO BE COMPLEX
-      DO I=1,NLOOPLINE
-        M2LGOLEM(I)=M2L(I)
-      ENDDO
-
-C     CONVERT THE MOMENTA FLOWING IN THE LOOP LINES TO CT CONVENTIONS
-      DO I=0,3
-        ABSPCT(I) = 0.D0
-        DO J=0,(NLOOPLINE-1)
-          PCT(I,J)=0.D0
-        ENDDO
-      ENDDO
-      DO I=0,3
-        DO J=1,NLOOPLINE
-          PCT(I,0)=PCT(I,0)+PL(I,J)
-          ABSPCT(I)=ABSPCT(I)+ABS(PL(I,J))
-        ENDDO
-      ENDDO
-      REF_P = MAX(ABSPCT(0), ABSPCT(1),ABSPCT(2),ABSPCT(3))
-      DO I=0,3
-        ABSPCT(I) = MAX(REF_P*1E-6, ABSPCT(I))
-      ENDDO
-      IF (CHECKPCONSERVATION.AND.REF_P.GT.1D-8) THEN
-        IF ((PCT(0,0)/ABSPCT(0)).GT.1.D-6) THEN
-          WRITE(*,*) 'energy is not conserved ',PCT(0,0)
-          STOP 'energy is not conserved'
-        ELSEIF ((PCT(1,0)/ABSPCT(1)).GT.1.D-6) THEN
-          WRITE(*,*) 'px is not conserved ',PCT(1,0)
-          STOP 'px is not conserved'
-        ELSEIF ((PCT(2,0)/ABSPCT(2)).GT.1.D-6) THEN
-          WRITE(*,*) 'py is not conserved ',PCT(2,0)
-          STOP 'py is not conserved'
-        ELSEIF ((PCT(3,0)/ABSPCT(3)).GT.1.D-6) THEN
-          WRITE(*,*) 'pz is not conserved ',PCT(3,0)
-          STOP 'pz is not conserved'
-        ENDIF
-      ENDIF
-      DO I=0,3
-        DO J=1,(NLOOPLINE-1)
-          DO K=1,J
-            PCT(I,J)=PCT(I,J)+PL(I,K)
-          ENDDO
-        ENDDO
-      ENDDO
-
-C     Now convert the loop momenta to Golem95 conventions
-      DO I=0,3
-        PGOLEM(1,I)=0.0E0_KI
-        DO J=2,NLOOPLINE
-          PGOLEM(J,I)=PCT(I,J-1)
-        ENDDO
-      ENDDO
-
-C     Fill in the kinematic s-matrix while taking care of on-shell
-C      limits.
-      CALL SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2LGOLEM)
-C     Construct the golem internal matrices derived from the kinetic
-C      one.
-      CALL PREPARESMATRIX()
-
-C     Fill in the golem coefficents and compute the loop
-      IF(GOLEM_RUN_MODE.EQ.2)THEN
-        RES_GOLEM = EVALUATE_B(GOLEM_LOOPNUM,PGOLEM,0,RANK)
-      ELSE
-        PASS_COEF_CHECK=.TRUE.
-        SELECT CASE(RANK)
-        CASE(0)
-        CALL FILL_GOLEM_COEFFS_0(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS0)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          COEFFS0_REC = GOLEM_LOOPNUM(NULL_VEC,0.0_KI)
-          PASS_COEF_CHECK=COMPARE_COEFS_0(COEFFS0,COEFFS0_REC)
-        ENDIF
-        CASE(1)
-        CALL FILL_GOLEM_COEFFS_1(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS1)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT1(GOLEM_LOOPNUM,COEFFS1_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_1(COEFFS1,COEFFS1_REC)
-        ENDIF
-        CASE(2)
-        CALL FILL_GOLEM_COEFFS_2(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS2)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT2(GOLEM_LOOPNUM,COEFFS2_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_2(COEFFS2,COEFFS2_REC)
-        ENDIF
-        CASE(3)
-        CALL FILL_GOLEM_COEFFS_3(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS3)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT3(GOLEM_LOOPNUM,COEFFS3_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_3(COEFFS3,COEFFS3_REC)
-        ENDIF
-        CASE(4)
-        CALL FILL_GOLEM_COEFFS_4(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS4)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT4(GOLEM_LOOPNUM,COEFFS4_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_4(COEFFS4,COEFFS4_REC)
-        ENDIF
-        CASE(5)
-        CALL FILL_GOLEM_COEFFS_5(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS5)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT5(GOLEM_LOOPNUM,COEFFS5_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_5(COEFFS5,COEFFS5_REC)
-        ENDIF
-        CASE(6)
-        CALL FILL_GOLEM_COEFFS_6(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS6)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT6(GOLEM_LOOPNUM,COEFFS6_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_6(COEFFS6,COEFFS6_REC)
-        ENDIF
-        CASE DEFAULT
-        WRITE(*,*)'Not yet implemented in Golem95 for rank= ',RANK
-        STOP
-        END SELECT
-
-        IF(.NOT.PASS_COEF_CHECK)THEN
-          WRITE(*,*)'Coefs mismatch for ID ',ID,' and rank ',RANK
-          WRITE(*,*)'Coefs form MadLoop5:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6)
-          END SELECT
-          WRITE(*,*)'Coefs reconstructed by Golem95:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0_REC
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1_REC)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2_REC)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3_REC)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4_REC)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5_REC)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6_REC)
-          END SELECT
-          STOP
-        ENDIF
-
-        SELECT CASE(NLOOPLINE)
-        CASE(1)
-        WRITE(*,*)'Golem95 cannot handle with tadpole yet'
-        STOP
-        CASE(2)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A20(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT2_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT2_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT2_3(COEFFS3,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',2,RANK
-        STOP
-        END SELECT
-        CASE(3)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A30(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT3_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT3_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT3_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT3_4(COEFFS4,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',3,RANK
-        STOP
-        END SELECT
-        CASE(4)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A40(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT4_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT4_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT4_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT4_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT4_5(COEFFS5,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',4,RANK
-        STOP
-        END SELECT
-        CASE(5)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A50(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT5_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT5_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT5_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT5_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT5_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT5_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',5,RANK
-        STOP
-        END SELECT
-        CASE(6)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A60(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT6_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT6_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT6_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT6_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT6_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT6_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',6,RANK
-        STOP
-        END SELECT
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N = ',NLOOPLINE
-        STOP
-        END SELECT
-      ENDIF
-
-      RES(1)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%C+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%B+2.0*LOG(MU_R)**2*RES_GOLEM%%A)
-      RES(2)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%B+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%A)
-      RES(3)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%A)
-C     WRITE(*,*) 'Loop ID',ID,' =',RES(1),RES(2),RES(3)
-
-C     Finally free golem memory and cache
-      CALL EXITGOLEM95()
-
-      END
-
-      FUNCTION GOLEM_LOOPNUM(Q,MU2)
-      USE PRECISION_GOLEM, ONLY: KI
-      REAL(KI),DIMENSION(0:3),INTENT(IN)::Q
-      REAL(KI),INTENT(IN)::MU2
-      COMPLEX(KI)::GOLEM_LOOPNUM
-
-      COMPLEX*16 QQ(0:3),NUM
-      INTEGER I
-
-      DO I=0,3
-        QQ(I)=CMPLX(Q(I),0.0D0,KIND=16)
-      ENDDO
-
-      CALL LOOPNUM(QQ,NUM)
-      GOLEM_LOOPNUM=NUM
-      RETURN
-      END FUNCTION
-
-      SUBROUTINE INITGOLEM()
-C     
-C     INITIALISATION OF GOLEM
-C     
-C     
-C     MODULE
-C     
-      USE PARAMETRE
-C     
-C     LOCAL VARIABLES 
-C     
-      REAL*8 THRS
-      LOGICAL EXT_NUM_FOR_R1
-C     
-C     GLOBAL VARIABLES 
-C     
-      INCLUDE 'MadLoopParams.inc'
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     DEFAULT PARAMETERS FOR GOLEM
-C     -------------------------------  
-C     One can chose here to have either just the rational R1 piece 
-C     or everything but the R2
-      RAT_OR_TOT_PAR = TOT
-
-      END
-
-      SUBROUTINE SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2L)
-C     
-C     MODULE
-C     
-      USE MATRICE_S
-      USE PRECISION_GOLEM, ONLY: KI
-C     
-C     ARGUMENTS
-C     
-      INTEGER NLOOPLINE
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX(KI) M2L(NLOOPLINE)
-C     
-C     LOCAL VARIABLES
-C     
-      INTEGER I,J
-      COMPLEX*16 S_MAT_FROM_MG(NLOOPLINE,NLOOPLINE)
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-      CALL BUILD_KINEMATIC_MATRIX(NLOOPLINE,PGOLEM,M2L,S_MAT_FROM_MG)
-
-      DO I=1,NLOOPLINE
-        DO J=1,NLOOPLINE
-          S_MAT(I,J)=S_MAT_FROM_MG(I,J)
-        ENDDO
-      ENDDO
-
-      END
-
-      FUNCTION COMPARE_COEFS_0(COEFS_A,COEFS_B)
-
-      USE PRECISION_GOLEM, ONLY: KI
-      COMPLEX(KI) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_0
-
-      NUM = ABS(COEFS_A-COEFS_B)
-      DENOM = ABS(COEFS_A+COEFS_B)
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_0=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_0=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_1(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_1)COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_1
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )
-
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_1=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_1=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_2(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_2) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_2
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_2=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_2=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_3(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_3) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_3
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_3=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_3=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_4(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_4) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_4
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_4=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_4=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_5(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_5) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_5
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_5=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_5=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_6(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_6) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_6
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_6=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_6=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-
-      SUBROUTINE FILL_GOLEM_COEFFS_0(ML_COEFS,GOLEM_COEFS)
-      USE PRECISION_GOLEM, ONLY: KI
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      COMPLEX(KI) GOLEM_COEFS
-      GOLEM_COEFS=ML_COEFS(0)
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_1(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_1
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_1) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_2(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_2
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_2) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(0)^2
-      GOLEM_COEFS%%C1(1,2)= ML_COEFS(5)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(1)^2
-      GOLEM_COEFS%%C1(2,2)= ML_COEFS(7)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(2)^2
-      GOLEM_COEFS%%C1(3,2)= ML_COEFS(10)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-C     Coefficient q(3)^2
-      GOLEM_COEFS%%C1(4,2)= ML_COEFS(14)
-C     Coefficient q(0)*q(1)
-      GOLEM_COEFS%%C2(1,1)= ML_COEFS(6)
-C     Coefficient q(0)*q(2)
-      GOLEM_COEFS%%C2(2,1)= ML_COEFS(8)
-C     Coefficient q(0)*q(3)
-      GOLEM_COEFS%%C2(3,1)= ML_COEFS(11)
-C     Coefficient q(1)*q(2)
-      GOLEM_COEFS%%C2(4,1)= ML_COEFS(9)
-C     Coefficient q(1)*q(3)
-      GOLEM_COEFS%%C2(5,1)= ML_COEFS(12)
-C     Coefficient q(2)*q(3)
-      GOLEM_COEFS%%C2(6,1)= ML_COEFS(13)
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_3(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_3
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_3) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_3
-      STOP 'ERROR: 3 > 2'
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_4(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_4
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_4) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_4
-      STOP 'ERROR: 4 > 2'
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_5(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_5
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_5) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_5
-      STOP 'ERROR: 5 > 2'
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_6(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_6
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_6) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_6
-      STOP 'ERROR: 6 > 2'
-      END
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%TIR_interface.f'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%TIR_interface.f	2016-08-07 07:16:07 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%TIR_interface.f	2018-11-09 08:24:41 +0000
@@ -96,17 +96,9 @@
       ENDIF
 
       IF (MLREDUCTIONLIB(I_LIB).EQ.4) THEN
-C       Using Golem95
-C       PDEN is dummy for Golem95 so we just initialize it to zero
-C        here so as to use it for the function SWITCHORDER
-        DO I=0,3
-          DO J=1,NLOOPLINE-1
-            PDEN(I,J)=0.0D0
-          ENDDO
-        ENDDO
-        CALL SWITCH_ORDER(CTMODE,NLOOPLINE,PL,PDEN,M2L)
-        CALL GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-        RETURN
+C       Golem95 not available
+        WRITE(*,*) 'ERROR:: Golem95 is not interfaced.'
+        STOP
       ENDIF
 
 C     INITIALIZE TIR IF NEEDED
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%loop_matrix.f'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%loop_matrix.f	2017-08-14 07:04:15 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%loop_matrix.f	2018-11-09 08:24:41 +0000
@@ -244,7 +244,7 @@
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING LOOPLIBS ARE
 C      AVAILABLE OR NOT
       LOGICAL LOOPLIBS_AVAILABLE(NLOOPLIB)
-      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.TRUE.,.FALSE.
+      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.FALSE.,.FALSE.
      $ ,.TRUE.,.TRUE./
       COMMON/LOOPLIBS_AV/ LOOPLIBS_AVAILABLE
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING DIRECTION TESTS
 
=== removed file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%GOLEM_interface.f'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%GOLEM_interface.f	2016-08-07 07:16:07 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%GOLEM_interface.f	1970-01-01 00:00:00 +0000
@@ -1,748 +0,0 @@
-      SUBROUTINE GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-C     
-C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
-C     By the MadGraph5_aMC@NLO Development Team
-C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
-C     
-C     Interface between MG5 and Golem95.
-C     The Golem95 version should be higher than 1.3.0.
-C     It supports RANK = NLOOPLINE + 1 tensor integrals when 1 <
-C      NLOOPLINE < 6.
-C     
-C     Process: u d~ > w+ QED<=1 WEIGHTED<=2 [ all = QCD ]
-C     Process: c s~ > w+ QED<=1 WEIGHTED<=2 [ all = QCD ]
-C     
-C     
-C     MODULES
-C     
-      USE MATRICE_S
-      USE FORM_FACTOR_TYPE, ONLY: FORM_FACTOR
-      USE PRECISION_GOLEM, ONLY: KI
-      USE TENS_COMB
-      USE TENS_REC
-      USE FORM_FACTOR_1P, ONLY: A10
-      USE FORM_FACTOR_2P, ONLY: A20
-      USE FORM_FACTOR_3P, ONLY: A30
-      USE FORM_FACTOR_4P, ONLY: A40
-      USE FORM_FACTOR_5P, ONLY: A50
-      USE FORM_FACTOR_6P, ONLY: A60
-      IMPLICIT NONE
-C     
-C     CONSTANTS 
-C     
-      INTEGER    NEXTERNAL
-      PARAMETER (NEXTERNAL=3)
-      LOGICAL CHECKPCONSERVATION
-      PARAMETER (CHECKPCONSERVATION=.TRUE.)
-      REAL*8 NORMALIZATION
-      PARAMETER (NORMALIZATION = 1.D0/(16.D0*3.14159265358979323846D0*
-     $ *2))
-      REAL(KI),DIMENSION(0:3),PARAMETER::NULL_VEC = (/0.0_KI,0.0_KI
-     $ ,0.0_KI,0.0_KI/)
-C     GOLEM_RUN_MODE = 1: Use directly MadLoop tensorial coefficients
-C     GOLEM_RUN_MODE = 2: Reconstruct the tensorial coefficeints
-C      directly from 
-C     numerator using golem internal reconstruction routine
-C     GOLEM_RUN_MODE = 3: Cross-checked reconstructed coefficients
-C      against
-C     MadLoop internal ones.
-      INTEGER GOLEM_RUN_MODE
-      PARAMETER (GOLEM_RUN_MODE=1)
-C     The following is the acceptance threshold used for
-C      GOLEM_RUN_MODE = 3
-      REAL*8 COEF_CHECK_THRS
-      DATA COEF_CHECK_THRS/1.0D-13/
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-
-      LOGICAL PASS_COEF_CHECK
-C     
-C     ARGUMENTS 
-C     
-      INTEGER NLOOPLINE, RANK
-      REAL*8 PL(0:3,NLOOPLINE)
-      REAL*8 PCT(0:3,0:NLOOPLINE-1), ABSPCT(0:3)
-      REAL*8 REF_P
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX*16 M2L(NLOOPLINE)
-      COMPLEX(KI) M2LGOLEM(NLOOPLINE)
-      COMPLEX*16 RES(3)
-      LOGICAL STABLE
-C     
-C     LOCAL VARIABLES 
-C     
-      INTEGER I, J, K
-      TYPE(FORM_FACTOR)::RES_GOLEM
-
-      COMPLEX(KI)::COEFFS0,COEFFS0_REC
-      TYPE(COEFF_TYPE_1)::COEFFS1,COEFFS1_REC
-      TYPE(COEFF_TYPE_2)::COEFFS2,COEFFS2_REC
-      TYPE(COEFF_TYPE_3)::COEFFS3,COEFFS3_REC
-      TYPE(COEFF_TYPE_4)::COEFFS4,COEFFS4_REC
-      TYPE(COEFF_TYPE_5)::COEFFS5,COEFFS5_REC
-      TYPE(COEFF_TYPE_6)::COEFFS6,COEFFS6_REC
-
-C     The pinch propagator optimization is not used, so for now it is
-C     always 0.
-      INTEGER PINCH
-C     
-C     EXTERNAL FUNCTIONS
-C     
-      COMPLEX(KI) GOLEM_LOOPNUM
-      EXTERNAL GOLEM_LOOPNUM
-      LOGICAL COMPARE_COEFS_0
-      LOGICAL COMPARE_COEFS_1
-      LOGICAL COMPARE_COEFS_2
-      LOGICAL COMPARE_COEFS_3
-      LOGICAL COMPARE_COEFS_4
-      LOGICAL COMPARE_COEFS_5
-      LOGICAL COMPARE_COEFS_6
-C     
-C     GLOBAL VARIABLES
-C     
-      INCLUDE 'coupl.inc'
-      INTEGER CTMODE
-      REAL*8 LSCALE
-      COMMON/CT/LSCALE,CTMODE
-
-      INTEGER ID,SQSOINDEX,R
-      COMMON/LOOP/ID,SQSOINDEX,R
-
-      LOGICAL CTINIT, TIRINIT, GOLEMINIT, SAMURAIINIT, NINJAINIT
-     $ ,COLLIERINIT
-      COMMON/REDUCTIONCODEINIT/CTINIT, TIRINIT,GOLEMINIT,SAMURAIINIT
-     $ ,NINJAINIT,COLLIERINIT
-
-      INTEGER NLOOPGROUPS
-      PARAMETER (NLOOPGROUPS=1)
-      INTEGER NSQUAREDSO
-      PARAMETER (NSQUAREDSO=1)
-      INCLUDE 'loop_max_coefs.inc'
-
-      COMPLEX*16 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
-      COMMON/LCOEFS/LOOPCOEFS
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     The CT initialization is also performed here if not done already
-C      because it calls MPINIT of OneLOop which is necessary on some
-C      system
-      IF (CTINIT) THEN
-        CTINIT=.FALSE.
-        CALL INITCT()
-      ENDIF
-
-C     INITIALIZE GOLEM IF NEEDED
-      IF (GOLEMINIT) THEN
-        GOLEMINIT=.FALSE.
-        CALL INITGOLEM()
-      ENDIF
-
-C     No stability test intrisic to Golem95 now
-      STABLE=.TRUE.
-
-C     This initialization must be done for each reduction because we
-C     have not setup anyoptimization using pinched propagators yet.
-      CALL INITGOLEM95(NLOOPLINE)
-      PINCH = 0
-
-C     YOU CAN FIND THE DETAILS ABOUT THE DIFFERENT CTMODE AT THE
-C      BEGINNING OF THE FILE CTS_CUTS.F90 IN THE CUTTOOLS DISTRIBUTION
-
-C     CONVERT THE MASSES TO BE COMPLEX
-      DO I=1,NLOOPLINE
-        M2LGOLEM(I)=M2L(I)
-      ENDDO
-
-C     CONVERT THE MOMENTA FLOWING IN THE LOOP LINES TO CT CONVENTIONS
-      DO I=0,3
-        ABSPCT(I) = 0.D0
-        DO J=0,(NLOOPLINE-1)
-          PCT(I,J)=0.D0
-        ENDDO
-      ENDDO
-      DO I=0,3
-        DO J=1,NLOOPLINE
-          PCT(I,0)=PCT(I,0)+PL(I,J)
-          ABSPCT(I)=ABSPCT(I)+ABS(PL(I,J))
-        ENDDO
-      ENDDO
-      REF_P = MAX(ABSPCT(0), ABSPCT(1),ABSPCT(2),ABSPCT(3))
-      DO I=0,3
-        ABSPCT(I) = MAX(REF_P*1E-6, ABSPCT(I))
-      ENDDO
-      IF (CHECKPCONSERVATION.AND.REF_P.GT.1D-8) THEN
-        IF ((PCT(0,0)/ABSPCT(0)).GT.1.D-6) THEN
-          WRITE(*,*) 'energy is not conserved ',PCT(0,0)
-          STOP 'energy is not conserved'
-        ELSEIF ((PCT(1,0)/ABSPCT(1)).GT.1.D-6) THEN
-          WRITE(*,*) 'px is not conserved ',PCT(1,0)
-          STOP 'px is not conserved'
-        ELSEIF ((PCT(2,0)/ABSPCT(2)).GT.1.D-6) THEN
-          WRITE(*,*) 'py is not conserved ',PCT(2,0)
-          STOP 'py is not conserved'
-        ELSEIF ((PCT(3,0)/ABSPCT(3)).GT.1.D-6) THEN
-          WRITE(*,*) 'pz is not conserved ',PCT(3,0)
-          STOP 'pz is not conserved'
-        ENDIF
-      ENDIF
-      DO I=0,3
-        DO J=1,(NLOOPLINE-1)
-          DO K=1,J
-            PCT(I,J)=PCT(I,J)+PL(I,K)
-          ENDDO
-        ENDDO
-      ENDDO
-
-C     Now convert the loop momenta to Golem95 conventions
-      DO I=0,3
-        PGOLEM(1,I)=0.0E0_KI
-        DO J=2,NLOOPLINE
-          PGOLEM(J,I)=PCT(I,J-1)
-        ENDDO
-      ENDDO
-
-C     Fill in the kinematic s-matrix while taking care of on-shell
-C      limits.
-      CALL SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2LGOLEM)
-C     Construct the golem internal matrices derived from the kinetic
-C      one.
-      CALL PREPARESMATRIX()
-
-C     Fill in the golem coefficents and compute the loop
-      IF(GOLEM_RUN_MODE.EQ.2)THEN
-        RES_GOLEM = EVALUATE_B(GOLEM_LOOPNUM,PGOLEM,0,RANK)
-      ELSE
-        PASS_COEF_CHECK=.TRUE.
-        SELECT CASE(RANK)
-        CASE(0)
-        CALL FILL_GOLEM_COEFFS_0(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS0)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          COEFFS0_REC = GOLEM_LOOPNUM(NULL_VEC,0.0_KI)
-          PASS_COEF_CHECK=COMPARE_COEFS_0(COEFFS0,COEFFS0_REC)
-        ENDIF
-        CASE(1)
-        CALL FILL_GOLEM_COEFFS_1(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS1)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT1(GOLEM_LOOPNUM,COEFFS1_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_1(COEFFS1,COEFFS1_REC)
-        ENDIF
-        CASE(2)
-        CALL FILL_GOLEM_COEFFS_2(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS2)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT2(GOLEM_LOOPNUM,COEFFS2_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_2(COEFFS2,COEFFS2_REC)
-        ENDIF
-        CASE(3)
-        CALL FILL_GOLEM_COEFFS_3(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS3)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT3(GOLEM_LOOPNUM,COEFFS3_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_3(COEFFS3,COEFFS3_REC)
-        ENDIF
-        CASE(4)
-        CALL FILL_GOLEM_COEFFS_4(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS4)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT4(GOLEM_LOOPNUM,COEFFS4_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_4(COEFFS4,COEFFS4_REC)
-        ENDIF
-        CASE(5)
-        CALL FILL_GOLEM_COEFFS_5(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS5)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT5(GOLEM_LOOPNUM,COEFFS5_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_5(COEFFS5,COEFFS5_REC)
-        ENDIF
-        CASE(6)
-        CALL FILL_GOLEM_COEFFS_6(LOOPCOEFS(0,SQSOINDEX,ID),COEFFS6)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT6(GOLEM_LOOPNUM,COEFFS6_REC)
-          PASS_COEF_CHECK=COMPARE_COEFS_6(COEFFS6,COEFFS6_REC)
-        ENDIF
-        CASE DEFAULT
-        WRITE(*,*)'Not yet implemented in Golem95 for rank= ',RANK
-        STOP
-        END SELECT
-
-        IF(.NOT.PASS_COEF_CHECK)THEN
-          WRITE(*,*)'Coefs mismatch for ID ',ID,' and rank ',RANK
-          WRITE(*,*)'Coefs form MadLoop5:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6)
-          END SELECT
-          WRITE(*,*)'Coefs reconstructed by Golem95:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0_REC
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1_REC)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2_REC)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3_REC)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4_REC)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5_REC)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6_REC)
-          END SELECT
-          STOP
-        ENDIF
-
-        SELECT CASE(NLOOPLINE)
-        CASE(1)
-        WRITE(*,*)'Golem95 cannot handle with tadpole yet'
-        STOP
-        CASE(2)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A20(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT2_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT2_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT2_3(COEFFS3,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',2,RANK
-        STOP
-        END SELECT
-        CASE(3)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A30(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT3_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT3_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT3_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT3_4(COEFFS4,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',3,RANK
-        STOP
-        END SELECT
-        CASE(4)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A40(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT4_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT4_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT4_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT4_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT4_5(COEFFS5,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',4,RANK
-        STOP
-        END SELECT
-        CASE(5)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A50(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT5_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT5_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT5_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT5_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT5_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT5_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',5,RANK
-        STOP
-        END SELECT
-        CASE(6)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A60(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT6_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT6_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT6_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT6_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT6_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT6_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',6,RANK
-        STOP
-        END SELECT
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N = ',NLOOPLINE
-        STOP
-        END SELECT
-      ENDIF
-
-      RES(1)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%C+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%B+2.0*LOG(MU_R)**2*RES_GOLEM%%A)
-      RES(2)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%B+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%A)
-      RES(3)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%A)
-C     WRITE(*,*) 'Loop ID',ID,' =',RES(1),RES(2),RES(3)
-
-C     Finally free golem memory and cache
-      CALL EXITGOLEM95()
-
-      END
-
-      FUNCTION GOLEM_LOOPNUM(Q,MU2)
-      USE PRECISION_GOLEM, ONLY: KI
-      REAL(KI),DIMENSION(0:3),INTENT(IN)::Q
-      REAL(KI),INTENT(IN)::MU2
-      COMPLEX(KI)::GOLEM_LOOPNUM
-
-      COMPLEX*16 QQ(0:3),NUM
-      INTEGER I
-
-      DO I=0,3
-        QQ(I)=CMPLX(Q(I),0.0D0,KIND=16)
-      ENDDO
-
-      CALL LOOPNUM(QQ,NUM)
-      GOLEM_LOOPNUM=NUM
-      RETURN
-      END FUNCTION
-
-      SUBROUTINE INITGOLEM()
-C     
-C     INITIALISATION OF GOLEM
-C     
-C     
-C     MODULE
-C     
-      USE PARAMETRE
-C     
-C     LOCAL VARIABLES 
-C     
-      REAL*8 THRS
-      LOGICAL EXT_NUM_FOR_R1
-C     
-C     GLOBAL VARIABLES 
-C     
-      INCLUDE 'MadLoopParams.inc'
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     DEFAULT PARAMETERS FOR GOLEM
-C     -------------------------------  
-C     One can chose here to have either just the rational R1 piece 
-C     or everything but the R2
-      RAT_OR_TOT_PAR = TOT
-
-      END
-
-      SUBROUTINE SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2L)
-C     
-C     MODULE
-C     
-      USE MATRICE_S
-      USE PRECISION_GOLEM, ONLY: KI
-C     
-C     ARGUMENTS
-C     
-      INTEGER NLOOPLINE
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX(KI) M2L(NLOOPLINE)
-C     
-C     LOCAL VARIABLES
-C     
-      INTEGER I,J
-      COMPLEX*16 S_MAT_FROM_MG(NLOOPLINE,NLOOPLINE)
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-      CALL BUILD_KINEMATIC_MATRIX(NLOOPLINE,PGOLEM,M2L,S_MAT_FROM_MG)
-
-      DO I=1,NLOOPLINE
-        DO J=1,NLOOPLINE
-          S_MAT(I,J)=S_MAT_FROM_MG(I,J)
-        ENDDO
-      ENDDO
-
-      END
-
-      FUNCTION COMPARE_COEFS_0(COEFS_A,COEFS_B)
-
-      USE PRECISION_GOLEM, ONLY: KI
-      COMPLEX(KI) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_0
-
-      NUM = ABS(COEFS_A-COEFS_B)
-      DENOM = ABS(COEFS_A+COEFS_B)
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_0=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_0=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_1(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_1)COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_1
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )
-
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_1=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_1=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_2(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_2) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_2
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_2=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_2=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_3(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_3) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_3
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_3=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_3=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_4(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_4) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_4
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_4=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_4=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_5(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_5) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_5
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_5=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_5=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION COMPARE_COEFS_6(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_6) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL COMPARE_COEFS_6
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        COMPARE_COEFS_6=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        COMPARE_COEFS_6=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-
-      SUBROUTINE FILL_GOLEM_COEFFS_0(ML_COEFS,GOLEM_COEFS)
-      USE PRECISION_GOLEM, ONLY: KI
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      COMPLEX(KI) GOLEM_COEFS
-      GOLEM_COEFS=ML_COEFS(0)
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_1(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_1
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_1) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_2(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_2
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_2) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(0)^2
-      GOLEM_COEFS%%C1(1,2)= ML_COEFS(5)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(1)^2
-      GOLEM_COEFS%%C1(2,2)= ML_COEFS(7)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(2)^2
-      GOLEM_COEFS%%C1(3,2)= ML_COEFS(10)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-C     Coefficient q(3)^2
-      GOLEM_COEFS%%C1(4,2)= ML_COEFS(14)
-C     Coefficient q(0)*q(1)
-      GOLEM_COEFS%%C2(1,1)= ML_COEFS(6)
-C     Coefficient q(0)*q(2)
-      GOLEM_COEFS%%C2(2,1)= ML_COEFS(8)
-C     Coefficient q(0)*q(3)
-      GOLEM_COEFS%%C2(3,1)= ML_COEFS(11)
-C     Coefficient q(1)*q(2)
-      GOLEM_COEFS%%C2(4,1)= ML_COEFS(9)
-C     Coefficient q(1)*q(3)
-      GOLEM_COEFS%%C2(5,1)= ML_COEFS(12)
-C     Coefficient q(2)*q(3)
-      GOLEM_COEFS%%C2(6,1)= ML_COEFS(13)
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_3(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_3
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_3) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_3
-      STOP 'ERROR: 3 > 2'
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_4(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_4
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_4) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_4
-      STOP 'ERROR: 4 > 2'
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_5(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_5
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_5) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_5
-      STOP 'ERROR: 5 > 2'
-      END
-
-      SUBROUTINE FILL_GOLEM_COEFFS_6(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_6
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_6) GOLEM_COEFS
-C     Dummy routine for FILL_GOLEM_COEFS_6
-      STOP 'ERROR: 6 > 2'
-      END
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%TIR_interface.f'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%TIR_interface.f	2016-08-07 07:16:07 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%TIR_interface.f	2018-11-09 08:24:41 +0000
@@ -96,17 +96,9 @@
       ENDIF
 
       IF (MLREDUCTIONLIB(I_LIB).EQ.4) THEN
-C       Using Golem95
-C       PDEN is dummy for Golem95 so we just initialize it to zero
-C        here so as to use it for the function SWITCHORDER
-        DO I=0,3
-          DO J=1,NLOOPLINE-1
-            PDEN(I,J)=0.0D0
-          ENDDO
-        ENDDO
-        CALL SWITCH_ORDER(CTMODE,NLOOPLINE,PL,PDEN,M2L)
-        CALL GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-        RETURN
+C       Golem95 not available
+        WRITE(*,*) 'ERROR:: Golem95 is not interfaced.'
+        STOP
       ENDIF
 
 C     INITIALIZE TIR IF NEEDED
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%loop_matrix.f'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%loop_matrix.f	2017-08-14 07:04:15 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%loop_matrix.f	2018-11-09 08:24:41 +0000
@@ -244,7 +244,7 @@
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING LOOPLIBS ARE
 C      AVAILABLE OR NOT
       LOGICAL LOOPLIBS_AVAILABLE(NLOOPLIB)
-      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.TRUE.,.FALSE.
+      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.FALSE.,.FALSE.
      $ ,.TRUE.,.TRUE./
       COMMON/LOOPLIBS_AV/ LOOPLIBS_AVAILABLE
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING DIRECTION TESTS
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_group/matrix1.f'
--- tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_group/matrix1.f	2017-11-24 12:17:04 +0000
+++ tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_group/matrix1.f	2018-11-09 08:24:41 +0000
@@ -73,6 +73,9 @@
       REAL*8 POL(2)
       COMMON/TO_POLARIZATION/ POL
 
+      DOUBLE PRECISION SMALL_WIDTH_TREATMENT
+      COMMON/NARROW_WIDTH/SMALL_WIDTH_TREATMENT
+
       INTEGER          ISUM_HEL
       LOGICAL                    MULTI_CHANNEL
       COMMON/TO_MATRIX/ISUM_HEL, MULTI_CHANNEL
@@ -110,6 +113,7 @@
 C     ----------
 C     BEGIN CODE
 C     ----------
+
       NTRY(IMIRROR)=NTRY(IMIRROR)+1
       THIS_NTRY(IMIRROR) = THIS_NTRY(IMIRROR)+1
       DO I=1,NEXTERNAL
@@ -290,6 +294,15 @@
 C     Needed for v4 models
       COMPLEX*16 DUM0,DUM1
       DATA DUM0, DUM1/(0D0, 0D0), (1D0, 0D0)/
+
+      DOUBLE PRECISION FK_ZERO
+      DOUBLE PRECISION FK_WZ
+      SAVE FK_ZERO
+      SAVE FK_WZ
+
+      LOGICAL FIRST
+      DATA FIRST /.TRUE./
+      SAVE FIRST
 C     
 C     FUNCTION
 C     
@@ -300,6 +313,9 @@
       DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
       COMMON/TO_AMPS/  AMP2,       JAMP2
       INCLUDE 'coupl.inc'
+
+      DOUBLE PRECISION SMALL_WIDTH_TREATMENT
+      COMMON/NARROW_WIDTH/SMALL_WIDTH_TREATMENT
 C     
 C     COLOR DATA
 C     
@@ -312,26 +328,33 @@
 C     ----------
 C     BEGIN CODE
 C     ----------
+      IF (FIRST) THEN
+        FIRST=.FALSE.
+        FK_WZ = SIGN(MAX(ABS(WZ), ABS(MZ*SMALL_WIDTH_TREATMENT)), WZ)
+        FK_ZERO = 0D0
+      ENDIF
+
+
       CALL IXXXXX(P(0,1),ZERO,NHEL(1),+1*IC(1),W(1,1))
       CALL OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
       CALL OXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
       CALL IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
-      CALL FFV1_3(W(1,1),W(1,2),GQQ,ZERO,ZERO,W(1,5))
+      CALL FFV1_3(W(1,1),W(1,2),GQQ,ZERO, FK_ZERO,W(1,5))
 C     Amplitude(s) for diagram number 1
       CALL FFV1_0(W(1,4),W(1,3),W(1,5),GQQ,AMP(1))
-      CALL FFV1_3(W(1,1),W(1,2),GQED,ZERO,ZERO,W(1,5))
+      CALL FFV1_3(W(1,1),W(1,2),GQED,ZERO, FK_ZERO,W(1,5))
 C     Amplitude(s) for diagram number 2
       CALL FFV1_0(W(1,4),W(1,3),W(1,5),GQED,AMP(2))
-      CALL FFV1_2_3(W(1,1),W(1,2),GUZ1,GUZ2,MZ,WZ,W(1,5))
+      CALL FFV1_2_3(W(1,1),W(1,2),GUZ1,GUZ2,MZ, FK_WZ,W(1,5))
 C     Amplitude(s) for diagram number 3
       CALL FFV1_2_0(W(1,4),W(1,3),W(1,5),GUZ1,GUZ2,AMP(3))
-      CALL FFV1_3(W(1,1),W(1,3),GQQ,ZERO,ZERO,W(1,5))
+      CALL FFV1_3(W(1,1),W(1,3),GQQ,ZERO, FK_ZERO,W(1,5))
 C     Amplitude(s) for diagram number 4
       CALL FFV1_0(W(1,4),W(1,2),W(1,5),GQQ,AMP(4))
-      CALL FFV1_3(W(1,1),W(1,3),GQED,ZERO,ZERO,W(1,5))
+      CALL FFV1_3(W(1,1),W(1,3),GQED,ZERO, FK_ZERO,W(1,5))
 C     Amplitude(s) for diagram number 5
       CALL FFV1_0(W(1,4),W(1,2),W(1,5),GQED,AMP(5))
-      CALL FFV1_2_3(W(1,1),W(1,3),GUZ1,GUZ2,MZ,WZ,W(1,5))
+      CALL FFV1_2_3(W(1,1),W(1,3),GUZ1,GUZ2,MZ, FK_WZ,W(1,5))
 C     Amplitude(s) for diagram number 6
       CALL FFV1_2_0(W(1,4),W(1,2),W(1,5),GUZ1,GUZ2,AMP(6))
 C     JAMPs contributing to orders ALL_ORDERS=1
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_nogroup/matrix.f'
--- tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_nogroup/matrix.f	2017-11-26 00:21:49 +0000
+++ tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_nogroup/matrix.f	2018-11-09 08:24:41 +0000
@@ -51,6 +51,7 @@
       INTEGER IDUM, NGOOD, IGOOD(NCOMB), JHEL, J, JJ
       REAL     XRAN1
       EXTERNAL XRAN1
+
 C     
 C     GLOBAL VARIABLES
 C     
@@ -71,6 +72,7 @@
       DATA IDUM /-1/
       DATA XTRY, XREJ, NGOOD /0,0,0/
       SAVE YFRAC, IGOOD, JHEL
+
       DATA (NHEL(I,   1),I=1,4) / 1,-1,-1,-1/
       DATA (NHEL(I,   2),I=1,4) / 1,-1,-1, 1/
       DATA (NHEL(I,   3),I=1,4) / 1,-1, 1,-1/
@@ -261,6 +263,13 @@
 C     Needed for v4 models
       COMPLEX*16 DUM0,DUM1
       DATA DUM0, DUM1/(0D0, 0D0), (1D0, 0D0)/
+
+      DOUBLE PRECISION FK_ZERO
+      SAVE FK_ZERO
+
+      LOGICAL FIRST
+      DATA FIRST /.TRUE./
+      SAVE FIRST
 C     
 C     FUNCTION
 C     
@@ -271,6 +280,9 @@
       DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
       COMMON/TO_AMPS/  AMP2,       JAMP2
       INCLUDE 'coupl.inc'
+
+      DOUBLE PRECISION SMALL_WIDTH_TREATMENT
+      COMMON/NARROW_WIDTH/SMALL_WIDTH_TREATMENT
 C     
 C     COLOR DATA
 C     
@@ -283,17 +295,22 @@
 C     ----------
 C     BEGIN CODE
 C     ----------
+      IF (FIRST) THEN
+        FIRST=.FALSE.
+        FK_ZERO = 0D0
+      ENDIF
+
       CALL IXXXXX(P(0,1),ZERO,NHEL(1),+1*IC(1),W(1,1))
       CALL OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
       CALL VXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
       CALL VXXXXX(P(0,4),ZERO,NHEL(4),+1*IC(4),W(1,4))
-      CALL FFV1_3(W(1,1),W(1,2),GQQ,ZERO,ZERO,W(1,5))
+      CALL FFV1_3(W(1,1),W(1,2),GQQ,ZERO, FK_ZERO,W(1,5))
 C     Amplitude(s) for diagram number 1
       CALL VVV1_0(W(1,3),W(1,4),W(1,5),G,AMP(1))
-      CALL FFV1_2(W(1,1),W(1,3),GQQ,ZERO,ZERO,W(1,5))
+      CALL FFV1_2(W(1,1),W(1,3),GQQ,ZERO, FK_ZERO,W(1,5))
 C     Amplitude(s) for diagram number 2
       CALL FFV1_0(W(1,5),W(1,2),W(1,4),GQQ,AMP(2))
-      CALL FFV1_2(W(1,1),W(1,4),GQQ,ZERO,ZERO,W(1,5))
+      CALL FFV1_2(W(1,1),W(1,4),GQQ,ZERO, FK_ZERO,W(1,5))
 C     Amplitude(s) for diagram number 3
       CALL FFV1_0(W(1,5),W(1,2),W(1,3),GQQ,AMP(3))
 C     JAMPs contributing to orders ALL_ORDERS=1
 
=== removed file 'tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%GOLEM_interface.f'
--- tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%GOLEM_interface.f	2016-05-26 03:24:09 +0000
+++ tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%GOLEM_interface.f	1970-01-01 00:00:00 +0000
@@ -1,823 +0,0 @@
-      SUBROUTINE ML5_0_GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-C     
-C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
-C     By the MadGraph5_aMC@NLO Development Team
-C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
-C     
-C     Interface between MG5 and Golem95.
-C     The Golem95 version should be higher than 1.3.0.
-C     It supports RANK = NLOOPLINE + 1 tensor integrals when 1 <
-C      NLOOPLINE < 6.
-C     
-C     Process: g g > t t~ [ virt = QCD ]
-C     
-C     
-C     MODULES
-C     
-      USE MATRICE_S
-      USE FORM_FACTOR_TYPE, ONLY: FORM_FACTOR
-      USE PRECISION_GOLEM, ONLY: KI
-      USE TENS_COMB
-      USE TENS_REC
-      USE FORM_FACTOR_1P, ONLY: A10
-      USE FORM_FACTOR_2P, ONLY: A20
-      USE FORM_FACTOR_3P, ONLY: A30
-      USE FORM_FACTOR_4P, ONLY: A40
-      USE FORM_FACTOR_5P, ONLY: A50
-      USE FORM_FACTOR_6P, ONLY: A60
-      IMPLICIT NONE
-C     
-C     CONSTANTS 
-C     
-      INTEGER    NEXTERNAL
-      PARAMETER (NEXTERNAL=4)
-      LOGICAL CHECKPCONSERVATION
-      PARAMETER (CHECKPCONSERVATION=.TRUE.)
-      REAL*8 NORMALIZATION
-      PARAMETER (NORMALIZATION = 1.D0/(16.D0*3.14159265358979323846D0*
-     $ *2))
-      REAL(KI),DIMENSION(0:3),PARAMETER::NULL_VEC = (/0.0_KI,0.0_KI
-     $ ,0.0_KI,0.0_KI/)
-C     GOLEM_RUN_MODE = 1: Use directly MadLoop tensorial coefficients
-C     GOLEM_RUN_MODE = 2: Reconstruct the tensorial coefficeints
-C      directly from 
-C     numerator using golem internal reconstruction routine
-C     GOLEM_RUN_MODE = 3: Cross-checked reconstructed coefficients
-C      against
-C     MadLoop internal ones.
-      INTEGER GOLEM_RUN_MODE
-      PARAMETER (GOLEM_RUN_MODE=1)
-C     The following is the acceptance threshold used for
-C      GOLEM_RUN_MODE = 3
-      REAL*8 COEF_CHECK_THRS
-      DATA COEF_CHECK_THRS/1.0D-13/
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-
-      LOGICAL PASS_COEF_CHECK
-C     
-C     ARGUMENTS 
-C     
-      INTEGER NLOOPLINE, RANK
-      REAL*8 PL(0:3,NLOOPLINE)
-      REAL*8 PCT(0:3,0:NLOOPLINE-1), ABSPCT(0:3)
-      REAL*8 REF_P
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX*16 M2L(NLOOPLINE)
-      COMPLEX(KI) M2LGOLEM(NLOOPLINE)
-      COMPLEX*16 RES(3)
-      LOGICAL STABLE
-C     
-C     LOCAL VARIABLES 
-C     
-      INTEGER I, J, K
-      TYPE(FORM_FACTOR)::RES_GOLEM
-
-      COMPLEX(KI)::COEFFS0,COEFFS0_REC
-      TYPE(COEFF_TYPE_1)::COEFFS1,COEFFS1_REC
-      TYPE(COEFF_TYPE_2)::COEFFS2,COEFFS2_REC
-      TYPE(COEFF_TYPE_3)::COEFFS3,COEFFS3_REC
-      TYPE(COEFF_TYPE_4)::COEFFS4,COEFFS4_REC
-      TYPE(COEFF_TYPE_5)::COEFFS5,COEFFS5_REC
-      TYPE(COEFF_TYPE_6)::COEFFS6,COEFFS6_REC
-
-C     The pinch propagator optimization is not used, so for now it is
-C     always 0.
-      INTEGER PINCH
-C     
-C     EXTERNAL FUNCTIONS
-C     
-      COMPLEX(KI) ML5_0_GOLEM_LOOPNUM
-      EXTERNAL ML5_0_GOLEM_LOOPNUM
-      LOGICAL ML5_0_COMPARE_COEFS_0
-      LOGICAL ML5_0_COMPARE_COEFS_1
-      LOGICAL ML5_0_COMPARE_COEFS_2
-      LOGICAL ML5_0_COMPARE_COEFS_3
-      LOGICAL ML5_0_COMPARE_COEFS_4
-      LOGICAL ML5_0_COMPARE_COEFS_5
-      LOGICAL ML5_0_COMPARE_COEFS_6
-C     
-C     GLOBAL VARIABLES
-C     
-      INCLUDE 'coupl.inc'
-      INTEGER CTMODE
-      REAL*8 LSCALE
-      COMMON/ML5_0_CT/LSCALE,CTMODE
-
-      INTEGER ID,SQSOINDEX,R
-      COMMON/ML5_0_LOOP/ID,SQSOINDEX,R
-
-      LOGICAL CTINIT, TIRINIT, GOLEMINIT, SAMURAIINIT, NINJAINIT
-     $ ,COLLIERINIT
-      COMMON/REDUCTIONCODEINIT/CTINIT, TIRINIT,GOLEMINIT,SAMURAIINIT
-     $ ,NINJAINIT,COLLIERINIT
-
-      INTEGER NLOOPGROUPS
-      PARAMETER (NLOOPGROUPS=26)
-      INTEGER NSQUAREDSO
-      PARAMETER (NSQUAREDSO=1)
-      INCLUDE 'loop_max_coefs.inc'
-
-      COMPLEX*16 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
-      COMMON/ML5_0_LCOEFS/LOOPCOEFS
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     The CT initialization is also performed here if not done already
-C      because it calls MPINIT of OneLOop which is necessary on some
-C      system
-      IF (CTINIT) THEN
-        CTINIT=.FALSE.
-        CALL ML5_0_INITCT()
-      ENDIF
-
-C     INITIALIZE GOLEM IF NEEDED
-      IF (GOLEMINIT) THEN
-        GOLEMINIT=.FALSE.
-        CALL ML5_0_INITGOLEM()
-      ENDIF
-
-C     No stability test intrisic to Golem95 now
-      STABLE=.TRUE.
-
-C     This initialization must be done for each reduction because we
-C     have not setup anyoptimization using pinched propagators yet.
-      CALL INITGOLEM95(NLOOPLINE)
-      PINCH = 0
-
-C     YOU CAN FIND THE DETAILS ABOUT THE DIFFERENT CTMODE AT THE
-C      BEGINNING OF THE FILE CTS_CUTS.F90 IN THE CUTTOOLS DISTRIBUTION
-
-C     CONVERT THE MASSES TO BE COMPLEX
-      DO I=1,NLOOPLINE
-        M2LGOLEM(I)=M2L(I)
-      ENDDO
-
-C     CONVERT THE MOMENTA FLOWING IN THE LOOP LINES TO CT CONVENTIONS
-      DO I=0,3
-        ABSPCT(I) = 0.D0
-        DO J=0,(NLOOPLINE-1)
-          PCT(I,J)=0.D0
-        ENDDO
-      ENDDO
-      DO I=0,3
-        DO J=1,NLOOPLINE
-          PCT(I,0)=PCT(I,0)+PL(I,J)
-          ABSPCT(I)=ABSPCT(I)+ABS(PL(I,J))
-        ENDDO
-      ENDDO
-      REF_P = MAX(ABSPCT(0), ABSPCT(1),ABSPCT(2),ABSPCT(3))
-      DO I=0,3
-        ABSPCT(I) = MAX(REF_P*1E-6, ABSPCT(I))
-      ENDDO
-      IF (CHECKPCONSERVATION.AND.REF_P.GT.1D-8) THEN
-        IF ((PCT(0,0)/ABSPCT(0)).GT.1.D-6) THEN
-          WRITE(*,*) 'energy is not conserved ',PCT(0,0)
-          STOP 'energy is not conserved'
-        ELSEIF ((PCT(1,0)/ABSPCT(1)).GT.1.D-6) THEN
-          WRITE(*,*) 'px is not conserved ',PCT(1,0)
-          STOP 'px is not conserved'
-        ELSEIF ((PCT(2,0)/ABSPCT(2)).GT.1.D-6) THEN
-          WRITE(*,*) 'py is not conserved ',PCT(2,0)
-          STOP 'py is not conserved'
-        ELSEIF ((PCT(3,0)/ABSPCT(3)).GT.1.D-6) THEN
-          WRITE(*,*) 'pz is not conserved ',PCT(3,0)
-          STOP 'pz is not conserved'
-        ENDIF
-      ENDIF
-      DO I=0,3
-        DO J=1,(NLOOPLINE-1)
-          DO K=1,J
-            PCT(I,J)=PCT(I,J)+PL(I,K)
-          ENDDO
-        ENDDO
-      ENDDO
-
-C     Now convert the loop momenta to Golem95 conventions
-      DO I=0,3
-        PGOLEM(1,I)=0.0E0_KI
-        DO J=2,NLOOPLINE
-          PGOLEM(J,I)=PCT(I,J-1)
-        ENDDO
-      ENDDO
-
-C     Fill in the kinematic s-matrix while taking care of on-shell
-C      limits.
-      CALL ML5_0_SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2LGOLEM)
-C     Construct the golem internal matrices derived from the kinetic
-C      one.
-      CALL PREPARESMATRIX()
-
-C     Fill in the golem coefficents and compute the loop
-      IF(GOLEM_RUN_MODE.EQ.2)THEN
-        RES_GOLEM = EVALUATE_B(ML5_0_GOLEM_LOOPNUM,PGOLEM,0,RANK)
-      ELSE
-        PASS_COEF_CHECK=.TRUE.
-        SELECT CASE(RANK)
-        CASE(0)
-        CALL ML5_0_FILL_GOLEM_COEFFS_0(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS0)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          COEFFS0_REC = ML5_0_GOLEM_LOOPNUM(NULL_VEC,0.0_KI)
-          PASS_COEF_CHECK=ML5_0_COMPARE_COEFS_0(COEFFS0,COEFFS0_REC)
-        ENDIF
-        CASE(1)
-        CALL ML5_0_FILL_GOLEM_COEFFS_1(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS1)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT1(ML5_0_GOLEM_LOOPNUM,COEFFS1_REC)
-          PASS_COEF_CHECK=ML5_0_COMPARE_COEFS_1(COEFFS1,COEFFS1_REC)
-        ENDIF
-        CASE(2)
-        CALL ML5_0_FILL_GOLEM_COEFFS_2(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS2)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT2(ML5_0_GOLEM_LOOPNUM,COEFFS2_REC)
-          PASS_COEF_CHECK=ML5_0_COMPARE_COEFS_2(COEFFS2,COEFFS2_REC)
-        ENDIF
-        CASE(3)
-        CALL ML5_0_FILL_GOLEM_COEFFS_3(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS3)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT3(ML5_0_GOLEM_LOOPNUM,COEFFS3_REC)
-          PASS_COEF_CHECK=ML5_0_COMPARE_COEFS_3(COEFFS3,COEFFS3_REC)
-        ENDIF
-        CASE(4)
-        CALL ML5_0_FILL_GOLEM_COEFFS_4(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS4)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT4(ML5_0_GOLEM_LOOPNUM,COEFFS4_REC)
-          PASS_COEF_CHECK=ML5_0_COMPARE_COEFS_4(COEFFS4,COEFFS4_REC)
-        ENDIF
-        CASE(5)
-        CALL ML5_0_FILL_GOLEM_COEFFS_5(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS5)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT5(ML5_0_GOLEM_LOOPNUM,COEFFS5_REC)
-          PASS_COEF_CHECK=ML5_0_COMPARE_COEFS_5(COEFFS5,COEFFS5_REC)
-        ENDIF
-        CASE(6)
-        CALL ML5_0_FILL_GOLEM_COEFFS_6(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS6)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT6(ML5_0_GOLEM_LOOPNUM,COEFFS6_REC)
-          PASS_COEF_CHECK=ML5_0_COMPARE_COEFS_6(COEFFS6,COEFFS6_REC)
-        ENDIF
-        CASE DEFAULT
-        WRITE(*,*)'Not yet implemented in Golem95 for rank= ',RANK
-        STOP
-        END SELECT
-
-        IF(.NOT.PASS_COEF_CHECK)THEN
-          WRITE(*,*)'Coefs mismatch for ID ',ID,' and rank ',RANK
-          WRITE(*,*)'Coefs form MadLoop5:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6)
-          END SELECT
-          WRITE(*,*)'Coefs reconstructed by Golem95:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0_REC
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1_REC)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2_REC)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3_REC)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4_REC)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5_REC)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6_REC)
-          END SELECT
-          STOP
-        ENDIF
-
-        SELECT CASE(NLOOPLINE)
-        CASE(1)
-        WRITE(*,*)'Golem95 cannot handle with tadpole yet'
-        STOP
-        CASE(2)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A20(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT2_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT2_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT2_3(COEFFS3,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',2,RANK
-        STOP
-        END SELECT
-        CASE(3)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A30(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT3_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT3_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT3_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT3_4(COEFFS4,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',3,RANK
-        STOP
-        END SELECT
-        CASE(4)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A40(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT4_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT4_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT4_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT4_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT4_5(COEFFS5,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',4,RANK
-        STOP
-        END SELECT
-        CASE(5)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A50(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT5_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT5_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT5_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT5_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT5_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT5_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',5,RANK
-        STOP
-        END SELECT
-        CASE(6)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A60(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT6_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT6_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT6_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT6_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT6_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT6_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',6,RANK
-        STOP
-        END SELECT
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N = ',NLOOPLINE
-        STOP
-        END SELECT
-      ENDIF
-
-      RES(1)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%C+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%B+2.0*LOG(MU_R)**2*RES_GOLEM%%A)
-      RES(2)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%B+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%A)
-      RES(3)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%A)
-C     WRITE(*,*) 'Loop ID',ID,' =',RES(1),RES(2),RES(3)
-
-C     Finally free golem memory and cache
-      CALL EXITGOLEM95()
-
-      END
-
-      FUNCTION ML5_0_GOLEM_LOOPNUM(Q,MU2)
-      USE PRECISION_GOLEM, ONLY: KI
-      REAL(KI),DIMENSION(0:3),INTENT(IN)::Q
-      REAL(KI),INTENT(IN)::MU2
-      COMPLEX(KI)::ML5_0_GOLEM_LOOPNUM
-
-      COMPLEX*16 QQ(0:3),NUM
-      INTEGER I
-
-      DO I=0,3
-        QQ(I)=CMPLX(Q(I),0.0D0,KIND=16)
-      ENDDO
-
-      CALL ML5_0_LOOPNUM(QQ,NUM)
-      ML5_0_GOLEM_LOOPNUM=NUM
-      RETURN
-      END FUNCTION
-
-      SUBROUTINE ML5_0_INITGOLEM()
-C     
-C     INITIALISATION OF GOLEM
-C     
-C     
-C     MODULE
-C     
-      USE PARAMETRE
-C     
-C     LOCAL VARIABLES 
-C     
-      REAL*8 THRS
-      LOGICAL EXT_NUM_FOR_R1
-C     
-C     GLOBAL VARIABLES 
-C     
-      INCLUDE 'MadLoopParams.inc'
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     DEFAULT PARAMETERS FOR GOLEM
-C     -------------------------------  
-C     One can chose here to have either just the rational R1 piece 
-C     or everything but the R2
-      RAT_OR_TOT_PAR = TOT
-
-      END
-
-      SUBROUTINE ML5_0_SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2L)
-C     
-C     MODULE
-C     
-      USE MATRICE_S
-      USE PRECISION_GOLEM, ONLY: KI
-C     
-C     ARGUMENTS
-C     
-      INTEGER NLOOPLINE
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX(KI) M2L(NLOOPLINE)
-C     
-C     LOCAL VARIABLES
-C     
-      INTEGER I,J
-      COMPLEX*16 S_MAT_FROM_MG(NLOOPLINE,NLOOPLINE)
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-      CALL ML5_0_BUILD_KINEMATIC_MATRIX(NLOOPLINE,PGOLEM,M2L
-     $ ,S_MAT_FROM_MG)
-
-      DO I=1,NLOOPLINE
-        DO J=1,NLOOPLINE
-          S_MAT(I,J)=S_MAT_FROM_MG(I,J)
-        ENDDO
-      ENDDO
-
-      END
-
-      FUNCTION ML5_0_COMPARE_COEFS_0(COEFS_A,COEFS_B)
-
-      USE PRECISION_GOLEM, ONLY: KI
-      COMPLEX(KI) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL ML5_0_COMPARE_COEFS_0
-
-      NUM = ABS(COEFS_A-COEFS_B)
-      DENOM = ABS(COEFS_A+COEFS_B)
-      IF(DENOM.GT.0D0)THEN
-        ML5_0_COMPARE_COEFS_0=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        ML5_0_COMPARE_COEFS_0=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION ML5_0_COMPARE_COEFS_1(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_1)COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL ML5_0_COMPARE_COEFS_1
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )
-
-      IF(DENOM.GT.0D0)THEN
-        ML5_0_COMPARE_COEFS_1=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        ML5_0_COMPARE_COEFS_1=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION ML5_0_COMPARE_COEFS_2(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_2) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL ML5_0_COMPARE_COEFS_2
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))
-      IF(DENOM.GT.0D0)THEN
-        ML5_0_COMPARE_COEFS_2=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        ML5_0_COMPARE_COEFS_2=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION ML5_0_COMPARE_COEFS_3(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_3) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL ML5_0_COMPARE_COEFS_3
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )
-      IF(DENOM.GT.0D0)THEN
-        ML5_0_COMPARE_COEFS_3=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        ML5_0_COMPARE_COEFS_3=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION ML5_0_COMPARE_COEFS_4(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_4) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL ML5_0_COMPARE_COEFS_4
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        ML5_0_COMPARE_COEFS_4=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        ML5_0_COMPARE_COEFS_4=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION ML5_0_COMPARE_COEFS_5(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_5) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL ML5_0_COMPARE_COEFS_5
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        ML5_0_COMPARE_COEFS_5=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        ML5_0_COMPARE_COEFS_5=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION ML5_0_COMPARE_COEFS_6(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_6) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/ML5_0_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL ML5_0_COMPARE_COEFS_6
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        ML5_0_COMPARE_COEFS_6=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        ML5_0_COMPARE_COEFS_6=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-
-      SUBROUTINE ML5_0_FILL_GOLEM_COEFFS_0(ML_COEFS,GOLEM_COEFS)
-      USE PRECISION_GOLEM, ONLY: KI
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      COMPLEX(KI) GOLEM_COEFS
-      GOLEM_COEFS=ML_COEFS(0)
-      END
-
-      SUBROUTINE ML5_0_FILL_GOLEM_COEFFS_1(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_1
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_1) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-      END
-
-      SUBROUTINE ML5_0_FILL_GOLEM_COEFFS_2(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_2
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_2) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(0)^2
-      GOLEM_COEFS%%C1(1,2)= ML_COEFS(5)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(1)^2
-      GOLEM_COEFS%%C1(2,2)= ML_COEFS(7)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(2)^2
-      GOLEM_COEFS%%C1(3,2)= ML_COEFS(10)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-C     Coefficient q(3)^2
-      GOLEM_COEFS%%C1(4,2)= ML_COEFS(14)
-C     Coefficient q(0)*q(1)
-      GOLEM_COEFS%%C2(1,1)= ML_COEFS(6)
-C     Coefficient q(0)*q(2)
-      GOLEM_COEFS%%C2(2,1)= ML_COEFS(8)
-C     Coefficient q(0)*q(3)
-      GOLEM_COEFS%%C2(3,1)= ML_COEFS(11)
-C     Coefficient q(1)*q(2)
-      GOLEM_COEFS%%C2(4,1)= ML_COEFS(9)
-C     Coefficient q(1)*q(3)
-      GOLEM_COEFS%%C2(5,1)= ML_COEFS(12)
-C     Coefficient q(2)*q(3)
-      GOLEM_COEFS%%C2(6,1)= ML_COEFS(13)
-      END
-
-      SUBROUTINE ML5_0_FILL_GOLEM_COEFFS_3(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_3
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_3) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(0)^2
-      GOLEM_COEFS%%C1(1,2)= ML_COEFS(5)
-C     Coefficient q(0)^3
-      GOLEM_COEFS%%C1(1,3)=-ML_COEFS(15)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(1)^2
-      GOLEM_COEFS%%C1(2,2)= ML_COEFS(7)
-C     Coefficient q(1)^3
-      GOLEM_COEFS%%C1(2,3)=-ML_COEFS(18)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(2)^2
-      GOLEM_COEFS%%C1(3,2)= ML_COEFS(10)
-C     Coefficient q(2)^3
-      GOLEM_COEFS%%C1(3,3)=-ML_COEFS(24)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-C     Coefficient q(3)^2
-      GOLEM_COEFS%%C1(4,2)= ML_COEFS(14)
-C     Coefficient q(3)^3
-      GOLEM_COEFS%%C1(4,3)=-ML_COEFS(34)
-C     Coefficient q(0)*q(1)
-      GOLEM_COEFS%%C2(1,1)= ML_COEFS(6)
-C     Coefficient q(0)*q(1)^2
-      GOLEM_COEFS%%C2(1,2)=-ML_COEFS(17)
-C     Coefficient q(0)^2*q(1)
-      GOLEM_COEFS%%C2(1,3)=-ML_COEFS(16)
-C     Coefficient q(0)*q(2)
-      GOLEM_COEFS%%C2(2,1)= ML_COEFS(8)
-C     Coefficient q(0)*q(2)^2
-      GOLEM_COEFS%%C2(2,2)=-ML_COEFS(22)
-C     Coefficient q(0)^2*q(2)
-      GOLEM_COEFS%%C2(2,3)=-ML_COEFS(19)
-C     Coefficient q(0)*q(3)
-      GOLEM_COEFS%%C2(3,1)= ML_COEFS(11)
-C     Coefficient q(0)*q(3)^2
-      GOLEM_COEFS%%C2(3,2)=-ML_COEFS(31)
-C     Coefficient q(0)^2*q(3)
-      GOLEM_COEFS%%C2(3,3)=-ML_COEFS(25)
-C     Coefficient q(1)*q(2)
-      GOLEM_COEFS%%C2(4,1)= ML_COEFS(9)
-C     Coefficient q(1)*q(2)^2
-      GOLEM_COEFS%%C2(4,2)=-ML_COEFS(23)
-C     Coefficient q(1)^2*q(2)
-      GOLEM_COEFS%%C2(4,3)=-ML_COEFS(21)
-C     Coefficient q(1)*q(3)
-      GOLEM_COEFS%%C2(5,1)= ML_COEFS(12)
-C     Coefficient q(1)*q(3)^2
-      GOLEM_COEFS%%C2(5,2)=-ML_COEFS(32)
-C     Coefficient q(1)^2*q(3)
-      GOLEM_COEFS%%C2(5,3)=-ML_COEFS(27)
-C     Coefficient q(2)*q(3)
-      GOLEM_COEFS%%C2(6,1)= ML_COEFS(13)
-C     Coefficient q(2)*q(3)^2
-      GOLEM_COEFS%%C2(6,2)=-ML_COEFS(33)
-C     Coefficient q(2)^2*q(3)
-      GOLEM_COEFS%%C2(6,3)=-ML_COEFS(30)
-C     Coefficient q(0)*q(1)*q(2)
-      GOLEM_COEFS%%C3(1,1)=-ML_COEFS(20)
-C     Coefficient q(0)*q(1)*q(3)
-      GOLEM_COEFS%%C3(2,1)=-ML_COEFS(26)
-C     Coefficient q(0)*q(2)*q(3)
-      GOLEM_COEFS%%C3(3,1)=-ML_COEFS(28)
-C     Coefficient q(1)*q(2)*q(3)
-      GOLEM_COEFS%%C3(4,1)=-ML_COEFS(29)
-      END
-
-      SUBROUTINE ML5_0_FILL_GOLEM_COEFFS_4(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_4
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_4) GOLEM_COEFS
-C     Dummy routine for ML5_0_FILL_GOLEM_COEFS_4
-      STOP 'ERROR: 4 > 3'
-      END
-
-      SUBROUTINE ML5_0_FILL_GOLEM_COEFFS_5(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_5
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_5) GOLEM_COEFS
-C     Dummy routine for ML5_0_FILL_GOLEM_COEFS_5
-      STOP 'ERROR: 5 > 3'
-      END
-
-      SUBROUTINE ML5_0_FILL_GOLEM_COEFFS_6(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_6
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_6) GOLEM_COEFS
-C     Dummy routine for ML5_0_FILL_GOLEM_COEFS_6
-      STOP 'ERROR: 6 > 3'
-      END
 
=== modified file 'tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%TIR_interface.f'
--- tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%TIR_interface.f	2016-06-28 23:55:35 +0000
+++ tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%TIR_interface.f	2018-11-09 08:24:41 +0000
@@ -95,17 +95,9 @@
       ENDIF
 
       IF (MLREDUCTIONLIB(I_LIB).EQ.4) THEN
-C       Using Golem95
-C       PDEN is dummy for Golem95 so we just initialize it to zero
-C        here so as to use it for the function SWITCHORDER
-        DO I=0,3
-          DO J=1,NLOOPLINE-1
-            PDEN(I,J)=0.0D0
-          ENDDO
-        ENDDO
-        CALL ML5_0_SWITCH_ORDER(CTMODE,NLOOPLINE,PL,PDEN,M2L)
-        CALL ML5_0_GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-        RETURN
+C       Golem95 not available
+        WRITE(*,*) 'ERROR:: Golem95 is not interfaced.'
+        STOP
       ENDIF
 
 C     INITIALIZE TIR IF NEEDED
 
=== modified file 'tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%loop_matrix.f'
--- tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%loop_matrix.f	2017-08-14 07:04:15 +0000
+++ tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%loop_matrix.f	2018-11-09 08:24:41 +0000
@@ -243,7 +243,7 @@
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING LOOPLIBS ARE
 C      AVAILABLE OR NOT
       LOGICAL LOOPLIBS_AVAILABLE(NLOOPLIB)
-      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.TRUE.,.FALSE.
+      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.FALSE.,.FALSE.
      $ ,.TRUE.,.TRUE./
       COMMON/ML5_0_LOOPLIBS_AV/ LOOPLIBS_AVAILABLE
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING DIRECTION TESTS
 
=== removed file 'tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%GOLEM_interface.f'
--- tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%GOLEM_interface.f	2016-05-26 03:24:09 +0000
+++ tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%GOLEM_interface.f	1970-01-01 00:00:00 +0000
@@ -1,755 +0,0 @@
-      SUBROUTINE MG5_1_GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-C     
-C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
-C     By the MadGraph5_aMC@NLO Development Team
-C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
-C     
-C     Interface between MG5 and Golem95.
-C     The Golem95 version should be higher than 1.3.0.
-C     It supports RANK = NLOOPLINE + 1 tensor integrals when 1 <
-C      NLOOPLINE < 6.
-C     
-C     Process: u u~ > u u~ [ virt = QCD ] @1
-C     
-C     
-C     MODULES
-C     
-      USE MATRICE_S
-      USE FORM_FACTOR_TYPE, ONLY: FORM_FACTOR
-      USE PRECISION_GOLEM, ONLY: KI
-      USE TENS_COMB
-      USE TENS_REC
-      USE FORM_FACTOR_1P, ONLY: A10
-      USE FORM_FACTOR_2P, ONLY: A20
-      USE FORM_FACTOR_3P, ONLY: A30
-      USE FORM_FACTOR_4P, ONLY: A40
-      USE FORM_FACTOR_5P, ONLY: A50
-      USE FORM_FACTOR_6P, ONLY: A60
-      IMPLICIT NONE
-C     
-C     CONSTANTS 
-C     
-      INTEGER    NEXTERNAL
-      PARAMETER (NEXTERNAL=4)
-      LOGICAL CHECKPCONSERVATION
-      PARAMETER (CHECKPCONSERVATION=.TRUE.)
-      REAL*8 NORMALIZATION
-      PARAMETER (NORMALIZATION = 1.D0/(16.D0*3.14159265358979323846D0*
-     $ *2))
-      REAL(KI),DIMENSION(0:3),PARAMETER::NULL_VEC = (/0.0_KI,0.0_KI
-     $ ,0.0_KI,0.0_KI/)
-C     GOLEM_RUN_MODE = 1: Use directly MadLoop tensorial coefficients
-C     GOLEM_RUN_MODE = 2: Reconstruct the tensorial coefficeints
-C      directly from 
-C     numerator using golem internal reconstruction routine
-C     GOLEM_RUN_MODE = 3: Cross-checked reconstructed coefficients
-C      against
-C     MadLoop internal ones.
-      INTEGER GOLEM_RUN_MODE
-      PARAMETER (GOLEM_RUN_MODE=1)
-C     The following is the acceptance threshold used for
-C      GOLEM_RUN_MODE = 3
-      REAL*8 COEF_CHECK_THRS
-      DATA COEF_CHECK_THRS/1.0D-13/
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-
-      LOGICAL PASS_COEF_CHECK
-C     
-C     ARGUMENTS 
-C     
-      INTEGER NLOOPLINE, RANK
-      REAL*8 PL(0:3,NLOOPLINE)
-      REAL*8 PCT(0:3,0:NLOOPLINE-1), ABSPCT(0:3)
-      REAL*8 REF_P
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX*16 M2L(NLOOPLINE)
-      COMPLEX(KI) M2LGOLEM(NLOOPLINE)
-      COMPLEX*16 RES(3)
-      LOGICAL STABLE
-C     
-C     LOCAL VARIABLES 
-C     
-      INTEGER I, J, K
-      TYPE(FORM_FACTOR)::RES_GOLEM
-
-      COMPLEX(KI)::COEFFS0,COEFFS0_REC
-      TYPE(COEFF_TYPE_1)::COEFFS1,COEFFS1_REC
-      TYPE(COEFF_TYPE_2)::COEFFS2,COEFFS2_REC
-      TYPE(COEFF_TYPE_3)::COEFFS3,COEFFS3_REC
-      TYPE(COEFF_TYPE_4)::COEFFS4,COEFFS4_REC
-      TYPE(COEFF_TYPE_5)::COEFFS5,COEFFS5_REC
-      TYPE(COEFF_TYPE_6)::COEFFS6,COEFFS6_REC
-
-C     The pinch propagator optimization is not used, so for now it is
-C     always 0.
-      INTEGER PINCH
-C     
-C     EXTERNAL FUNCTIONS
-C     
-      COMPLEX(KI) MG5_1_GOLEM_LOOPNUM
-      EXTERNAL MG5_1_GOLEM_LOOPNUM
-      LOGICAL MG5_1_COMPARE_COEFS_0
-      LOGICAL MG5_1_COMPARE_COEFS_1
-      LOGICAL MG5_1_COMPARE_COEFS_2
-      LOGICAL MG5_1_COMPARE_COEFS_3
-      LOGICAL MG5_1_COMPARE_COEFS_4
-      LOGICAL MG5_1_COMPARE_COEFS_5
-      LOGICAL MG5_1_COMPARE_COEFS_6
-C     
-C     GLOBAL VARIABLES
-C     
-      INCLUDE 'coupl.inc'
-      INTEGER CTMODE
-      REAL*8 LSCALE
-      COMMON/MG5_1_CT/LSCALE,CTMODE
-
-      INTEGER ID,SQSOINDEX,R
-      COMMON/MG5_1_LOOP/ID,SQSOINDEX,R
-
-      LOGICAL CTINIT, TIRINIT, GOLEMINIT, SAMURAIINIT, NINJAINIT
-     $ ,COLLIERINIT
-      COMMON/REDUCTIONCODEINIT/CTINIT, TIRINIT,GOLEMINIT,SAMURAIINIT
-     $ ,NINJAINIT,COLLIERINIT
-
-      INTEGER NLOOPGROUPS
-      PARAMETER (NLOOPGROUPS=13)
-      INTEGER NSQUAREDSO
-      PARAMETER (NSQUAREDSO=1)
-      INCLUDE 'loop_max_coefs.inc'
-
-      COMPLEX*16 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
-      COMMON/MG5_1_LCOEFS/LOOPCOEFS
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     The CT initialization is also performed here if not done already
-C      because it calls MPINIT of OneLOop which is necessary on some
-C      system
-      IF (CTINIT) THEN
-        CTINIT=.FALSE.
-        CALL MG5_1_INITCT()
-      ENDIF
-
-C     INITIALIZE GOLEM IF NEEDED
-      IF (GOLEMINIT) THEN
-        GOLEMINIT=.FALSE.
-        CALL MG5_1_INITGOLEM()
-      ENDIF
-
-C     No stability test intrisic to Golem95 now
-      STABLE=.TRUE.
-
-C     This initialization must be done for each reduction because we
-C     have not setup anyoptimization using pinched propagators yet.
-      CALL INITGOLEM95(NLOOPLINE)
-      PINCH = 0
-
-C     YOU CAN FIND THE DETAILS ABOUT THE DIFFERENT CTMODE AT THE
-C      BEGINNING OF THE FILE CTS_CUTS.F90 IN THE CUTTOOLS DISTRIBUTION
-
-C     CONVERT THE MASSES TO BE COMPLEX
-      DO I=1,NLOOPLINE
-        M2LGOLEM(I)=M2L(I)
-      ENDDO
-
-C     CONVERT THE MOMENTA FLOWING IN THE LOOP LINES TO CT CONVENTIONS
-      DO I=0,3
-        ABSPCT(I) = 0.D0
-        DO J=0,(NLOOPLINE-1)
-          PCT(I,J)=0.D0
-        ENDDO
-      ENDDO
-      DO I=0,3
-        DO J=1,NLOOPLINE
-          PCT(I,0)=PCT(I,0)+PL(I,J)
-          ABSPCT(I)=ABSPCT(I)+ABS(PL(I,J))
-        ENDDO
-      ENDDO
-      REF_P = MAX(ABSPCT(0), ABSPCT(1),ABSPCT(2),ABSPCT(3))
-      DO I=0,3
-        ABSPCT(I) = MAX(REF_P*1E-6, ABSPCT(I))
-      ENDDO
-      IF (CHECKPCONSERVATION.AND.REF_P.GT.1D-8) THEN
-        IF ((PCT(0,0)/ABSPCT(0)).GT.1.D-6) THEN
-          WRITE(*,*) 'energy is not conserved ',PCT(0,0)
-          STOP 'energy is not conserved'
-        ELSEIF ((PCT(1,0)/ABSPCT(1)).GT.1.D-6) THEN
-          WRITE(*,*) 'px is not conserved ',PCT(1,0)
-          STOP 'px is not conserved'
-        ELSEIF ((PCT(2,0)/ABSPCT(2)).GT.1.D-6) THEN
-          WRITE(*,*) 'py is not conserved ',PCT(2,0)
-          STOP 'py is not conserved'
-        ELSEIF ((PCT(3,0)/ABSPCT(3)).GT.1.D-6) THEN
-          WRITE(*,*) 'pz is not conserved ',PCT(3,0)
-          STOP 'pz is not conserved'
-        ENDIF
-      ENDIF
-      DO I=0,3
-        DO J=1,(NLOOPLINE-1)
-          DO K=1,J
-            PCT(I,J)=PCT(I,J)+PL(I,K)
-          ENDDO
-        ENDDO
-      ENDDO
-
-C     Now convert the loop momenta to Golem95 conventions
-      DO I=0,3
-        PGOLEM(1,I)=0.0E0_KI
-        DO J=2,NLOOPLINE
-          PGOLEM(J,I)=PCT(I,J-1)
-        ENDDO
-      ENDDO
-
-C     Fill in the kinematic s-matrix while taking care of on-shell
-C      limits.
-      CALL MG5_1_SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2LGOLEM)
-C     Construct the golem internal matrices derived from the kinetic
-C      one.
-      CALL PREPARESMATRIX()
-
-C     Fill in the golem coefficents and compute the loop
-      IF(GOLEM_RUN_MODE.EQ.2)THEN
-        RES_GOLEM = EVALUATE_B(MG5_1_GOLEM_LOOPNUM,PGOLEM,0,RANK)
-      ELSE
-        PASS_COEF_CHECK=.TRUE.
-        SELECT CASE(RANK)
-        CASE(0)
-        CALL MG5_1_FILL_GOLEM_COEFFS_0(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS0)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          COEFFS0_REC = MG5_1_GOLEM_LOOPNUM(NULL_VEC,0.0_KI)
-          PASS_COEF_CHECK=MG5_1_COMPARE_COEFS_0(COEFFS0,COEFFS0_REC)
-        ENDIF
-        CASE(1)
-        CALL MG5_1_FILL_GOLEM_COEFFS_1(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS1)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT1(MG5_1_GOLEM_LOOPNUM,COEFFS1_REC)
-          PASS_COEF_CHECK=MG5_1_COMPARE_COEFS_1(COEFFS1,COEFFS1_REC)
-        ENDIF
-        CASE(2)
-        CALL MG5_1_FILL_GOLEM_COEFFS_2(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS2)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT2(MG5_1_GOLEM_LOOPNUM,COEFFS2_REC)
-          PASS_COEF_CHECK=MG5_1_COMPARE_COEFS_2(COEFFS2,COEFFS2_REC)
-        ENDIF
-        CASE(3)
-        CALL MG5_1_FILL_GOLEM_COEFFS_3(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS3)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT3(MG5_1_GOLEM_LOOPNUM,COEFFS3_REC)
-          PASS_COEF_CHECK=MG5_1_COMPARE_COEFS_3(COEFFS3,COEFFS3_REC)
-        ENDIF
-        CASE(4)
-        CALL MG5_1_FILL_GOLEM_COEFFS_4(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS4)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT4(MG5_1_GOLEM_LOOPNUM,COEFFS4_REC)
-          PASS_COEF_CHECK=MG5_1_COMPARE_COEFS_4(COEFFS4,COEFFS4_REC)
-        ENDIF
-        CASE(5)
-        CALL MG5_1_FILL_GOLEM_COEFFS_5(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS5)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT5(MG5_1_GOLEM_LOOPNUM,COEFFS5_REC)
-          PASS_COEF_CHECK=MG5_1_COMPARE_COEFS_5(COEFFS5,COEFFS5_REC)
-        ENDIF
-        CASE(6)
-        CALL MG5_1_FILL_GOLEM_COEFFS_6(LOOPCOEFS(0,SQSOINDEX,ID)
-     $   ,COEFFS6)
-        IF(GOLEM_RUN_MODE.EQ.3)THEN
-          CALL RECONSTRUCT6(MG5_1_GOLEM_LOOPNUM,COEFFS6_REC)
-          PASS_COEF_CHECK=MG5_1_COMPARE_COEFS_6(COEFFS6,COEFFS6_REC)
-        ENDIF
-        CASE DEFAULT
-        WRITE(*,*)'Not yet implemented in Golem95 for rank= ',RANK
-        STOP
-        END SELECT
-
-        IF(.NOT.PASS_COEF_CHECK)THEN
-          WRITE(*,*)'Coefs mismatch for ID ',ID,' and rank ',RANK
-          WRITE(*,*)'Coefs form MadLoop5:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6)
-          END SELECT
-          WRITE(*,*)'Coefs reconstructed by Golem95:'
-          SELECT CASE(RANK)
-          CASE(0)
-          WRITE(*,*)'Constant coef = ',COEFFS0_REC
-          CASE(1)
-          CALL PRINT_COEFFS(COEFFS1_REC)
-          CASE(2)
-          CALL PRINT_COEFFS(COEFFS2_REC)
-          CASE(3)
-          CALL PRINT_COEFFS(COEFFS3_REC)
-          CASE(4)
-          CALL PRINT_COEFFS(COEFFS4_REC)
-          CASE(5)
-          CALL PRINT_COEFFS(COEFFS5_REC)
-          CASE(6)
-          CALL PRINT_COEFFS(COEFFS6_REC)
-          END SELECT
-          STOP
-        ENDIF
-
-        SELECT CASE(NLOOPLINE)
-        CASE(1)
-        WRITE(*,*)'Golem95 cannot handle with tadpole yet'
-        STOP
-        CASE(2)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A20(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT2_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT2_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT2_3(COEFFS3,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',2,RANK
-        STOP
-        END SELECT
-        CASE(3)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A30(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT3_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT3_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT3_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT3_4(COEFFS4,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',3,RANK
-        STOP
-        END SELECT
-        CASE(4)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A40(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT4_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT4_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT4_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT4_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT4_5(COEFFS5,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',4,RANK
-        STOP
-        END SELECT
-        CASE(5)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A50(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT5_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT5_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT5_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT5_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT5_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT5_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',5,RANK
-        STOP
-        END SELECT
-        CASE(6)
-        SELECT CASE(RANK)
-        CASE(0)
-        RES_GOLEM = COEFFS0*A60(PINCH)
-        CASE(1)
-        RES_GOLEM = CONTRACT6_1(COEFFS1,PGOLEM,PINCH)
-        CASE(2)
-        RES_GOLEM = CONTRACT6_2(COEFFS2,PGOLEM,PINCH)
-        CASE(3)
-        RES_GOLEM = CONTRACT6_3(COEFFS3,PGOLEM,PINCH)
-        CASE(4)
-        RES_GOLEM = CONTRACT6_4(COEFFS4,PGOLEM,PINCH)
-        CASE(5)
-        RES_GOLEM = CONTRACT6_5(COEFFS5,PGOLEM,PINCH)
-        CASE(6)
-        RES_GOLEM = CONTRACT6_6(COEFFS6,PGOLEM,PINCH)
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N,r = ',6,RANK
-        STOP
-        END SELECT
-        CASE DEFAULT
-        WRITE(*,*)'Golem95 cannot handle with: N = ',NLOOPLINE
-        STOP
-        END SELECT
-      ENDIF
-
-      RES(1)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%C+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%B+2.0*LOG(MU_R)**2*RES_GOLEM%%A)
-      RES(2)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%B+2.0*LOG(MU_R)
-     $ *RES_GOLEM%%A)
-      RES(3)=NORMALIZATION*2.0D0*DBLE(RES_GOLEM%%A)
-C     WRITE(*,*) 'Loop ID',ID,' =',RES(1),RES(2),RES(3)
-
-C     Finally free golem memory and cache
-      CALL EXITGOLEM95()
-
-      END
-
-      FUNCTION MG5_1_GOLEM_LOOPNUM(Q,MU2)
-      USE PRECISION_GOLEM, ONLY: KI
-      REAL(KI),DIMENSION(0:3),INTENT(IN)::Q
-      REAL(KI),INTENT(IN)::MU2
-      COMPLEX(KI)::MG5_1_GOLEM_LOOPNUM
-
-      COMPLEX*16 QQ(0:3),NUM
-      INTEGER I
-
-      DO I=0,3
-        QQ(I)=CMPLX(Q(I),0.0D0,KIND=16)
-      ENDDO
-
-      CALL MG5_1_LOOPNUM(QQ,NUM)
-      MG5_1_GOLEM_LOOPNUM=NUM
-      RETURN
-      END FUNCTION
-
-      SUBROUTINE MG5_1_INITGOLEM()
-C     
-C     INITIALISATION OF GOLEM
-C     
-C     
-C     MODULE
-C     
-      USE PARAMETRE
-C     
-C     LOCAL VARIABLES 
-C     
-      REAL*8 THRS
-      LOGICAL EXT_NUM_FOR_R1
-C     
-C     GLOBAL VARIABLES 
-C     
-      INCLUDE 'MadLoopParams.inc'
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-C     DEFAULT PARAMETERS FOR GOLEM
-C     -------------------------------  
-C     One can chose here to have either just the rational R1 piece 
-C     or everything but the R2
-      RAT_OR_TOT_PAR = TOT
-
-      END
-
-      SUBROUTINE MG5_1_SETUP_KIN_MATRIX(NLOOPLINE,PGOLEM,M2L)
-C     
-C     MODULE
-C     
-      USE MATRICE_S
-      USE PRECISION_GOLEM, ONLY: KI
-C     
-C     ARGUMENTS
-C     
-      INTEGER NLOOPLINE
-      REAL(KI) PGOLEM(NLOOPLINE,0:3)
-      COMPLEX(KI) M2L(NLOOPLINE)
-C     
-C     LOCAL VARIABLES
-C     
-      INTEGER I,J
-      COMPLEX*16 S_MAT_FROM_MG(NLOOPLINE,NLOOPLINE)
-C     ----------
-C     BEGIN CODE
-C     ----------
-
-      CALL MG5_1_BUILD_KINEMATIC_MATRIX(NLOOPLINE,PGOLEM,M2L
-     $ ,S_MAT_FROM_MG)
-
-      DO I=1,NLOOPLINE
-        DO J=1,NLOOPLINE
-          S_MAT(I,J)=S_MAT_FROM_MG(I,J)
-        ENDDO
-      ENDDO
-
-      END
-
-      FUNCTION MG5_1_COMPARE_COEFS_0(COEFS_A,COEFS_B)
-
-      USE PRECISION_GOLEM, ONLY: KI
-      COMPLEX(KI) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL MG5_1_COMPARE_COEFS_0
-
-      NUM = ABS(COEFS_A-COEFS_B)
-      DENOM = ABS(COEFS_A+COEFS_B)
-      IF(DENOM.GT.0D0)THEN
-        MG5_1_COMPARE_COEFS_0=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        MG5_1_COMPARE_COEFS_0=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION MG5_1_COMPARE_COEFS_1(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_1)COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL MG5_1_COMPARE_COEFS_1
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )
-
-      IF(DENOM.GT.0D0)THEN
-        MG5_1_COMPARE_COEFS_1=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        MG5_1_COMPARE_COEFS_1=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION MG5_1_COMPARE_COEFS_2(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_2) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL MG5_1_COMPARE_COEFS_2
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))
-      IF(DENOM.GT.0D0)THEN
-        MG5_1_COMPARE_COEFS_2=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        MG5_1_COMPARE_COEFS_2=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION MG5_1_COMPARE_COEFS_3(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_3) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL MG5_1_COMPARE_COEFS_3
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )
-      IF(DENOM.GT.0D0)THEN
-        MG5_1_COMPARE_COEFS_3=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        MG5_1_COMPARE_COEFS_3=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION MG5_1_COMPARE_COEFS_4(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_4) COEFS_A, COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL MG5_1_COMPARE_COEFS_4
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        MG5_1_COMPARE_COEFS_4=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        MG5_1_COMPARE_COEFS_4=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION MG5_1_COMPARE_COEFS_5(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_5) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL MG5_1_COMPARE_COEFS_5
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        MG5_1_COMPARE_COEFS_5=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        MG5_1_COMPARE_COEFS_5=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-      FUNCTION MG5_1_COMPARE_COEFS_6(COEFS_A,COEFS_B)
-
-      USE TENS_REC
-      TYPE(COEFF_TYPE_6) COEFS_A,COEFS_B
-      REAL*8 COEF_CHECK_THRS
-      COMMON/MG5_1_COEF_CHECK_THRS/COEF_CHECK_THRS
-      REAL*8 DENOM,NUM
-      LOGICAL MG5_1_COMPARE_COEFS_6
-
-      NUM = ABS(COEFS_A%%C0-COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1-COEFS_B%%C1))
-     $ +SUM(ABS(COEFS_A%%C2-COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3-COEFS_B%%C3))
-     $ +SUM(ABS(COEFS_A%%C4-COEFS_B%%C4))
-      DENOM = ABS(COEFS_A%%C0+COEFS_B%%C0)+SUM(ABS(COEFS_A%%C1+COEFS_B%%C1)
-     $ )+SUM(ABS(COEFS_A%%C2+COEFS_B%%C2))+SUM(ABS(COEFS_A%%C3+COEFS_B%%C3)
-     $ )+SUM(ABS(COEFS_A%%C4+COEFS_B%%C4))
-      IF(DENOM.GT.0D0)THEN
-        MG5_1_COMPARE_COEFS_6=((NUM/DENOM).LT.COEF_CHECK_THRS)
-      ELSE
-        MG5_1_COMPARE_COEFS_6=(NUM.LT.COEF_CHECK_THRS)
-      ENDIF
-
-      END
-
-
-      SUBROUTINE MG5_1_FILL_GOLEM_COEFFS_0(ML_COEFS,GOLEM_COEFS)
-      USE PRECISION_GOLEM, ONLY: KI
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      COMPLEX(KI) GOLEM_COEFS
-      GOLEM_COEFS=ML_COEFS(0)
-      END
-
-      SUBROUTINE MG5_1_FILL_GOLEM_COEFFS_1(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_1
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_1) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-      END
-
-      SUBROUTINE MG5_1_FILL_GOLEM_COEFFS_2(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_2
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_2) GOLEM_COEFS
-C     Constant coefficient 
-      GOLEM_COEFS%%C0=ML_COEFS(0)
-C     Coefficient q(0)
-      GOLEM_COEFS%%C1(1,1)=-ML_COEFS(1)
-C     Coefficient q(0)^2
-      GOLEM_COEFS%%C1(1,2)= ML_COEFS(5)
-C     Coefficient q(1)
-      GOLEM_COEFS%%C1(2,1)=-ML_COEFS(2)
-C     Coefficient q(1)^2
-      GOLEM_COEFS%%C1(2,2)= ML_COEFS(7)
-C     Coefficient q(2)
-      GOLEM_COEFS%%C1(3,1)=-ML_COEFS(3)
-C     Coefficient q(2)^2
-      GOLEM_COEFS%%C1(3,2)= ML_COEFS(10)
-C     Coefficient q(3)
-      GOLEM_COEFS%%C1(4,1)=-ML_COEFS(4)
-C     Coefficient q(3)^2
-      GOLEM_COEFS%%C1(4,2)= ML_COEFS(14)
-C     Coefficient q(0)*q(1)
-      GOLEM_COEFS%%C2(1,1)= ML_COEFS(6)
-C     Coefficient q(0)*q(2)
-      GOLEM_COEFS%%C2(2,1)= ML_COEFS(8)
-C     Coefficient q(0)*q(3)
-      GOLEM_COEFS%%C2(3,1)= ML_COEFS(11)
-C     Coefficient q(1)*q(2)
-      GOLEM_COEFS%%C2(4,1)= ML_COEFS(9)
-C     Coefficient q(1)*q(3)
-      GOLEM_COEFS%%C2(5,1)= ML_COEFS(12)
-C     Coefficient q(2)*q(3)
-      GOLEM_COEFS%%C2(6,1)= ML_COEFS(13)
-      END
-
-      SUBROUTINE MG5_1_FILL_GOLEM_COEFFS_3(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_3
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_3) GOLEM_COEFS
-C     Dummy routine for MG5_1_FILL_GOLEM_COEFS_3
-      STOP 'ERROR: 3 > 2'
-      END
-
-      SUBROUTINE MG5_1_FILL_GOLEM_COEFFS_4(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_4
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_4) GOLEM_COEFS
-C     Dummy routine for MG5_1_FILL_GOLEM_COEFS_4
-      STOP 'ERROR: 4 > 2'
-      END
-
-      SUBROUTINE MG5_1_FILL_GOLEM_COEFFS_5(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_5
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_5) GOLEM_COEFS
-C     Dummy routine for MG5_1_FILL_GOLEM_COEFS_5
-      STOP 'ERROR: 5 > 2'
-      END
-
-      SUBROUTINE MG5_1_FILL_GOLEM_COEFFS_6(ML_COEFS,GOLEM_COEFS)
-      USE TENS_REC, ONLY: COEFF_TYPE_6
-      INCLUDE 'coef_specs.inc'
-      INCLUDE 'loop_max_coefs.inc'
-      COMPLEX*16 ML_COEFS(0:LOOPMAXCOEFS-1)
-      TYPE(COEFF_TYPE_6) GOLEM_COEFS
-C     Dummy routine for MG5_1_FILL_GOLEM_COEFS_6
-      STOP 'ERROR: 6 > 2'
-      END
 
=== modified file 'tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%TIR_interface.f'
--- tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%TIR_interface.f	2016-06-28 23:55:35 +0000
+++ tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%TIR_interface.f	2018-11-09 08:24:41 +0000
@@ -95,17 +95,9 @@
       ENDIF
 
       IF (MLREDUCTIONLIB(I_LIB).EQ.4) THEN
-C       Using Golem95
-C       PDEN is dummy for Golem95 so we just initialize it to zero
-C        here so as to use it for the function SWITCHORDER
-        DO I=0,3
-          DO J=1,NLOOPLINE-1
-            PDEN(I,J)=0.0D0
-          ENDDO
-        ENDDO
-        CALL MG5_1_SWITCH_ORDER(CTMODE,NLOOPLINE,PL,PDEN,M2L)
-        CALL MG5_1_GOLEMLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
-        RETURN
+C       Golem95 not available
+        WRITE(*,*) 'ERROR:: Golem95 is not interfaced.'
+        STOP
       ENDIF
 
 C     INITIALIZE TIR IF NEEDED
 
=== modified file 'tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%loop_matrix.f'
--- tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%loop_matrix.f	2017-08-15 19:34:40 +0000
+++ tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%loop_matrix.f	2018-11-09 08:24:41 +0000
@@ -243,7 +243,7 @@
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING LOOPLIBS ARE
 C      AVAILABLE OR NOT
       LOGICAL LOOPLIBS_AVAILABLE(NLOOPLIB)
-      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.TRUE.,.FALSE.
+      DATA LOOPLIBS_AVAILABLE/.TRUE.,.FALSE.,.TRUE.,.FALSE.,.FALSE.
      $ ,.TRUE.,.TRUE./
       COMMON/MG5_1_LOOPLIBS_AV/ LOOPLIBS_AVAILABLE
 C     A FLAG TO DENOTE WHETHER THE CORRESPONDING DIRECTION TESTS
 
=== modified file 'tests/unit_tests/iolibs/test_export_v4.py'
--- tests/unit_tests/iolibs/test_export_v4.py	2018-06-14 20:20:22 +0000
+++ tests/unit_tests/iolibs/test_export_v4.py	2018-11-09 08:24:41 +0000
@@ -2455,6 +2455,13 @@
     created_files = ['leshouche'
                     ]
 
+
+    def setUp(self):
+        
+        
+        helas_call_writers.HelasCallWriter.customize_argument_for_all_other_helas_object =\
+            staticmethod(helas_call_writers.HelasCallWriter.default_customize_argument_for_all_other_helas_object)
+        
     tearDown = test_file_writers.CheckFileCreate.clean_files
 
     def test_generate_helas_diagrams_ea_ae(self):