jit(TensorElement("CG", V.cell(), 1, shape=(1,1,3)))

Asked by David M. Rogers

I am getting a strange error when instantiating tensor function spaces.
TensorFunctionSpace(UnitInterval(4), "CG", 1, (1,2,3))
works, but
TensorFunctionSpace(UnitInterval(4), "CG", 1, (1,1,3))
does not.

It is coming from a call to jit of the tensor element:
>>> V = FunctionSpace(UnitInterval(4), "CG", 1)
>>> jit(TensorElement("CG", V.cell(), 1, shape=(1,2,3)))
=> works fine
>>> jit(TensorElement("CG", V.cell(), 1, shape=(1,1,3)))
Fixed index out of range!
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.6/site-packages/dolfin/compilemodules/jit.py", line 66, in mpi_jit
    return local_jit(*args, **kwargs)
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.6/site-packages/dolfin/compilemodules/jit.py", line 154, in jit
    return jit_compile(form, parameters=p, common_cell=common_cell)
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.6/site-packages/ffc/jitcompiler.py", line 71, in jit
    return jit_element(ufl_object, parameters)
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.6/site-packages/ffc/jitcompiler.py", line 174, in jit_element
    v = v[i]
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.6/site-packages/ufl/exproperators.py", line 310, in _getitem
    a = Indexed(self, indices)
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.6/site-packages/ufl/indexed.py", line 51, in __init__
    error("Fixed index out of range!")
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.6/site-packages/ufl/log.py", line 148, in error
    raise self._exception_type(self._format_raw(*message))
ufl.log.UFLException: Fixed index out of range!

I'm using the pre-compiled fenics dmg package for mac osx 10.6 fyi.

~ David M. Rogers.

Question information

Language:
English Edit question
Status:
Answered
For:
DOLFIN Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Marie Rognes (meg-simula) said :
#1

Sorry for the lag in response. This error is not present in the development version of FEniCS, please consider upgrading.

Can you help with this problem?

Provide an answer of your own, or ask David M. Rogers for more information if necessary.

To post a message you must log in.