handle no coverage inputs in coverage dependence

Author: mjohnson541

src/Calculators/RateCoverageDependence.jl

@@ -28,13 +28,19 @@ ms corresponds to m (float)
     indms::Dict{N,L} = Dict()
 end
 @inline function getcovdepactivationbarriercorrection(covdep::PolynomialRateCoverageDependence,T,coverages)
+    if length(coverages) == 0
+        return zero(typeof(coverages).parameters[1])
+    end
     E = 0.0
     for (ind,Epoly) in covdep.indEpolys
         E += evalpoly(coverages[ind],indEpoly)
     end
     return E
 end
 @inline function getcovdepfactorcorrection(covdep::PolynomialRateCoverageDependence,T,coverages)
+    if length(coverages) == 0
+        return zero(typeof(coverages).parameters[1])
+    end
     av = 0.0
     for (ind,a) in covdep.indavals
        av += coverages[ind]*a

src/Calculators/ThermoCoverageDependence.jl

@@ -23,6 +23,9 @@ end
     return PolynomialThermoEnergyCoverageDependence(polys,Dict{Int64,Vector{Float64}}())
 end
 @inline function getcovdepenthalpycorrection(covdep::PolynomialThermoEnergyCoverageDependence,T,coverages)
+    if length(coverages) == 0
+        return zero(typeof(coverages).parameters[1])
+    end
     E = zero(coverages[1])
     for (ind,poly) in covdep.indpolys
         E += evalpoly(coverages[ind],poly)