- No, MATLAB Coder doesn't yet support sparse. http://www.mathworks.com/help/releases/R2015a/coder/ug/matlab-language-features-supported-for-code-generation.html
- Perhaps gmres
Backslash for solving very big sparse matrix in Ax=b
5 次查看(过去 30 天)
显示 更早的评论
In my code at each time step sparse matrix A is built as follows:
1< DX <120000
A = sparse(i,j,s,DX,DX,6*DX)
b = (1,DX)
The problem that I am dealing is a sort of discretization problem. I have maximum 120000 nodes. Each of these nodes are having special charachters and I choose only the ones that meet my criteria. the number of these chosen ones is DX and is completely dependent on the physical process.
I am using backslash in x = A\b. But as the size of A could become quite big, the computational time rises drastically (more than 10e5 time steps are having DX > 6e4). As far as I know, backslash operation is already well optimized in MATLAB but I would like to know:
1. Would it make sense to use codegen and convert the code to C?
2. Does any one know an alternative method instead of backslash, so that the compuattional time decreases (maybe an iterative method?)?
0 个评论
回答(3 个)
Sean de Wolski
2015-7-16
0 个评论
John D'Errico
2015-7-16
NO. The solver in backslash is ALREADY compiled and highly optimized. Why do you think that converting the call to C would help? NOT.
You can try an iterative method. They are designed for problems that are too large to solve in memory. Don't be surprised if you need to learn about preconditioning matrices.
Royi Avital
2015-8-9
You can use it without a preconditioner to have classic Conjugate Gradient.
The beauty of this method is you never multiply the whole matrix by the vector. You only use do products as needed.
Its speeds depends on the condition number of the Matrix, hence with a good preconditioner things gets much better.
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!