Sir.
I want to using BDDC solve a linear system. I already have a matrix A and
vector b. I want to solve Ax=b.
I have no idea how BDDC works. I just want to experience the performance of
BDDC preconditioner.
I wonder does the BDDC precondtioner can be used as a black box as other
Sir, I rewrite the code.
The code do not include any domain concept.
The code is simple allocate space, insert value, compress, set attribute, then
convert the matrix type.
It can help to reproduce the problem.
Please take a look.
As you suggested, I found that the problem comes from the line
```
Sir, I already dump the matrix into binary file and attached in the mail, and
please have a look.
Thanks for your time.
The follow code can be used to reproduce problem.
---
std::string path = "..";
PetscViewer viewer;
Sir, I can not understand why my convert from MPIAIJ to MATIS failed.
I already have a matrix A, create by follow commands
-
Mat A;
MatCreateAIJ(...)
MatSetValues(...)
MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY);
MatAssemblyEnd(A,
Thanks sir.
I now modify my code into following.
Everything works good.
-
// stage 1:
Vec x1, b2;
Vec x1, b2;
Mat A, P, F;
PC pc;
// solve first system
MatCreateAIJ(A, ...)
MatSetVaules(A, ...)
MatAssembleBegin(A, ...)
MatAssembleBegin(A, ...)
Sir, I have a problem when using petsc.
I want to solve a series of linear equations.
A1*x1=b1, A2*x2=b2, A3*x3=b3 ...
The A1,A2,A3 have the same sparstiy pattern.
I want to use MUMPS to solve the system.
In order to enhance performance, I want to reuse the symbolic factorization.
Here