Add a workspace struct

[vchuravy]

• Add NewtonKrylovWorkspace to pre-allocate all Newton-Krylov buffers
• Add evaluate! and refactor linesearch to operate on workspace

[github-actions]

Once the build has completed, you can preview your PR at this URL: https://NumericalMathematics.github.io/Ariadne.jl/previews/PR130/ in a couple of minutes.

[ranocha]

Thanks!

[github-actions]

Your PR requires formatting changes to meet the project's style guidelines.
Please consider running Runic (git runic main) to apply these changes.

Click here to view the suggested changes.

diff --git a/src/Ariadne.jl b/src/Ariadne.jl
index 14628d2..40d9dbd 100644
--- a/src/Ariadne.jl
+++ b/src/Ariadne.jl
@@ -338,7 +338,7 @@ struct NewtonKrylovWorkspace{F, A, P, JOp <: AbstractJacobianOperator, KW}
 end
 
 function NewtonKrylovWorkspace(
-        F!, u::AbstractArray, p, res::AbstractArray, ::Val{Algo}=Val(:gmres);
+        F!, u::AbstractArray, p, res::AbstractArray, ::Val{Algo} = Val(:gmres);
         assume_p_const::Bool = false
     ) where {Algo}
     # res .= 0 might ignore ghost cells
diff --git a/test/workspace.jl b/test/workspace.jl
index 3ade878..b592a75 100644
--- a/test/workspace.jl
+++ b/test/workspace.jl
@@ -96,5 +96,5 @@ Base.zero(v::WrappedVector) = WrappedVector(zero(v.data))
     res = zero(x₀)
     x, stats = newton_krylov!(F!, x₀, nothing, res)
     @test stats.solved
-    @test x ≈ [1.0, 1.0] atol = 1e-5
+    @test x ≈ [1.0, 1.0] atol = 1.0e-5
 end