Comment 1 for bug 495602

Not sure. I've been having these bugs crop up.
I have company this weekend and might not get to this.
The quick attempt to work-around is to replace numpy.sqrt with
math.sqrt.

Rob

On Fri, Dec 11, 2009 at 12:58 PM, Garth Wells <email address hidden> wrote:

> Public bug reported:
>
> When trying to compile an elasticity problem using the FIAT_NEW and
> experimental FFC, I get
>
> File "/usr/local/lib/python2.6/dist-packages/FIAT_NEW/expansions.py",
> line 400, in tabulate
> results[idx(p,q,r)] *= numpy.sqrt((p+0.5)*(p+q+1.0)*(p+q+r+1.5))
> TypeError: unsupported operand type(s) for *=: 'numpy.ndarray' and
> 'numpy.float64'
>
> Is this a FIAT or a NumPy bug?
>
> ** Affects: fiat
> Importance: Undecided
> Assignee: Robert Kirby (robert-c-kirby)
> Status: New
>
> ** Changed in: fiat
> Assignee: (unassigned) => Robert Kirby (robert-c-kirby)
>
> ** Branch linked: lp:~fiat-core/fiat/fiat-new
>
> --
> FIAT_NEW + NumPy 64bit problem
> https://bugs.launchpad.net/bugs/495602
> You received this bug notification because you are a bug assignee.
>
> Status in FInite element Automatic Tabulator: New
>
> Bug description:
> When trying to compile an elasticity problem using the FIAT_NEW and
> experimental FFC, I get
>
> File "/usr/local/lib/python2.6/dist-packages/FIAT_NEW/expansions.py",
> line 400, in tabulate
> results[idx(p,q,r)] *= numpy.sqrt((p+0.5)*(p+q+1.0)*(p+q+r+1.5))
> TypeError: unsupported operand type(s) for *=: 'numpy.ndarray' and
> 'numpy.float64'
>
> Is this a FIAT or a NumPy bug?
>
>
>