I see, thank you!
Yuyun
From: Matthew Knepley
Sent: Friday, March 8, 2019 4:51 PM
To: Yuyun Yang
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] Append vector to existing file
On Fri, Mar 8, 2019 at 5:53 PM Yuyun Yang via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Hello team,
Hello team,
This is a very simple question, but I just want to be sure I understand how the
viewer works. I have a file with some vectors already written into it. Now I'm
calling PetscViewerASCIIOpen, setting the format and then the file mode into
append:
PetscViewer viewer;
ierr =
Just seeing this now. It is hard to imagine how bad GAMG could be on a
coarse grid, but you can run with -info and grep on GAMG and send that. You
will see listing of levels, number of equations and number of non-zeros
(nnz). You can send that and I can get some sense of GAMG is going nuts.
Mark
MatCreate is collective so you want to check that it is not seeing load
imbalance from earlier code.
And duplicating communicators can be expensive on some systems.
On Fri, Mar 8, 2019 at 10:21 AM Ale Foggia via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hello all,
>
> I have a problem
https://www.mcs.anl.gov/petsc/documentation/faq.html#efficient-assembly
> On Mar 8, 2019, at 9:19 AM, Ale Foggia via petsc-users
> wrote:
>
> Hello all,
>
> I have a problem with the scaling of the MatCreate() function. I wrote a code
> to diagonalize sparse matrices and I'm running it in
It may not address the memory issue, but can you build 3.10 with the
same options you used for 3.6? It is currently a debugging build:
##
##
#
Dear petsc4py experts,
I'd like to ask why several PETSc functions are not wrapped in petsc4py. I'd
need to use DMPlexCreateFromDAG from python. Could you explain with this
function as an example why there is no python wrapper available? Do I have to
expect severe difficulties when I try this
This is very unusual. MatCreate() does no work, merely dup'ing a
communicator (or referencing an inner communicator if this is not the
first PetscObject on the provided communicator). What size matrices are
you working with? Can you send some performance data and (if feasible)
a reproducer?
Dear Matt,
Thank you for your answer but I get an error even before compiling ex6. The
error is:
error: too many arguments to function ‘DMSetField’
ierr = DMSetField(dm, 0, NULL, (PetscObject) fe);CHKERRQ(ierr);
Eda
Matthew Knepley , 7 Mar 2019 Per, 15:07 tarihinde şunu
yazdı:
> On Thu,