Hi everyone,
I’m very new to FEniCS. Assemble does not work properly getting some error in line writing likes this given below.
So, I modified the code, but I can’t solve this error. Help me~~
from dolfin import *
from ufl import *
import sys
class Forms(object):
def V0constrainedE(self):
mesh = self.parameters["mesh"]
u = self.parameters["displacement_variable"]
ds = Measure('ds',subdomain_data = self.parameters["facetboundaries"])
dsendo = ds(self.parameters["endoid"], domain = self.parameters["mesh"], subdomain_data = self.parameters["facetboundaries"])
pendo = self.parameters["volconst_variable"]
V0= self.parameters["constrained_vol"]
X = SpatialCoordinate(mesh)
x = u + X
F = self.Fmat()
N = self.parameters["facet_normal"]
n = cofac(F)*N
area = assemble(Constant(1.0) * dsendo)
V_u = -Constant(1.0/3.0) * inner(x, n)
Wvol = (Constant(1.0/area) * pendo * V0 * dsendo) - (pendo * V_u *dsendo)
return (Wvol)
rm: cannot remove '*.pvtu': No such file or directory
Invalid cell 1.0.
Traceback (most recent call last):
File "LV_closed.py", line 172, in <module>
Wvol = uflforms.V0constrainedE()
File "/home/dmryu/FEniCS/forms.py", line 94, in V0constrainedE
area = assemble(Constant(1.0) * dsendo)
File "/usr/lib/python3/dist-packages/ufl/coefficient.py", line 125, in Constant
domain = as_domain(domain)
File "/usr/lib/python3/dist-packages/ufl/domain.py", line 296, in as_domain
cell = as_cell(domain)
File "/usr/lib/python3/dist-packages/ufl/cell.py", line 329, in as_cell
error("Invalid cell %s." % cell)
File "/usr/lib/python3/dist-packages/ufl/log.py", line 172, in error
raise self._exception_type(self._format_raw(*message))
ufl.log.UFLException: Invalid cell 1.0.