I am trying to solve a sparse stiffness matrix (positive definite) of size 1986360x1986360 SparseMatrixCSC{Float64, Int64} with 18327153 stored entries. Here’s a sample code. using SparseArrays using LinearAlgebra A = sprandn(1986360, 1986360, 2e-6).*20; A = A+A'+100I b = zeros(size(a,1)); b[[1:…

How to solve this Ax=b faster?

Jae-Mo_Lihm May 26, 2022, 8:19am 2

Have you tried iterative solvers?

There is Krylov.jl and many others.

4 Likes

Topic		Replies	Views
Linear solver \(A, B) performance vs Matlab A\b General Usage	32	7842	May 21, 2017
From toy to production code with sparse matrix and sparse-direct factorizations Performance sparse , cholmod	17	1234	November 1, 2022
Sparse Cholesky of Gram Matrix (SuiteSparse?) Numerics question , linearalgebra	25	2717	April 11, 2020
Solving Sparse Linear Systems fast Performance sparse , linearsolve	11	5375	June 23, 2022
Cannot solve A\b = x with sparse matrix A - Julia V 0.6.2 Specific Domains	11	3141	May 20, 2018