IndexError and Mesh Artifacts about femwell HOT 3 CLOSED

Hey Tristan :)

to me it seems that the code you posted works properly as the mesh looks fine.

The problem seems to be in the part where you define the epsilon. Could you also post that part? I'd guess you use a wrong Basis or something like that
(if so, I'll convert this to a discussion, as it doesn't seem to be a problem with femwell)

from femwell.

@st4eve we were following the docs, but it's possible something was off in the element definition

from femwell.

basis0 = Basis(mesh, ElementTriP0()) epsilon = basis0.zeros() for subdomain, n in {"sin": 1.9963, "box": 1.444, "clad": 1}.items(): epsilon[basis0.get_dofs(elements=subdomain)] = n**2 basis0.plot(epsilon, colorbar=True).show() plt.show()

This should be all of it. Apologies

from femwell.