diff --git a/ext/IntegralsArblibExt.jl b/ext/IntegralsArblibExt.jl
index 6e9d5ed..0eac48a 100644
--- a/ext/IntegralsArblibExt.jl
+++ b/ext/IntegralsArblibExt.jl
@@ -21,16 +21,13 @@ function Integrals.__solvebp_call(
         val = Arblib.integrate!(f_, res, lb, ub, atol = abstol, rtol = reltol,
             check_analytic = alg.check_analytic, take_prec = alg.take_prec,
             warn_on_no_convergence = alg.warn_on_no_convergence, opts = alg.opts)
-        SciMLBase.build_solution(
-            prob, alg, val, get_radius(val), retcode = ReturnCode.Success)
     else
         f_ = (x; kws...) -> only(prob.f(x, p; kws...))
         val = Arblib.integrate(f_, lb, ub, atol = abstol, rtol = reltol,
             check_analytic = alg.check_analytic, take_prec = alg.take_prec,
             warn_on_no_convergence = alg.warn_on_no_convergence, opts = alg.opts)
-        SciMLBase.build_solution(
-            prob, alg, val, get_radius(val), retcode = ReturnCode.Success)
     end
+    SciMLBase.build_solution(prob, alg, val, get_radius(val), retcode = ReturnCode.Success)
 end
 
 function get_radius(ball)
diff --git a/src/Integrals.jl b/src/Integrals.jl
index 844eaf1..1627304 100644
--- a/src/Integrals.jl
+++ b/src/Integrals.jl
@@ -156,7 +156,6 @@ function __solvebp_call(cache::IntegralCache, alg::QuadGKJL, sensealg, domain, p
             end
             val, err = quadgk(_f, lb, ub, segbuf = cache.cacheval, maxevals = maxiters,
                 rtol = reltol, atol = abstol, order = alg.order, norm = alg.norm)
-            SciMLBase.build_solution(prob, alg, val, err, retcode = ReturnCode.Success)
         else
             prototype = f(typeof(mid)[], p)
             _f = if prototype isa AbstractVector
@@ -170,7 +169,6 @@ function __solvebp_call(cache::IntegralCache, alg::QuadGKJL, sensealg, domain, p
             end
             val, err = quadgk(_f, lb, ub, segbuf = cache.cacheval, maxevals = maxiters,
                 rtol = reltol, atol = abstol, order = alg.order, norm = alg.norm)
-            SciMLBase.build_solution(prob, alg, val, err, retcode = ReturnCode.Success)
         end
     else
         if isinplace(f)
@@ -179,14 +177,13 @@ function __solvebp_call(cache::IntegralCache, alg::QuadGKJL, sensealg, domain, p
             val, err = quadgk!(_f, result, lb, ub, segbuf = cache.cacheval,
                 maxevals = maxiters,
                 rtol = reltol, atol = abstol, order = alg.order, norm = alg.norm)
-            SciMLBase.build_solution(prob, alg, val, err, retcode = ReturnCode.Success)
         else
             _f = u -> f(u, p)
             val, err = quadgk(_f, lb, ub, segbuf = cache.cacheval, maxevals = maxiters,
                 rtol = reltol, atol = abstol, order = alg.order, norm = alg.norm)
-            SciMLBase.build_solution(prob, alg, val, err, retcode = ReturnCode.Success)
         end
     end
+    SciMLBase.build_solution(prob, alg, val, err, retcode = ReturnCode.Success)
 end
 
 function init_cacheval(alg::HCubatureJL, prob::IntegralProblem)