Quick update for the list: Matt and I were emailing back and forth a bit
and I at least have a workaround for now.
It turns out that MVAPICH2 does include their own implementations of
malloc/calloc . Matt believes (and I agree) that they should be private to
the library though. It looks like so
On Mon, Dec 5, 2016 at 10:56 PM, Derek Gaston wrote:
> Please excuse the slightly off-topic post: but I'm pulling my hair out
> here and I'm hoping someone else has seen this before.
>
> I'm calling PETSc from Julia and it's working great on my Mac with MPICH
> but I'm seeing a segfault on my Lin
Please excuse the slightly off-topic post: but I'm pulling my hair out here
and I'm hoping someone else has seen this before.
I'm calling PETSc from Julia and it's working great on my Mac with MPICH
but I'm seeing a segfault on my Linux cluster using MVAPICH2. I get the
same segfault both with th
> On Dec 5, 2016, at 2:47 PM, K. N. Ramachandran wrote:
>
> Hello PETSc-Users,
>
> I am working on an application where we capture the A matrix in a linear
> system Ax=b, which is solved using Petsc. Let us also say that the matrix A
> can change after a few iterations. We want to capture onl
"K. N. Ramachandran" writes:
> Hello PETSc-Users,
>
> I am working on an application where we capture the A matrix in a linear
> system Ax=b, which is solved using Petsc. Let us also say that the matrix A
> can change after a few iterations. We want to capture only the changed
> matrices and simp
Hello PETSc-Users,
I am working on an application where we capture the A matrix in a linear
system Ax=b, which is solved using Petsc. Let us also say that the matrix A
can change after a few iterations. We want to capture only the changed
matrices and simply avoid the duplicate ones.
I was consid
On 5 December 2016 at 16:49, Massoud Rezavand wrote:
> Dear Petsc team,
>
> In order to create a parallel matrix and solve by KSP, is it possible to
> directly use MatSetValues() in runtime when each matrix entry is just
> created without MatMPIAIJSetPreallocation()?
>
Yes, but performance will
Dear Petsc team,
In order to create a parallel matrix and solve by KSP, is it possible to
directly use MatSetValues() in runtime when each matrix entry is just
created without MatMPIAIJSetPreallocation()?
I mean, when you only know the global size of Mat, and the number of
nonzeros per row is no
Hi Matt:
Thank you for your kind reply. I am aware of this problem from my test
case. I simulate the lid driven cavity
by the code, and the grid is a 100x100 2D domain. I use the routine
DMPlexReconstructGradientsFVM to compute
the gradients and limiters. The limiter which I used in the code
