Facet-spaces and hybrid methods iTutorial

Hi,

I noticed an issue with the Facet-spaces and hybrid methods iTutorial. When input block 4 executes, an error occurs:

[code]sigma,u,uhat = X.TrialFunction()
tau,v,vhat = X.TestFunction()

a = BilinearForm(X)
a += SymbolicBFI(1/lam * sigmatau + div(sigma)v + div(tau)u)
n = specialcf.normal(mesh.dim)
a += SymbolicBFI(sigma
n
vhat+tau
n*uhat, element_boundary=True)
a.Assemble()

f = LinearForm(X)
f += SymbolicLFI(-sourcev)
f += SymbolicLFI(g
vhat, BND)
f.Assemble()
[/code]


RuntimeError Traceback (most recent call last)
in ()
11 f += SymbolicLFI(-sourcev)
12 f += SymbolicLFI(g
vhat, BND)
—> 13 f.Assemble()

RuntimeError: cannot evaluate facet-fe inside element, add trans simdin Assemble LinearForm

Best,
Dow

There is no longer an issue with the current version of the iTutorial. The current code in that block (shown below) runs without any errors.

[code]sigma,u,uhat = X.TrialFunction()
tau,v,vhat = X.TestFunction()

condense=True

a = BilinearForm(X, eliminate_internal=condense)
a += SymbolicBFI(1/lam * sigmatau + div(sigma)v + div(tau)u)
n = specialcf.normal(mesh.dim)
a += SymbolicBFI(-sigma
n
vhat-tau
n*uhat, element_boundary=True)

c = Preconditioner(a, “bddc”)

f = LinearForm(X)
f += SymbolicLFI(-sourcev)
f += SymbolicLFI(-g
vhat.Trace(), BND)

a.Assemble()
print (“A non-zero”, a.mat.nze)

gfu = GridFunction(X)[/code]

Thanks!
Dow