Just a quick question. Does anyone know if it is possible to supply a seed for the iterative methods in LinearSolve?
I have a large sparse system, and even using some Krylov iterative solver, my estimate is Mathematica should take around 12 hours to find a solution. However, I have an idea of what the solution should look like, and I think I could speed up the process by supplying Mathematica with an accurate initial guess. Does anyone knows if there is an option for that? My current code looks like this:
SparseArray`KrylovLinearSolve[thisLHSOp, thisRHS, Method -> "GMRES", Tolerance -> 10^ (-10), MaxIterations -> Automatic]
Any advice would be appreciated.