Dear sir,
I am solving the 3D problem with 60000 elements and the degrees of freedom are 2 at each node. And, the global stiffness matrix is large but sparse in our case. please suggest me fast and accurate linear solver for sparse matrix. OR GPU(paralell ) based solver.
Why not use KSP: Linear System Solvers — PETSc v3.17.2-605-g145e647609 documentation, which are the default solver backend of DOLFIN?
Scipy, as far as Im aware, does not do parallel distributed matrices, see for instance: Parallelize Scipy iterative methods for linear equation systems(bicgstab) in Python - Computational Science Stack Exchange