[petsc-users] Euclid or Boomeramg vs ILU: questions.

ory costs with *ILU, Euclid, > Boomeramg*? > > 5. At what stage are memory leaks most likely? > > > > In any case, thank you so much for your attention! Will be grateful for > any response. > > Kind regards, > Viktor Nazdrachev > R senior researcher > Geosteering Technologies LLC -- Ed Bueler Dept of Mathematics and Statistics University of Alaska Fairbanks Fairbanks, AK 99775-6660 306C Chapman

[petsc-users] re Example finite volume silver in PETSc

. > > Thank you. > > Regards, -- Ed Bueler Dept of Mathematics and Statistics University of Alaska Fairbanks Fairbanks, AK 99775-6660 306C Chapman

[petsc-users] new book introducing PETSc for PDEs

k with too; send them your book idea! Ed -- Ed Bueler Dept of Mathematics and Statistics University of Alaska Fairbanks Fairbanks, AK 99775-6660 306C Chapman

Re: [petsc-users] recommended error-checking for PetscInitialize() and PetscFinalize()

Thanks Matt and Barry. Clear. Ed On Fri, Sep 11, 2020 at 2:05 PM Barry Smith wrote: > > > On Sep 11, 2020, at 2:58 PM, Matthew Knepley wrote: > > On Fri, Sep 11, 2020 at 3:57 PM Ed Bueler wrote: > >> Dear PETSc -- >> >> I notice in the users manual

[petsc-users] recommended error-checking for PetscInitialize() and PetscFinalize()

of these choices is that if PetscInitialize() fails then CHKERRQ(ierr) may not do the right thing, while if PetscFinalize() fails then that should be the result of main() (without any fiddling by CHKERRQ etc.). Is this the correct understanding? Thanks, Ed -- Ed Bueler Dept of Mathematics and Statistics

Re: [petsc-users] ARKIMEX produces incorrect values

/ts/tutorials/ examples seem to unconditionally supply all four, and I can't tell (e.g. from -ts_view) which parts are seen and called by the various methods. Ed On Mon, Aug 31, 2020 at 2:27 PM Constantinescu, Emil M. wrote: > > > > > > > > > > > On

Re: [petsc-users] ARKIMEX produces incorrect values

ges, for example, and probably others. > 5- other Not sure. Thanks, Ed On Mon, Aug 31, 2020 at 6:09 AM Constantinescu, Emil M. wrote: > > On 8/30/20 6:04 PM, Ed Bueler wrote: > > Actually, ARKIMEX is not off the hook. It still gets the wrong answer if > told the

Re: [petsc-users] ARKIMEX produces incorrect values

is the identity? Ed On Sun, Aug 30, 2020 at 3:24 PM Barry Smith wrote: > > > > On Aug 30, 2020, at 6:04 PM, Jed Brown wrote: > > > > Ed Bueler writes: > > > >> Darn, sorry. > >> > >> I realize the ARKIMEX page does say "Methods with a

Re: [petsc-users] ARKIMEX produces incorrect values

used with ODE in which the stiff part G(t,X,Xdot) has the form Xdot + Ghat(t,X)." means. Ed On Sun, Aug 30, 2020 at 2:57 PM Ed Bueler wrote: > Darn, sorry. > > I realize the ARKIMEX page does say "Methods with an explicit stage can > only be used with ODE in which the stiff pa

Re: [petsc-users] ARKIMEX produces incorrect values

2020 at 2:44 PM Ed Bueler wrote: > Dear PETSc -- > > I tried twice to make this an issue at the gitlab.com host site, but both > times got "something went wrong (500)". So this is a bug report by > old-fashioned means. > > I created a TS example, > https://git

[petsc-users] ARKIMEX produces incorrect values

It should appear SIAM Press in a couple of months. In all the examples in my book, only my diffusion-reaction system example using F(t,u,u') = G(t,u) is broken. Thus the motivation for a trivial ODE example as above. -- Ed Bueler Dept of Mathematics and Statistics University of Alaska

[petsc-users] what is the restriction for PCMG, and how to control it?

s/examples/tutorials/, ./ex5 -da_refine 3 -pc_type mg -help | grep XX returns nothing at all for XX=injection,restriction,interpolation. Would a new option -pc_mg_restriction_type transpose|injection make sense? Thanks, Ed -- Ed Bueler Dept of Math and Stat and Geophysical Institute University

Re: [petsc-users] Fwd: nondeterministic behavior with ./program -help | head

to know what help to print. > > > > On Apr 26, 2017, at 5:34 PM, Ed Bueler <elbue...@alaska.edu> wrote: > > > > [Missed reply all.] > > > > -- Forwarded message -- > > From: Ed Bueler <elbue...@alaska.edu> > > Date: Wed, Ap

[petsc-users] Fwd: nondeterministic behavior with ./program -help | head

[Missed reply all.] -- Forwarded message -- From: Ed Bueler <elbue...@alaska.edu> Date: Wed, Apr 26, 2017 at 2:33 PM Subject: Re: [petsc-users] nondeterministic behavior with ./program -help | head To: Barry Smith <bsm...@mcs.anl.gov> Yes, I thought of something

Re: [petsc-users] nondeterministic behavior with ./program -help | head

t 12:55 PM, Barry Smith <bsm...@mcs.anl.gov> wrote: > > Is > > ./ex5 -help -no_signal_handler | head > > good enough? > > Normally we won't want to turn off catching of PIPE errors since it might > hide other errors. > > Barry > > > > > On Apr

[petsc-users] nondeterministic behavior with ./program -help | head

re indicates an integer that MMS3 will be evaluated with 2^m_par, 2^n_par-- Petsc Development GIT revision: v3.7.6-3453-ge45481d470 GIT Date: 2017-04-26 13:00:15 -0500 ~/petsc/src/snes/examples/tutorials[master*]$

Re: [petsc-users] -snes_mf_operator yields "No support for this operation for this object type" in TS codes?

-snes_fd or -snes_fd_color or -snes_mf. > > > > Ed > > > > > > > > On Sun, Apr 2, 2017 at 6:23 PM, Jed Brown <j...@jedbrown.org> wrote: > > The issue is that we need to create > > > > a*I - Jrhs > > > > and this is currently done by

Re: [petsc-users] -snes_mf_operator yields "No support for this operation for this object type" in TS codes?

to fix that logic. Does this branch work for > you? > > https://bitbucket.org/petsc/petsc/pull-requests/655/fix- > flaw-with-tssetrhsjacobian-and-no/diff > > Ed Bueler <elbue...@alaska.edu> writes: > > > Dear PETSc -- > > > > I have a TS-using and DMDA-usi

[petsc-users] -snes_mf_operator yields "No support for this operation for this object type" in TS codes?

ch (commit 0c7851c55cba8e40da5083f79ba1ff846acd45b2). Thanks for your help and awesome library! Ed -- Ed Bueler Dept of Math and Stat and Geophysical Institute University of Alaska Fairbanks Fairbanks, AK 99775-6660 301C Chapman

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

if you look in p>2 cases, add additional degeneracy to the diffusivity, make the zero-order term signed, use SNESVI (because the source term is signed and so a VI/NCP must be solved), and you add a sliding flux. I was trying to work by analogy ... Ed On Wed, Feb 15, 2017 at 9:37 AM, Jed Br

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

user has decided with ts_type to use for production if it is > fully implicit or explicit then they can depending on the type selected, > write just a left hand side, just a right hand side for higher efficiency > (less update of ghost points, fewer iterations over loops etc). > >> &g

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

something I want, then indeed I will suck it up and learn it and write it. Recall my original request was for the ability to know if the TS type I have chosen will actually use the SNES I have set ... so long ago it now seems. Ed On Tue, Feb 14, 2017 at 6:56 PM, Jed Brown <j...@jedbrown.o

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

user code testing what method is running etc.. No such thing for now; my current code is not fragile. (Is an API based on "side" of an equation inherently suspect?) Ed On Tue, Feb 14, 2017 at 10:41 AM, Jed Brown <j...@jedbrown.org> wrote: > Ed Bueler <elbue...@alaska.edu

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

@mcs.anl.gov> wrote: > > > On Feb 13, 2017, at 1:53 PM, Ed Bueler <elbue...@alaska.edu> wrote: > > > > Dear Petsc -- > > > > This is the first of two short TS usage questions. > > > > My problem is both stiff (diffusive PDE) and constrained, so I req

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

norm supported". I was surprised because I figured if you could implement one (vector) norm you could do any. Why? If there is a good reason, *which I expect*, is there a citation on the adaptive fundamentals that would explain this surprise? Thanks! Ed -- Ed Bueler Dept of Mat

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

, but it is not a DAE and it *is* constrained. Thanks! Ed PS I'd prefer not to enumerate the existing TS types and error on the bad ones. It is not nicely-maintainable. -- Ed Bueler Dept of Math and Stat and Geophysical Institute University of Alaska Fairbanks Fairbanks, AK 99775-6660 301C

Re: [petsc-users] meaning of constants in classical iterations

ick Sanan <patrick.sa...@gmail.com> wrote: > On Fri, Jan 20, 2017 at 2:13 AM, Ed Bueler <elbue...@alaska.edu> wrote: > > Dear PETSc -- > > > > In my humble opinion, the answers to the following rather basic questions > > are missing from the petsc users manual

[petsc-users] meaning of constants in classical iterations

ve two questions by telling me that classical iterations are lame? Answer: Not the point! These iterations are used blockwise (e.g. ASM) and as smoothers in many proper PETSc applications. Clarity is helpful! Thanks for the great tool, as usual! Ed -- Ed Bueler Dept of Math and Stat and Ge

Re: [petsc-users] SetVariableBounds vs ComputeVariableBounds

Justin -- Yeah, good point. SNESVISetVariableBounds() works fine, at least in ex9.c (see attached patch). The reason for the other choice, which I found in my 5 year old email, was some bug in petsc3.2. Ed Date: Wed, 22 Jun 2016 08:42:33 +0100 > From: Justin Chang > To:

Re: [petsc-users] petsc-users Digest, Vol 90, Issue 34

ive function at all. (That is a question for the TAO developers?) In the SNESVINEWTONXSLS solvers an L^2 merit function (i.e. square of the residual) replaces the not-even-present objective function when doing line search, just as happens in solving nonlinear equations. Ed -- Ed Bueler Dept of Math and S

Re: [petsc-users] use of VecGetArray() with dof>1 (block size > 1)

hout a DM, you probably have some > analogous data structure and I would suggest making a Vec accessor in > the language of that discretization. If you really think the flat array > of structs that you're getting from VecGetArray (with an awkward cast) > is the best interface for your

Re: [petsc-users] use of VecGetArray() with dof>1 (block size > 1)

yhe...@gmail.com> wrote: > > > On 27 May 2016 at 21:24, Ed Bueler <elbue...@alaska.edu> wrote: > >> Dave -- >> >> Perhaps you should re-read my questions. >> > > Actually - maybe we got our wires crossed from the beginning. > I

Re: [petsc-users] use of VecGetArray() with dof>1 (block size > 1)

inside DMDAVecGetArray(). The same cast is simply exposed to the user if they want to get a struct-valued pointer from VecGetArray(). Ed On Fri, May 27, 2016 at 12:06 PM, Dave May <dave.mayhe...@gmail.com> wrote: > > > On 27 May 2016 at 20:34, Ed Bueler <elbue...@alaska.edu&

Re: [petsc-users] -snes_mf and -snes_fd

> We should have -snes_jacobian <analytic,mf,mf_operator,fd,fd_color> > or similar. That would be nice! Ed -- Ed Bueler Dept of Math and Stat and Geophysical Institute University of Alaska Fairbanks Fairbanks, AK 99775-6660 301C Chapman and 410D Elvey 907 474-7693 and 907 474-71

Re: [petsc-users] how to get full trajectory from TS into petsc binary

s in the same file, one file per timestep.) > you could then use a variation of the above python code to read from those > files (I think the vector comes first in the trajectory file). > > Barry > > > On Mar 1, 2016, at 1:26 AM, Ed Bueler <elbue...@alaska.edu> wrote: &g

Re: [petsc-users] how to get full trajectory from TS into petsc binary

tes)[0] else: break f.close() However, was there a more elegant intended method? I am disturbed by the apparent need to specify big-endian-ness (= '>d') in the struct.unpack() call. Ed On Mon, Feb 29, 2016 at 9:39 PM, Ed Bueler <elbue...@alaska.edu> wrote:

Re: [petsc-users] how to get full trajectory from TS into petsc binary

that > PetscBinaryIO.py could read them in currently. I don't know how far we want > to go in "spiffing up" the PETSc binary format to do more elaborate things > since, presumably, other more power IO tools exist that would be used for > "real" problems? > >

[petsc-users] how to get full trajectory from TS into petsc binary

ss. (I am not sure what the boolean option -ts_save_trajectory does, by the way.) Thanks! Ed PS Sorry if this is a "RTFM" question, but so far I can't find the documentation. -- Ed Bueler Dept of Math and Stat and Geophysical Institute University of Alaska Fairbanks Fairbanks, AK 997

Re: [petsc-users] master branch option "-snes_monitor_solution"

> Further work is needed to unify and simplified the various monitor options. > > > On Nov 22, 2015, at 8:55 PM, Ed Bueler <elbue...@alaska.edu> wrote: > > > > Barry -- > > > > That is reassuring, actually. That is, knowing that occasionally ya

[petsc-users] master branch option "-snes_monitor_solution"

ETSC ERROR: -snes_monitor_solution [0]PETSC ERROR: End of Error Message ---send entire error message to petsc-ma...@mcs.anl.gov-- application called MPI_Abort(MPI_COMM_WORLD, 85) - process 0 [unset]: aborting job: application called MPI_Abort(MPI_COMM_WORLD, 85) - process 0

[petsc-users] portability question for reduce operations

but no error in serial. On my two other machines there is no error with option 1. All three machines work fine with option 3. I have not tested option 2. -- Ed Bueler Dept of Math and Stat and Geophysical Institute University of Alaska Fairbanks Fairbanks, AK 99775-6660 301C Chapman and 410D Elvey 907

[petsc-users] how to catch zero pivot errors underneath SNES

*is* a DIVERGED_LINEAR_SOLVE error anyway. If not, can I at least catch the zero pivot error so it is not an error (i.e. so the program continues), without ignoring all errors? I think I don't want to push PetscIgnoreErrorHandler, because I am happy with how other errors are handled. Thanks, Ed -- Ed

Re: [petsc-users] how to catch zero pivot errors underneath SNES

, 2015 at 10:08 AM, Ed Bueler elbue...@alaska.edu wrote: Dear PETSc -- Context: I am solving a not very smooth and very nonlinear problem using SNES. Robustness is an issue but I think there is no theory to help me, and I am stuck mollifying the problem a bit when divergence happens. Thus I

[petsc-users] evidence that VINEWTONRSLS is broken (in next, but not in maint)

= CONVERGED_FNORM_RELATIVE errors:av |u-uexact| = 7.156e-03 |u-uexact|_inf = 2.550e-02 -- Ed Bueler Dept of Math and Stat and Geophysical Institute University of Alaska Fairbanks Fairbanks, AK 99775-6660 301C Chapman and 410D Elvey 907 474-7693 and 907 474-7199 (fax 907 474-5394)

[petsc-users] how to read netcdf into petsc

, this is a taunt. Thanks for help! Ed -- Ed Bueler Dept of Math and Stat and Geophysical Institute University of Alaska Fairbanks Fairbanks, AK 99775-6660 301C Chapman and 410D Elvey 907 474-7693 and 907 474-7199 (fax 907 474-5394)