Hi,
I tried to run the following code from MiroK
from dolfin import *
mesh = UnitSquareMesh(2, 2)
V = VectorFunctionSpace(mesh, 'CG', 1)
v = TestFunction(V)
b = assemble(inner(Constant((0, 0)), v)*dx)
x = Point(0.25, 0.25)
# Apply to components
p0 = PointSource(V.sub(0), x, 1) # Assuming that V.sub(0) is a scalar-valued space!
p1 = PointSource(V.sub(1), x, 2) # Assuming that V.sub(1) is a scalar-valued space!
p0.apply(b)
print b.array()
p1.apply(b)
print b.array()
and I get the error: Segmentation fault: 11
I have Fenics 1.3.0 on Mac OSX 10.7.5, but there is the same problem on another computer with Linux (also with Fenics 1.3.0)
Also I get the following error message:
Python quit unexpectedly while using the libdolfin.1.3.dylib plug-in.
To be precise, the error happens exactly at the line
p0.apply(b)
Do you think it is because the version is too old?
Do you have the same problem?
thanks