Hi all Following on from a post a while ago:
I was wondering if there was any interest in allowing deal.II FullMatrix and Tensors classes to be instantiated (via templates) on Sacado types (e.g. Sacado::Fad::DFad). The deal.II FullMatrix class is already templatized for a variety of floating point types whereas the Tensor classes are fixed for double precision. The main objective would be to construct tangent matrices from residual expressions that involve products of Tensors. Currently, you cano use Table to do this. This is bit of a pain for products of fourth and second-order tensors, but possible. 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. I chatted to Luca about this and he reckons it would be easy to implement. Any suggestions? Any complications that would make this a bad idea? Regards Andrew _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
