@@ -18,7 +18,7 @@ function my_linsolve(A,b,u,p,newA,Pl,Pr,solverdata;verbose=true, kwargs...)
18
18
end
19
19
20
20
prob = LinearProblem (Diagonal (rand (4 )), rand (4 ))
21
- alg = LinearSolveFunction (my_linsolve),
21
+ alg = LinearSolveFunction (my_linsolve)
22
22
sol = solve (prob, alg)
23
23
```
24
24
The inputs to the function are as follows:
@@ -30,8 +30,8 @@ The inputs to the function are as follows:
30
30
- ` Pl ` , left-preconditioner
31
31
- ` Pr ` , right-preconditioner
32
32
- ` 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
+
35
35
The function ` my_linsolve ` must accept the above specified arguments, and return
36
36
the solution, ` u ` . As memory for ` u ` is already allocated, the user may choose
37
37
to modify ` u ` in place as follows:
47
47
alg = LinearSolveFunction (my_linsolve!)
48
48
sol = solve (prob, alg)
49
49
```
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