eval_jacobi_batch crashes on interval

[Kristian B. Ølgaard]

The following simple for crashes FIAT:

element = FiniteElement("Lagrange", interval, 3)

v = TestFunction(element)
u = TrialFunction(element)

a = inner(v, u)*dx

  File "/home/oelgaard/fenics/branches/dev/scripts/ffc", line 167, in <module>
    sys.exit(main(sys.argv[1:]))
  File "/home/oelgaard/fenics/branches/dev/scripts/ffc", line 155, in main
    compile_form(ufd.forms, prefix, options)
  File "/home/oelgaard/fenics/branches/dev/ffc/compiler.py", line 133, in compile_form
    ir = compute_ir(preprocessed_form, form_data, options)
  File "/home/oelgaard/fenics/branches/dev/ffc/representation.py", line 50, in compute_ir
    ir_elements = [compute_element_ir(e, form_data) for e in form_data.unique_sub_elements]
  File "/home/oelgaard/fenics/branches/dev/ffc/representation.py", line 65, in compute_element_ir
    element = create_element(ufl_element)
  File "/home/oelgaard/fenics/branches/dev/ffc/fiatinterface.py", line 67, in create_element
    element = create_fiat_element(ufl_element)
  File "/home/oelgaard/fenics/branches/dev/ffc/fiatinterface.py", line 79, in create_fiat_element
    element = ElementClass(cell, ufl_element.degree())
  File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/lagrange.py", line 35, in __init__
    poly_set = polynomial_set.ONPolynomialSet( ref_el , degree )
  File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/polynomial_set.py", line 166, in __init__
    dv = expansion_set.tabulate( degree , dpts )
  File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/expansions.py", line 117, in tabulate
    psitilde_as = jacobi.eval_jacobi_batch(0,0,n,ref_pts)
  File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/jacobi.py", line 53, in eval_jacobi_batch
    result[1,:] = 0.5 * ( a - b + ( a + b + 2.0 ) * xsnew )
ValueError: setting an array element with a sequence.

[Garth Wells (garth-wells)]

On Jan 11 2010, Robert Kirby wrote:

>** Changed in: fiat
> Status: New => Fix Committed
>
>

I don't see a change in the commit log. Has it been pushed?

Garth

[Robert Kirby (robert-c-kirby)]

oops...committed but did not push. Should be there now.

On Mon, Jan 11, 2010 at 11:49 AM, Garth Wells <email address hidden> wrote:

> On Jan 11 2010, Robert Kirby wrote:
>
> >** Changed in: fiat
> > Status: New => Fix Committed
> >
> >
>
> I don't see a change in the commit log. Has it been pushed?
>
> Garth
>
> --
> eval_jacobi_batch crashes on interval
> https://bugs.launchpad.net/bugs/504747
> You received this bug notification because you are a bug assignee.
>
> Status in FInite element Automatic Tabulator: Fix Committed
>
> Bug description:
> The following simple for crashes FIAT:
>
> element = FiniteElement("Lagrange", interval, 3)
>
> v = TestFunction(element)
> u = TrialFunction(element)
>
> a = inner(v, u)*dx
>
> File "/home/oelgaard/fenics/branches/dev/scripts/ffc", line 167, in
> <module>
> sys.exit(main(sys.argv[1:]))
> File "/home/oelgaard/fenics/branches/dev/scripts/ffc", line 155, in main
> compile_form(ufd.forms, prefix, options)
> File "/home/oelgaard/fenics/branches/dev/ffc/compiler.py", line 133, in
> compile_form
> ir = compute_ir(preprocessed_form, form_data, options)
> File "/home/oelgaard/fenics/branches/dev/ffc/representation.py", line 50,
> in compute_ir
> ir_elements = [compute_element_ir(e, form_data) for e in
> form_data.unique_sub_elements]
> File "/home/oelgaard/fenics/branches/dev/ffc/representation.py", line 65,
> in compute_element_ir
> element = create_element(ufl_element)
> File "/home/oelgaard/fenics/branches/dev/ffc/fiatinterface.py", line 67,
> in create_element
> element = create_fiat_element(ufl_element)
> File "/home/oelgaard/fenics/branches/dev/ffc/fiatinterface.py", line 79,
> in create_fiat_element
> element = ElementClass(cell, ufl_element.degree())
> File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/lagrange.py", line
> 35, in __init__
> poly_set = polynomial_set.ONPolynomialSet( ref_el , degree )
> File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/polynomial_set.py",
> line 166, in __init__
> dv = expansion_set.tabulate( degree , dpts )
> File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/expansions.py",
> line 117, in tabulate
> psitilde_as = jacobi.eval_jacobi_batch(0,0,n,ref_pts)
> File "/home/oelgaard/fenics/branches/fiat-new/FIAT_NEW/jacobi.py", line
> 53, in eval_jacobi_batch
> result[1,:] = 0.5 * ( a - b + ( a + b + 2.0 ) * xsnew )
> ValueError: setting an array element with a sequence.
>
>
>