Facet-spaces and hybrid methods iTutorial

ddrake
Topic Author
Offline
Premium Member

7 years 2 months ago #587 by ddrake

Facet-spaces and hybrid methods iTutorial was created by ddrake

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 * sigma*tau + 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(-source*v)
f += SymbolicLFI(g*vhat, BND)
f.Assemble()

RuntimeError Traceback (most recent call last)
<ipython-input-4-44aa5d55b0e0> in <module>()
11 f += SymbolicLFI(-source*v)
12 f += SymbolicLFI(g*vhat, BND)
---> 13 f.Assemble()

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

Best,
Dow

ddrake
Topic Author
Offline
Premium Member

7 years 1 month ago #604 by ddrake

Replied by ddrake on topic Facet-spaces and hybrid methods iTutorial

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 * sigma*tau + 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(-source*v)
f += SymbolicLFI(-g*vhat.Trace(), BND)

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

gfu = GridFunction(X)

Thanks!
Dow

Time to create page: 0.095 seconds