Hello all,
It souds that the best way to introduce petsc in a code is not to introduce
it, but develop the code over the petsc structure.
It is probably true but my problem is that my existing code already is
equipped with a domain decomposition based on MPI (a typical themal
hydraulic with carte
Without having seen your code, it sounds to me like the best strategy here
is to:
1) Produce a mirror of your mesh using DMStag
2) Use that DM to construct the linear block problems, which can then be
solved by PETSc
Since the PETSc grid matches your own, you can share the solution
vec
On Sun, Oct 13, 2019 at 5:25 AM Pierre Gubernatis via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hello all,
>
> It souds that the best way to introduce petsc in a code is not to
> introduce it, but develop the code over the petsc structure.
>
All things being equal, yes, but few users start
> On Oct 13, 2019, at 4:24 AM, Pierre Gubernatis via petsc-users
> wrote:
>
> Hello all,
>
> It souds that the best way to introduce petsc in a code is not to introduce
> it, but develop the code over the petsc structure.
>
> It is probably true but my problem is that my existing code alrea
omain in sub-domains and construct a linear problem
> by block: each sub-domain assembles its part of the operator and its part
> of the RHS.
>
> I am wondering what is the best way now to introduce petsc (considering
> that I don?t want to assemble a global operator on a given pro
Thanks. It seems that the iterative method is solved the entire linear system
very accurately. For the last linear solve the initial true residual norm is
around 10^6 (i.e. the norm of b) and the final one around 10^-10
The initial true residual norm on each block (part of a slab on one p