|Re: [eigen] news on the refactoring of the expression template mechanism|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: Eigen mailinglist <eigen@xxxxxxxxxxxxxxxxxxx>
- Subject: Re: [eigen] news on the refactoring of the expression template mechanism
- From: Jitse Niesen <jitseniesen@xxxxxxxxx>
- Date: Sat, 12 Apr 2014 13:56:11 +0100 (BST)
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=yahoo.com; s=s1024; t=1397307372; bh=oQv8/Gm/qhAkc5l+eFhkwy58Z8i1S9GZLlDpTd9mitk=; h=X-Yahoo-Newman-Id:X-Yahoo-Newman-Property:X-YMail-OSG:X-Yahoo-SMTP:X-Rocket-Received:Date:From:X-X-Sender:To:Subject:In-Reply-To:Message-ID:References:User-Agent:MIME-Version:Content-Type; b=BI81gWRllUKaN6vV+fZQEvdYWjbgg46Ue9x1ycZjHDES9yVlUADvvUHttoMHXKAJB+mxzCQ5NgW+M531QEhPM5bkGY1qs+QYsBEXTvpt7xpN2Z6jqf2luP2Bn+cXV9BhJQjRyIarrQqaG8SqiNbPK0x5egeayzU0P+STL0qC0lI=
Some initial reactions on the open questions on the wiki page:
This is more technical, so I put it in a separate message.
What about DirectAccessBit?
Currently it is still handled by the expressions because DenseCoeffBase
need it to enable data/outerStride/innerStride and the likes. However,
it not clear that the propagation of DirectAccess is generic and it
seems to rather be an evaluation details.
Yes, it is more naturally part of evaluators. One possibility is to have
expr.data() be defined for all expression objects and have it call
evaluator.data() giving a compile-time error if the expression is not
direct access and the evaluator objects does not have a .data() method.
This will not give a very nice compile-time error though so the other
possibility is to have the expression object depend on DirectAccessBit as
defined in the evaluator flag in order to enable .data(). This introduces
a dependency, but that seems to me not too bad.
What to do with TriangularView::solve<OnTheRight>(lhs)?
We want to deprecate this API and use lhs * tri.inverse().
So TriangularView::solve<OnTheRight> could return a
Product<Lhs,Inverse<TriangularView>> instead of a Solve<> expression for
the "OnTheLeft" case.
I feel a bit uncomfortable telling people to write inverse to solve
because I (and others) teach people not to use inverse. Of course I know
(or at least hope!) that Eigen does not actually compute the inverse but
that may be a subtlety lost on some.
On the other hand, solve on the right is not a very clear terminology; I
still do not know immediately what's on the left and what's on the right.
Clarify StorageKind versus XprShape
I am a bit lost here with the terminology. My understanding is the