- Nov 20, 2015
-
-
Christian Engwer authored
In the current state the scaling parameter for matrix-vector products (e.g. the alpha parameter for DenseMatrix<...>::usmv) is the field_type of the matrix. I you consider mixed precision, like in the xblas library, you would actually want the scaling to be of the same accuracy as the output vector; an other reason is vectorization for one matrix and multiple vectors. The patch contains two changes: a) the FieldTraits are extended, so that the contain specializations for C and C++ vectors b) the thre matrix implementations are changed, such that they deduce the type as FieldTraits<Y>::field_type, where Y is the type of the result vector. In the case of identity matrix the old implementation was laso completely wrong, as it used Y::axpy to implement usmv and Y::axpy expects alpha to be Y::value_type, so that we already had a type clash, which just went unnoticed, as we never used mixed-types.
-
- Oct 06, 2014
-
-
Ansgar Burchardt authored
-
- May 31, 2012
-
-
Oliver Sander authored
With the exception of mpicollectivecommunication.hh and propertymap.hh, which I didn't understand. [[Imported from SVN: r6784]]
-
- Nov 29, 2010
-
-
Andreas Dedner authored
revert r6278 and r6277 - if possible we should not keep the trunk dysfunctional for more than a few hours (I need alugrid) [[Imported from SVN: r6280]]
-
Markus Blatt authored
specify types as typedefs. Someghow I missed this file. [[Imported from SVN: r6278]]
-
- Oct 11, 2010
-
-
Christian Engwer authored
[[Imported from SVN: r6163]]
-
- Sep 08, 2010
-
-
Oliver Sander authored
[[Imported from SVN: r6137]]
-
Christian Engwer authored
fix typo [[Imported from SVN: r6136]]
-
Christian Engwer authored
[[Imported from SVN: r6135]]
-
- Oct 15, 2009
-
-
Christian Engwer authored
type * make *norm* in fvector and fmatrix use this real type [[Imported from SVN: r5656]]
-