To me, it looks fine.
Matt
On Tue, Sep 5, 2017 at 7:26 PM, Tobin Isaac wrote:
> I just came across Random123 [1]: portable, parallel, high-quality
> pseudorandom number generators. It's developed by DE Shaw, but the license
> [2] looks to me like we should be able to make a repo for it and
In my opinion, the software engineering problem with ADIFOR(C) is that they
are designed to run on whole programs. What
we really want is for small pieces of code, like routines for element
vectors, to be differentiated. Once you get whole programs,
you need to parse system headers, etc. which is a
On Fri, Sep 22, 2017 at 12:06 AM, Richard Tran Mills
wrote:
> Thanks for sharing this, Barry. I haven't had time to read their paper,
> but it looks worth a read.
>
> Hong, since many machine-learning or data-mining problems can be cast as
> linear algebra problems (several examples involving eig
On Thu, Sep 21, 2017 at 3:51 PM, Zhang, Hong wrote:
> Great news! According to their papers, MLSVM works only in serial. I am
> not sure what is stopping them using PETSc in parallel.
>
I think its because they use FLANN (http://www.cs.ubc.ca/research/flann/)
which appears to be serial.
Richard
I did not even make the Top 10. Sad.
Matt
On Fri, Sep 22, 2017 at 9:01 PM, Lisandro Dalcin wrote:
> 32 CharacteristicPackageInitialized
> 32 CharacteristicSendCoordinatesEnd
> 32 DMForestGetBaseCoordinateMapping
> 32 DMForestSetBaseCoordinateMapping
> 32 DMGlobalToLocalBeginDefaultShell
> 32
On Sun, Sep 24, 2017 at 8:57 PM, Adrian Croucher
wrote:
>
>
> On 25/09/17 13:37, Adrian Croucher wrote:
>
>>
>> Actually it looks like DMGetDefaultSF() does already work in Fortran
>> though, and could be used to do the same thing.
>>
>> Do you think that would be a reasonable way to do it?
>>
>>
On Mon, Sep 25, 2017 at 11:39 AM, Kong, Fande wrote:
> Hi All,
>
> A matrix is created with the right preallocation, and then MatAssembly is
> called. The preallocation info will be removed. We insert any values then,
> and will encounter an malloc error.
>
> My question is that we was intending
On Tue, Sep 26, 2017 at 11:49 PM, Adrian Croucher wrote:
> hi Matt,
>
> On 25/09/17 23:12, Matthew Knepley wrote:
>
>
> If you truly need the exact same SF for your grid, you should be able to
> use DMGet/SetPointSF() since it will just reference count it for you. The
On Sat, Oct 7, 2017 at 6:26 AM, Barry Smith wrote:
>
>We'll be preparing another release in the next couple of months. It
> would be nice to transition more examples to the new test harness.
>
>How can we organize this? My original plan was to try to do a
> test/tutorial directory at a ti
On Sat, Oct 7, 2017 at 9:08 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Sat, Oct 7, 2017 at 6:26 AM, Barry Smith wrote:
> >
> >>
> >>We'll be preparing another release in the next couple of months. It
> >> would be nice t
On Sat, Oct 7, 2017 at 9:41 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > This seems masochistic. I want to run a test with a different number
> > of processes. So my best option is print out the options, cut & paste
> > that soewhere, alter it to what I want
On Mon, Oct 9, 2017 at 3:36 PM, Scott Kruger wrote:
>
> It's already there -- I fixed your pull request:
>
Great! Thanks Scott.
Matt
> gabrielle 70: git blame gmakefile | grep NP
> b95da926fad (Matthew Knepley 2017-06-14 17:33:01 -0600 103) ifdef NP
> b95da926fad (
On Tue, Oct 10, 2017 at 8:25 AM, Patrick Sanan
wrote:
> I am reading some of the DMPlex docs and would like to take the
> opportunity to clean up things I notice. There are several places in the
> man pages and source that "sieve" is mentioned. As far as I know this is
>
- the name of the project
On Tue, Oct 10, 2017 at 9:48 AM, Satish Balay wrote:
> I never understood --host option of autoconf - and it would rarely
> work for me.
>
I do not understand this option either. Jed, what does it do?
Matt
> For packages where '--with-batch --download-package' did not work -
> our suggstion
On Tue, Oct 10, 2017 at 10:36 AM, Kong, Fande wrote:
> Hi Developers,
>
> I am trying to use the native LU in PETSc to compute a matrix (arising
> from a fluid-structure interaction problem). I encounter a zero pivot. The
> matrix gives the wrong Newton direction when using "-pc_factor_shift_type
On Wed, Oct 11, 2017 at 4:54 AM, Patrick Sanan
wrote:
> Another terminology question to help with the docs.
>
> What's the origin of the term "chart" in DMPlex? I'm only previously
> familiar with the term in the context of manifolds (where chart =
> homeomorphism from an open set in the manifold
On Tue, Oct 10, 2017 at 9:07 PM, Adrian Croucher
wrote:
> hi
>
> On 03/08/17 13:10, Matthew Knepley wrote:
>
> On Wed, Aug 2, 2017 at 8:06 PM, Adrian Croucher > wrote:
>
>> I haven't tried this myself, but two of my colleagues are reporting
>> problems acce
On Wed, Oct 11, 2017 at 10:55 AM, Patrick Sanan
wrote:
> On Wed, Oct 11, 2017 at 10:58 AM, Matthew Knepley
> wrote:
>
>> On Wed, Oct 11, 2017 at 4:54 AM, Patrick Sanan
>> wrote:
>>
>>> Another terminology question to help with the docs.
>>>
>&g
you want the wrapper to take in F90 pointer arguments, then you have to
declare it or you get an SEGV. These
get autogenerated in include/petsc/finclude/ftn-auto/*.h90 when you run
'make allfortranstubs'. Did this happen
for your function?
Thanks,
Matt
>
> - Adrian
>
>
On Fri, Oct 13, 2017 at 10:26 AM, Vaclav Hapla
wrote:
> Hello
>
> In DMPlexDistribute, when it is run on 1 process and sf != NULL, the
> output *sf is undefined. This in turn leads to crash of PetscSFDestroy,
> e.g. in src/dm/impls/plex/examples/tutorials/ex5.c:45.
>
> I created the corresponding
Someone had to do it.
I will not try to frame the entire discussion. Barry has already thrown
down the "show me your interface" gauntlet. However, I want to emphasize
one point that may have been lost in the prior discussion. Every example I
have looked at so far is focused on the reduced space fo
> the tangent linear model), reuses the callbacks set by TSSetGradientDAE and
> TSSetGradientIC. Hong, do you also need to integrate some quadrature
> variable in your TSForward code?
>
For people like me who do not have all the terminology down, is the Tangent
Linear Model (TLM) just
On Sun, Oct 15, 2017 at 2:03 PM, Jed Brown wrote:
> Barry Smith writes:
>
> >> I don't see how any of the present interfaces work for waveform
> >> relaxation. I also think that is rarely a desirable technique -- too
> >> many awkward limitations. My recollection is that Borzi only uses it
> >
On Sun, Oct 15, 2017 at 2:42 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > I don't think the full space method is of questionable practically in the
> > time-independent case.
> > This is still important in this discussion because the way we talk about
On Sun, Oct 15, 2017 at 6:54 PM, Adrian Croucher
wrote:
> hi
>
> On 14/10/17 00:47, Matthew Knepley wrote:
>
>
> If you want the wrapper to take in F90 pointer arguments, then you have to
> declare it or you get an SEGV. These
> get autogenerated in include/petsc/finclude
On Sun, Oct 15, 2017 at 9:13 PM, Adrian Croucher
wrote:
> On 16/10/17 13:12, Matthew Knepley wrote:
>
>
> Okay, I was wrong about what gets produced. Barry rewrote the Fortran
> bindings last, and I did not understand how everything
> was put together. You will need to make
We had a previous error with pdgssvx in SuperLU I think. Maybe searching
petsc-maint would get it?
Matt
On Mon, Oct 16, 2017 at 12:21 PM, Mark Adams wrote:
> I just ran this and have a little bit of a stack trace. This is on my
> laptop and MPI can be a little flaky here (eg, IBarrier does n
On Tue, Oct 17, 2017 at 8:16 AM, Jed Brown wrote:
> Stefano Zampini writes:
>
> >> >
> >> >
> >> >> In case of multiple objectives, there may be a performance reason to
> >> >> amortize evaluation of several at once, though the list interface is
> >> >> convenient. Consider common objectives be
On Tue, Oct 17, 2017 at 8:51 AM, Jed Brown wrote:
> Matthew Knepley writes:
>
> >> It's a recipe for confusion. Either the parameters are never passed
> >> explicitly or they are always passed explicitly and should not be stored
> >> redundantly in the co
I already pushed a fix this morning. What did you push?
Matt
On Tue, Oct 17, 2017 at 3:59 PM, Satish Balay wrote:
> This has been broken for a couple of days. I pushed a fix [hoping
> builds don't abort anymore] - but there are more things to fix..
>
> Satish
>
> On Tue, 17 Oct 2017, PETSc ch
ng
>
> >>>>>>>>
> /sandbox/petsc/petsc.clone/src/snes/utils/ftn-auto/dmadaptf.c: In
> function ‘petscglobalminmax_’:
> /sandbox/petsc/petsc.clone/src/snes/utils/ftn-auto/dmadaptf.c:125:1:
> warning: implicit declaration of function ‘PetscGlobalMinMax’
On Tue, Oct 17, 2017 at 11:35 PM, Adrian Croucher wrote:
> hi
>
> On 16/10/17 15:16, Jed Brown wrote:
>
>> Adrian Croucher writes:
>>
>> So do you think the SF stuff in petsc/finclude/petscis.h should be taken
>>> out, and put into a new petsc/finclude/petscsf.h ?
>>>
>> I think that's desirable
I am confused by some of the terminology in TS. At the top level, IFunction
appears to mean the entire equation
F(u, u_t, x) = 0
However, this appears to mean something different than ifunction at the
function pointer level. Inside TSCompiteIFunction(), it uses both ifunction
and rhsfunction. T
On Fri, Oct 20, 2017 at 9:23 AM, Emil Constantinescu
wrote:
> On 10/20/17 7:57 AM, Matthew Knepley wrote:
>
>> I am confused by some of the terminology in TS. At the top level,
>> IFunction appears to mean the entire equation
>>
>>F(u, u_t, x) = 0
>>
>
On Fri, Oct 20, 2017 at 11:22 AM, Emil Constantinescu
wrote:
> On 10/20/17 9:11 AM, Matthew Knepley wrote:
>
>> On Fri, Oct 20, 2017 at 9:23 AM, Emil Constantinescu <
>> emcon...@mcs.anl.gov <mailto:emcon...@mcs.anl.gov>> wrote:
>>
>>
On Fri, Oct 20, 2017 at 11:58 AM, Emil Constantinescu
wrote:
> On 10/20/17 10:43 AM, Matthew Knepley wrote:
>
>> On Fri, Oct 20, 2017 at 11:22 AM, Emil Constantinescu <
>> emcon...@mcs.anl.gov <mailto:emcon...@mcs.anl.gov>> wrote:
>>
>>
On Fri, Oct 20, 2017 at 12:34 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Fri, Oct 20, 2017 at 11:58 AM, Emil Constantinescu <
> emcon...@mcs.anl.gov>
> > wrote:
> >
> >> On 10/20/17 10:43 AM, Matthew Knepley wrote:
> >>
> >
On Fri, Oct 20, 2017 at 1:20 PM, Jed Brown wrote:
> Barry Smith writes:
>
> >> On Oct 20, 2017, at 11:53 AM, Jed Brown wrote:
> >>
> >> Barry Smith writes:
> >>
> >>> The name absolutely has to be changed. But to what? And the manual
> page is WRONG! You cannot justify that no matter how muc
On Fri, Oct 20, 2017 at 4:32 PM, Blaise A Bourdin wrote:
> Hi,
>
> Recent versions of exodus are implemented with netcdf4, which is based on
> hdf5.
> Instead of going through the process of rewriting exodus viewers, I would
> like to use PETSc’s HDF5 I/O operations.
> Problem is that as far as I
On Fri, Oct 20, 2017 at 6:36 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > I think you are right that the SNES thing needs fixed as well. The
> > problem has not come up because no one ever uses that interface. It is
> > only there because it makes FAS a lot ea
On Fri, Oct 20, 2017 at 6:44 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Fri, Oct 20, 2017 at 6:36 PM, Jed Brown wrote:
> >
> >> Matthew Knepley writes:
> >>
> >> > I think you are right that the SNES thing needs fixed as well. T
On Fri, Oct 20, 2017 at 4:21 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> >> Note that TSComputeIFunction is very much like SNESComputeFunction,
> >> which includes
> >>
> >> if (snes->vec_rhs) {
> >> ierr = VecAXPY(y,-1.0,snes-&
On Sat, Oct 21, 2017 at 9:28 AM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Fri, Oct 20, 2017 at 4:21 PM, Jed Brown wrote:
> >
> >> Matthew Knepley writes:
> >>
> >> >> Note that TSComputeIFunction is very much like SNESComputeFu
On Tue, Oct 24, 2017 at 4:22 PM, Adrian Croucher
wrote:
> hi,
>
> On 11/10/17 22:18, Matthew Knepley wrote:
>
> On Tue, Oct 10, 2017 at 9:07 PM, Adrian Croucher <
> a.crouc...@auckland.ac.nz> wrote:
>
>> hi
>>
>> On 03/08/17 13:10, Matthew Kneple
TSGetSNES()
SNESGetIterationNumber()
Or do you want something else?
Matt
On Wed, Oct 25, 2017 at 7:59 AM, Mark Adams wrote:
> I want to modify the TS time step, in a post-step function, and would
> like to get the number of Newton iterations that were used in the time
> step. I am not seei
On Sat, Oct 28, 2017 at 7:11 AM, Franck Houssen
wrote:
> I just added that before EPSSetOperators:
> PetscViewer viewerA; PetscViewerBinaryOpen(PETSC_
> COMM_WORLD,"Atau.out",FILE_MODE_WRITE,&viewerA); MatView(A,viewerA);
> PetscViewer viewerB; PetscViewerBinaryOpen(PETSC_
> COMM_WORLD,"Btau.out"
On Tue, Oct 31, 2017 at 11:36 AM, Blaise A Bourdin wrote:
> Hi,
>
> I am losing my mind over segfaults in DMPlexVecGetClosure.
>
> In the example below, I get a segfault whenever cval is not set to NULL
> before calling DMPlexVecGetClosure. Is this expected?
>
Yes. If cval is NULL, it allocates
On Wed, Nov 1, 2017 at 6:50 PM, Blaise A Bourdin wrote:
> Hi,
>
> I have just spent 2h helping a student debugging a code, and I think that
> the problem is not ours…
> See the attached example: 1 create 1 DM and 2 SNES.
> If I assign the same DM to both SNES, the function and Jacobean of the
> s
On Thu, Nov 2, 2017 at 12:08 PM, Vaclav Hapla
wrote:
> It seems that DMPlexCreateMedFromFile leaves out some mesh elements. I
> found it out when investigating why ParMetis redistribution crashes.
>
> I attach the datafile
> $PETSC_DIR/share/petsc/datafiles/meshes/blockcylinder-50.exo
> converte
On Sun, Nov 5, 2017 at 8:57 AM, Franck Houssen
wrote:
> Hello,
>
> Does PETSc provide support for CMake ? I would say this is an on-going
> feature that is not yet fully working, right ?
>
> My understanding is that "./configure" generates the CMakeLists.txt. But
> then, when using cmake, make fa
On Sun, Nov 5, 2017 at 9:37 AM, Franck Houssen
wrote:
> --
>
> *De: *"Matthew Knepley"
> *À: *"Franck Houssen"
> *Cc: *"For users of the development version of PETSc" <
> petsc-dev@mcs.anl.gov>
> *Envoyé: *Di
On Thu, Nov 2, 2017 at 12:08 PM, Vaclav Hapla
wrote:
> It seems that DMPlexCreateMedFromFile leaves out some mesh elements. I
> found it out when investigating why ParMetis redistribution crashes.
>
> I attach the datafile
> $PETSC_DIR/share/petsc/datafiles/meshes/blockcylinder-50.exo
> converte
On Mon, Nov 6, 2017 at 8:09 AM, Vaclav Hapla
wrote:
> Hello
>
> The whole PetscPartitioner class sources are part of
> src/dm/impls/plex/plexpartition.c,
> mixed together with some DMPlex* functions.
> If you don't mind, I would move the PetscPartitioner* stuff into the
> separate file petscpart
On Mon, Nov 6, 2017 at 8:24 AM, Smith, Barry F. wrote:
>
>Vaclav,
>
> Actually you should not just do this! PETSc already has a full class
> for managing partitioning (that Matt ignored for no good reason)
Obviously, a good reason existed. All sorts of horrible Mat-specific crap
was i
On Mon, Nov 6, 2017 at 8:34 AM, Smith, Barry F. wrote:
> > On Nov 6, 2017, at 7:27 AM, Matthew Knepley wrote:
> >
> > On Mon, Nov 6, 2017 at 8:24 AM, Smith, Barry F.
> wrote:
> >
> >Vaclav,
> >
> > Actually you should not just do this!
On Mon, Nov 6, 2017 at 12:49 PM, Franck Houssen
wrote:
> - Mail original -
> > De: "Satish Balay"
> > À: "Franck Houssen"
> > Cc: "petsc-dev"
> > Envoyé: Lundi 6 Novembre 2017 16:08:55
> > Objet: Re: [petsc-dev] CMake: make, install, find_package ?
> >
> > On Mon, 6 Nov 2017, Franck Ho
oken?
>
It appears to be an MED incompatibility
https://bitbucket.org/petsc/petsc/src/17bd883d72f40a596f2d89b5afda5a233b621464/src/dm/impls/plex/plexmed.c?at=master&fileviewer=file-view-default#plexmed.c-86
We should ask Nico about it.
Thanks,
Matt
> Thanks
>
> Vaclav
>
>
> 5. 1
On Wed, Nov 8, 2017 at 4:52 AM, Vaclav Hapla
wrote:
>
> > 8. 11. 2017 v 9:06, Lisandro Dalcin :
> >
> > On 8 November 2017 at 05:51, Smith, Barry F. wrote:
> >>
> >>> On Nov 7, 2017, at 1:33 AM, Lisandro Dalcin wrote:
> >>>
> >>> The only concern I have about PetscPartitioner is that the API de
On Wed, Nov 8, 2017 at 4:52 AM, Vaclav Hapla
wrote:
>
> > 8. 11. 2017 v 9:06, Lisandro Dalcin :
> >
> > On 8 November 2017 at 05:51, Smith, Barry F. wrote:
> >>
> >>> On Nov 7, 2017, at 1:33 AM, Lisandro Dalcin wrote:
> >>>
> >>> The only concern I have about PetscPartitioner is that the API de
On Wed, Nov 8, 2017 at 9:14 AM, Vaclav Hapla
wrote:
>
> 8. 11. 2017 v 14:52, Jed Brown :
>
> Matthew Knepley writes:
>
> On Wed, Nov 8, 2017 at 4:52 AM, Vaclav Hapla
> wrote:
>
>
> 8. 11. 2017 v 9:06, Lisandro Dalcin :
>
> On 8 November 2017 at 05:51, Sm
On Wed, Nov 8, 2017 at 10:04 AM, Vaclav Hapla
wrote:
>
> 8. 11. 2017 v 15:59, Matthew Knepley :
>
> On Wed, Nov 8, 2017 at 9:14 AM, Vaclav Hapla
> wrote:
>
>>
>> 8. 11. 2017 v 14:52, Jed Brown :
>>
>> Matthew Knepley writes:
>>
>>
On Wed, Nov 8, 2017 at 11:21 AM, Munson, Todd wrote:
>
> I had a question. There is VecGetArrayRead and VecRestoreArrayRead.
> However, the
> sameconcept does not appear with IS. Is there a reason for not having
> ISGetIndicesRead and ISRestoreIndicesRead?
>
Actually, ISGetIndices() is always
On Wed, Nov 8, 2017 at 1:02 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Wed, Nov 8, 2017 at 10:04 AM, Vaclav Hapla >
> > wrote:
> >
> >>
> >> 8. 11. 2017 v 15:59, Matthew Knepley :
> >>
> >> On Wed, Nov 8, 2017 at 9:14
On Wed, Nov 8, 2017 at 1:49 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> >> > No, this is the right structure.
> >>
> >> Oh come on. You're defending a quadratic algorithm.
> >>
> >> ierr = ParMETIS_V3_PartKway(vtxdist,
On Wed, Nov 8, 2017 at 2:09 PM, Munson, Todd wrote:
>
> > On Nov 8, 2017, at 10:32 AM, Matthew Knepley wrote:
> >
> > On Wed, Nov 8, 2017 at 11:21 AM, Munson, Todd
> wrote:
> >
> > I had a question. There is VecGetArrayRead and VecRestoreArrayRead.
>
On Wed, Nov 8, 2017 at 2:27 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Wed, Nov 8, 2017 at 1:49 PM, Jed Brown wrote:
> >
> >> Matthew Knepley writes:
> >>
> >> >> > No, this is the right structure.
> >>
On Thu, Nov 9, 2017 at 5:53 AM, Vaclav Hapla
wrote:
> Thanks for the discussion. I think we got further.
>
> 8. 11. 2017 v 20:42, Matthew Knepley :
>
> On Wed, Nov 8, 2017 at 2:27 PM, Jed Brown wrote:
>
>> Matthew Knepley writes:
>>
>> > On Wed, N
On Thu, Nov 9, 2017 at 11:15 AM, Vaclav Hapla
wrote:
>
> 7. 11. 2017 v 19:51, Matthew Knepley :
>
> On Tue, Nov 7, 2017 at 11:44 AM, Vaclav Hapla
> wrote:
>
>> Nico Schlömer just introduced MED support into his meshio tool, based on
>> my issue request. When
On Fri, Nov 10, 2017 at 4:50 PM, Satish Balay wrote:
>
> On Fri, 10 Nov 2017, Richard Tran Mills wrote:
>
> > Hi Satish,
> >
> > Thanks for taking the initiative to switch to testing next-tmp to help
> > clear up the constipation with moving things into master. It looks like
> > there hasn't been
On Fri, Nov 10, 2017 at 8:56 AM, Jed Brown wrote:
> Vaclav Hapla writes:
>
> >> 10. 11. 2017 v 5:09, Smith, Barry F. :
> >>
> >>
> >>
> >>> On Nov 8, 2017, at 3:52 AM, Vaclav Hapla
> wrote:
> >>>
> >>>
> 8. 11. 2017 v 9:06, Lisandro Dalcin :
>
> On 8 November 2017 at 05:51, Smith
On Sat, Nov 11, 2017 at 12:37 PM, Satish Balay wrote:
> On Sat, 11 Nov 2017, Matthew Knepley wrote:
>
> > > In the long term - Barry wants to get rid of next..
> >
> >
> > 1) I think next really prevents master from getting screwed up (witness
> > next)
>
On Sat, Nov 11, 2017 at 12:53 PM, Satish Balay wrote:
> On Sat, 11 Nov 2017, Matthew Knepley wrote:
>
> > On Sat, Nov 11, 2017 at 12:37 PM, Satish Balay
> wrote:
> >
> > > On Sat, 11 Nov 2017, Matthew Knepley wrote:
> > >
> > > >
On Sat, Nov 11, 2017 at 1:12 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> >> Matrix and graph are equivalent concepts.
> >
> >
> > This is clearly wrong. A matrix is the coordinate representation of a
> > linear operator, and thus has a
On Sat, Nov 11, 2017 at 1:15 PM, Satish Balay wrote:
> On Sat, 11 Nov 2017, Matthew Knepley wrote:
>
> > >
> > > BTW: Ultimlately if you want to improve current next model - everyone
> > > has to do a 'make alltests DIFF=$PETSC_DIR/bin/petscdiff' fo
On Sat, Nov 11, 2017 at 1:33 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Sat, Nov 11, 2017 at 1:12 PM, Jed Brown wrote:
> >
> >> Matthew Knepley writes:
> >>
> >> >> Matrix and graph are equivalent concepts.
> >> >
On Sat, Nov 11, 2017 at 9:40 PM, Satish Balay wrote:
> [Starting a new thread for this]
>
> Tonight's builds will be on 'next-tmp' with:
>
How easy would it be to setup a system that throttles (or in some way
intelligently selects) the branches
for 'next'? And when a branch graduates to master,
On Sun, Nov 12, 2017 at 10:36 AM, Satish Balay wrote:
> On Sun, 12 Nov 2017, Satish Balay wrote:
>
> > [Starting a new thread for this]
> >
> > Tonight's builds will be on 'next-tmp' with:
> >
> > $ git fetch -p && comm -12 <(git branch -r --merged origin/next-tmp |
> sort) <(git branch -r --no-m
On Sun, Nov 12, 2017 at 11:45 AM, Jed Brown wrote:
> Satish Balay writes:
>
> > On Sun, 12 Nov 2017, Smith, Barry F. wrote:
> >
> >> If the testing/fixing would be faster if we bought more machines
> then decide what machines we need and we order them now. Buying machines is
> far better th
On Sun, Nov 12, 2017 at 12:17 PM, Satish Balay wrote:
> On Sun, 12 Nov 2017, Matthew Knepley wrote:
>
> >
> > Have we tried histogramming test times? It would be nice to know how much
> > cumulative
> > time it takes to run 37%, 67%, 95%, etc.
>
> I'
Do you have any idea what they changed the name of CVSpgmr to in 3.0?
Matt
On Wed, Nov 15, 2017 at 2:53 PM, Antonio Trande
wrote:
> Hi all.
>
> I can't build 'PETSc 3.8.*' with MPI against new 'Sundials-3.0.0':
>
>
> **
On Wed, Nov 15, 2017 at 8:56 PM, Richard Tran Mills wrote:
> On Sun, Nov 12, 2017 at 9:54 AM, Smith, Barry F.
> wrote:
>
>>
>>
>> > On Nov 12, 2017, at 11:21 AM, Matthew Knepley
>> wrote:
>> >
>> > On Sun, Nov 12, 2017 at 12:17 PM, Sat
On Wed, Nov 15, 2017 at 10:48 PM, Smith, Barry F.
wrote:
>
>
> > On Nov 15, 2017, at 8:30 PM, Matthew Knepley wrote:
> >
> > On Wed, Nov 15, 2017 at 8:56 PM, Richard Tran Mills
> wrote:
> > On Sun, Nov 12, 2017 at 9:54 AM, Smith, Barry F.
> wrote:
> &
Satish is the God of Configure Debugging (consult whole thread, all ye who
enter here)
Matt
-- Forwarded message --
From: Satish Balay
Date: Fri, Nov 17, 2017 at 1:06 PM
Subject: Re: [petsc-maint] Could not find a suitable archiver
To: Patrick Welche
Cc: petsc-maint
On Fri
Why does this petsc-maint message have a messed up Reply-To: header?
Matt
-- Forwarded message --
From: Yann Jobic
Date: Wed, Nov 22, 2017 at 12:39 PM
Subject: [petsc-users] 2 Dirichlet conditions for one Element in PetscFE
To: PETSc
Hello,
I've found a strange behavior wh
On Mon, Nov 27, 2017 at 8:08 PM, Jed Brown wrote:
> I don't know the answer to your question (Matt?), but do you really need
> to reorder the entire mesh or would it be sufficient to label your
> points with their original numbering?
Maybe I am wrong, but I think it amounts to the same thing. I
On Mon, Nov 27, 2017 at 9:24 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Mon, Nov 27, 2017 at 8:08 PM, Jed Brown wrote:
> >
> >> I don't know the answer to your question (Matt?), but do you really need
> >> to reorder the entire mes
On Wed, Nov 29, 2017 at 10:38 AM, Scott Kruger wrote:
> rm $PETSC_ARCH/tests//
>
>
> I assume you really mean something else, specifically related to the
> generated test script, but there isn't enough information to answer that.
>
> Or more specifically, the test language was not written to gene
y
> directories are not removed automatically by the test harness as it hasn't
> been a problem yet.
>
Okay, then it would be nice to give a list of files to be removed after the
test.
Thanks,
Matt
> Scott
>
>
> On 11/29/17 9:43 AM, Matthew Knepley wrote:
>
On Fri, Dec 1, 2017 at 11:10 AM, Vaclav Hapla
wrote:
> Hello
>
> I noticed DMPlexCreateFromFile ignores the interpolate flag for HDF5. It's
> the only format for which PetscViewer API is used for loading, and this API
> has no means to specify interpolation. I guess there should be at least
>
>
On Mon, Dec 4, 2017 at 1:06 PM, Matthew Knepley wrote:
> On Fri, Dec 1, 2017 at 11:10 AM, Vaclav Hapla
> wrote:
>
>> Hello
>>
>> I noticed DMPlexCreateFromFile ignores the interpolate flag for HDF5.
>> It's the only format for which PetscViewer API is us
On Thu, Dec 7, 2017 at 10:26 AM, Vaclav Hapla
wrote:
>
>
> > 9. 11. 2017 v 12:53, Matthew Knepley :
> >
> > I think I need to create a proof-of-concept. I would start by employing
> MatPartitioning in PetscPartitionerPartition with anything outside of this
>
On Tue, Dec 5, 2017 at 7:47 PM, Blaise A Bourdin wrote:
> Hi,
>
> I am building torture test for DMPlex with somewhat complicated dof
> layout. In some situations, DMCreateSubDM fails.
>
> DMUAS has 3 fields corresponding to
> U: corresponding to vector-valued P1 or P2 Lagrange FE, Sigma,
> vecto
On Thu, Nov 23, 2017 at 11:15 AM, Vaclav Hapla
wrote:
> Hello
>
> The DMPlex example src/dm/impls/plex/examples/tests/ex18.c fails with
> -cell_simplex 0 :
>
> [0]PETSC ERROR: Null argument, when expecting valid pointer
> [0]PETSC ERROR: Null Object: Parameter # 1
> [0]PETSC ERROR: See http://www
On Mon, Dec 18, 2017 at 7:42 AM, Karin&NiKo wrote:
> Dear PETSc team,
>
> I would like to implement and possibly commit to PETSc the Golub-Kahan
> bidiagonalization algorithm (GK) describe in Arioli's paper :
> http://epubs.siam.org/doi/pdf/10.1137/120866543.
> In this work, Mario Arioli uses GK
that is how we are structuring it. You could imagine that we
break it down further, so
that all the block solvers are separate KSP and just require a
PCFIELDSPLIT, but that seems like
overkill to me.
Thanks,
Matt
> Thanks,
> Nicolas
>
> 2017-12-18 14:03 GMT+01:00 Matthew Knepley :
&g
On Mon, Dec 18, 2017 at 12:30 PM, Jose E. Roman wrote:
> I find the following definitions in petscconf.h, which are wrong because
> the corresponding subroutines are present.
>
> #define PETSC_MISSING_LAPACK_UNGQR 1
> #define PETSC_MISSING_LAPACK_HETRS 1
> #define PETSC_MISSING_LAPACK_HETRF 1
> #
On Mon, Dec 18, 2017 at 1:24 PM, Jose E. Roman wrote:
> > El 18 dic 2017, a las 18:58, Matthew Knepley
> escribió:
> >
> > On Mon, Dec 18, 2017 at 12:30 PM, Jose E. Roman
> wrote:
> > I find the following definitions in petscconf.h, which are wrong because
>
On Mon, Dec 18, 2017 at 1:29 PM, Jose E. Roman wrote:
> > El 18 dic 2017, a las 19:24, Matthew Knepley
> escribió:
> >
> > On Mon, Dec 18, 2017 at 1:24 PM, Jose E. Roman
> wrote:
> > > El 18 dic 2017, a las 18:58, Matthew Knepley
> escribió:
> > >
The configure is broken, perhaps beyond fixing. I need to give --with-mpi
to get anything
to work because the way it checks for MPI is screwed up. I will see if I
can throw away all its idiotic configure stuff in favor of just telling it
everything, but I am not eager to depend on something so rick
On Tue, Dec 19, 2017 at 1:54 PM, Matthew Knepley wrote:
> The configure is broken, perhaps beyond fixing. I need to give --with-mpi
> to get anything
> to work because the way it checks for MPI is screwed up. I will see if I
> can throw away all its idiotic configure stuff in f
501 - 600 of 3111 matches
Mail list logo