Hi Barry,
I've followed your suggestion to use PetscMallocDump after the call to KSPSolve
and I've looked at the output.
What I find is that for a smaller sized model that I can run on my laptop, I
don't see any increase in the memory after it is allocated the first pass
through.
However,
Randy,
Just email me the PetscMallocDump output file directly.
Barry
On Mar 6, 2015, at 12:10 PM, Randall Mackie rlmackie...@gmail.com wrote:
Hi Barry,
I've followed your suggestion to use PetscMallocDump after the call to
KSPSolve and I've looked at the output.
What I find
Barry Smith bsm...@mcs.anl.gov writes:
This should not happen, probably you are working with a MPIAIJ
matrix
Randall, it sounds like you are creating an MPIAIJ matrix even in
serial. Meanwhile, GAMG creates a SeqAIJ prolongation operator (when
running in serial). Why are you creating an
Please keep the list Cc'd.
Randall Mackie rlmackie...@gmail.com writes:
On Mar 5, 2015, at 11:17 AM, Jed Brown j...@jedbrown.org wrote:
Barry Smith bsm...@mcs.anl.gov writes:
This should not happen, probably you are working with a MPIAIJ
matrix
Randall, it sounds like you are
I have fixed the bug of PCGAMG crashing on one process when used with MPIAIJ
matrix in the branch barry/clean-gamg and merged to next for testing.
As Jed notes you can create your matrix as AIJ (which is recommended
approach anyways) instead of MPIAIJ to work around the bug.
Barry
Randy,
I've not been able to reproduce this; let us know if get to the point of
having something we can run and debug.
Barry
On Mar 4, 2015, at 7:45 PM, Randall Mackie rlmackie...@gmail.com wrote:
In my application, I am repeatedly calling KSPSolve with the following
options:
On Mar 4, 2015, at 7:30 PM, Barry Smith bsm...@mcs.anl.gov wrote:
On Mar 4, 2015, at 7:45 PM, Randall Mackie rlmackie...@gmail.com wrote:
In my application, I am repeatedly calling KSPSolve with the following
options:
-ksp_type gmres \
-pc_type gamg \
-pc_gamg_type agg \
This should not happen, probably you are working with a MPIAIJ matrix and
GAMG incorrectly generates a SeqAIJ matrix internally and tries to use them
together. I'll try to reproduce
Barry
On Mar 5, 2015, at 12:07 PM, Randall Mackie rlmackie...@gmail.com wrote:
On Mar 4, 2015, at
On Mar 4, 2015, at 7:45 PM, Randall Mackie rlmackie...@gmail.com wrote:
In my application, I am repeatedly calling KSPSolve with the following
options:
-ksp_type gmres \
-pc_type gamg \
-pc_gamg_type agg \
-pc_gamg_agg_nsmooths 1\
each call is after the matrix and right hand side
In my application, I am repeatedly calling KSPSolve with the following options:
-ksp_type gmres \
-pc_type gamg \
-pc_gamg_type agg \
-pc_gamg_agg_nsmooths 1\
each call is after the matrix and right hand side have been updated.
This works well in the sense that it solves the system in a
10 matches
Mail list logo