[petsc-users] Scientific Software Days Conference, April, 2017

Apologies if this is not the appropriate place. I thought folks here would be interested in this. Registration is open for the 8th Annual Scientific Software Days conference. Register here: http://scisoftdays.org/ The 8th Annual Scientific Software Days Conference (SSD) targets users and

[petsc-users] lock matrices in EPS

Hi ALL, I am solving a generalized eigenvalue problem Ax=\lambda Bx. I want to retrieve matrices A and B from EPS, but got the following errors: [0]PETSC ERROR: No support for this operation for this object type [0]PETSC ERROR: Cannot retrieve original matrices (have been modified) [0]PETSC

Re: [petsc-users] TS question 1: how to stop explicit methods because they do not use SNES(VI)?

On Feb 13, 2017, at 3:08 PM, Ed Bueler > wrote: Barry -- > Sounds like a bug to me. The methods should be checking if an > IFunction is being provided and error out in that case. I don't think it is that simple. I.e. having an IFunction and

Re: [petsc-users] TS question 1: how to stop explicit methods because they do not use SNES(VI)?

> On Feb 13, 2017, at 3:08 PM, Ed Bueler wrote: > > Barry -- > > > Sounds like a bug to me. The methods should be checking if an > > IFunction is being provided and error out in that case. > > I don't think it is that simple. I.e. having an IFunction and an explicit >

Re: [petsc-users] TS question 1: how to stop explicit methods because they do not use SNES(VI)?

Barry -- > Sounds like a bug to me. The methods should be checking if an > IFunction is being provided and error out in that case. I don't think it is that simple. I.e. having an IFunction and an explicit scheme is not, by itself, a bug. I think. If the user has provided IFunction

Re: [petsc-users] TS question 2: why not -ts_adapt_wnormtype 1 ?

I think it just wasn't implemented. Note that these have custom code because of the weighting vectors -- it isn't just a call to VecNorm. Patches welcome. Ed Bueler writes: > Dear Petsc -- > > This is the second of two short TS usage questions. > > I would like my TS, for

Re: [petsc-users] TS question 1: how to stop explicit methods because they do not use SNES(VI)?

> On Feb 13, 2017, at 1:53 PM, Ed Bueler wrote: > > Dear Petsc -- > > This is the first of two short TS usage questions. > > My problem is both stiff (diffusive PDE) and constrained, so I require > >-snes_type vinewton{rs|ss}ls > > *and* I split my ODE system into

[petsc-users] TS question 2: why not -ts_adapt_wnormtype 1 ?

Dear Petsc -- This is the second of two short TS usage questions. I would like my TS, for MOL on a PDE based on a spatial grid (DMDA) of course, to be tolerant of crap at a few locations. Thus I set -ts_adapt_wnormtype 1 and I got a very clear error saying "Only 2-norm and infinite norm

[petsc-users] TS question 1: how to stop explicit methods because they do not use SNES(VI)?

Dear Petsc -- This is the first of two short TS usage questions. My problem is both stiff (diffusive PDE) and constrained, so I require -snes_type vinewton{rs|ss}ls *and* I split my ODE system into IFunction and RHSFunction F(t,u,u_t) = G(t,u) where F(t,u,u_t) = u_t + f(t,u) in my case

Re: [petsc-users] TS routine calling order

Gideon Simpson writes: > In terms of the TS PreStep, PostStep, Monitor, etc., what is the calling > order? I’m specifically curious as to whether the Monitor routine is run > before or after the PostStep, but I’d be interested to know what the whole > workflow is.

[petsc-users] TS routine calling order

In terms of the TS PreStep, PostStep, Monitor, etc., what is the calling order? I’m specifically curious as to whether the Monitor routine is run before or after the PostStep, but I’d be interested to know what the whole workflow is. -gideon

Re: [petsc-users] TSSetPostStep for projected methods

Yup, implemented that with DM. Works as desired. Thanks. -gideon > On Feb 13, 2017, at 8:03 AM, Matthew Knepley wrote: > > On Sun, Feb 12, 2017 at 3:21 PM, Barry Smith > wrote: > > > On Feb 12, 2017, at 11:30 AM, Gideon

Re: [petsc-users] VecScatter betweeen Vectors with different parallel layouts

"Barletta, Ivano" writes: > Hello PETSc users > > I was wondering if you've got any suggestion for the issue > mentioned in the topic. > > I'm solving a linear system on my unstructured grid that has its > own partitioning. When I set up the linear system in PETSc I > use the

Re: [petsc-users] VecScatter betweeen Vectors with different parallel layouts

On Mon, Feb 13, 2017 at 5:51 AM, Barletta, Ivano wrote: > Hello PETSc users > > I was wondering if you've got any suggestion for the issue > mentioned in the topic. > > I'm solving a linear system on my unstructured grid that has its > own partitioning. When I set up the

Re: [petsc-users] TSSetPostStep for projected methods

On Sun, Feb 12, 2017 at 3:21 PM, Barry Smith wrote: > > > On Feb 12, 2017, at 11:30 AM, Gideon Simpson > wrote: > > > > I’ve begun working on implementing a projected time stepping method for > solving y’ = f(y), where g(y(t)) is an invariant of the

[petsc-users] VecScatter betweeen Vectors with different parallel layouts

Hello PETSc users I was wondering if you've got any suggestion for the issue mentioned in the topic. I'm solving a linear system on my unstructured grid that has its own partitioning. When I set up the linear system in PETSc I use the model partitioning whose grid points might not be evenly