Re: [petsc-users] Scatter error

2013-10-01 Thread Jed Brown
Mohammad Mirzadeh mirza...@gmail.com writes:
 I assume I'm getting this since they somehow share a VecScatter
 internally because of VecDuplicate? If so how can I fix this?

Hmm, there is no VecScatterDuplicate so I think you'd have to make the
vectors independently.

 Also as a side question, if I stored fxx and fyy inside a block
 vector, would I gain any performance in terms of MPI traffic? What
 about in 3D?

It would be half the messages.  Why are these vectors separate in the
first place?


pgp3uxLLJzYPu.pgp
Description: PGP signature


Re: [petsc-users] Scatter error

2013-10-01 Thread Mohammad Mirzadeh
 It would be half the messages.  Why are these vectors separate in the
 first place?

I had not worked with block vectors before so just convenience i
guess! When working with blocked vectors, how are items stored? Is it
like fxx[0], fyy[0], fxx[1], fyy[1], ... or is it like
fxx[0],fxx[1],...,fyy[0],fyy[1],...


Re: [petsc-users] Scatter error

2013-10-01 Thread Matthew Knepley
On Tue, Oct 1, 2013 at 7:02 PM, Mohammad Mirzadeh mirza...@gmail.comwrote:

  It would be half the messages.  Why are these vectors separate in the
  first place?

 I had not worked with block vectors before so just convenience i
 guess! When working with blocked vectors, how are items stored? Is it
 like fxx[0], fyy[0], fxx[1], fyy[1], ... or is it like


This one.

   Matt


 fxx[0],fxx[1],...,fyy[0],fyy[1],...




-- 
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


Re: [petsc-users] Scatter error

2013-10-01 Thread Mohammad Mirzadeh
Thanks Matt.

On Tue, Oct 1, 2013 at 5:26 PM, Matthew Knepley knep...@gmail.com wrote:
 On Tue, Oct 1, 2013 at 7:02 PM, Mohammad Mirzadeh mirza...@gmail.com
 wrote:

  It would be half the messages.  Why are these vectors separate in the
  first place?

 I had not worked with block vectors before so just convenience i
 guess! When working with blocked vectors, how are items stored? Is it
 like fxx[0], fyy[0], fxx[1], fyy[1], ... or is it like


 This one.

Matt


 fxx[0],fxx[1],...,fyy[0],fyy[1],...




 --
 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