Re: [petsc-dev] petsc4py doc problem

2023-07-02 Thread Lisandro Dalcin
=260s, a very well respected person in the Python community. PETSc can very well increase the limit to anything you consider sensible, as long as the agreed limit is somehow enforced, which is what Stefano just implemented. -- Lisandro Dalcin Senior Research Scientist Extreme Computing Res

Re: [petsc-dev] PETSc future starting as a new design layer that runs on top of PETSc 3?

2022-07-31 Thread Lisandro Dalcin
On Sun, 31 Jul 2022 at 17:07, Matthew Knepley wrote: > On Sun, Jul 31, 2022 at 9:06 AM Lisandro Dalcin wrote: > >> On Sun, 31 Jul 2022 at 16:41, Jacob Faibussowitsch >> wrote: >> >>> >>> > Please don't take my words as advocacy for C++ &g

Re: [petsc-dev] PETSc future starting as a new design layer that runs on top of PETSc 3?

2022-07-31 Thread Lisandro Dalcin
dern C++ bindings implemented by wrapping the C interface. And then map these C++ bindings to Python. -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] PETSc future starting as a new design layer that runs on top of PETSc 3?

2022-07-31 Thread Lisandro Dalcin
C). I'm just pointing out that mapping C -> Python is not so easy to automate, unless you are OK with coding Python with a wrapper that is C-like and thus very unpythonic. -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah Universi

[petsc-dev] PetscUse/TryMethod

2022-04-03 Thread Lisandro Dalcin
The recent PetscUse/TryMethod changes are backward incompatible. Third-party codes cannot compile without modification. Our users deserve better. -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology

Re: [petsc-dev] PETSc-3.16.3 and Python-3.11

2022-01-16 Thread Lisandro Dalcin
at Bugzilla's bug: > https://bugzilla.redhat.com/show_bug.cgi?id=2039365#c3 > > -- > --- > Antonio Trande > Fedora Project > mailto: sagit...@fedoraproject.org > GPG key: 0xCC1CFEF30920C8AE > GPG key server: https://keyserver1.pgp.com/ > -- Lisandro Dalcin ==

Re: [petsc-dev] Problem with the reorganization of complex with clang

2021-05-16 Thread Lisandro Dalcin
://gitlab.com/petsc/petsc/-/jobs/1260342204 >> > >>> > > >> > >>> > > The compile fails building a C++ file, demangle.cxx. It fails >> at the >> > >>> first >> > >>> > > line, >> > >>> > > including , down in petscsystypes.h. >> It >> > >>> bombs >> > >>> > > during the definition of complex because it looks like the >> compiler >> > >>> is not >> > >>> > > defining __cplusplus, and thus takes the wrong branch. Is this >> what >> > >>> is >> > >>> > > happening? >> > >>> > > I cannot access this machine. >> > >>> > > >> > >>> > > 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/ >> > >> <http://www.cse.buffalo.edu/~knepley/> >> > >> >> > > >> > > >> > > -- >> > > 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/ >> > > <http://www.cse.buffalo.edu/~knepley/> >> > > >> > >> >> -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] Tag name v3.15

2021-03-31 Thread Lisandro Dalcin
get https://ftp.mcs.anl.gov/pub/petsc/release-snapshots/petsc-$x.$y.$z.tar.gz tar zxf petsc-$x.$y.$z.tar.gz cd petsc-$x.$y.$z ... But I cannot do git clone ... git checkout v$x.$y.$z because there is no tag v3.15.0, just v3.15 -- Lisandro Dalcin Senior Research Scientist Extreme

[petsc-dev] Tag name v3.15

2021-03-31 Thread Lisandro Dalcin
Satish, if we use semantic versioning for the PETSc version itself and the tarball, why don't we use the same format for tags? In short, why aren't we using v3.15.0 instead of v3.15 ? -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah

Re: [petsc-dev] integrate petsc4py tarball generation with petsc tarball generation

2021-03-15 Thread Lisandro Dalcin
solution is to get rid of Epydoc and replace it for something else, most likely Sphinx-based. If generating this epydoc stuff is too much of an annoyance for you, just kill it for good and let's forget about it. -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] integrate petsc4py tarball generation with petsc tarball generation

2021-03-15 Thread Lisandro Dalcin
ist". But this way the tarball will miss the API reference. PS: All this could be converted to a script that installs petsc4py and epydoc in a Python 2 virtual environment, and next users another venv for the Python 3 Sphinx stuff. -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] [petsc-users] headsup: switch git default branch from 'master' to 'main'

2021-02-23 Thread Lisandro Dalcin
witch on friday the 26th] > > We've previously switched 'maint' branch to 'release' before 3.14 > release - and this change (to 'main') is the next step in this direction. > > Satish > > -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Cent

Re: [petsc-dev] obscure changes in TSGetStages_Theta

2021-01-24 Thread Lisandro Dalcin
ges` array to the caller. Easy and clean, no API changes, no impact in third party code. -- Lisandro Dalcin Senior Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] Building petc4py inside PETSc source archive

2020-12-08 Thread Lisandro Dalcin
d like to get this going - but due to "git subtree" errors - we are > unable to make progress here [i.e keep this in sync with changes in petsc > repo] > > Satish, I think at this point we should kill the petsc/petsc4py repo. This will keep confusing people. -- Lisandro Dalcin =

Re: [petsc-dev] C++ error with just PETSc includes

2020-07-13 Thread Lisandro Dalcin
gt; > On Mon, 13 Jul 2020, Mark Adams wrote: > > > I have a wip C++ test that just calls PetscInitialize and Finalize, and I > > am getting some sort of error in the include file: > > > > https://gitlab.com/petsc/petsc/-/jobs/635768809 > > > > All I do is: &g

Re: [petsc-dev] Sphinx error

2020-07-09 Thread Lisandro Dalcin
should point to /usr/local/bin. PS: 99% of Python packages, and particularly the pure-Python ones, can be installed with pip under Homebrew Python. Go this way ASAP and forget about system Apple Python. -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] Sphinx error

2020-07-07 Thread Lisandro Dalcin
nd forget about the system Python install. >From brew you get `/usr/local/bin/python3`, and then you can freely `python3 -m pip install ` I further symlink `/usr/local/bin/{python|pip}3` to `~/bin/{python|pip}`, so I don't have to deal ever again with the EOLed system Apple Python 2. -

Re: [petsc-dev] Meaning of PETSc matrices with zero rows but nonzero columns?

2020-06-01 Thread Lisandro Dalcin
t; > Just FYI, this is how NumPy handles it: In [1]: import numpy as np In [2]: A = np.empty([3,0]) In [3]: B = np.empty([0,8]) In [4]: C = np.dot(A,B) In [5]: C.shape Out[5]: (3, 8) If we do not follow these rules, it will bite us in the future. -- Lisandro Dalcin Research

Re: [petsc-dev] Configure won't find python

2020-04-17 Thread Lisandro Dalcin
r code in case things go awry > > Echo “PYTHON NOT FOUND” >2 > > "Echo" ? Are you really running this on Ubuntu? Or you just hacked in on macOS and assumed that it will work everywhere? BTW, Are you aware of the behavior of `command -v` with shell aliases? $ alias p

Re: [petsc-dev] snes_tutorials-ex19_cuda_1

2020-04-01 Thread Lisandro Dalcin
Well, my request will not fix the problem: https://gitlab.com/petsc/petsc/-/jobs/495147366#L5231 On Thu, 2 Apr 2020 at 03:26, Lisandro Dalcin wrote: > Can anyone messing with CPUs please update test snes_tutorials-ex19_cuda_1 > to use -ksp_monitor_short and update its output with REP

[petsc-dev] snes_tutorials-ex19_cuda_1

2020-04-01 Thread Lisandro Dalcin
Can anyone messing with CPUs please update test snes_tutorials-ex19_cuda_1 to use -ksp_monitor_short and update its output with REPLACE=1 ? Please do it in maint, or cherry-pick if already fixed in master. Regards, -- Lisandro Dalcin Research Scientist Extreme Computing Research

[petsc-dev] New std CPP in MSVC

2020-03-28 Thread Lisandro Dalcin
https://devblogs.microsoft.com/cppblog/announcing-full-support-for-a-c-c-conformant-preprocessor-in-msvc/ -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] Globsearch fails for me when running tests

2020-03-27 Thread Lisandro Dalcin
print(' '.join(list(itertools.chain.from_iterable(m" '$(globsearch)' > '$(alltesttargets)') > > For this feature, probably put them in an argsfile > > $(file >$(TESTDIR)/globsearch.args,$(alltesttargets)) > > IIRC, the $(file ...) function does not wo

Re: [petsc-dev] Suggestions for MatProductCreate()

2020-03-27 Thread Lisandro Dalcin
outine, right? In that >> case, you are right, no changes/refactors needed here. >> > Yes, MatProductSetFromOptions() only > dispatches SetFromOptions_typeA_typeB_producttype() which handles local > options and polymorphism. The organization of new API is described at the > top of petsc/src/mat/interface/matproduct.c. > > OK, all good with that one. Regards, -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] Suggestions for MatProductCreate()

2020-03-24 Thread Lisandro Dalcin
ns() is logically equivalent to MatSetFromOptions(), > MatProductSymbolic() is logically equivalent to MatSetUp(), > and > MatProductNumeric() assembles the matproduct D repeatedly when A or B are > updated. > > OK, so... your point it that what looks like "setup"

Re: [petsc-dev] Suggestions for MatProductCreate()

2020-03-23 Thread Lisandro Dalcin
e proper way to do it, IMHO. This way is the one that gives maximum flexibility to power users/developers, because some times developers may want to configure objects with some parameters, and do a SetUp() without calling SetFromOptions() to prevent users to mess up with the parameters the d

Re: [petsc-dev] Suggestions for MatProductCreate()

2020-03-23 Thread Lisandro Dalcin
r of the previous API. > I only see one or two PETSc subroutines call them. I do not think we > need provide backward compatibility PETSC_DEPRECATED calls for 6 pairs of > such routines. > > OK, you have a good point here, the effort is not worth it. Forget about this one. -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] Suggestions for MatProductCreate()

2020-03-22 Thread Lisandro Dalcin
MatProductSetFromOptions() is not really setting any options. It looks like the correct name for this API is MatProductSetUp(). In any case, it looks like both SetFromOptions() and SetUp() should be provided. On Mon, 23 Mar 2020 at 01:48, Lisandro Dalcin wrote: > * Please consider fix

[petsc-dev] Suggestions for MatProductCreate()

2020-03-22 Thread Lisandro Dalcin
nzero pattern of A^2. * It should be also considered to provide backward compatibility PETSC_DEPRECATED calls to the previous MatMatMultSymbolic() and MatMatMultNumeric(). It looks like it would be trivial to do, though I may be getting it wrong because I have not looked at all the details. Regards,

[petsc-dev] lazygit

2020-03-20 Thread Lisandro Dalcin
a statically-linked binary (a bit heavy, ~ 15MB), just drop it in your ~/bin or ~/.local/bin, and you are ready to go. Profit! -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] Question about Binary-IO in READ mode with POSIX APIs

2020-03-16 Thread Lisandro Dalcin
ncerns about it. > But of course, when not using MPI-IO, only process zero should open the > file. > > That's what I was asking about. So it seems we all are agree about this. I'll go ahead and submit a MR. -- Lisandro Dalcin Research Scientist Extreme Computing Research Cen

Re: [petsc-dev] Question about Binary-IO in READ mode with POSIX APIs

2020-03-16 Thread Lisandro Dalcin
On Mon, 16 Mar 2020 at 16:35, Jed Brown wrote: > Lisandro Dalcin writes: > > > Currently, binary viewers using POSIX file descriptors with READ mode > open > > the file in ALL processes in the communicator. For WRITE mode, only > process > > zero opens

[petsc-dev] Question about Binary-IO in READ mode with POSIX APIs

2020-03-16 Thread Lisandro Dalcin
the enhancements to binary viewers, we should consider whether we should turn MPI-IO on by default. -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

[petsc-dev] Add -malloc_dump to Stage-1 ?

2020-03-11 Thread Lisandro Dalcin
Dear Satish, I pushed a test that did not deallocate a matrix. I got a failure in Stage-2 out of -malloc_dump https://gitlab.com/petsc/petsc/pipelines/125455265 Would't it make sense to have -malloc_dump on any of the Stage-1 builds? -- Lisandro Dalcin Research Scientist Extreme

Re: [petsc-dev] Request for comments: allow C99 internally

2020-03-11 Thread Lisandro Dalcin
(*a)[n] = (double (*)[n])p; > > a[i][j] = 42; > Oh, boy... PS: BTW, if you use typeof() [*], you can simplify these VLA declarations: double (*a)[m][n] = (typeof(a)) p; a[i][j][k] = 42; [*] Don't remember if it is part of C99 or up -- Lisandro Dalcin Research Scientist Ext

Re: [petsc-dev] Request for comments: allow C99 internally

2020-03-09 Thread Lisandro Dalcin
o the public headers, so a user of PETSc could still build > with -std=c89 and the like. > > Definitely +1 for variadic macros and for-loop declarations, but not VLAs. -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University

Re: [petsc-dev] Signal handling

2020-03-09 Thread Lisandro Dalcin
the current API is a bit weird... PS: Do you really think this may be actually useful in production to be worth to push into PETSc? -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

[petsc-dev] Signal handling

2020-03-05 Thread Lisandro Dalcin
may be useful, but I'm not sure. Do you think this may be of some value for core PETSc? I'm asking before submitting a MR because that would require writing some docs, I don't want to do the doc work before knowing your opinion first :-). Regards, -- Lisandro Dalcin Research Scient

Re: [petsc-dev] [petsc-users] Matrix-free method in PETSc

2020-02-20 Thread Lisandro Dalcin
> Extending to 2 and 3d is straight forward. > >> > > >> > Any questions let us know. > >> > > >> >Barry > >> > > >> >If you like this would make a great merge request with your > code to improve our examples. > >> > > >> > > >> > > On Feb 15, 2020, at 9:42 PM, Yuyun Yang > wrote: > >> > > > >> > > Hello team, > >> > > > >> > > I wanted to apply the Krylov subspace method to a matrix-free > implementation of a stencil, such that the iterative method acts on the > operation without ever constructing the matrix explicitly (for example, > when doing backward Euler). > >> > > > >> > > I'm not sure whether there is already an example for that > somewhere. If so, could you point me to a relevant example? > >> > > > >> > > Thank you! > >> > > > >> > > Best regards, > >> > > Yuyun > >> > >> > >> > >> > >> > >> -- > >> 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/ > > > > -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] complex fix and -Wfloat-equal

2020-01-29 Thread Lisandro Dalcin
On Thu, 30 Jan 2020 at 08:25, Smith, Barry F. wrote: > > > > On Jan 27, 2020, at 3:42 PM, Lisandro Dalcin wrote: > > > > +#if defined(__GNUC__) && (__GNUC__ >= 4) > > +#pragma GCC system_header > > +#endif > > + > > #define PETSC_CXX_C

[petsc-dev] complex fix and -Wfloat-equal

2020-01-27 Thread Lisandro Dalcin
&& const_cast(lhs).real() == PetscReal(rhs); } \ | ^ /home/devel/petsc/dev/include/petsccxxcomplexfix.h:60:1: note: in expansion of macro ‘PETSC_CXX_COMPLEX_FIX’ 60 | PETSC_CXX_COMPLEX_FIX(PetscInt) -- Lisa

[petsc-dev] PETSC_HAVE_PRAGMA_DISJOINT

2020-01-27 Thread Lisandro Dalcin
) src/mat/impls/sell/seq/sell.c:#if defined(PETSC_HAVE_PRAGMA_DISJOINT) src/mat/impls/sell/seq/sell.c:#if defined(PETSC_HAVE_PRAGMA_DISJOINT) src/mat/impls/sell/seq/sell.c:#if defined(PETSC_HAVE_PRAGMA_DISJOINT) src/vec/vec/impls/seq/dvec2.c:#if defined(PETSC_HAVE_PRAGMA_DISJOINT) -- Lisandro

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

2020-01-20 Thread Lisandro Dalcin
his > behavior for me. > > > > Satish > > > > - > > > > wget > https://kojipkgs.fedoraproject.org//packages/ccache/3.4.3/2.fc30/x86_64/ccache-3.4.3-2.fc30.x86_64.rpm > > rpm -Uvh --oldpackage ccache-3.4.3-2.fc30.x86_64.rpm > > > > -

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

2020-01-20 Thread Lisandro Dalcin
achine. Would be interesting to check a vanilla build or gcc-10 > branch. I'm still learning toward it being unintentional. > > On Jan 20, 2020 07:57, Lisandro Dalcin wrote: > > > > On Mon, 20 Jan 2020 at 17:07, Jed Brown wrote: > > From my man page (contradicti

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

2020-01-20 Thread Lisandro Dalcin
MT, but it is cumbersome to do, unless you use "-MT $@", though no idea how to update configure for that. -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

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

2020-01-20 Thread Lisandro Dalcin
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 >> dependency generation, and things are now broken for PETSc. I do not have a >> c

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

2020-01-20 Thread Lisandro Dalcin
multiple '-MT' options. For example, '-MT '$(objpfx)foo.o'' might give $(objpfx)foo.o: foo.c -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

[petsc-dev] make: print path to mpi.h

2019-12-16 Thread Lisandro Dalcin
flags somewhere. -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] AVX kernels, old gcc, still broken

2019-10-25 Thread Lisandro Dalcin via petsc-dev
ifort', 'PETSC_ARCH=arch-gnu-opt', ] configure.petsc_configure(configure_options) -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/

Re: [petsc-dev] petsc-dev post from tabrez...@gmail.com requires approval

2019-08-28 Thread Lisandro Dalcin via petsc-dev
On Wed, 28 Aug 2019 at 06:45, Smith, Barry F. wrote: > > > > On Aug 27, 2019, at 3:27 AM, Matthew Knepley wrote: > > > > On Tue, Aug 27, 2019 at 4:11 AM Lisandro Dalcin > wrote: > > On Mon, 26 Aug 2019 at 21:26, Smith, Barry F. > wrote: > > > &

[petsc-dev] Scotch 6.0.8

2019-08-27 Thread Lisandro Dalcin via petsc-dev
There is a new scotch release. We should update configure to this version. I'm a bit confused about the current status PTscotch.py in BuildSystem, it seems to be using the scotch/master from the official repo. -- Lisandro Dalcin Research Scientist Extreme Computing Research Center

Re: [petsc-dev] petsc-dev post from tabrez...@gmail.com requires approval

2019-08-27 Thread Lisandro Dalcin via petsc-dev
On Mon, 26 Aug 2019 at 21:26, Smith, Barry F. wrote: > > > > On Aug 26, 2019, at 10:11 AM, Lisandro Dalcin wrote: > > > > > > > > On Sun, 25 Aug 2019 at 18:37, Smith, Barry F. via petsc-dev < > petsc-dev@mcs.anl.gov> wrote: > > > >

[petsc-dev] Bad use of defined(MPI_XXX)

2019-05-24 Thread Lisandro Dalcin via petsc-dev
These two are definitely wrong, we need PETSC_HAVE_MPI_XXX instead. include/petscsf.h:#if defined(MPI_REPLACE) src/sys/objects/init.c:#if defined(PETSC_USE_64BIT_INDICES) || !defined(MPI_2INT) -- Lisandro Dalcin Research Scientist Extreme Computing Research Center (ECRC) King

[petsc-dev] Bogus tag

2019-05-03 Thread Lisandro Dalcin via petsc-dev
$ git branch * master $ git describe WIP-112-g872c329e82 $ git show WIP tag WIP Tagger: Jed Brown Date: Wed Apr 24 22:56:11 2019 + Work In Progress (do not merge) ... Is there any way we can get rid of that? -- Lisandro Dalcin Research Scientist Extreme Computing Research

[petsc-dev] DMLabel struct

2019-03-27 Thread Lisandro Dalcin via petsc-dev
DMLabel is now a PetscObject, but the struct name does not match the usual conventions, it should be _p_DMLabel, but instead: include/petsc/private/dmlabelimpl.h:struct _n_DMLabel { include/petscdmlabel.h:typedef struct _n_DMLabel *DMLabel; -- Lisandro Dalcin Research Scientist

Re: [petsc-dev] MATLMVM

2018-09-11 Thread Lisandro Dalcin
the Dark Lord would do it??"... and the answer should always revolve around reducing top level API complexity. Then you do it that way, and it's going to be the right way. -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CE

[petsc-dev] MATLMVM

2018-09-11 Thread Lisandro Dalcin
brdn" #define MATLMVMBADBRDN "lmvmbadbrdn" #define MATLMVMSYMBRDN "lmvmsymbrdn" #define MATLMVMSYMBADBRDN "lmvmsymbadbrdn" #define MATLMVMDIAGBRDN"lmvmdiagbrdn" -- Lisandro Dalcin Research Scientist Computer, Electrical and

[petsc-dev] TSBASICSYMPLECTIC

2018-09-11 Thread Lisandro Dalcin
. Or, if there are no plans to have a family, then why to name it BASIC in the first place? PS: Not an expert in the field, feel free to hammer me about my ignorance. -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Ext

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

2018-09-04 Thread Lisandro Dalcin
PlexTopologyInterpolate?) - > would still refer to "interpolation" if somebody is used to that and the > same time more clearly distinguish from DMInterpolate. And as Matt says, > it's usually used up to once per program so the longer name is not a pain. > > Vaclav > -- L

Re: [petsc-dev] Adding p4est to nightly tests

2018-07-01 Thread Lisandro Dalcin
un, 1 Jul 2018 at 01:07, Stefano Zampini wrote: > > Satish, > > when you have time can you add p4est to some of the archs we tests nightly? > > Matt and/or Toby, is there any reason why p4est is not under the testsuite? > > Thanks, > Stefano -- Lisandro Dalcin =

Re: [petsc-dev] PetscObjectViewFromOptions() pull request 1005 1006

2018-06-29 Thread Lisandro Dalcin
he > configure-time flag turning off the viewers dropped that down to ~1%. I don't > know what PCPATCH is, but it sounds like it may be doing something similar. > That sounds like a more realistic number. -- Lisandro Dalcin Research Scientist Computer, Electrical a

Re: [petsc-dev] PetscObjectViewFromOptions() pull request 1005 1006

2018-06-29 Thread Lisandro Dalcin
bject (and trivial matrix, and ILU preconditioner) is used over and over again (all memory will reside in cache?), and all the solves are just one iteration. The fact that you get an improvement in Matt's branch respect to turning the viewers off smells to just C function call over

Re: [petsc-dev] PetscObjectViewFromOptions() pull request 1005 1006

2018-06-28 Thread Lisandro Dalcin
I didn't do > further exhaustive benchmarking to check for further problems. > There you have! Thanks for confirming my previous guess. I would really like that PR being split, and leave the changes to XXXSetFromOptions() on hold in a separate PR. -- Lisandro Dalcin Research

Re: [petsc-dev] PetscObjectViewFromOptions() pull request 1005 1006

2018-06-28 Thread Lisandro Dalcin
nges the core KSP infrastructure, and I know it is going to be merged as it is, and I know in advance I'm already defeated. Let it be, Let it be... -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing R

[petsc-dev] Opening XML performance log files in browsers

2018-06-13 Thread Lisandro Dalcin
nd ADDRESS Specify bind address -p PORT, --port PORT Specify port $ ./petsc-performance-view log.xml # should open browser with log.xml # and the script should exit immediately # to not seize your command line -- Lisandro Dalcin Research Scientist Computer,

Re: [petsc-dev] Nested XML logging and stages

2018-05-17 Thread Lisandro Dalcin
https://bitbucket.org/petsc/petsc/pull-requests/964/logxml-various-enhancements/diff On Mon, 14 May 2018 at 23:51, Smith, Barry F. <bsm...@mcs.anl.gov> wrote: >I am fine with the suggested changes. > > On May 13, 2018, at 1:06 PM, Lisandro Dalcin <dalc...@gmail.com> w

[petsc-dev] Nested XML logging and stages

2018-05-13 Thread Lisandro Dalcin
generated. So I guess we need something like PetscMax(threshold,1e-10). Why PetscLogView_Nested() calls PetscLogNestedEnd() ? While convenient and quick to implement, it smells like bad API usage. -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematica

Re: [petsc-dev] Proposed changes to TS API

2018-05-11 Thread Lisandro Dalcin
ts (which is definitely not efficient and is what Hong does also). Yes, that's basically the proposal. User could even provide both IJacobian and SplitJacobians, and TS can then make the best use of these callbacks depending on what is more efficient. -- Lisandro Dalcin Research Scientist

Re: [petsc-dev] GCC8 Fortran length changes from int to size_t

2018-05-07 Thread Lisandro Dalcin
On Mon, 7 May 2018 at 03:11, Jed Brown <j...@jedbrown.org> wrote: > Lisandro Dalcin <dalc...@gmail.com> writes: > >> > How do you exactly want to implement that? Totally replace these > >> > special functions with the BIND(C) interface that calls di

Re: [petsc-dev] GCC8 Fortran length changes from int to size_t

2018-05-06 Thread Lisandro Dalcin
On Thu, 3 May 2018 at 20:21, Jed Brown <j...@jedbrown.org> wrote: > Lisandro Dalcin <dalc...@gmail.com> writes: > > On 3 May 2018 at 18:50, Smith, Barry F. <bsm...@mcs.anl.gov> wrote: > >> > >>Jeff, (and others), > >> >

Re: [petsc-dev] GCC8 Fortran length changes from int to size_t

2018-05-03 Thread Lisandro Dalcin
t calls the C function through a BIND(C) interface? -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc

Re: [petsc-dev] GCC8 Fortran length changes from int to size_t

2018-05-02 Thread Lisandro Dalcin
On Wed, 2 May 2018 at 17:29, Satish Balay <ba...@mcs.anl.gov> wrote: > So we need a 64bit arm with gcc8 - for this testcase failure? Or a big-endian machine/OS ? -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (C

[petsc-dev] Add PetscFree1() ?

2018-04-04 Thread Lisandro Dalcin
ters) I guess we should just rewrite PetscFree() to be safer, but too close to release right now to play these games... -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King

Re: [petsc-dev] futurize

2018-04-02 Thread Lisandro Dalcin
port > configuration changes from 'maint' to 'master' after the release, under > the assumption that we'll want to support Python-3 before a 3.10 > release. Well, while I'm a bit worried about dependencies, I think your goal makes total sense and I would support it. -- Lisandro Dalcin ==

[petsc-dev] Fix VecScatterDestroy() and issue with VecScatterCreateEmpty()

2018-04-02 Thread Lisandro Dalcin
. -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/ 4700 King Abdullah University of Sci

Re: [petsc-dev] No MatGetOperation() ?

2018-04-01 Thread Lisandro Dalcin
Changed my mind a little about Mat{Set|Get}Operation() and fixed some mess in MATSHELL: https://bitbucket.org/petsc/petsc/pull-requests/915/add-matgetoperation-and-fixes-to-matshell/diff On 1 April 2018 at 12:30, Lisandro Dalcin <dalc...@gmail.com> wrote: > Moreover, after a second look

Re: [petsc-dev] No MatGetOperation() ?

2018-04-01 Thread Lisandro Dalcin
Moreover, after a second look, I think MatSetOperation() is wrong. It should special-case MATSHELL and dispatch to MatShellSetOperation() if the type match, otherwise bad things would happen. On 1 April 2018 at 12:08, Lisandro Dalcin <dalc...@gmail.com> wrote: > Barry, you made the di

[petsc-dev] No MatGetOperation() ?

2018-04-01 Thread Lisandro Dalcin
Barry, you made the distinction MatSetOperation() vs. MatShellSetOperation(). 1) MatSetOperation() is misplaced, it should go somewhere in src/mat/interface/matrix.c 2) We should add MatGetOperation(), right? Should it just return the routine set in mat-ops? Are you ok with this? -- Lisandro

[petsc-dev] TS: -info_exclude and -log_exclude

2018-03-29 Thread Lisandro Dalcin
add `-xxx_exclude <dmts,tsadapt,tstrajectory>` for fine-grained deactivation. Option 3: We implement both option 1 and option 2. So, Barry, what's your preference? The quickest one is just (1) but the others are also rather trivial to implement. -- Lisandro Dalcin Re

Re: [petsc-dev] MPI I/O tests failing in master after merging my branch

2018-03-22 Thread Lisandro Dalcin
those that are missing (as we do for MPIU_INT). Is > there more > to it than that? > We have that (I mean, each PETSc typedef). What we are somehow missing is the mapping plain C 'int' -> MPI_INT. That's important for binary I/O with non-PETSc formats (e.g. GMSH). -- Lisandro Da

Re: [petsc-dev] MPI I/O tests failing in master after merging my branch

2018-03-22 Thread Lisandro Dalcin
On 22 March 2018 at 12:51, Matthew Knepley <knep...@gmail.com> wrote: > On Thu, Mar 22, 2018 at 3:44 AM, Lisandro Dalcin <dalc...@gmail.com> wrote: >> >> >> ftp://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2018/03/21/examples_master_arch-linux-cxx-cmplx-pkg

Re: [petsc-dev] MPI I/O tests failing in master after merging my branch

2018-03-22 Thread Lisandro Dalcin
On 22 March 2018 at 11:10, Lisandro Dalcin <dalc...@gmail.com> wrote: > On 22 March 2018 at 10:44, Lisandro Dalcin <dalc...@gmail.com> wrote: >> >> @Barry or anyone: the only issue for me is to decide on the name, I'll >> call it PETSC_C_INT (and maybe al

Re: [petsc-dev] MPI I/O tests failing in master after merging my branch

2018-03-22 Thread Lisandro Dalcin
On 22 March 2018 at 10:44, Lisandro Dalcin <dalc...@gmail.com> wrote: > > @Barry or anyone: the only issue for me is to decide on the name, I'll > call it PETSC_C_INT (and maybe also add PETSC_C_LONG along the way). > If you don't like and have a better naming convention to sug

[petsc-dev] MPI I/O tests failing in master after merging my branch

2018-03-22 Thread Lisandro Dalcin
: the only issue for me is to decide on the name, I'll call it PETSC_C_INT (and maybe also add PETSC_C_LONG along the way). If you don't like and have a better naming convention to suggest, please let me know. -- Lisandro Dalcin Research Scientist Computer, Electrical

Re: [petsc-dev] [Mike McQuaid] Re: [Homebrew/homebrew-core] PETSc: import from homebrew-science (#23598)

2018-02-25 Thread Lisandro Dalcin
scripts. Why is "fixing" the installer so difficult that you guys > just leave it in a broken state for years? (The reason I don't "fix" it > is that I don't know what needs to be fixed since I don't interact with the > package managers. > Well, maybe

Re: [petsc-dev] [Mike McQuaid] Re: [Homebrew/homebrew-core] PETSc: import from homebrew-science (#23598)

2018-02-24 Thread Lisandro Dalcin
at to go unchallenged. >> >> — >> You are receiving this because you were mentioned. >> Reply to this email directly, view it on GitHub >> <https://github.com/Homebrew/homebrew-core/pull/23598#issuecomment-368122755>, >> or mute the thread >> <https://gi

Re: [petsc-dev] SETERRQ in fortran

2018-01-31 Thread Lisandro Dalcin
>>> will make it easier to understand and avoid the circus antics when used >>>> in any context other than a positive conditional with no else clause. >>> >>> I'll take this under advisement. Of course in our examples the endif will >>> ALWAYS be on t

Re: [petsc-dev] Possible bug in PCFIELDSPLIT

2018-01-10 Thread Lisandro Dalcin
your diagnostic is correct, the code is just missing the inner "atype ==" part you wrote. -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King Abdullah Univer

[petsc-dev] Possible bug in PCFIELDSPLIT

2018-01-10 Thread Lisandro Dalcin
A00's %sdiagonal's inverse\n",atype == MAT_SCHUR_COMPLEMENT_AINV_DIAG ? "" : (MAT_SCHUR_COMPLEMENT_AINV_BLOCK_DIAG ? "block " : "lumped "));CHKERRQ(ierr);break; -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & En

Re: [petsc-dev] some change in DM plex output breaks next?

2017-11-25 Thread Lisandro Dalcin
t; Vertex Sets: 1 strata with value/size (1 (1)) > # 10a10 > # > Vertex Sets: 1 strata with value/size (1 (1)) > -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC)

Re: [petsc-dev] proposed minor PetscPartitioner changes

2017-11-08 Thread Lisandro Dalcin
tor... (just kidding, actually, Barry is way too much benevolent with all us damn beasts). -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King Abdullah University of S

Re: [petsc-dev] proposed minor PetscPartitioner changes

2017-11-08 Thread Lisandro Dalcin
On 8 November 2017 at 05:51, Smith, Barry F. <bsm...@mcs.anl.gov> wrote: > >> On Nov 7, 2017, at 1:33 AM, Lisandro Dalcin <dalc...@gmail.com> wrote: >> >> The only concern I have about PetscPartitioner is that the API depends >> on DM (Pets

Re: [petsc-dev] proposed minor PetscPartitioner changes

2017-11-06 Thread Lisandro Dalcin
ere is some sort of deficiency there. > However, I do not believe that there is. > The only concern I have about PetscPartitioner is that the API depends on DM (PetscPartitionerPartition_ routines). Maybe PetscPartitioner should eventually move to became more agnostic, and that way it can be used to partition mat

Re: [petsc-dev] is this wrong?

2017-11-06 Thread Lisandro Dalcin
onst char *v = PetscBools[currentvalue]; > ierr = (*PetscHelpPrintf)(PetscOptionsObject->comm," -%s%s: <%s> %s > (%s)\n",PetscOptionsObject->prefix?PetscOptionsObject->prefix:"",opt+1,v,text,ManSection(man));CHKERRQ(ierr); > } > PetscFunctionRet

Re: [petsc-dev] Improving nightly builds for maint

2017-10-21 Thread Lisandro Dalcin
ke` does fails with non-zero exit status), are a bit useless. Let's keep Travis-CI as a quick way to get fast feedback on commits that would make the build fail, a no more. Anyway, we all depend on merging into next and checking the nightly builds. -- Lisandro Dalcin Resear

Re: [petsc-dev] Improving nightly builds for maint

2017-10-21 Thread Lisandro Dalcin
On 20 October 2017 at 21:12, Jed Brown <j...@jedbrown.org> wrote: > Barry Smith <bsm...@mcs.anl.gov> writes: > >>> On Oct 18, 2017, at 3:10 PM, Lisandro Dalcin <dalc...@gmail.com> wrote: >>> >>> On 18 October 2017 at 22:36, Satish Balay <ba...@

Re: [petsc-dev] Improving nightly builds for maint

2017-10-18 Thread Lisandro Dalcin
y build matrix, right? -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/ 4700 King Abdul

Re: [petsc-dev] Improving nightly builds for maint

2017-10-18 Thread Lisandro Dalcin
On 18 October 2017 at 22:36, Satish Balay <ba...@mcs.anl.gov> wrote: > On Wed, 18 Oct 2017, Lisandro Dalcin wrote: > > Alternate is to remove/increase the timout. It doesn't work for some > builds anyway [esp valgrind runs] > > This feature relies on sending 'kill

[petsc-dev] Improving nightly builds for maint

2017-10-18 Thread Lisandro Dalcin
valgrind? -- Lisandro Dalcin Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/ 4700 King Abdullah Univer

Re: [petsc-dev] (no subject)

2017-10-09 Thread Lisandro Dalcin
Please ignore the previous email. The INSTALL instructions in the tarball explicitly state that both libraries should be linked. On 9 October 2017 at 15:29, Lisandro Dalcin <dalc...@gmail.com> wrote: > Satish, in PTScotch.py, you check for the parallel (libpt*.a) as well > as sequenti

  1   2   3   4   5   6   7   8   >