Facet-spaces and hybrid methods iTutorial

Main Forum Old-Forum - do not post here
1

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(sigmanvhat+taun*uhat, element_boundary=True)
a.Assemble()

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

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

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

Best,
Dow

2

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(-sigmanvhat-taun*uhat, element_boundary=True)

c = Preconditioner(a, “bddc”)

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

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

gfu = GridFunction(X)[/code]

Thanks!
Dow