| [eigen] Array as the Scalar type |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen <eigen@xxxxxxxxxxxxxxxxxxx>
- Subject: [eigen] Array as the Scalar type
- From: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>
- Date: Fri, 7 May 2010 16:18:21 -0400
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:received:date:message-id :subject:from:to:content-type; bh=gi1zwlIwFneeR8XRID/pxzfvKngCPxsqJCuceGgB0jU=; b=KLG9ax93vNSS2pdKgzEDTeok/Drd/Y2SzS0dCtdDqa1WBXELweiqIq9++aR9EJkkeI HKLYdxRnEZ3pg1syy8rDjD/LURfizMk7XkN8k36EnJkgV3Ua1uVMM4C/Y149wU/pAXqv QNNF9PyBrhWdPPYewkpYVKY5tEdRcnL62mXgI=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:date:message-id:subject:from:to:content-type; b=ZzMdsSe6fl/TTfGkbUjIu8rM8kkCradPdVeZL7rlB1rlwfqAfHF4sWXBo2P6JQvNoB gf+TejeOxCeW92fGN39BxrWpNdigKjGNG9t47FMLh1vStU1oC3KtWLYVD7xAQoaqOBd/ ga+rJWWe1MUJ8lJSxkVbknge+eP5mVBry2OJQ=
Hi,
Mostly a question to Gael:
We wanted to enable support for Array as the scalar type, i.e.
Matrix<ArrayType, ...>.
But all our design is based on the assumption that it's OK to return
Scalars by value. E.g. the coeff() const methods in rvalue expressions
all return a Scalar by value.
So? Are we saying that using Array as Scalar type only is efficient
when the Array is small enough to fit entirely in registers? Should we
then discourage using large Array types as Scalar ? Note that would be
a bit disappointing as it would be the end of our "horizontal
vectorization" dreams.
Or do you want us to refactor Eigen so that we support expression
templates for the Scalar type?
In my rework of MathFunctions.h i have made sure to allow supporting
expression templates for the Scalar type, and have implemented that
for Array expressions, but if it's the only part of Eigen that does
that, it makes no sense! Should I drop xpr template Scalar support,
which would simplify MathFunctions.h quite a bit (e.g. no need for
SFINAE anymore).
Benoit