Solving the sparse system is returning a Matrix while converting to dense and then solving returns a Vector as expected. Is this the expected behavior? Or am I doing something incorrectly? using SparseArrays # Constructing sparse array function hess_sparse(x) return [-sin(x[1] + x[2]) + 1, -s…

You are artificially sparsifying the rhs. Your function grad_sparse already returns the rhs as is. You could then even save the additional memory allocation due to gradinds, if you left the rhs dense. Also, you would save the allocation of the index vector for the result, and, if you don’t need grad…

Solving sparse linear system returns a Matrix

Specific Domains Numerics

stevengj January 5, 2023, 2:05pm 11

See also Sparse solve with sparse rhs

1 Like

Topic		Replies	Views
Solving sparse, linear equations? Numerics linearalgebra , sparse	10	1561	November 15, 2022
Cannot solve A\b = x with sparse matrix A - Julia V 0.6.2 Specific Domains	11	3103	May 20, 2018
Missing sparse methods Numerics linearalgebra	4	2277	May 24, 2022
Julia cannot solve Ax=b with LU where A is sparse, General Usage question , linearalgebra , sparse	12	1167	May 5, 2023
Solving sparse matrix using backslash, got MethodError: no method matching lu! Numerics linearalgebra , sparse	3	1036	July 11, 2020

Solving sparse linear system returns a Matrix

Related topics