Il giorno gio 10 apr 2025 alle ore 11:25 Peder Jørgensgaard Olesen via
petsc-users ha scritto:
> I would like to output the contents of a Vec to an ASCII file in which the
> entries are all on one line, as "x[0] x[1] x[2] ... x[N-1]". This can
> be done in a slightly roundabout way by putting
Hi, Steven
Thank you for the bug report and test example. You were right. The
MatCopy(A, B,..) implementation was wrong when B was on the device.
I have a fix at
https://urldefense.us/v3/__https://gitlab.com/petsc/petsc/-/merge_requests/8288__;!!G_uCfscf7eWS!e15LlybplgvD-3BTmMELbBy9EiU1biKGxn
Sounds correct.
> On Apr 10, 2025, at 2:04 PM, Anton Popov wrote:
>
>
> Hi guys,
> I have a custom multigrid preconditioner that uses Galerkin coarsening for
> the staggered grid finite difference. Now I want to replace a few top level
> operators with the matrix-free shell matrices.
>
Hi guys,
I have a custom multigrid preconditioner that uses Galerkin coarsening
for the staggered grid finite difference. Now I want to replace a few
top level operators with the matrix-free shell matrices.
Here is my plan:
1. Set PC_MG_GALERKIN_NONE
2. Call PCMGSetRestriction and PCMGSetIn
Hi Peder,
You can run with '-help -vec_view' (or whatever your view option is) and
that prints all options that are queried for.
This is very noisy so grep on 'vec_view' and you will see the option
information: what it was, what it is after the query, and what options it
accepts.
I don't know of
I would like to output the contents of a Vec to an ASCII file in which the
entries are all on one line, as "x[0] x[1] x[2] ... x[N-1]". This can be
done in a slightly roundabout way by putting the values in a 1xN dense Mat,
assembling the matrix, and then use MatView with a suitable format,
