-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Dear all,
until version 1.3 (where it still works) I have used the following
script to compute the discrete divergence of a given velocity field.
### ===============================================================
import dolfin
from dolfin import dx, div
import numpy as np
N = 10
mesh = dolfin.UnitSquareMesh(N, N)
# print 'we use Crouzieux-Raviart elements !'
# V = dolfin.VectorFunctionSpace(mesh, "CR", 1)
# Q = dolfin.FunctionSpace(mesh, "DG", 0)
print 'we use Taylor-Hood elements'
V = dolfin.VectorFunctionSpace(mesh, "CG", 2)
Q = dolfin.FunctionSpace(mesh, "CG", 1)
nv = V.dim()
# make up a test vector
velvec = np.ones((nv, 1))
v = dolfin.Function(V)
v.vector().set_local(velvec)
q = dolfin.TrialFunction(Q)
divv = dolfin.assemble(q*div(v)*dx)
### ==================================================
Since 1.4 it fails with exception (full traceback is attached)
**Exception: Unable to extract all indices.**
Last time I checked for this I thought it is related to this issue:
https://bitbucket.org/fenics-project/dolfin/issue/352
However, also in version 1.5 this error is still there.
Any ideas? Workarounds??
Best regards,
Jan
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1
iQEcBAEBAgAGBQJU1ei3AAoJED0q33mesbTn0QsH/RsgiUEUeWh/D6Y54LNBZPOR
d8si7yhtYQyu9+eYccknXpf9x4HhciUdR8RKw3bwGcvy2H/nkXT1EgIi5CokvIhq
OYkcZxWppPiANdXR5WI/qwRe14peikKcjzhnV2Y2SbiNkQxkqwakpaDEteKxIlLC
FiPRCVeVmTbE2sWdCJHWBvRIKNH/s3bnCydFolmH+RH8XOKk1sD8oEatoHKLqSQ1
/BaMk0oJ+YOfR1sk/lqB7JO4Z5o4HbujezVNHrkbmaPVoW0QyzWNq+L3xd9zt7qg
dqzguiIjsFkkjFXFHSM/qZr9hEDadiVzoy2gOtVb8H+ev8Z2wqSqwRKdDe1t67c=
=EHac
-----END PGP SIGNATURE-----
we use Taylor-Hood elements
Calling FFC just-in-time (JIT) compiler, this may take some time.
Unable to extract all indices.
Traceback (most recent call last):
File "test_div.py", line 23, in <module>
divv = dolfin.assemble(q*div(v)*dx)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/dolfin/fem/assembling.py",
line 193, in assemble
cell_domains, exterior_facet_domains, interior_facet_domains)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/dolfin/fem/assembling.py",
line 75, in _create_dolfin_form
form_compiler_parameters=form_compiler_parameters)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/dolfin/fem/form.py",
line 70, in __init__
mpi_comm=mesh.mpi_comm())
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/dolfin/compilemodules/jit.py",
line 64, in mpi_jit
return local_jit(*args, **kwargs)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/dolfin/compilemodules/jit.py",
line 128, in jit
return form_compiler.jit(form, parameters=p)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/jitcompiler.py",
line 74, in jit
return jit_form(ufl_object, parameters)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/jitcompiler.py",
line 130, in jit_form
parameters=parameters)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/compiler.py",
line 156, in compile_form
ir = compute_ir(analysis, parameters)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/representation.py",
line 94, in compute_ir
for (i, fd) in enumerate(form_datas)]
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/representation.py",
line 218, in _compute_integral_ir
parameters)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/tensor/tensorrepresentation.py",
line 80, in compute_integral_ir
terms = compute_terms(None, None)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/tensor/tensorrepresentation.py",
line 75, in <lambda>
facet_cellname)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/tensor/tensorrepresentation.py",
line 125, in _compute_terms
facet_cellname)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/tensor/referencetensor.py",
line 63, in __init__
self.primary_multi_index = create_multiindex(primary_indices)
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ffc/tensor/multiindex.py",
line 55, in create_multiindex
error("Unable to extract all indices.")
File "<string>", line 1, in <lambda>
File
"/afs/mpi-magdeburg.mpg.de/data/csc/software/64bit/fenics-1.5.0-ubuntu1204/lib/python2.7/site-packages/ufl/log.py",
line 151, in error
raise self._exception_type(self._format_raw(*message))
Exception: Unable to extract all indices.
_______________________________________________
fenics-support mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics-support
