### Is it possible to solve linear system generated from cbc.block?

56

views

0

Hello all,

I have assembled the matrix and applied boundary condition via cbc.block:

For example, I want to solve it directly using mumps.

Is there any possible way to do it?

Thanks in advance.

Best,

I have assembled the matrix and applied boundary condition via 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
'''
```

Now I want to solve the linear system without the solver provided by cbc.block.For example, I want to solve it directly using mumps.

Is there any possible way to do it?

Thanks in advance.

Best,

Community: FEniCS Project

Please login to add an answer/comment or follow this question.

https://www.allanswered.com/post/bevkr/

and

https://github.com/MiroK/cbc.block-eig/blob/master/stokes.py

It seems that I have to first transferred my cbc.block.mat to PETSc.mat and then solve the system?