Re: Gradient at ghost cells only?

["Matthew Knepley" <knepley@mcs.anl.gov>]

I think it would be really useful to look at the DA code first. This has
all the correct structures, but a simple implementation for structured
grids. For instance, the gradient is the action of the Jacobian on a
vector. That is really how it should be coded. You have a bunch of
local actions, and then a global assembly.

   Matt

On 3/4/06, billy@dem.uminho.pt <billy@dem.uminho.pt> wrote:

Hi,

I need the gradient of 5 scalars at ghost cells and I was thinking of creating
5x3 vectors with VecCreateGhost(). Right now, I calculate the gradient during
each iteration and I don't store it's value.
I would like to calculate, store and exchange the gradient at the ghost cells
only, to minimize memory requirements. How is the best way to do that?

Billy.

-- 
"Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness