I have a matrix system Ax = b, A is of type MatSeqAIJ or MatMPIAIJ, depending
on the number of cores.
I try to solve this problem by pc_type ilu and ksp_type bcgs, it does not
converge. The options I specify are:
-pc_type hypre -pc_hypre_type pilut -pc_hypre_pilut_maxiter 1000
I’m gathering from your suggestions that I would need, a priori, knowledge of
how many ghost points I would need, is that right?
-gideon
On Feb 17, 2015, at 9:10 AM, Matthew Knepley knep...@gmail.com wrote:
On Tue, Feb 17, 2015 at 7:46 AM, Gideon Simpson gideon.simp...@gmail.com
On Tue, Feb 17, 2015 at 8:15 AM, Gideon Simpson gideon.simp...@gmail.com
wrote:
I’m gathering from your suggestions that I would need, a priori, knowledge
of how many ghost points I would need, is that right?
We have to be more precise about a priori. You can certainly create a
VecScatter on
This warrants a deeper description. I was thinking about trying to use petsc
to solve a 1D PDE problem using the equidistribution principle for adaptive
meshing. Briefly, the idea is that given a function u(x), where x is the
physical coordinate, a monitor function, w[u(x)]0, is introduced
On Tue, Feb 17, 2015 at 7:46 AM, Gideon Simpson gideon.simp...@gmail.com
wrote:
Suppose I have data in Vec x and Vec y, and I want to interpolate this
onto Vec xx, storing the values in Vec yy. All vectors have the same
layout. The problem is that, for example, some of the values in xx on
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
On 10/02/15 01:31, Jed Brown wrote:
Lawrence Mitchell lawrence.mitch...@imperial.ac.uk writes:
Having just tried a build with --download-mpich, I notice this
problem does not occur. So should I shout at the OpenMPI team?
Open MPI has many
Suppose I have data in Vec x and Vec y, and I want to interpolate this onto Vec
xx, storing the values in Vec yy. All vectors have the same layout. The
problem is that, for example, some of the values in xx on processor 0 may need
the values of x and y on processor 1, and so on. Aside from