New question #150875 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/150875
I am trying to solve for a system of equations with two symmetric tensors as
unknowns, and encountering a number of errors with a MixedFunctionSpace
consisting of two TensorFunctionSpaces with symmetry when assembling.
I have written a short piece of code which creates the same error:
from dolfin import *
mesh = Rectangle(0,0,1,1,10,10)
dim = 2 #assume 2D
symm = dict(((i,j), (j,i))
for i in range(dim) for j in range(dim) if i > j )
T1 = TensorFunctionSpace(mesh, 'CG', 1, symmetry=symm)
T2 = TensorFunctionSpace(mesh, 'CG', 1, symmetry=symm)
TT = T1*T2
R, S = TrialFunctions(TT)
v_R, v_S = TestFunctions(TT)
P = Function(T1)
Q = Function(T2)
Fr = inner(R,v_R)*dx+ inner(P,v_R)*dx
Fs = inner( S, v_S )*dx + inner( Q, v_S )*dx
F = Fr + Fs
A = lhs(F)
B = rhs(F)
a = Matrix()
a = assemble(A, tensor=a)
b = Vector()
b = assemble(B, tensor=b)
I get the following error:
Calling FFC just-in-time (JIT) compiler, this may take some time.
Traceback (most recent call last):
File "tensortest2.py", line 28, in <module>
a = assemble(A, tensor=a)
File "/usr/lib/python2.6/dist-packages/dolfin/fem/assemble.py", line 100, in
assemble
common_cell=common_cell)
File "/usr/lib/python2.6/dist-packages/dolfin/fem/form.py", line 34, in
__init__
(self._compiled_form, module, self.form_data) = jit(form,
form_compiler_parameters, common_cell)
File "/usr/lib/python2.6/dist-packages/dolfin/compilemodules/jit.py", line
47, in mpi_jit
return local_jit(*args, **kwargs)
File "/usr/lib/python2.6/dist-packages/dolfin/compilemodules/jit.py", line
114, in jit
return jit_compile(form, parameters=p, common_cell=common_cell)
File "/usr/lib/python2.6/dist-packages/ffc/jitcompiler.py", line 64, in jit
return jit_form(object, parameters, common_cell)
File "/usr/lib/python2.6/dist-packages/ffc/jitcompiler.py", line 122, in
jit_form
compile_form(preprocessed_form, prefix=jit_object.signature(),
parameters=parameters)
File "/usr/lib/python2.6/dist-packages/ffc/compiler.py", line 140, in
compile_form
ir = compute_ir(analysis, parameters)
File "/usr/lib/python2.6/dist-packages/ffc/representation.py", line 66, in
compute_ir
irs = [_compute_integral_ir(f, i, parameters) for (i, f) in
enumerate(forms)]
File "/usr/lib/python2.6/dist-packages/ffc/representation.py", line 186, in
_compute_integral_ir
parameters)
File "/usr/lib/python2.6/dist-packages/ffc/tensor/tensorrepresentation.py",
line 59, in compute_integral_ir
ir["AK"] = _compute_terms(monomial_form, None, None, domain_type,
quadrature_degree)
File "/usr/lib/python2.6/dist-packages/ffc/tensor/tensorrepresentation.py",
line 98, in _compute_terms
quadrature_degree)
File "/usr/lib/python2.6/dist-packages/ffc/tensor/referencetensor.py", line
28, in __init__
self.A0 = integrate(monomial, domain_type, facet0, facet1, quadrature_order)
File "/usr/lib/python2.6/dist-packages/ffc/tensor/monomialintegration.py",
line 50, in integrate
psis = [_compute_psi(v, table, len(points), domain_type) for v in
monomial.arguments]
File "/usr/lib/python2.6/dist-packages/ffc/tensor/monomialintegration.py",
line 169, in _compute_psi
Psi[component][tuple(dlist)] = etable[dtuple][:,
cindex[0].index_range[component], :]
What's going wrong here?
--
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : dolfin@lists.launchpad.net
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
