At 03:56 PM 6/18/98 -0400, jeremic@sokocalo.cee.clarkson.edu wrote:
>The main operator* has to know what to contract or expand so parsing
>part must be there. You can make many small functions that will do
>specific contractions or expansions, like for example:
>
>tensor test = A("ijkl") * B("klpq");
>
>might be done with:
>
>tensor test = TensorMultiply_0011_1100(A,B);
The idea of putting a bunch of stuff in the types of the tensors is that
you could then do:
tensor test = TensorMultiply(A,B)
and let the overload resolution take care of the rest. This does require
some trickery in the types -- it probably can't be reasonably done without
an engine-like technique.
-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