Here is another idea on stencils which could be interesting but we haven't
fleshed it out very far yet.
The idea is to express a stencil in some way like the following:
Array<3> A, B;
A = apply_stencil(coef(j)*B(j),SixPoint(j));
The idea here is that apply_stencil(expression,footprint) would invoke some
sort of template metaprogram that would evaluate the expression at the
points specified by the footprint.
In a POOMA II, context, that means it would build an Array with an engine
which, when evaluated at a point, would evaluate a template metaprogram
which does the stencil. Like the previous, this could be wrapped so that
the user says
A = laplace(B).
This would take some trickery, but it could be done. Something that this
does is it separates the specification of the stencil footprint and the
calculation using the footprint. This may or may not be an advantage.
Right now, the other stencil technique looks much simpler to me, and I
don't see really major advantages to this one. What do others think?
-Steve Karmesin
======================================================================
Pooma Team Leader
karmesin@lanl.gov
phone: (505)665-6019
fax : (505)667-4939
For PGP public key, send message with subject line: SEND PUBLIC KEY
======================================================================
--------------------- blitz-dev list --------------------------------
* To subscribe/unsubscribe: mail to majordomo@oonumerics.org, with
"subscribe blitz-dev" or "unsubscribe blitz-dev" in the body of the message
* Blitz++ web page: http://oonumerics.org/blitz/
This archive was generated by hypermail 2b29 : Wed Feb 20 2002 - 04:30:05 EST