Andrew,
> I think a neat way to do this would be to follow the suggestions in
> the Template instantiations documentation that allow users to include
> there own template arguments. The only issue then would be to
> overload mathematical operators such as std::sqrt to allow for use of
> the corresponding Sacado operators. The changes to the Tensor class
> would be as those that allow FullMatrices to be templatized on type.
This looks like a great idea, I'd be willing to help implement it. One
of the weird-and-wonderful things I have wanted is a Tensor of integers;
which fits into your suggestion (but is slightly less complicated and
does not involve Sacado).
Instantiations for double, float (and int?) would be a good starting
point to get going later with Sacado. Let me know if I can help out as
far as getting the instantiation model setup here, ok?
> I chatted to Luca about this and he reckons it would be easy to
> implement.
A short concentrated effort is probably all this needs.
Best,
Toby
--
Toby D. Young
Assistant Professor Philosophy & Physics
Polish Academy of Sciences
www: http://www.ippt.gov.pl/~tyoung
skype: stenografia
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
