@@ -18,7 +18,7 @@ function my_linsolve(A,b,u,p,newA,Pl,Pr,solverdata;verbose=true, kwargs...)
1818end
1919
2020prob = LinearProblem (Diagonal (rand (4 )), rand (4 ))
21- alg = LinearSolveFunction (my_linsolve),
21+ alg = LinearSolveFunction (my_linsolve)
2222sol = solve (prob, alg)
2323```
2424The inputs to the function are as follows:
@@ -30,8 +30,8 @@ The inputs to the function are as follows:
3030- ` Pl ` , left-preconditioner
3131- ` Pr ` , right-preconditioner
3232- ` solverdata ` , solver cache set to ` nothing ` if solver hasn't been initialized
33- - ` kwargs ` , standard SciML keyword arguments such as ` verbose ` , ` maxiters ` ,
34- ` abstol ` , ` reltol `
33+ - ` kwargs ` , standard SciML keyword arguments such as ` verbose ` , ` maxiters ` , ` abstol ` , ` reltol `
34+
3535The function ` my_linsolve ` must accept the above specified arguments, and return
3636the solution, ` u ` . As memory for ` u ` is already allocated, the user may choose
3737to modify ` u ` in place as follows:
4747alg = LinearSolveFunction (my_linsolve!)
4848sol = solve (prob, alg)
4949```
50- Finally, note that ` LinearSolveFunction() ` dispatches to the default linear solve
51- algorithm handling if no arguments are passed in.
52- ``` julia
53- alg = LinearSolveFunction ()
54- sol = solve (prob, alg) # same as solve(prob, nothing)
55- ```
0 commit comments