"""
    solve_multigrad(problem, gradients)

Solve `problem` for multiple magnetic field gradients.
"""
function solve_multigrad(problem, gradients)
    itertimes = zeros(0)
    results = Vector{output_type(problem)}()
    for grad ∈ gradients
        @info "Solving" typeof(problem) grad
        starttime = time()
        push!(results, solve(problem, grad))
        push!(itertimes, time() - starttime)
    end
    results, itertimes
end
