Re: [petsc-users] Question about memory usage in Multigrid preconditioner

2016-09-09 Thread Hengjie Wang
Hi Barry, I checked. On the supercomputer, I had the option "-ksp_view_pre" but it is not in file I sent you. I am sorry for the confusion. Regards, Frank On Friday, September 9, 2016, Barry Smith wrote: > > > On Sep 9, 2016, at 3:11 PM, frank

Re: [petsc-users] Question about memory usage in Multigrid preconditioner

2016-09-09 Thread Barry Smith
> On Sep 9, 2016, at 3:11 PM, frank wrote: > > Hi Barry, > > I think the first KSP view output is from -ksp_view_pre. Before I submitted > the test, I was not sure whether there would be OOM error or not. So I added > both -ksp_view_pre and -ksp_view. But the options

Re: [petsc-users] Sorted CSR Matrix and Multigrid PC.

2016-09-09 Thread Manuel Valera
Thank you SO much for helping me out on this. Dumb error from my part not to notice. This means the common /mypcs/ elements are preconfigure internally in PETSc? Regards and happy weekend, On Fri, Sep 9, 2016 at 5:19 PM, Barry Smith wrote: > > The missing third argument

Re: [petsc-users] Sorted CSR Matrix and Multigrid PC.

2016-09-09 Thread Barry Smith
The missing third argument to PCApply means that you haven't allocated the vector work? > On Sep 9, 2016, at 3:34 PM, Manuel Valera wrote: > > Hello everyone, > > I'm having an error with my program that i cannot understand, the weird part > is that the same

Re: [petsc-users] DMPlex problem

2016-09-09 Thread Matthew Knepley
On Fri, Sep 9, 2016 at 9:49 AM, Mark Lohry wrote: > Regarding DMPlex, I'm unclear on the separation between Index Sets and > Plex for unstructured solvers. I've got an existing unstructured serial > solver using petsc's newton-krylov > Index Sets are lists of integers.

Re: [petsc-users] Question about memory usage in Multigrid preconditioner

2016-09-09 Thread frank
Hi Barry, I think the first KSP view output is from -ksp_view_pre. Before I submitted the test, I was not sure whether there would be OOM error or not. So I added both -ksp_view_pre and -ksp_view. Frank On 09/09/2016 12:38 PM, Barry Smith wrote: Why does ksp_view2.txt have two KSP

Re: [petsc-users] Diagnosing a difference between "unpreconditioned" and "true" residual norms

2016-09-09 Thread Barry Smith
Patrick, I have only seen this when the "linear" operator turned out to not actually be linear or at least not linear in double precision. Are you using differencing or anything in your MatShell that might make it not be a linear operator in full precision? Since your problem is so

Re: [petsc-users] Question about memory usage in Multigrid preconditioner

2016-09-09 Thread Barry Smith
Why does ksp_view2.txt have two KSP views in it while ksp_view1.txt has only one KSPView in it? Did you run two different solves in the 2 case but not the one? Barry > On Sep 9, 2016, at 10:56 AM, frank wrote: > > Hi, > > I want to continue digging into the memory

Re: [petsc-users] Question about memory usage in Multigrid preconditioner

2016-09-09 Thread frank
Hi, I want to continue digging into the memory problem here. I did find a work around in the past, which is to use less cores per node so that each core has 8G memory. However this is deficient and expensive. I hope to locate the place that uses the most memory. Here is a brief summary of

Re: [petsc-users] DMPlex problem

2016-09-09 Thread Mark Lohry
Regarding DMPlex, I'm unclear on the separation between Index Sets and Plex for unstructured solvers. I've got an existing unstructured serial solver using petsc's newton-krylov solvers. Should I be looking to parallelize this via IS or plex? Is there an interface for either that allows me to

[petsc-users] Diagnosing a difference between "unpreconditioned" and "true" residual norms

2016-09-09 Thread Patrick Sanan
I am debugging a linear solver which uses a custom operator and preconditioner, via MATSHELL and PCSHELL. Convergence seems to be fine, except that I unexpectedly see a difference between the "unpreconditioned" and "true" residual norms when I use -ksp_monitor_true_residual with a

Re: [petsc-users] DMPlex problem

2016-09-09 Thread Matthew Knepley
On Fri, Sep 9, 2016 at 4:04 AM, Morten Nobel-Jørgensen wrote: > Dear PETSc developers and users, > > Last week we posted a question regarding an error with DMPlex and multiple > dofs and have not gotten any feedback yet. This is uncharted waters for us, > since we have gotten used

[petsc-users] DMPlex problem

2016-09-09 Thread Morten Nobel-Jørgensen
Dear PETSc developers and users, Last week we posted a question regarding an error with DMPlex and multiple dofs and have not gotten any feedback yet. This is uncharted waters for us, since we have gotten used to an extremely fast feedback from the PETSc crew. So - with the chance of sounding