Re: [petsc-dev] GAMG and custom MatMults in smoothers

2018-07-02 Thread Matthew Knepley
On Mon, Jul 2, 2018 at 9:33 AM Vaclav Hapla wrote: > > > 2. 7. 2018 v 15:50, Matthew Knepley : > > On Mon, Jul 2, 2018 at 8:28 AM Vaclav Hapla > wrote: > >> 2. 7. 2018 v 15:05, Matthew Knepley : >> >> On Mon, Jul 2, 2018 at 7:54 AM Vaclav Hapla >

Re: [petsc-dev] GAMG and custom MatMults in smoothers

2018-07-02 Thread Matthew Knepley
On Mon, Jul 2, 2018 at 10:19 AM Vaclav Hapla wrote: > 2. 7. 2018 v 16:38, Matthew Knepley : > > On Mon, Jul 2, 2018 at 9:33 AM Vaclav Hapla > wrote: > >> >> >> 2. 7. 2018 v 15:50, Matthew Knepley : >> >> On Mon, Jul 2, 2018 at 8:28 AM Vaclav Hapla

Re: [petsc-dev] GAMG error with MKL

2018-07-04 Thread Matthew Knepley
On Tue, Jul 3, 2018 at 10:32 PM Jeff Hammond wrote: > > > On Tue, Jul 3, 2018 at 4:35 PM Mark Adams wrote: > >> On Tue, Jul 3, 2018 at 1:00 PM Richard Tran Mills >> wrote: >> >>> Hi Mark, >>> >>> I'm glad to see you trying out the AIJMKL stuff. I think you are the >>> first person trying to act

Re: [petsc-dev] snes/ex56 regression

2018-07-04 Thread Matthew Knepley
On Wed, Jul 4, 2018 at 8:47 AM Mark Adams wrote: > Sorry, this error was with MKL matrices. It is still an error inside of > Plex. My non-Plex tests (KSP) now work in a branch (pushed). > Ah, I might know what this is. When you ask the DM to create a matrix, it tries to convert it to the interna

Re: [petsc-dev] GAMG error with MKL

2018-07-04 Thread Matthew Knepley
On Wed, Jul 4, 2018 at 4:51 PM Jeff Hammond wrote: > On Wed, Jul 4, 2018 at 6:31 AM Matthew Knepley wrote: > >> On Tue, Jul 3, 2018 at 10:32 PM Jeff Hammond >> wrote: >> >>> >>> >>> On Tue, Jul 3, 2018 at 4:35 PM Mark Adams wrote: >>

Re: [petsc-dev] GAMG error with MKL

2018-07-05 Thread Matthew Knepley
On Thu, Jul 5, 2018 at 2:04 PM Mark Adams wrote: > > > On Thu, Jul 5, 2018 at 12:41 PM Tobin Isaac wrote: > >> On Thu, Jul 05, 2018 at 09:28:16AM -0400, Mark Adams wrote: >> > > >> > > >> > > Please share the results of your experiments that prove OpenMP does >> not >> > > improve performance fo

[petsc-dev] My configure is non-working again because of CMake

2018-07-05 Thread Matthew Knepley
Now we require a fully C++11 compliant compiler, or the CMake build dies. This is crazy. Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buf

Re: [petsc-dev] My configure is non-working again because of CMake

2018-07-05 Thread Matthew Knepley
does not require C++11 compliant compiler? We > could downgrade to that.. > I reverted back to 3.7 which worked, so somewhere in between 3.7 and 3.11 they made the change. Thanks, Matt > Satish > > On Thu, 5 Jul 2018, Matthew Knepley wrote: > > > Now we require a f

Re: [petsc-dev] My configure is non-working again because of CMake

2018-07-05 Thread Matthew Knepley
to: > https://bitbucket.org/petsc/petsc/commits/e5b329dd4ac > > No point in going back to having old issues.. > > Satish > > On Thu, 5 Jul 2018, Matthew Knepley wrote: > > > On Thu, Jul 5, 2018 at 5:18 PM Satish Balay wrote: > > > > > Logs? > > > &g

Re: [petsc-dev] My configure is non-working again because of CMake

2018-07-05 Thread Matthew Knepley
he version between 3.7 and 3.11 that works on your > box? > > > > The previous update was done due to: > > https://bitbucket.org/petsc/petsc/commits/e5b329dd4ac > > > > No point in going back to having old issues.. > > > > Satish > > > > On Thu

Re: [petsc-dev] Remove legacy tests?

2018-07-05 Thread Matthew Knepley
On Thu, Jul 5, 2018 at 5:37 PM Jed Brown wrote: > When can we delete the legacy test system? Are we currently using it > anywhere? > Not I said the fox. Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to

Re: [petsc-dev] Remove legacy tests?

2018-07-05 Thread Matthew Knepley
On Thu, Jul 5, 2018 at 6:07 PM Satish Balay wrote: > Ok - there is some confusion here. > > The test harness does not use targets in > examples/[tests,tutorials]/makefile. > > However we want to keep the functionality of us/users compiling > individual examples manually - without going through th

Re: [petsc-dev] My configure is non-working again because of CMake

2018-07-05 Thread Matthew Knepley
ed. We can move back to that one if we want. Maybe at the meeting we can issue a declaration of war on CMake (or at least get Zimmerman to send a telegram to Autotools). Thanks, Matt > > Satish > > > On Thu, 5 Jul 2018, Matthew Knepley wrote: > > > On Thu, Jul

Re: [petsc-dev] GAMG error with MKL

2018-07-06 Thread Matthew Knepley
good solutions. > > --Richard > > On Wed, Jul 4, 2018 at 11:38 PM, Smith, Barry F. > wrote: > >> >>Jed, >> >> You could use your same argument to argue PETSc should do >> "something" to help people who have (rightly or wron

Re: [petsc-dev] GAMG error with MKL

2018-07-09 Thread Matthew Knepley
gt; > Clever ideas from anyone on this list about how to use an adequate >> number of MPI ranks for PETSc while using only a subset of these ranks for >> the MPI+OpenMP application code will be appreciated, though I don't know if >> there are any good solutions. >&

Re: [petsc-dev] GAMG error with MKL

2018-07-09 Thread Matthew Knepley
On Mon, Jul 9, 2018 at 6:53 PM Jeff Hammond wrote: > On Mon, Jul 9, 2018 at 6:38 AM, Matthew Knepley wrote: > >> On Mon, Jul 9, 2018 at 9:34 AM Jeff Hammond >> wrote: >> >>> On Fri, Jul 6, 2018 at 4:28 PM, Smith, Barry F. >>> wrote: >>> >

Re: [petsc-dev] PetscSF and/or VecScatter with device pointers

2018-07-12 Thread Matthew Knepley
On Thu, Jul 12, 2018 at 6:47 AM Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > Dear petsc-dev, > > we're starting to explore (with Andreas cc'd) residual assembly on > GPUs. The question naturally arises: how to do GlobalToLocal and > LocalToGlobal. > There is not a lot of Mem Ba

Re: [petsc-dev] Reading periodic meshes from GMSH

2018-07-12 Thread Matthew Knepley
On Thu, Jul 12, 2018 at 6:20 AM Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > Dear petsc-dev, > > I am attempting to understand how Plex reads GMSH files with > periodicity. In particular, what happens to the coordinate localisation. > Stefano wrote this so he will have to help.

Re: [petsc-dev] Reading periodic meshes from GMSH

2018-07-12 Thread Matthew Knepley
On Thu, Jul 12, 2018 at 10:00 AM Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > > > On 12 Jul 2018, at 12:27, Stefano Zampini > wrote: > > > > Lawrence, > > > > I wrote the periodic GMSH reader and I can answer for it. > > > > With DMPlex, when we have periodicity, we currently ha

Re: [petsc-dev] PETSc blame digest (next) 2018-07-16

2018-07-16 Thread Matthew Knepley
On Mon, Jul 16, 2018 at 11:38 AM Fande Kong wrote: > Hi Satish and Barry, > > Any way to reproduce these warnings? I just tried the same configure > script on my desktop, and I did NOT see any warning. > You would need the same compiler/options. Is it not clear why they happen? Matt > Fand

Re: [petsc-dev] Should we replace hypre PILUT with hypre Euclid?

2018-07-17 Thread Matthew Knepley
On Tue, Jul 17, 2018 at 2:55 PM Jed Brown wrote: > PETSc had a Euclid interface until Barry removed it. > > > https://bitbucket.org/petsc/petsc/commits/22b6d1ca80e4f79123036fb6afd7eb4463b3692a > > See also > > https://www.mail-archive.com/petsc-users@mcs.anl.gov/msg33196.html I would also add t

Re: [petsc-dev] PETSc blame digest (next-tmp) 2018-07-20

2018-07-20 Thread Matthew Knepley
On Fri, Jul 20, 2018 at 4:28 PM PETSc checkBuilds < petsc-checkbui...@mcs.anl.gov> wrote: > > > Dear PETSc developer, > > This email contains listings of contributions attributed to you by > `git blame` that caused compiler errors or warnings in PETSc automated > testing. Follow the links to see

Re: [petsc-dev] Bad scaling of GAMG in FieldSplit

2018-07-26 Thread Matthew Knepley
On Thu, Jul 26, 2018 at 12:56 PM Fande Kong wrote: > > > On Thu, Jul 26, 2018 at 10:35 AM, Junchao Zhang > wrote: > >> On Thu, Jul 26, 2018 at 11:15 AM, Fande Kong wrote: >> >>> >>> >>> On Thu, Jul 26, 2018 at 9:51 AM, Junchao Zhang >>> wrote: >>> Hi, Pierre, From your log_view fil

Re: [petsc-dev] Bad scaling of GAMG in FieldSplit

2018-07-26 Thread Matthew Knepley
On Thu, Jul 26, 2018 at 2:43 PM Jed Brown wrote: > Matthew Knepley writes: > > > On Thu, Jul 26, 2018 at 12:56 PM Fande Kong wrote: > > > >> > >> > >> On Thu, Jul 26, 2018 at 10:35 AM, Junchao Zhang > >> wrote: > >&

Re: [petsc-dev] Petsc 19 conflict calendar

2018-08-07 Thread Matthew Knepley
On Tue, Aug 7, 2018 at 12:24 PM Tobin Isaac wrote: > > Hi all, > > If you subscribe to petsc-dev, you might be interested in getting in > on the ground floor of planning our Summer 2019 meeting in Atlanta. > > It would be good to get a list of conflicts for June-ish. We've just > found out, for

Re: [petsc-dev] any script to extract data from "-log_view"?

2018-08-08 Thread Matthew Knepley
On Wed, Aug 8, 2018 at 7:23 PM Fande Kong wrote: > Hi all, > > If we are going to do a scaling study, can we automatically make a plot > based on a bunch of output files (from -log_view)? It is a big pain to > manually copy data if we have a lot of output files. > I use -log_view :run1.py:asc

Re: [petsc-dev] How to enforce private data/methods in PETSc?

2018-08-10 Thread Matthew Knepley
On Fri, Aug 10, 2018 at 5:43 PM Junchao Zhang wrote: > I met several bugs that remind me to raise this question. In PETSc, > object of type A can arbitrarily access object of type B's data. But > designer of B may later change the meaning of its data (and of course, > update B's interfaces, whic

Re: [petsc-dev] Petsc 19 directly after SC 19 PC?

2018-08-10 Thread Matthew Knepley
On Fri, Aug 10, 2018 at 8:16 PM Tobin Isaac wrote: > In another thread it came up that the SC'19 planning committee meeting > June 3-4 is actually in Atlanta, too. I've checked and there are plenty of > options on campus for a meeting of our size that I can reserve for the > second half of the w

Re: [petsc-dev] How to enforce private data/methods in PETSc?

2018-08-10 Thread Matthew Knepley
ctions break the rules. Matt > --Junchao Zhang > > > On Fri, Aug 10, 2018 at 9:43 PM Matthew Knepley wrote: > >> On Fri, Aug 10, 2018 at 5:43 PM Junchao Zhang >> wrote: >> >>> I met several bugs that remind me to raise this question. In PETSc, >&

Re: [petsc-dev] How to enforce private data/methods in PETSc?

2018-08-11 Thread Matthew Knepley
pages/Sys/PetscMalloc4.html#PetscMalloc4>(B->cmap->N*from->starts[from->n],&contents->rvalues,412: >B->cmap->N*to->starts[to->n],&contents->svalues,413: > from->n,&contents->rwaits,414:

Re: [petsc-dev] IS block size meaning

2018-08-14 Thread Matthew Knepley
On Tue, Aug 14, 2018 at 5:32 AM Vaclav Hapla wrote: > Hello > > I would like to move the discussion from the PR > https://bitbucket.org/petsc/petsc/pull-requests/1076/fix-issetblocksize_general > because > it's quite a general concern. > > It's about the meaning of block size for IS. > > Barry t

Re: [petsc-dev] Undefined symbols for _kspfgmresmodifypcksp_ and _kspfgmresmodifypcnochange_ when rebuilding

2018-08-22 Thread Matthew Knepley
On Wed, Aug 22, 2018 at 6:35 AM Lawrence Mitchell wrote: > > > On 22 Aug 2018, at 10:04, Patrick Sanan wrote: > > > > This happens fairly frequently when I try to switch/update branches of > PETSc (here invoked by building my own code, but the error message looks > the same with "make check"): >

Re: [petsc-dev] public functions with no manual pages

2018-08-24 Thread Matthew Knepley
Fixed. Matt On Wed, Aug 22, 2018 at 4:15 PM Smith, Barry F. wrote: > > DMSNESCheckFromOptions() > > DMTSCheckFromOptions() > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- N

Re: [petsc-dev] Gathering information: types of interpolations in PETSc and which ones are supported?

2018-08-31 Thread Matthew Knepley
On Fri, Aug 31, 2018 at 2:20 PM Smith, Barry F. wrote: > > PETSc developers, > >There are a variety of "interpolation" modules in PETSc but the > documentation is scattered (mostly missing). Could everyone who knows > anything about the various modules provide a little information abo

Re: [petsc-dev] Gathering information: types of interpolations in PETSc and which ones are supported?

2018-09-01 Thread Matthew Knepley
On Fri, Aug 31, 2018 at 4:11 PM Smith, Barry F. wrote: > > On Aug 31, 2018, at 2:52 PM, Matthew Knepley wrote: > > > > On Fri, Aug 31, 2018 at 2:20 PM Smith, Barry F. > wrote: > > > > PETSc developers, > > > >There are a variet

Re: [petsc-dev] Compiling PETsc against openblas on Linux

2018-09-02 Thread Matthew Knepley
On Sun, Sep 2, 2018 at 11:11 AM Antonio wrote: > --with-blas-lapack-lib=-lopenblas --known-64-bit-blas-indices=0 > I am curious why you are giving the 64-bit option. This should be figured automatically unless you cannot executable programs on the platform on which you are configuring. What goe

Re: [petsc-dev] PETSc - MPI3 functionality

2018-09-08 Thread Matthew Knepley
On Sat, Sep 8, 2018 at 4:43 AM Tamara Dancheva wrote: > Hi Barry, > > I see the issue.. > > In the FEM library and solver that I am working on, PETSc is used all > throughout for the data distribution, synchronization of functions, > assembly. There is another UPC alternation of using the JANPAC

Re: [petsc-dev] MATLMVM

2018-09-11 Thread Matthew Knepley
On Tue, Sep 11, 2018 at 9:47 AM Dener, Alp wrote: > 1 base type and 8 subtypes. If there’s a better convention/structure to do > this in PETSc, I’d be happy to get a refactoring done ASAP this week so > that it’s cleaner in the release. > I guess what Lisandro is asking is, what backend differen

Re: [petsc-dev] TSBASICSYMPLECTIC

2018-09-11 Thread Matthew Knepley
On Tue, Sep 11, 2018 at 10:32 AM Zhang, Hong wrote: > A few related discussions can be found at > https://bitbucket.org/petsc/petsc/pull-requests/1108/rename-bsi-to-symplectic/diff > > In addition, what we have in PETSc now is "Basic Symplectic Integrators" > as introduced in Ernst Hairer's artic

Re: [petsc-dev] TSBASICSYMPLECTIC

2018-09-11 Thread Matthew Knepley
On Tue, Sep 11, 2018 at 11:53 AM Jed Brown wrote: > Matthew Knepley writes: > > > On Tue, Sep 11, 2018 at 10:32 AM Zhang, Hong wrote: > > > >> A few related discussions can be found at > >> > https://bitbucket.org/petsc/petsc/pull-requests/1108/ren

Re: [petsc-dev] Add dependency to configure

2018-09-19 Thread Matthew Knepley
On Wed, Sep 19, 2018 at 6:22 AM Marek Pecha wrote: > Hello, > > I work on clustering framework uses PETSc. Basically, I would like to use > this framework for segmentation images. So, I had idea - implement viewer > for loading images (JPEG,TIFF, etc) into PETSc. For this purpose, I would > like

Re: [petsc-dev] Can MATAIJ support block size?

2018-09-20 Thread Matthew Knepley
On Thu, Sep 20, 2018 at 1:26 PM Fande Kong wrote: > Hi Developers, > > MATBAIJ actually assumes that the point-block is dense. It is fine if the > block size is small for example less than 50, but we may suffer from a > memory issue if the block is large such as 1000. The block is not necessary >

[petsc-dev] Test output broken for test sets

2018-09-24 Thread Matthew Knepley
At least on my machine, this does not work correctly testset: nsize: {{1 2 4}} args: -use_generator -faces {{2,2,2 1,3,5 3,4,7}} args: -interpolate -distribute -interpolate_after_distribute {{0 1}} args: -dm_plex_check_pointsf test: #TODO some combinations generate wr

Re: [petsc-dev] DM not propagating into Fieldsplit when using external DMShell?

2018-09-24 Thread Matthew Knepley
On Mon, Sep 24, 2018 at 2:48 PM Boris Boutkov wrote: > Hello all, > > I am trying to solve a Stokes type problem using PETSc's Fieldsplit > infrastructure while providing external GMG information coming from libMesh > and am running into issues when using the two simultaneously. > > I have verifi

Re: [petsc-dev] Test output broken for test sets

2018-09-27 Thread Matthew Knepley
_6_tet.out > > not found. > Warning: > /Users/kruger/ptroot/upstream/petsc/src/sys/examples/tests/output/ex1_faces-2__2__2_6_tet.out > > not found. > > > Is this what you really want? > > Scott > > P.S. Jakub: Your issue is separate. I'll have a fix pu

Re: [petsc-dev] No LocalToGlobalMapping in DMPlex

2018-10-07 Thread Matthew Knepley
On Fri, Oct 5, 2018 at 2:58 PM Tristan Konolige < tristan.konol...@colorado.edu> wrote: > Hello PETSc Developers, > > I’d like to get an ISLocalToGlobalMapping from a DMPlex. But, > currently, DMGetLocalToGlobalMapping returns an empty > ISLocalToGlobalMapping. I can see that ISLocalToGlobalMappin

Re: [petsc-dev] MUMPS silent errors

2018-10-10 Thread Matthew Knepley
On Wed, Oct 10, 2018 at 3:12 PM Zhang, Junchao wrote: > I met a case where MUMPS numeric factorization returned an error code -9 in > mumps->id.INFOG(1) but A->erroriffailure was false in the following code in > mumps.c > > 1199: *PetscErrorCode >

Re: [petsc-dev] Question on extracting submatrixes from global projection when using FS + external DMShell + MG

2018-10-10 Thread Matthew Knepley
On Wed, Oct 10, 2018 at 12:23 PM Boris Boutkov wrote: > Hello again all, > > I've been trying to get GMG+FS to work with an external DMShell provided > by libMesh, see for reference > https://www.mail-archive.com/petsc-dev@mcs.anl.gov/msg23497.html . > > The situation is that while I've managed t

Re: [petsc-dev] Question on extracting submatrixes from global projection when using FS + external DMShell + MG

2018-10-11 Thread Matthew Knepley
Matt > Thanks again for the help, > - Boris > > > > On Wed, Oct 10, 2018 at 8:21 PM Matthew Knepley wrote: > >> On Wed, Oct 10, 2018 at 12:23 PM Boris Boutkov >> wrote: >> >>> Hello again all, >>> >>> I've been tryin

Re: [petsc-dev] PETSc OpenMPI library test interminable

2018-10-15 Thread Matthew Knepley
On Mon, Oct 15, 2018 at 6:21 AM Antonio Trande wrote: > Hi all. > > PETSc-3.10.2 compiled with OpenMPI-2.1.5 on Fedora 30 (devel branch) (or > with OpenMPI-2.1.1 on Fedora 28), on x86 architectures only: > This looks like an OpenMPI bug. Does it work with MPICH? Thanks, Matt > the 'sr

Re: [petsc-dev] PETSc OpenMPI library test interminable

2018-10-16 Thread Matthew Knepley
On Tue, Oct 16, 2018 at 5:30 AM Antonio Trande wrote: > > This looks like an OpenMPI bug. Does it work with MPICH? > > > > Thanks, > > > > Matt > > Yes, it works with MPICH. > Did you report this to OpenMPI? They are usually good with using PETSc to reproduce. Thanks, Matt > -- >

Re: [petsc-dev] PETSc OpenMPI library test interminable

2018-10-16 Thread Matthew Knepley
On Tue, Oct 16, 2018 at 12:49 PM Antonio wrote: > > > On 16/10/18 12:07, Matthew Knepley wrote: > > On Tue, Oct 16, 2018 at 5:30 AM Antonio Trande > <mailto:anto.tra...@gmail.com>> wrote: > > > > > This looks like an OpenMPI bug. Does it wo

Re: [petsc-dev] --download and spack

2018-10-21 Thread Matthew Knepley
On Sat, Oct 20, 2018 at 11:47 PM Jed Brown wrote: > We're currently maintaining a (somewhat clumsy and limited, albeit > functional) package manager in the form of the --download functionality. > What do people think of outsourcing this to Spack, which would handle > dependency resolution and rep

[petsc-dev] I just removed Fortran support from my development build

2018-10-22 Thread Matthew Knepley
Jenkins will now catch bad Fortran bindings. However, this is a larger problem. GFortran takes forever to build the bindings. Is it similar for other compilers? Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results

Re: [petsc-dev] I just removed Fortran support from my development build

2018-10-22 Thread Matthew Knepley
On Mon, Oct 22, 2018 at 10:58 AM Dener, Alp wrote: > > > On Oct 22, 2018, at 8:18 AM, Matthew Knepley wrote: > > Jenkins will now catch bad Fortran bindings. However, this is a larger > problem. GFortran takes forever to build the bindings. Is it similar for > other c

Re: [petsc-dev] Testing for deprecated MKL functions?

2018-10-23 Thread Matthew Knepley
On Tue, Oct 23, 2018 at 2:20 PM Mills, Richard Tran wrote: > Barry, Satish, and others, > > In MKL 2018 update 2, the old, non-inspector executor sparse BLAS > interfaces were deprecated in favor of the newer, inspector-executor > versions. To keep lots of warnings from being generated, I put in

Re: [petsc-dev] plex test 26 failing for complex

2018-10-24 Thread Matthew Knepley
On Wed, Oct 24, 2018 at 9:54 AM Hapla Vaclav wrote: > Tests dm_impls_plex_tests-ex26% fail for exodusii+complex build. > > E.g. > not ok dm_impls_plex_tests-ex26_15 > # [0]PETSC ERROR: - Error Message > -- > #

[petsc-dev] Problem with DMKSP/DMSNES and GMG

2018-10-26 Thread Matthew Knepley
I am having a problem, and every solution I can think of would seem to do violence to encapsulation or to DMKSP/DMSNES. I want to have GMG automatically work on the velocity part of Stokes. So I give -fieldsplit_velocity_pc_type mg in SNES ex62. This complains that the KSP has a DM but KSPSetC

Re: [petsc-dev] Problem with DMKSP/DMSNES and GMG

2018-10-26 Thread Matthew Knepley
On Fri, Oct 26, 2018 at 11:06 AM Jed Brown wrote: > Matthew Knepley writes: > > > I am having a problem, and every solution I can think of would seem to do > > violence to encapsulation or to DMKSP/DMSNES. > > > > I want to have GMG automatically work on the v

Re: [petsc-dev] Problem with DMKSP/DMSNES and GMG

2018-10-26 Thread Matthew Knepley
On Fri, Oct 26, 2018 at 11:24 AM Jed Brown wrote: > Matthew Knepley writes: > > > On Fri, Oct 26, 2018 at 11:06 AM Jed Brown wrote: > > > >> Matthew Knepley writes: > >> > >> > I am having a problem, and every solution I can think of would see

Re: [petsc-dev] Problem with DMKSP/DMSNES and GMG

2018-10-26 Thread Matthew Knepley
On Fri, Oct 26, 2018 at 12:13 PM Jed Brown wrote: > Matthew Knepley writes: > > > On Fri, Oct 26, 2018 at 11:24 AM Jed Brown wrote: > > > >> Matthew Knepley writes: > >> > >> > On Fri, Oct 26, 2018 at 11:06 AM Jed Brown wrote: > >>

Re: [petsc-dev] error with --download-amgx

2019-12-03 Thread Matthew Knepley
Could you have run out of disk space? /usr/bin/ranlib: libamgx_base.a: Input/output error gmake[2]: *** [base/libamgx_base.a] Error 1 gmake[2]: *** Deleting file `base/libamgx_base.a' gmake[1]: *** [base/CMakeFiles/amgx_base.dir/all] Error 2 gmake: *** [all] Error 2 Matt On Tue, Dec 3, 2019 a

[petsc-dev] Valgrind problems

2019-12-07 Thread Matthew Knepley
I am trying to clean up valgrind errors. However this one dm_impls_plex_tests-ex1_fluent_2 is valgrind clean on my machine. Does anyone get it to output something? Thanks, Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interest

Re: [petsc-dev] Valgrind problems

2019-12-07 Thread Matthew Knepley
e. > I wonder if it was a not updating master thing: https://gitlab.com/petsc/petsc/-/jobs/372455942 Thanks, Matt >Thanks > > > > On Dec 7, 2019, at 2:45 PM, Matthew Knepley wrote: > > > > I am trying to clean up valgrind errors. However this

Re: [petsc-dev] Valgrind problems

2019-12-07 Thread Matthew Knepley
gt; ==4990==by 0x588492F: DMPlexCreateFluentFromFile (plexfluent.c:30) > ==4990==by 0x577E2B7: DMPlexCreateFromFile (plexcreate.c:3254) > ==4990==by 0x10BEDF: CreateMesh (ex1.c:170) > ==4990==by 0x10A20B: main (ex1.c:430) > ==4990== Uninitialised value was created by a stack

Re: [petsc-dev] Valgrind problems

2019-12-07 Thread Matthew Knepley
; > ==4990==by 0x588147D: DMPlexCreateFluent_ReadSection > (plexfluent.c:105) > > ==4990==by 0x5882E3A: DMPlexCreateFluent (plexfluent.c:246) > > ==4990==by 0x588492F: DMPlexCreateFluentFromFile (plexfluent.c:30) > > ==4990==by 0x577E2B7: DMPlexCreateFromFile (plex

[petsc-dev] Testing a new configure addition

2019-12-19 Thread Matthew Knepley
How do I test a new package added to configure? Thanks, Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buffalo.edu/~knepley/

Re: [petsc-dev] Testing a new configure addition

2019-12-19 Thread Matthew Knepley
tests configure > scripts are best updated in the same branch/MR > Thanks! Matt > Satish > > On Thu, 19 Dec 2019, Matthew Knepley wrote: > > > How do I test a new package added to configure? > > > > Thanks, > > > > Matt > > > >

[petsc-dev] Failure in Master

2020-01-01 Thread Matthew Knepley
Does anyone else see this? not ok diff-ts_tutorials_power_grid-ex9opt_2 # Error code: 1 5631 TEST arch-ci-mswin-uni/tests/counts/ts_tutorials-ex18_adv_delta_yper_fim_0.counts 5632 # 4c4 5633 # < Testing hand-coded Gradient, if (for double precision runs) ||G - Gfd||_F/||G||_F is 5634 # -

Re: [petsc-dev] Bash tool for handling switching between multiple PETSc branches and configures

2020-01-03 Thread Matthew Knepley
On Fri, Jan 3, 2020 at 3:52 PM Smith, Barry F. wrote: > >A long time ago Oana suggested a tool that allowed switching between > PETSc configurations after pulls etc that didn't require waiting to > recompile code or rerun configure. Based on her ideas I finally got > something that has been b

Re: [petsc-dev] Bash tool for handling switching between multiple PETSc branches and configures

2020-01-03 Thread Matthew Knepley
> I didn't realize there were any active developers not using ccache. We > should add a note to the manual. > > > > On Jan 3, 2020 15:04, "Smith, Barry F." wrote: > > > > > > > On Jan 3, 2020, at 3:07 PM, Matthew Knepley > wrote: >

[petsc-dev] Weird errors with enum in 64-bit ints

2020-01-13 Thread Matthew Knepley
I have some unexplainable errors with 64-bit ints in this branch: https://gitlab.com/petsc/petsc/-/jobs/400852717 from https://gitlab.com/petsc/petsc/pipelines/108833203 I have built my own 64-bit arches, and the one compiled with new gcc works fine and is valgrind clean. With the bad old A

Re: [petsc-dev] Weird errors with enum in 64-bit ints

2020-01-13 Thread Matthew Knepley
NL disabled my login. Matt > Barry > > > > On Jan 13, 2020, at 3:52 AM, Matthew Knepley wrote: > > > > I have some unexplainable errors with 64-bit ints in this branch: > > > > https://gitlab.com/petsc/petsc/-/jobs/400852717 > > > > from >

Re: [petsc-dev] ccache tips?

2020-01-17 Thread Matthew Knepley
I configured MPI with CC="ccache gcc" knepley/feature-dm-remove-hybrid *$:/PETSc3/petsc/petsc-dev$ /PETSc3/petsc/bin/mpicc -show /Users/knepley/MacSoftware/bin/ccache gcc -Qunused-arguments -fstack-protector -Qunused-arguments -g3 -Wl,-flat_namespace -I/PETSc3/petsc/include -L/PETSc3/petsc/lib -lp

Re: [petsc-dev] how do see artifacts?

2020-01-20 Thread Matthew Knepley
Try Firefox or Chrome? Matt On Sun, Jan 19, 2020 at 11:41 PM Smith, Barry F. via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > > > > On Jan 19, 2020, at 10:39 PM, Jed Brown wrote: > > > > "Smith, Barry F. via petsc-dev" writes: > > > >> I have no left hand side. > > > > Right hand side? I c

Re: [petsc-dev] Issue with dependency generation in GCC 9

2020-01-20 Thread Matthew Knepley
On Mon, Jan 20, 2020 at 3:06 AM Lisandro Dalcin wrote: > I'm using GCC 9 from Fedora 31. A few things have changed regarding > dependency generation, and things are now broken for PETSc. I do not have a > clean solution to propose. > > The problem is the way targets are emitted : > > $ head -n 3

Re: [petsc-dev] Issue with dependency generation in GCC 9

2020-01-20 Thread Matthew Knepley
On Mon, Jan 20, 2020 at 7:06 AM Lisandro Dalcin wrote: > On Mon, 20 Jan 2020 at 14:14, Matthew Knepley wrote: > >> On Mon, Jan 20, 2020 at 3:06 AM Lisandro Dalcin >> wrote: >> >>> I'm using GCC 9 from Fedora 31. A few things have changed regarding >>

Re: [petsc-dev] Issue with dependency generation in GCC 9

2020-01-20 Thread Matthew Knepley
On Mon, Jan 20, 2020 at 7:22 AM Matthew Knepley wrote: > On Mon, Jan 20, 2020 at 7:06 AM Lisandro Dalcin wrote: > >> On Mon, 20 Jan 2020 at 14:14, Matthew Knepley wrote: >> >>> On Mon, Jan 20, 2020 at 3:06 AM Lisandro Dalcin >>> wrote: >>> >

Re: [petsc-dev] Issue with dependency generation in GCC 9

2020-01-20 Thread Matthew Knepley
On Mon, Jan 20, 2020 at 9:58 AM Lisandro Dalcin wrote: > On Mon, 20 Jan 2020 at 17:07, Jed Brown wrote: > >> From my man page (contradicting the behavior Lisandro observed): >> >>-MD -MD is equivalent to -M -MF file, except that -E is not implied. >> The driver >>determines file base

Re: [petsc-dev] Best DM For Structured Atom Layers

2020-01-28 Thread Matthew Knepley
On Tue, Jan 28, 2020 at 2:32 PM Jacob Faibussowitsch wrote: > Hello All, > > I am currently working on rewriting an existing code into PETSc. This > particular code performs tight binding and works on a tiled hexagonal > carbon lattice as a mesh. I am wondering what the best DM format is to > rep

Re: [petsc-dev] Fwd: [SIAM-CSE] Introducing hIPPYlib, a python-based inverse problems solver library

2020-02-05 Thread Matthew Knepley
Yes. I had a discussion with Noemi about using it with PETSc. The big requirement is formulating the discretization in FEniCS. I think it might be possible to peel back one layer of interface and use it directly. Once I find the right student, I will have them investigate. There are some very nice

Re: [petsc-dev] First call to cudaMalloc or cudaFree is very slow on summit

2020-02-07 Thread Matthew Knepley
On Fri, Feb 7, 2020 at 1:23 PM Zhang, Hong via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > Hi all, > > Previously I have noticed that the first call to a CUDA function such as > cudaMalloc and cudaFree in PETSc takes a long time (7.5 seconds) on summit. > Then I prepared a simple example as attach

Re: [petsc-dev] First call to cudaMalloc or cudaFree is very slow on summit

2020-02-08 Thread Matthew Knepley
On Sat, Feb 8, 2020 at 4:34 PM Zhang, Hong via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > I did some further investigation. The overhead persists for both the PETSc > shared library and the static library. In the previous example, it does not > call any PETSc function, the first CUDA function bec

Re: [petsc-dev] DMPlex + P_1 FE interpolation

2020-02-09 Thread Matthew Knepley
On Sun, Feb 9, 2020 at 3:23 PM Pierre Jolivet wrote: > Hello, > I’ve a hard time answering the following DMPlex questions by just looking > at some of the examples and manual. > Considering two DMPlex dm and dma, as in > petsc/src/dm/impls/plex/examples/tests/ex19.c, I’d like to interpolate a > s

Re: [petsc-dev] DMPlex + P_1 FE interpolation

2020-02-10 Thread Matthew Knepley
On Sun, Feb 9, 2020 at 11:11 PM Pierre Jolivet wrote: > > > On 10 Feb 2020, at 6:20 AM, Matthew Knepley wrote: > > On Sun, Feb 9, 2020 at 3:23 PM Pierre Jolivet > wrote: > >> Hello, >> I’ve a hard time answering the following DMPlex questions by just looki

Re: [petsc-dev] DMPlex + P_1 FE interpolation

2020-02-10 Thread Matthew Knepley
On Mon, Feb 10, 2020 at 5:08 AM Dave May wrote: > > > On Mon 10. Feb 2020 at 13:09, Matthew Knepley wrote: > >> On Sun, Feb 9, 2020 at 11:11 PM Pierre Jolivet < >> pierre.joli...@enseeiht.fr> wrote: >> >>> >>> >>> On 10 Feb 2020,

Re: [petsc-dev] DMPlex + P_1 FE interpolation

2020-02-10 Thread Matthew Knepley
On Mon, Feb 10, 2020 at 5:31 AM Dave May wrote: > > > On Mon 10. Feb 2020 at 14:17, Matthew Knepley wrote: > >> On Mon, Feb 10, 2020 at 5:08 AM Dave May wrote: >> >>> >>> >>> On Mon 10. Feb 2020 at 13:09, Matthew Knepley wrote: >

Re: [petsc-dev] DMPlex + P_1 FE interpolation

2020-02-10 Thread Matthew Knepley
On Mon, Feb 10, 2020 at 6:49 AM Pierre Jolivet wrote: > > > On 10 Feb 2020, at 1:08 PM, Matthew Knepley wrote: > > On Sun, Feb 9, 2020 at 11:11 PM Pierre Jolivet > wrote: > >> >> >> On 10 Feb 2020, at 6:20 AM, Matthew Knepley wrote: >> >

Re: [petsc-dev] First call to cudaMalloc or cudaFree is very slow on summit

2020-02-12 Thread Matthew Knepley
2020, at 11:14 AM, Smith, Barry F. wrote: > > > gprof or some similar tool? > > > On Feb 10, 2020, at 11:18 AM, Zhang, Hong via petsc-dev < > petsc-dev@mcs.anl.gov> wrote: > > -cuda_initialize 0 does not make any difference. Actually this issue has > nothin

Re: [petsc-dev] First call to cudaMalloc or cudaFree is very slow on summit

2020-02-12 Thread Matthew Knepley
On Wed, Feb 12, 2020 at 11:38 AM Zhang, Hong wrote: > > > On Feb 12, 2020, at 11:09 AM, Matthew Knepley wrote: > > On Wed, Feb 12, 2020 at 11:06 AM Zhang, Hong via petsc-dev < > petsc-dev@mcs.anl.gov> wrote: > >> Sorry for the long post. Here are replies I have

Re: [petsc-dev] DMPlex + P_1 FE interpolation

2020-02-17 Thread Matthew Knepley
oking at the > wrong place and overcomplicating all this? > I think you may be. I am still trying to understand what you want to do. Thanks, Matt > Thank you in advance for your help, > Pierre > > On 10 Feb 2020, at 4:30 PM, Matthew Knepley wrote: > > On Mon, Feb 1

Re: [petsc-dev] Proper matrix size to choose when evaluating MatMult?

2020-02-21 Thread Matthew Knepley
I think Karl goes into these issues here: https://arxiv.org/pdf/1410.4054.pdf Thanks, Matt On Fri, Feb 21, 2020 at 5:58 PM Junchao Zhang via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > > > On Fri, Feb 21, 2020 at 4:38 PM Mark Adams wrote: > >> >> >> On Fri, Feb 21, 2020 at 4:51 PM Juncha

Re: [petsc-dev] Exodus viewer

2020-02-25 Thread Matthew Knepley
On Tue, Feb 25, 2020 at 10:44 AM Blaise A Bourdin wrote: > Hi, > > I just noticed recent changes to the exodus I/O and a new exodus viewer, > which had been on my todo list forever… > Is there a test or example available? Is there anything that remains I > could help with? > I know! It was on mi

Re: [petsc-dev] Unions of IS and MatNest

2020-03-12 Thread Matthew Knepley
On Thu, Mar 12, 2020 at 5:59 PM Jed Brown wrote: > Pierre Jolivet writes: > > > Hello, > > Has there been any follow-up on this > https://lists.mcs.anl.gov/pipermail/petsc-users/2015-January/024020.html < > https://lists.mcs.anl.gov/pipermail/petsc-users/2015-January/024020.html>? > > Given a 3x

Re: [petsc-dev] Unions of IS and MatNest

2020-03-13 Thread Matthew Knepley
On Fri, Mar 13, 2020 at 3:19 AM Pierre Jolivet wrote: > > > On 12 Mar 2020, at 11:40 PM, Matthew Knepley wrote: > > On Thu, Mar 12, 2020 at 5:59 PM Jed Brown wrote: > >> Pierre Jolivet writes: >> >> > Hello, >> > Has there been any follow-u

Re: [petsc-dev] Unions of IS and MatNest

2020-03-13 Thread Matthew Knepley
On Fri, Mar 13, 2020 at 11:20 AM Pierre Jolivet wrote: > > > On 13 Mar 2020, at 2:36 PM, Matthew Knepley wrote: > > On Fri, Mar 13, 2020 at 3:19 AM Pierre Jolivet > wrote: > >> >> >> On 12 Mar 2020, at 11:40 PM, Matthew Knepley wrote: >> >&g

Re: [petsc-dev] print statements in fortran interface functions

2020-03-17 Thread Matthew Knepley
On different OSes, the buffering for printf is different. Do you mean that the program ends correctly and prints nothing, or do you mean that you execute a print but nothing immediately outputs to the screen. We do not explicitly flush, so it could be that this weird OS is buffering in a way you do

Re: [petsc-dev] MatCreateVecs from Fortran problem

2020-03-17 Thread Matthew Knepley
On Tue, Mar 17, 2020 at 1:34 PM Mark Adams wrote: > Our code broke in moving from v3.7 to current. The problem seems to be in > MatCreateVecs > > Our code has: > > Vec::XVec > Vec::BVec > > this%xVec2 = PETSC_NULL_VEC > this%bVec2 = PETSC_NULL_VEC > Passing NULL to MatCreateVecs() means th

Re: [petsc-dev] MatCreateVecs from Fortran problem

2020-03-17 Thread Matthew Knepley
On Tue, Mar 17, 2020 at 2:16 PM Mark Adams wrote: > > >> Passing NULL to MatCreateVecs() means that you do not want a vector out: >> >> > Yes, I want both vectors out and we pass it Vec that have been initialized > with PETSC_NULL_VEC > That is wrong. > > >> >> https://www.mcs.anl.gov/petsc/pe

Re: [petsc-dev] MatCreateVecs from Fortran problem

2020-03-17 Thread Matthew Knepley
On Tue, Mar 17, 2020 at 2:48 PM Mark Adams wrote: > > > On Tue, Mar 17, 2020 at 2:25 PM Matthew Knepley wrote: > >> On Tue, Mar 17, 2020 at 2:16 PM Mark Adams wrote: >> >>> >>> >>>> Passing NULL to MatCreateVecs() means that you do not wan

  1   2   3   4   5   6   7   8   9   10   >