Unreachable code

[schillic]

The current implementation (see below) contains code that should be unreachable: If Re_λ₁ = 0, then R = Inf and the check case should have triggered before reaching the case iszero(Re_λ₁). The only way to reach there is by deactivating check, but then the result is not valid anyway as far as I understand.

CarlemanLinearization.jl/src/error_bounds.jl

Lines 85 to 99 in 39b92e8

	function error_bound_specabs(x₀, F₁, F₂; N, check=true)
	(R, Re_λ₁) = _error_bound_specabs_R(x₀, F₁, F₂; check=check)
	if check
	@assert R < 1 "expected R < 1, got R = $R; try scaling the ODE"
	end
	nx₀ = norm(x₀, 2)
	if iszero(Re_λ₁)
	nF₂ = opnorm(F₂, 2)
	ε = t -> nx₀ * (nx₀ * nF₂ * t)^N
	else
	ε = t -> nx₀ * R^N * (1 - exp(Re_λ₁ * t))^N
	end
	return ε
	end