Re: [deal.II] sparse matrix multiplication

2017-12-04 Thread Wolfgang Bangerth
On 12/04/2017 12:43 PM, Wolfgang Bangerth wrote: Can I figure out the pattern of B*B or BA^{-1}B ? In general, A^{-1} is going to be a dense matrix even if A is a sparse matrix. So building this matrix is not efficient. Of course, you say that A is diagonal, in which case this doesn't

Re: [deal.II] sparse matrix multiplication

2017-12-04 Thread Wolfgang Bangerth
On 12/04/2017 12:04 PM, Peimeng Yin wrote: I want to solve a linear system, which has the following form (A+B A^{-1} B) u = f Here A is a diagonal black matrix, B is a block matrix that I can figure out the pattern. Are there any solver that I can use to solve the linear system? Yes. Have

[deal.II] sparse matrix multiplication

2017-12-04 Thread Peimeng Yin
I want to solve a linear system, which has the following form (A+B A^{-1} B) u = f Here A is a diagonal black matrix, B is a block matrix that I can figure out the pattern. Are there any solver that I can use to solve the linear system? or Can I figure out the pattern of B*B or BA^{-1}B ?