using Integrals

# Recommended integration algorithms are:
# For scalar functions: QuadGKJL()
# For ≤ 8 dimensional vector functions: HCubatureJL()
# For > 8 dimensional vector functions: MonteCarloIntegration.vegas(f, st, en, kwargs...)

function Integrals.solve(p::IntegralProblem; kwargs...)
	if length(p.lb) <= 1
		return solve(p, QuadGKJL(); kwargs...)
	elseif length(p.lb) <= 8
		return solve(p, HCubatureJL(); kwargs...)
	else # length(p.lb) > 8
		return solve(p, VEGAS(); kwargs...)
	end
end