3 ms·
The title is incorrect. It's not a general partial differential equation solver. The paper doesn't even claim it to be. The issue is that not all, not even most
by ChrisRackauckas 8y ago
The title is incorrect. It's not a general partial differential equation solver. The paper doesn't even claim it to be. The issue is that not all, not even most, PDEs can be discretized to a linear system of equations. Normally these need to be discretized to either a nonlinear rootfinding problem G(x)=0 find x, or a stiff ODE (normally solved by an implicit method, but not necessarily).
However, both of those problems do include a subproblem which is solving linear systems (when the ODE is solved implicitly), and so PDE solvers which take these routes will generally be able to be accelerated by an improved linear equation solver. Thus it's not an be all and end all algorithm, but when combined with the tools that exist for PDEs in something like PETSc or Julia then this can be used to give a more efficient algorithm for most PDEs (though there are still other PDEs where linear solving is entirely avoided, for example with exponential integrators or Runge-Kutta Chebyshev methods, which can be more efficient depending on the problem).
- tanderson92 8y agoThe link was changed from the original Nature article which has the original title. The new link is unrelated. You're correct that the innovation is really about sparse linear algebra (though not just solvers) done in hardware. However the advantages are exploited not just by way of solving linear systems for implicit time-stepping. They implemented a Jacobi iteration in hardware for solving a system arising from discretization of the Poisson equation, but they also used an efficient forward-map of the coefficient matrix to do explicit time-stepping for a wave solver. So they solve both a time-independent and time-dependent problem, but neither with implicit time-stepping (though there appears to be no technical limitation here, it's just what they focused on).