
	if( sqrtm .lt. 8.94427190999916D0 )
     &    Warning('Extrapolating bbhSM('//Digit(h)//') in MHiggs')
	if( sqrts .lt. 7.D0 )
     &    Warning('Extrapolating bbhSM('//Digit(h)//') in sqrts')

        bbhSM(h) = exp(9.828582156488078D0 + 
     &     sqrtm*(-0.7744489491850788D0 + 
     &        0.004979441798667225D0*sqrtm + 
     &        0.010661983992092014D0*sqrts) + 
     &     sqrts*(0.4972580586805361D0 - 0.02099276251143286D0*sqrts))

#ifdef DETAILED_DEBUG
	DPROD "bbhSM(h) =", bbhSM(h) ENDL
#endif


	if( sqrtm .lt. 10.D0 )
     &    Warning('Extrapolating gghSM('//Digit(h)//') in MHiggs')
	if( logsqrts .lt. 0.6931471805599453D0 )
     &    Warning('Extrapolating gghSM('//Digit(h)//') in logsqrts')

	if( sqrtm.lt.18.61719635176038D0 ) then

        gghSM(h) = exp(7.839021207947174D0 + 
     &     sqrtm*(0.4597789726533755D0 + 
     &        sqrtm*(-0.09520410304013488D0 + 
     &           0.002595171851553549D0*sqrtm)) + 
     &     logsqrts*(1.0470241902681292D0 + 
     &        sqrtm*(0.2178045072534087D0 + 
     &           0.0014030210301997678D0*sqrtm) + 
     &        logsqrts*(-0.32569456924069123D0 + 
     &           0.10441152853302572D0*logsqrts - 
     &           0.04061951581618903D0*sqrtm)))

	else if( sqrtm.lt.28.284271247461902D0 ) then

        gghSM(h) = exp(-817.3651235357215D0 + 
     &     sqrtm*(167.82468547305393D0 + 
     &        sqrtm*(-13.650814245866401D0 + 
     &           sqrtm*(0.5536166109007861D0 + 
     &              sqrtm*(-0.011282053936726252D0 + 
     &                 0.00009266949162361316D0*sqrtm)))) + 
     &     logsqrts*(10.399254523856342D0 + 
     &        sqrtm*(-1.2219996437212877D0 + 
     &           sqrtm*(0.044525398408260586D0 + 
     &              sqrtm*(0.0015423282033737537D0 - 
     &                 0.00006062990881567466D0*sqrtm))) + 
     &        logsqrts*(-0.8392009160482815D0 + 
     &           sqrtm*(0.48169526275775426D0 + 
     &              sqrtm*(-0.04350564413455493D0 + 
     &                 0.0009781671647075918D0*sqrtm)) + 
     &           logsqrts*(-1.683260726356996D0 + 
     &              sqrtm*(0.14096784028684323D0 - 
     &                 0.0041717583993185064D0*sqrtm) + 
     &              logsqrts*
     &               (0.10338152541874138D0 - 
     &                 0.03780175829659603D0*logsqrts + 
     &                 0.008545416556998925D0*sqrtm)))))

	else

        gghSM(h) = exp(16.04925833582041D0 + 
     &     sqrtm*(-0.9108576388708007D0 + 
     &        0.003469185796746001D0*sqrtm) + 
     &     logsqrts*(3.995988922886652D0 - 
     &        0.9484622101428902D0*logsqrts + 
     &        0.12423688318064019D0*sqrtm))

	endif

#ifdef DETAILED_DEBUG
	DPROD "gghSM(h) =", gghSM(h) ENDL
#endif


	if( sqrtm .lt. 10.D0 )
     &    Warning('Extrapolating gghSMNLO('//Digit(h)//') in MHiggs')
	if( logsqrts .lt. 0.6931471805599453D0 )
     &    Warning('Extrapolating gghSMNLO('//Digit(h)//') in logsqrts')

	if( sqrtm.lt.18.61719635176038D0 ) then

        gghSMNLO(h) = exp(-453.9834306013304D0 + 
     &     sqrtm*(247.98469346922593D0 + 
     &        sqrtm*(-56.4198826236009D0 + 
     &           sqrtm*(7.0622502142848D0 + 
     &              sqrtm*(-0.5270068485931979D0 + 
     &                 sqrtm*
     &                  (0.023451257710193404D0 + 
     &                    sqrtm*
     &                     (-0.000576331165235695D0 + 
     &                       6.035687694001616D-6*sqrtm)))))) + 
     &     logsqrts*(-2.1352058896077684D0 + 
     &        sqrtm*(1.2485372212218022D0 + 
     &           sqrtm*(-0.19925936999165453D0 + 
     &              sqrtm*(0.02170451792690496D0 + 
     &                 sqrtm*
     &                  (-0.0012633605934095481D0 + 
     &                    sqrtm*
     &                     (0.000039277268531300396D0 - 
     &                       4.869386167860268D-7*sqrtm))))) + 
     &        logsqrts*(1.214796665301791D0 + 
     &           sqrtm*(0.005295108193192817D0 + 
     &              sqrtm*(-0.021777401752199714D0 + 
     &                 sqrtm*
     &                  (0.0013101589789698913D0 + 
     &                    sqrtm*
     &                     (-0.0000546603602490181D0 + 
     &                       4.4027583992401005D-7*sqrtm)))) + 
     &           logsqrts*(-1.4034574401087694D0 + 
     &              sqrtm*(0.08010651581071794D0 + 
     &                 sqrtm*
     &                  (0.0005776828983918536D0 + 
     &                    sqrtm*
     &                     (0.0000736749930537802D0 + 
     &                       4.015883540659236D-6*sqrtm))) + 
     &              logsqrts*
     &               (0.7134266598125986D0 + 
     &                 sqrtm*
     &                  (-0.032494833017523476D0 + 
     &                    sqrtm*
     &                     (-0.00038203133194332695D0 - 
     &                       0.00003820078563333181D0*sqrtm)) + 
     &                 logsqrts*
     &                  (-0.20811295823534173D0 + 
     &                    sqrtm*
     &                     (0.008041196780004973D0 + 
     &                       0.00019265588057849022D0*sqrtm) + 
     &                    logsqrts*
     &                     (0.030907792238600698D0 - 
     &                       0.001380462041748478D0*logsqrts - 
     &                       0.0011677233980619324D0*sqrtm)))))))

	else if( sqrtm.lt.28.284271247461902D0 ) then

        gghSMNLO(h) = exp(-19052.774939750147D0 + 
     &     sqrtm*(5593.608871380419D0 + 
     &        sqrtm*(-702.4771163574136D0 + 
     &           sqrtm*(48.9262225205997D0 + 
     &              sqrtm*(-2.0408982232951787D0 + 
     &                 sqrtm*
     &                  (0.05098005090602945D0 + 
     &                    sqrtm*
     &                     (-0.0007060198082261215D0 + 
     &                       4.181410399005631D-6*sqrtm)))))) + 
     &     logsqrts*(28.39756542570317D0 + 
     &        sqrtm*(-6.565102919750903D0 + 
     &           sqrtm*(0.7141153997488958D0 + 
     &              sqrtm*(-0.04036930356519674D0 + 
     &                 sqrtm*
     &                  (0.0013133613144338657D0 + 
     &                    sqrtm*
     &                     (-0.000021914025250901773D0 + 
     &                       1.6086790622741766D-7*sqrtm))))) + 
     &        logsqrts*(0.31632243258938314D0 + 
     &           sqrtm*(-0.4304111664312427D0 + 
     &              sqrtm*(0.04447372552299553D0 + 
     &                 sqrtm*
     &                  (-0.0022469646685933666D0 + 
     &                    sqrtm*
     &                     (0.000024969177638614434D0 - 
     &                       5.077183964679362D-7*sqrtm)))) + 
     &           logsqrts*(1.6199222053572735D0 + 
     &              sqrtm*(-0.12361988560896664D0 + 
     &                 sqrtm*
     &                  (0.00248946092870488D0 + 
     &                    sqrtm*
     &                     (0.0003245909873292568D0 + 
     &                       5.682560456962786D-6*sqrtm))) + 
     &              logsqrts*
     &               (-0.5398090279863337D0 + 
     &                 sqrtm*
     &                  (0.048737163729617825D0 + 
     &                    sqrtm*
     &                     (-0.003113497540786245D0 - 
     &                       0.00010201610490014155D0*sqrtm)) + 
     &                 logsqrts*
     &                  (0.05407903078325591D0 + 
     &                    sqrtm*
     &                     (0.0016326917197586107D0 + 
     &                       0.0009704919615231642D0*sqrtm) + 
     &                    logsqrts*
     &                     (-0.009114819861799376D0 + 
     &                       0.005769904537050648D0*logsqrts - 
     &                       0.0035247968760733226D0*sqrtm)))))))

	else

        gghSMNLO(h) = exp(7.073360187841716D0 + 
     &     sqrtm*(0.4001551467525027D0 + 
     &        sqrtm*(-0.05792447603512322D0 + 
     &           sqrtm*(0.0010870212999651242D0 - 
     &              0.000017010828214181483D0*sqrtm))) + 
     &     logsqrts*(0.6467820374631277D0 + 
     &        sqrtm*(0.0865334744075569D0 + 
     &           sqrtm*(0.010350263320617323D0 + 
     &              0.0004393527104899614D0*sqrtm)) + 
     &        logsqrts*(1.4798454680050523D0 + 
     &           sqrtm*(-0.12153731710854525D0 - 
     &              0.011434855281518611D0*sqrtm) + 
     &           logsqrts*(-0.14998519434441773D0 - 
     &              0.39038482059871255D0*logsqrts + 
     &              0.12122387264875743D0*sqrtm))))

	endif

#ifdef DETAILED_DEBUG
	DPROD "gghSMNLO(h) =", gghSMNLO(h) ENDL
#endif

