Re: [petsc-users] hypre / hip usage

2022-01-24 Thread Paul T. Bauman
On Mon, Jan 24, 2022 at 9:53 AM Jed Brown wrote: > "Paul T. Bauman" writes: > > > 1. `rocgdb` will be in your PATH when the `rocm` module is loaded. This > is > > gdb, but with some extra AMDGPU goodies. AFAIK, you cannot, yet, do > > stepping through a ker

Re: [petsc-users] hypre / hip usage

2022-01-24 Thread Paul T. Bauman
le through `rocgdb` as one would do with `gdb` (`rocgdb` is literally `gdb` built with extra AMD stuff (that stuff is either upstreamed or being upstreamed to gdb BTW)). You can do it in batch mode as well so you can dump the logs from each MPI process. > > On Mon, Jan 24, 2022 at 10:16 AM

Re: [petsc-users] hypre / hip usage

2022-01-24 Thread Paul T. Bauman
On Mon, Jan 24, 2022 at 8:53 AM Matthew Knepley wrote: > On Mon, Jan 24, 2022 at 9:24 AM Mark Adams wrote: > >> What is the fastest way to rebuild hypre? reconfiguring did not work and >> is slow. >> >> I am printf debugging to find this HSA_STATUS_ERROR_MEMORY_FAULT (no >> debuggers other than

Re: [petsc-users] hypre / hip usage

2022-01-21 Thread Paul T. Bauman
Dammit, didn't reply-all. Sorry. On Fri, Jan 21, 2022 at 8:52 AM Paul T. Bauman wrote: > > > On Fri, Jan 21, 2022 at 8:37 AM Jed Brown wrote: > >> "Paul T. Bauman" writes: >> >> > On Fri, Jan 21, 2022 at 8:19 AM Jed Brown wrote: >> &

Re: [petsc-users] hypre / hip usage

2022-01-21 Thread Paul T. Bauman
On Fri, Jan 21, 2022 at 8:19 AM Jed Brown wrote: > Mark Adams writes: > > > Two questions about hypre on HIP: > > > > * I am doing this now. Is this correct? > > > > '--download-hypre', > > '--download-hypre-configure-arguments=--enable-unified-memory', > Apologies for interjecting, but

Re: [petsc-users] DMPlex higher order elements

2016-08-24 Thread Paul T. Bauman
On Wed, Aug 24, 2016 at 1:22 PM, Matthew Knepley wrote: > Notice that Paraview does exactly the wrong thing here in that it has > straight lines connecting the midpoints and corners of the triangles. > Everytime I meet someone from Kitware, I complain about this and their representation of quadr

Re: [petsc-users] Providing Interpolation/Injections to Multigrid infrastructure

2016-02-12 Thread Paul T. Bauman
Hi Barry, On Thu, Feb 11, 2016 at 2:45 PM, Barry Smith wrote: > > > On Feb 11, 2016, at 1:36 PM, Boris Boutkov wrote: > > > > Hello All, > > > > I'm currently working on hooking into the PETSc multi-grid > infrastructure by creating my own DMShell and providing my own > interpolation and inject

Re: [petsc-users] 3.5 -> 3.6 Change in MatOrdering Behavior

2016-01-21 Thread Paul T. Bauman
On Thu, Jan 21, 2016 at 12:51 PM, Hong wrote: > Paul: > It might be caused by our changes in default shift strategy. > We previously used '-pc_factor_shift_type NONZERO' for ilu, then changed > to '-pc_factor_shift_type NONE'. > For your test, I get > ./ex10 -f0 test.mat -rhs 0 -pc_type asm -pc_a

Re: [petsc-users] 3.5 -> 3.6 Change in MatOrdering Behavior

2016-01-21 Thread Paul T. Bauman
Thanks Hong, On Thu, Jan 21, 2016 at 12:16 PM, Hong wrote: > Paul : > Using petsc-dev (we recently added feature for better displaying > convergence behavior), > OK, good to know, thanks. > I found that '-sub_pc_factor_mat_ordering_type 1wd' causes zero pivot: > I figured it was something al

[petsc-users] 3.5 -> 3.6 Change in MatOrdering Behavior

2016-01-21 Thread Paul T. Bauman
Greetings, We have a test that has started failing upon switching from 3.5.4 to 3.6.0 (actually went straight to 3.6.3 but checked this is repeatable with 3.6.0). I've attached the matrix generated with -mat_view binary and a small PETSc program that runs in serial that reproduces the behavior by

Re: [petsc-users] PETSc / GAMG user counts

2015-07-15 Thread Paul T. Bauman
On Wed, Jul 15, 2015 at 9:23 PM, Paul T. Bauman wrote: > (Shame on me, but is there a relevant GAMG publication?) > Arg, this came off bad, sorry. I meant, is there a particular publication for which GAMG should be cited?

Re: [petsc-users] PETSc / GAMG user counts

2015-07-15 Thread Paul T. Bauman
Maybe relevant pubs/bibtex could be added alongside the PETSc ones on the PETSc webpage to help this? (Shame on me, but is there a relevant GAMG publication?) It doesn't help with the disingenuous researcher not willing to go the extra step of proper citation, but would certainly help someone like

Re: [petsc-users] Crash when trying to use FD Jacobian

2015-06-05 Thread Paul T. Bauman
On Fri, Jun 5, 2015 at 3:37 PM, Barry Smith wrote: > > > On Jun 5, 2015, at 1:49 PM, Paul T. Bauman wrote: > > > > > > > > On Fri, Jun 5, 2015 at 2:42 PM, Harshad Sahasrabudhe < > hsaha...@purdue.edu> wrote: > > It is not the coloring you need

Re: [petsc-users] Crash when trying to use FD Jacobian

2015-06-05 Thread Paul T. Bauman
On Fri, Jun 5, 2015 at 2:42 PM, Harshad Sahasrabudhe wrote: > It is not the coloring you need to generate, just the nonzero structure. >> (From that PETSc computes the coloring) so do as I suggest in my other >> email. > > > Thanks, I'll try that out. > FYI, this information is contained in th

[petsc-users] KSPSetOperators/sub-preconditioner type

2012-07-11 Thread Paul T. Bauman
On Wed, Jul 11, 2012 at 8:47 AM, Jed Brown wrote: > > Thanks for your test case. Yes, this is a PETSc bug associated with > handling of DIFFERENT_NONZERO_PATTERN. Fixed here and will be in the next > patch release. > > http://petsc.cs.iit.edu/petsc/releases/petsc-3.3/rev/09d701958d66 > Awesome,

[petsc-users] KSPSetOperators/sub-preconditioner type

2012-07-10 Thread Paul T. Bauman
On Tue, Jul 10, 2012 at 1:34 PM, Paul T. Bauman wrote: > Greetings, > > This question is somewhat related to this thread: > http://lists.mcs.anl.gov/pipermail/petsc-users/2011-August/009583.html but > instead of resurrecting that thread, I decided to start a new one. > > Su

[petsc-users] KSPSetOperators/sub-preconditioner type

2012-07-10 Thread Paul T. Bauman
Greetings, This question is somewhat related to this thread: http://lists.mcs.anl.gov/pipermail/petsc-users/2011-August/009583.html but instead of resurrecting that thread, I decided to start a new one. Summary: In the libMesh wrapper of PETSc's KSP, a call is made to KSPSetOperators; I gather th

compiling PETSc with Intel MKL 10.0.1.14

2009-02-03 Thread Paul T. Bauman
pilation > doesn't work. > > Randy > > > > > Paul T. Bauman wrote: >> Sorry this took so long to get around to doing. So it turns out that >> there's a newer version of 2.3.3p13 posted at the PETSc ftp server. >> This worked flawlessly with the

compiling PETSc with Intel MKL 10.0.1.14

2008-08-30 Thread Paul T. Bauman
Please send the configure.log to petsc-maint at mcs.anl.gov > >Barry > > On Aug 15, 2008, at 1:39 PM, Paul T. Bauman wrote: > >> Was there ever a fix/workaround introduced for this? I'm using >> 2.3.3p13 and I'm having trouble getting the config to recogn

compiling PETSc with Intel MKL 10.0.1.14

2008-08-15 Thread Paul T. Bauman
Was there ever a fix/workaround introduced for this? I'm using 2.3.3p13 and I'm having trouble getting the config to recognize mkl 10.0.3.020. Thanks, Paul Barry Smith wrote: > > Could you email to petsc-maint at mcs.anl.gov ALL the messages as to > what goes wrong with > our current linking

CHKMEMQ

2008-07-23 Thread Paul T. Bauman
t; >Barry > > > On Jul 23, 2008, at 1:04 PM, Satish Balay wrote: > >> On Wed, 23 Jul 2008, Paul T. Bauman wrote: >> >>> Hello, >>> >>> Quick (and probably stupid) question: Does CHKMEMQ check *all* of >>> memory used >>&

CHKMEMQ

2008-07-23 Thread Paul T. Bauman
Hello, Quick (and probably stupid) question: Does CHKMEMQ check *all* of memory used by the program for corruption or just the memory PETSc allocated? I have a code that I'm assigned to work on that doesn't use PETSc (yet), but I was hoping to use this macro to track down what I think is a mem

PetscPrintf/ifort 9.1

2007-08-29 Thread Paul T. Bauman
esolves this problem. If you are > not using petsc-dev you can simply replace the file > src/sys/fileio/ftn-custom/zmprintf.c with the attached file then run > make lib shared > in that directory then relink your program. > >Barry > > > On Tue, 28 Aug 2007, Paul

PetscPrintf/ifort 9.1

2007-08-28 Thread Paul T. Bauman
shared > in that directory then relink your program. > >Barry > > > On Tue, 28 Aug 2007, Paul T. Bauman wrote: > > >> Hello, >> >> Kind of a non-numerical question - sorry. I'm using PetscPrintf from Fortran >> using the following call

PetscPrintf/ifort 9.1

2007-08-28 Thread Paul T. Bauman
(Since I don't > really know Fortran, I don't know how to tell a Fortran compiler to > embed a newline...) > > --Leif > > > Matthew Knepley wrote: >> On 8/28/07, Paul T. Bauman wrote: >>> Hello, >>> >>> Kind of a non-numer

PetscPrintf/ifort 9.1

2007-08-28 Thread Paul T. Bauman
Hello, Kind of a non-numerical question - sorry. I'm using PetscPrintf from Fortran using the following calling sequence (for example): call PetscPrintf(PETSC_COMM_WORLD, "=== \n", ierr) call PetscPrintf(PETSC_COMM_WORLD, " Be

FORTRAN usage of SNESSetConvergenceHistory

2007-08-07 Thread Paul T. Bauman
ory() it does not return the arrays, just the > count. You must keep around the arrays you passed into > SNESSetConvergenceHistory() and read from them. > >Barry > > > On Mon, 6 Aug 2007, Paul T. Bauman wrote: > > >> Hello, >> >> Is it possi

Programming in *.f90 free format with PETSc

2007-08-07 Thread Paul T. Bauman
Hi Ben, Here is the format I use in my .f90 modules. Note I do not know what the standard is, but this "style" has worked with ifort, gfortran, g95 and Visual Studio (so I'm told). module IAmAModule implicit none private !(you don't need this if you don't want it) #PREPROCESSING STATEMENTS

FORTRAN usage of SNESSetConvergenceHistory

2007-08-06 Thread Paul T. Bauman
Hello, Is it possible to use SNESSetConvergenceHistory/SNESGetConvergenceHistory with FORTRAN? Thanks, Paul

Preallocation for external direct solvers

2007-04-14 Thread Paul T. Bauman
t;Barry > > The bug is that MatSetType() does not know that AIJMUMPS is inherited from > AIJ and hence destroys the AIJ data structure (including the preallocation > information) before creating a whole new AIJMUMPS matrix. > > On Sat, 14 Apr 2007, Paul T. Bauman wrote: > &

Preallocation for external direct solvers

2007-04-14 Thread Paul T. Bauman
Hello, I'm developing some code and am using direct solvers (for the moment), in particular MUMPS and SuperLU. I noticed that PETSc is not preallocating the memory for it even though I specify maximum number of entries per row. call MatCreateSeqAIJ(PETSC_COMM_WORLD,total_num_dofs,total

need help with makefile

2007-03-03 Thread Paul T. Bauman
Aaron, Don't want to be rude and interject, but I looked at the unsteady.F file - I've always had trouble compiling when 'use statements' come after variable declarations. Give a try putting things in this order: program use modules implicit none include statements variable declarations c

PETSc and gfortran

2007-01-10 Thread Paul T. Bauman
rompt help. Paul Satish Balay wrote: > BUS error is a bit suspicious. What can you send the output from all > of the following commands.. > > cd src/ksp/ksp/examples/tutorials > make ex2f > ./ex2f > gdb ex2f > (gdb) run > > > Satish > > On Wed, 10 Jan 2007,

PETSc and gfortran

2007-01-10 Thread Paul T. Bauman
Yes, example ex5f broke in the exact same way. Note the preceding c/c++ examples worked perfectly. How would you suggest I proceed? Thanks, Paul Satish Balay wrote: > Can you reproduce this with a PETSc example? > > make test > > Satish > > On Wed, 10 Jan 2007,

PETSc and gfortran

2007-01-10 Thread Paul T. Bauman
Hello, Has anyone had any experience using PETSc and gfortran together? My code compiles, but when I run it, it crashes with the following error (mac Tiger 10.4.8, power pc, latest build of gfortran, MPICH2, petsc 2.3.2): [0]PETSC ERROR: ---

Two questions regarding SNES/KSP

2006-05-03 Thread Paul T. Bauman
Hello, I'm a new PETSc user, so if the questions I ask are redundant, I sincerely apologize. I was, however, unavailable to find the following information in the manual pages. I have two questions. First: When using the SNES package, is the initial guess for the /iterative linear solver/ (e.g