| Re: [eigen] overriding scalar_X_op for auto-dif? |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] overriding scalar_X_op for auto-dif?
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Thu, 19 May 2011 23:42:09 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:from:date :message-id:subject:to:content-type:content-transfer-encoding; bh=2brkyegIvfcD9kfGTHQetmUwHSu1r9Hzj3H6cr9TgJ4=; b=nSRbGKaD3aymcIgLC2NPgyzaIjT3HI+EqbjA5dAdP/0+a/KpNzjMeFp75+9NgrJr1d v+Y2dJPI+uPw9PZphBwbDdUCRx0WGCC/VXP9ENyw2TUNzCfUywZ+zXftXE9uRRkgymKT fHKUD8kxV3RHsshgrsSQt9P+yykidednCws2I=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type:content-transfer-encoding; b=Xv8kfRGfLYVrv9rjccDZB+wuCfpZAA3Ur1d6+9+p1NBYF5ylL9lq4IpL3B1ADW/Whk aXxkyAlEchS685wg0oEE7zIX5kq4nVK4JVnEZ1LvGoHJ28qn/VTDBqpp/SqT9M8upQpn 3jjiVs/j+ClXWJ2TA2Xs+mA1Cy9lzx4Qc18cA=
ok, actually it seems this is something we have to solve in Eigen
itself. Nothing difficult though.
gael
On Thu, May 19, 2011 at 10:48 PM, Bob Carpenter <carp@xxxxxxxxxxx> wrote:
> Thanks so much for getting back to me so quickly
> on this. I tried overloading scalar_product_traits
> and NumTraits, but still no luck getting multiplication
> of an auto-dif matrix times a double matrix to compile.
>
> I'd MUCH rather just overload the basic ops for now.
> Later, I can optimize with expressions and hand-coded
> gradients for some cases without disturbing the API.
> I'm afraid your auto-dif scalar code's a bit beyond
> my current understanding of Eigen and C++ templates.
>
> I've attached a small tar.gz ball. Check out the
> README.txt, which shows how to run our reverse-mode
> auto-dif, has an example using Eigen auto-dif var
> multiplication which works fine, and then a demo that
> doesn't compile using double * auto-dif variables.
> I include the specialization of scalar_product_traits
> you suggested (with both orders), as well as NumTraits,
> but still no luck.
>
> - Bob Carpenter
> Columbia Uni, Dept of Statistics
>