Andrew, > 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.
It should be easy enough for the Tensors if you add the type as last template argument (defaulted to 'double'). I'm not sure what's missing for FullMatrix but that should be simple to find out :-) Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
