flaport / fdtd

A 3D electromagnetic FDTD simulator written in Python with optional GPU support

Home Page:https://fdtd.readthedocs.io

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

outputing to VTK issue with interal grid structure

jstjep00 opened this issue · comments

Hey Floris,
First off love the repo, thank you for writing it.

I wanted to jump on implementing imageToVTK function from library pyevtk that is part of the visualization and then potentially creating a PR. I jumped on few issues regarding difference between points/cells within the VTK class references and the grid that is implemented within the code. As I understood both the flux and electric field are represented as cells within the grid, but imageToVTK function requires electric field to be in terms of points which there are the actual points of the grid cells (which there are nx + 1, ny + 1, nz + 1)

My question is whether or not this means the discretization algorithm needs to be changed or could this be sidestepped by expanding the shape that permittivity is expecting for.

King regards,
Josip