Issue on page /julia/waveguide_anisotropic.html

[muchalito]

Hello,

Can the Python/ScikitFem version of the code model anisotropic materials as well, or is this only available in Julia?

I'm looking at how to modify the statement in Python to assign the material parameters as 3x3 matrices but it flags an error as it's expecting only scalars:

basis0 = Basis(mesh, ElementTriP0())
epsilon = basis0.zeros()
if True:
for subdomain, n in {"core": index_core, "box": index_box, "clad": index_cladding}.items():
epsilon[basis0.get_dofs(elements=subdomain)] = (n**2)*np.eye(3)
basis0.plot(epsilon, colorbar=True).show()

[jetremblay]

Hello, I don't think anisotropic materials are supported in the Python version. I did start working on the feature, would there be some interest in a pull request?

[muchalito]

Yes, there would be interest in a pull request from me, at least. I wanted to use anisotropy to build a PML for accurate bend loss simulations.

For this we'd need at least "biaxial" anisotropy, or all diagonal terms of the tensors can be different with off-diagonals still zero. Full anisotropy would work too.

Thanks!

[simbilod]

Happy to review PRs!