New question #153352 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/153352
I don't know if it makes sense but I'd like to do:
m = UnitInterval(11)
Q= FunctionSpace(m,'CG',1)
W = VectorFunctionSpace(m,'CG',1,dim=3)
f = Function(Q)
df = grad(f)
g = project(df,W)
but I get the errpr below. Any help?
Thanks,
Chaffra
/usr/local/lib/python2.6/dist-packages/dolfin/fem/project.py in project(v, V,
mesh, solver_type, preconditioner_type, form_compiler_parameters)
50 Pv = TrialFunction(V)
51 a = ufl.inner(w, Pv)*ufl.dx
---> 52 L = ufl.inner(w, v)*ufl.dx
53
54 # Assemble linear system
/usr/local/lib/python2.6/dist-packages/ufl/operators.pyc in inner(a, b)
56 if a.shape() == () and b.shape() == ():
57 return a*b
---> 58 return Inner(a, b)
59 #return contraction(a, range(a.rank()), b, range(b.rank()))
60
/usr/local/lib/python2.6/dist-packages/ufl/tensoralgebra.pyc in __new__(cls, a,
b)
160
161 def __new__(cls, a, b):
--> 162 ufl_assert(a.shape() == b.shape(), "Shape mismatch.")
163 if isinstance(a, Zero) or isinstance(b, Zero):
164 free_indices, index_dimensions = merge_indices(a, b)
/usr/local/lib/python2.6/dist-packages/ufl/assertions.pyc in
ufl_assert(condition, *message)
18 def ufl_assert(condition, *message):
19 "Assert that condition is true and otherwise issue an error with
given message."
---> 20 if not condition: error(*message)
21
22
/usr/local/lib/python2.6/dist-packages/ufl/log.pyc in error(self, *message)
122 "Write error message and raise an exception."
123 self._log.error(*message)
--> 124 raise UFLException(self._format_raw(*message))
125
126 def begin(self, *message):
UFLException: Shape mismatch.
--
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
