Re: [petsc-dev] Valgrind MPI-Related Errors

2020-06-01 Thread Karl Rupp
Hi Jacob, the recommendation in the past was to use MPICH as it is (was?) valgrind-clean. Which MPI do you use? OpenMPI used to have these kinds of issues. (My information might be outdated) Best regards, Karli On 6/2/20 2:43 AM, Jacob Faibussowitsch wrote: Hello All, TL;DR: valgrind

Re: [petsc-dev] MatCreateTranspose semantics

2020-06-01 Thread Jeff Hammond
Jeff dumb. Make copy paste error. Sorry. Jeff On Mon, Jun 1, 2020 at 5:02 PM Jeff Hammond wrote: > I'm still unable to get a basic matrix transpose working. I may be > stupid, but I cannot figure out why the object is in the wrong state, no > matter what I do. > > This is the full code: > >

[petsc-dev] Valgrind MPI-Related Errors

2020-06-01 Thread Jacob Faibussowitsch
Hello All,TL;DR: valgrind always complains about "Syscall param write(buf) points to uninitialised byte(s)” for a LOT of MPI operations in petsc code, making debugging using valgrind fairly annoying since I have to sort through a ton of unrelated stuff. I have built valgrind from source, used apt

Re: [petsc-dev] MatCreateTranspose semantics

2020-06-01 Thread Jeff Hammond
I'm still unable to get a basic matrix transpose working. I may be stupid, but I cannot figure out why the object is in the wrong state, no matter what I do. This is the full code: https://github.com/jeffhammond/PRK/commit/617973dfbe07d64cc2c0418a1702c418d51802c5

Re: [petsc-dev] MatCreateTranspose semantics

2020-06-01 Thread Jed Brown
I would just use MatTranspose here. We've had limited demand for more sophisticated operations on objects of type MATTRANSPOSE and there probably isn't much benefit in fusing the parallel version here anyway. Jeff Hammond writes: > I am trying to understand how to use a transposed matrix view

[petsc-dev] MatCreateTranspose semantics

2020-06-01 Thread Jeff Hammond
I am trying to understand how to use a transposed matrix view along the lines of Numpy ( https://github.com/ParRes/Kernels/blob/master/PYTHON/transpose-numpy.py#L99 ). https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateTranspose.html says "Creates a new matrix object that

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

2020-06-01 Thread Stefano Zampini
MKL fails on gemms calls too. At least in the past, not sure about the latest version. > On Jun 1, 2020, at 10:48 PM, Jeff Hammond wrote: > > I assume it's MKL sparse that fails for empty matrices, because the BLAS and > LAPACK should follow the Netlib convention. > > I'm sorry that MKL

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

2020-06-01 Thread Jeff Hammond
I assume it's MKL sparse that fails for empty matrices, because the BLAS and LAPACK should follow the Netlib convention. I'm sorry that MKL is unwilling to do nothing for you :-) Jeff On Sat, May 30, 2020 at 4:26 PM Mills, Richard Tran via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > Thanks for

Re: [petsc-dev] why does the Ubuntu PETSc package install coarray Fortran?

2020-06-01 Thread Jeff Hammond
Thanks everybody for your help. I'll use this in the future, but for now, building from source with minimal features is working for me. Jeff On Sun, May 31, 2020 at 9:06 PM Drew Parsons wrote: > On 2020-06-01 10:48, Satish Balay wrote: > > This is likely a dependency of a dependency. > > > >

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

2020-06-01 Thread Lisandro Dalcin
On Sun, 31 May 2020 at 02:26, Mills, Richard Tran via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > Thanks for the replies, everyone. It suppose is not actually that hard for > me to handle these dimensions properly -- I just hadn't personally > encountered or thought much about when such