Hi all,
When I run my code in my lab's server, which is on Debian, I got an error as following,
"RuntimeError: '*** Error: All dofmaps must have same block size (for now)'",
which is at the line,
n3_express = inner(grad(p),v)*dx
n3_mat = assemble(n3_express)
The whole codes related are posted as following,
import numpy as np
import math
from dolfin import *
## Build the fluid dynamic model
# Create mesh
nx = 6
ny = 6
mesh = UnitSquareMesh(nx, ny)
Lend = 0
Rend = 1
# Define function spaces
V = VectorFunctionSpace(mesh, "CG", 2)
Q = FunctionSpace(mesh, "CG", 1)
T = FunctionSpace(mesh, "CG", 1)
# Define test and trail functions
u = TrialFunction(V)
v = TestFunction(V)
p = TrialFunction(Q)
q = TestFunction(Q)
te = TrialFunction(T)
de = TestFunction(T)
n3_express = inner(grad(p),v)*dx
n3_mat = assemble(n3_express)
n3 = n3_mat.array()
It seems to be a very problem with the OS, since when I run my code in Ubuntu, there is no such error. How can I fix it? Thanks a lot for your help!
