grads

[Roman Krylov]

Hi! Concerning GRADs - I've made only 2d vert2cell(attached) only to 
make PIC2D example available.
First, I did
for f in `ls Div*.h`;do cat $f|sed s/Div/Grad/g|sed s/div/grad/g>${f};done
:-)
well...
then I had converted that FieldStencil's functor GradVertToCell to be 
SCALAR->VECTOR from VECTOR->SCALAR operator.
and finally, implying OutputElement_t should b Vector<...>, curly braces 
member:

>   template<class F>
>   inline OutputElement_t
>   operator()(const F &f, int i1, int i2) const
>   {
>   return OutputElement_t(
>         0.5*(fact_m(0)*(
>         f(i1 + 1, i2    ) - f(i1    , i2    ) +
>         f(i1 + 1, i2 + 1) - f(i1    , i2 + 1))),
>         0.5*(fact_m(1)*(
>         f(i1    , i2 + 1) - f(i1    , i2    ) +
>         f(i1 + 1, i2 + 1) - f(i1 + 1, i2    )))
>         );

So, it exposed not to be unresolvable.
wbw,
Roman.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Grad.UR.h
Type: text/x-c-header
Size: 11151 bytes
Desc: not available
URL: <http://sourcerytools.com/pipermail/pooma-dev/attachments/20040712/a081e8c9/attachment.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Grad.h
Type: text/x-c-header
Size: 5899 bytes
Desc: not available
URL: <http://sourcerytools.com/pipermail/pooma-dev/attachments/20040712/a081e8c9/attachment-0001.bin>