fs/95504034.gadi-pbs/petsc/include/petsc/private/kspimpl.h:383
#25 0x15378d4d1097 in KSPSolve_CG (ksp=0x2cc7500) at
/jobfs/95504034.gadi-pbs/petsc/src/ksp/ksp/impls/cg/cg.c:162
Let me know if there is anything further we can try to debug this issue
Kind regards
Stephan Kramer
On 02/09/2023 01:58
re is anything further we can try to debug this issue
Kind regards
Stephan Kramer
On 02/09/2023 01:58, Mark Adams wrote:
Fantastic!
I fixed a memory free problem. You should be OK now.
I am pretty sure you are good but I would like to wait to get any feedback
from you.
We should have a release at
ee whether the
memory errors are still there. Will let you know
Best wishes
Stephan Kramer
Coarsening details with various options for Level 6 of the test case:
In our original setup (using "old" petsc), we had:
rows=516, cols=516, bs=6
rows=12660, cols=12660, bs=6
ace,
hence the bs=3 to bs=6.
We have tried different values for the gamg_threshold but it doesn't
really seem to significantly alter the coarsening amount in that first
step.
Do you have any suggestions for further things we should try/look at?
Any feedback would be much appreciated
Best wishes
Step
hings we should try/look at?
Any feedback would be much appreciated
Best wishes
Stephan Kramer
Full logs including log_view timings available from
https://github.com/stephankramer/petsc-scaling/
In particular:
https://github.com/stephankramer/petsc-scaling/blob/main/before/Level_5/output_2.dat
Stephan,
Sorry for the long delay. I have updated the branch
barry/maint/mv-matnullspace-to-mat to use the Amat argument for passing in the
null space instead of the pmat argument. Could you please try it with your
configuration of solvers and let us know if it resolves the problem?
Sorry if we're talking completely cross-purposes here: but the pcksp solve
(that doesn't actually have a nullspace) is inside a solve that does have a
nullspace. If what you are saying is that applying the nullspace inside the
pcksp solve does not affect the outer solve, I can only see that if
If we move the location of the nullspace used in removal from the pmat to
the mat would that completely resolve the problem for you?
Barry
That would indeed resolve the issue (and to me would also make the most sense)
Stephan
On Mon, Jun 22, 2015 at 1:13 PM, Stephan Kramer s.kramer at
imperial.ac.uk
wrote:
Dear petsc devs
I've been trying to move our code from using KSPSetNullSpace to use
MatSetNullSpace instead. Although I appreciate the convenience of the
nullspace being propagated automatically through
On Mon, Jun 22, 2015 at 1:13 PM, Stephan Kramer s.kramer at imperial.ac.uk
wrote:
Dear petsc devs
I've been trying to move our code from using KSPSetNullSpace to use
MatSetNullSpace instead. Although I appreciate the convenience of the
nullspace being propagated automatically through
Dear petsc devs
I've been trying to move our code from using KSPSetNullSpace to use
MatSetNullSpace instead. Although I appreciate the convenience of the
nullspace being propagated automatically through the solver hierarchy,
I'm still a little confused on how to deal with the case that
On 01/04/14 16:07, Mark Adams wrote:
Stephan, I have pushed a pull request to fix this but for now you can just
use -mg_levels_ksp_type chebyshev -mg_levels_pc_type jacobi. This used to
be the default be we move to SOR recently.
Mark
Ah, that's great news. Thanks a lot for the effort. You're
On 21/03/14 11:34, Stephan Kramer wrote:
On 21/03/14 04:24, Jed Brown wrote:
Stephan Kramer s.kra...@imperial.ac.uk writes:
We have been having some problems with GAMG on petsc-dev (master) for
cases that worked fine on petsc 3.4. We're solving a Stokes equation
(just the velocity block
On 21/03/14 04:24, Jed Brown wrote:
Stephan Kramer s.kra...@imperial.ac.uk writes:
We have been having some problems with GAMG on petsc-dev (master) for
cases that worked fine on petsc 3.4. We're solving a Stokes equation
(just the velocity block) for a simple convection in a square box
Hi guys,
We have been having some problems with GAMG on petsc-dev (master) for cases that worked fine on petsc 3.4. We're solving a Stokes equation (just the velocity block) for a simple convection in a square
box (isoviscous). The problem only occurs if we supply a near null space (via
From: Stephan Kramer s.kra...@imperial.ac.uk
Subject: Fwd: [petsc-users] MatZeroRowsColumns for mpibaij
Date: 31 October 2013 13:16:05 GMT+11:00
To: Rhodri Davies rhodri.dav...@anu.edu.au
Original Message
Subject: [petsc-users] MatZeroRowsColumns for mpibaij
Date: Wed, 19 Jun
Dear all,
We have found the MatZeroRowsColumns() routine to be very useful for lifting boundary conditions. Unfortunately we also found that it isn't implemented for matrix type mpibaij and we're quite keen on
keeping the baij structure as it seems to give us significant savings for DG problems
On 15/11/10 09:54, TAY wee-beng wrote:
Hi,
I have a modified makefile originally given to me by the PETSc team. It
worked well but there's a small problem. Sometimes, when I have edited a
fortran file and tries to rebuild, it gives this error:
/opt/openmpi-1.4.1/bin/mpif90 -r8 -w95 -c -O3
should *not* be declared
external. You do
not need that line.
Cheers
Stephan
2010/9/30 Stephan Kramer s.kramer at imperial.ac.uk
mailto:s.kramer at imperial.ac.uk
On 30/09/10 15:31, Leo van Kampenhout wrote:
Hi all,
since it is mandatory to declare all subroutines
On 30/09/10 15:31, Leo van Kampenhout wrote:
Hi all,
since it is mandatory to declare all subroutines as external in
Fortran, is it possible for Modules to have subroutines? I'm unable to
declare the subroutine external inside the module itself, nor in the
program which is using it. Not
On 20/06/10 16:21, Stephan Kramer wrote:
On 20/06/10 05:07, Barry Smith wrote:
Hmm, this is a good question. There are tons of places where some sort
of inner solve is imbedded in an outer solver, in fact many levels of
nesting. We should handle this in a systematic way which likely
Dear all,
Is there a way in petsc, when performing inner solves like PCKSP or
MatSchurComplementGetKSP, to make the outer solve stop immediately and report
back a negative convergence reason? I find that often
when such inner solves fail, the outer solve happily continues and sometimes
Hi all,
I'm trying to write a portable configure script for our software that uses
petsc, and needs to deal with petsc installations on a number of different
platforms. In petsc 3.0 I could do
make -f $(PETSC_DIR)/conf/base getlinklibs
and
make -f $(PETSC_DIR)/conf/base getincludedirs
Satish Balay wrote:
On Fri, 23 Apr 2010, Stephan Kramer wrote:
Hi all,
I'm trying to write a portable configure script for our software that uses
petsc, and needs to deal with petsc installations on a number of different
platforms. In petsc 3.0 I could do
make -f $(PETSC_DIR)/conf/base
Hi all,
I'm slightly confused by the fact that if I use my BAIJ matrix in parallel and
select the sor preconditioner. it complains about mat-ops-relax and
mat-ops-pbrelax not being defined. Also I can't
find a MatRelax_MPIBAIJ in src/mat/impls/baij/mpi/mpibaij.c as I would have
expected. Is
of the
matrix, mat-A, with its=lits and lits=PETSC_NULL (?). However the first
line of MatRelax_SeqAIJ then says: its = its*lits. Is that right?
Please tell me if I'm totally misunderstanding how the routine works,
thanks for any help.
Cheers
Stephan
--
Stephan Kramer s.kramer
Thanks for your answers
Barry Smith wrote:
On Sep 22, 2009, at 8:47 AM, Stephan Kramer wrote:
Hi all,
I have some questions basically about the MatRelax_SeqAIJ routine:
If I understand correctly there are 2 versions of the sor routine
depending on whether or not there is a zero guess
what you expect for the sparse matrices, and return a 1d array of length
nrows*ncolumns in the case of MATDENSE?
Cheers
Stephan Kramer
--
Stephan Kramer
Applied Modelling and Computation Group,
Department of Earth Science and Engineering,
Imperial College London
Satish Balay wrote:
On Sat, 30 May 2009, Stephan Kramer wrote:
Satish Balay wrote:
On Sat, 30 May 2009, Stephan Kramer wrote:
Thanks a lot for looking into this. The explicit fortran interfaces are in
general very useful. The problem occurred for me with petsc-3.0.0-p1. I'm
happy to try
debating how to handle it.
Barry
On May 23, 2009, at 9:00 AM, Stephan Kramer wrote:
Hi all,
First of all thanks of a lot for providing explicit fortran interfaces for
most functions in Petsc 3. This is of great help. I do however run into a
problem using MatGetInfo. The calling sequence
Satish Balay wrote:
On Sat, 30 May 2009, Stephan Kramer wrote:
Thanks a lot for looking into this. The explicit fortran interfaces are in
general very useful. The problem occurred for me with petsc-3.0.0-p1. I'm
happy to try it out with a more recent patch-level or with petsc-dev.
Did you
Hi all,
First of all thanks of a lot for providing explicit fortran interfaces
for most functions in Petsc 3. This is of great help. I do however run
into a problem using MatGetInfo. The calling sequence for fortran
(according to the manual) is:
double precision info(MAT_INFO_SIZE)
Mat
Hi,
I'm trying to pick up a linear system (matrix and rhs) that I've written
out in a previous parallel run using MatView and VecView in binary mode.
Now when I've read the matrix and rhs vector the rows are evenly
distributed by PETSc over the processes. The partioning previously used
in the
Matthew Knepley wrote:
On Tue, Feb 10, 2009 at 2:11 PM, Stephan Kramer
stephan.kramer at imperial.ac.uk wrote:
Hi,
I'm trying to pick up a linear system (matrix and rhs) that I've written out
in a previous parallel run using MatView and VecView in binary mode. Now
when I've read
Barry Smith wrote:
On Feb 10, 2009, at 5:01 PM, Stephan Kramer wrote:
You can certainly use a VecScatter to redistribute a vector. You
already have the
indices from MatGetSubMatrix() so it should be easy.
However, you can also use VecLoadIntoVector() and provide an initial
Shao-Ching Huang wrote:
Hi,
I am trying to use CG to solve a singular systems (i.e. Poisson
equation with periodic conditions on all boundaries).
The code works when I use GMRES, but it diverges when I switch CG.
Since the null space is a constant vector, in the code I have:
I might
Paul T. Bauman wrote:
Hi Barry,
To preface, I sincerely appreciate your prompt response and
willingness to produce quick fixes.
However, it was pointed out to me in my office today (by someone who
knows more about interfacing C and FORTRAN than I and who watches the
petsc-users list)
object,
and you thus expect the value of context to be stored.
Cheers
Stephan Kramer
38 matches
Mail list logo