KSP solver and Distributed arrays
Hi guys, A simple question. Can I solve a linear system [A]{x} = {b} using KSP solvers using Matrix and rhs, solution vectors with different number of ghost nodes (width)? As an example, can it be possible to solve for vector {x} DA, width=3, STAR stencil, [A] DA, width=1, STAR stencil {b} DA, width=1, STAR stencil I am suspecting that it is not possible. Thanks, Mohamad
The linear solvers only know about global dimensions and local sizes (they know nothing about local representations) so if the vectors and matrix are compatible the answer is yes. Barry On Dec 13, 2010, at 10:42 PM, Mohamad M. Nasr-Azadani wrote:
Hi guys,
A simple question. Can I solve a linear system [A]{x} = {b} using KSP solvers using Matrix and rhs, solution vectors with different number of ghost nodes (width)? As an example, can it be possible to solve for vector {x} DA, width=3, STAR stencil, [A] DA, width=1, STAR stencil {b} DA, width=1, STAR stencil
I am suspecting that it is not possible.
Thanks, Mohamad
Yes that is possible. The solver does not know about ghost nodes. Matt On Mon, Dec 13, 2010 at 8:42 PM, Mohamad M. Nasr-Azadani <[email protected]>wrote:
Hi guys,
A simple question. Can I solve a linear system [A]{x} = {b} using KSP solvers using Matrix and rhs, solution vectors with different number of ghost nodes (width)? As an example, can it be possible to solve for vector {x} DA, width=3, STAR stencil, [A] DA, width=1, STAR stencil {b} DA, width=1, STAR stencil
I am suspecting that it is not possible.
Thanks, Mohamad
-- 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
Thank you Barry and Matthew. Mohamad On Tue, Dec 14, 2010 at 7:14 AM, Matthew Knepley <[email protected]> wrote:
Yes that is possible. The solver does not know about ghost nodes.
Matt
On Mon, Dec 13, 2010 at 8:42 PM, Mohamad M. Nasr-Azadani <[email protected]
wrote:
Hi guys,
A simple question. Can I solve a linear system [A]{x} = {b} using KSP solvers using Matrix and rhs, solution vectors with different number of ghost nodes (width)? As an example, can it be possible to solve for vector {x} DA, width=3, STAR stencil, [A] DA, width=1, STAR stencil {b} DA, width=1, STAR stencil
I am suspecting that it is not possible.
Thanks, Mohamad
-- 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
participants (3)
-
Barry Smith -
Matthew Knepley -
Mohamad M. Nasr-Azadani