************************************************************************ * * initIntegralspDIS.f: * * This routine initializes the integrals of coefficient functions * and interpolation functions. * ************************************************************************ subroutine initIntegralspDIS * implicit none * include "../commons/grid.h" include "../commons/DynScVar.h" ** * Internal Variables * integer alpha,beta * * Check that the number of intervals does not exceed the maximum number * if(nin(igrid)+inter_degree(igrid).gt.nint_max_DIS)then write(6,*) "In initIntegralsSIA.f:" write(6,*) "Number of grid points too large:" write(6,*) "Maximum value allowed =",nint_max_DIS write(6,*) "You should reduce it." write(6,*) " " call exit(-10) endif * * Initialize integrals * if(IsExt(igrid))then * * If this is an external grid, compute the integrals for * the entire splitting matrix ... * do alpha=0,nin(igrid)-1 do beta=alpha,nin(igrid)-1 call RSLintegralspDIS(alpha,beta) enddo enddo else * * ... otherwise only for the first line * do alpha=0,nin(igrid)-1 call RSLintegralspDIS(0,alpha) enddo endif * * Scale variations * (only if the dynamical scale variation is disabled) * if(.not.DynScVar) call IncludeScaleVariation * return end