Is it possible to solve linear system generated from cbc.block?
19 days ago by
I have assembled the matrix and applied boundary condition via cbc.block:
Now I want to solve the linear system without the solver provided by cbc.block.
A = block_assemble([[a11, a12] [a21, a22]], bcs=bcs) b = block_assemble([L1, 0], bcs=bcs) ''' I want to solve: A(x,y)^T = b '''
For example, I want to solve it directly using mumps.
Is there any possible way to do it?
Thanks in advance.
Community: FEniCS Project
Please login to add an answer/comment or follow this question.