Define a non-cartesian mesh in fipy

[Henri Colaux]

Dear Fipy users,

I would like to draw your attention on the question I asked on stackoverflow :

https://stackoverflow.com/questions/69196642/define-a-non-cartesian-mesh-in-fipy

I am trying to simulate the elementary unit of a 2D system that have a P6mm symmetry in fipy and I would like to define a non-cartesian mesh that describes the system described below. Yet,

mesh = fipy.Grid2D(nx = 10, ny = 10, dx = 1., dy = 1.)

only returns uniform meshes. I thought of changing the FaceVariable, but it seems it only accepts Boolean variables. I could also simulate a cartesian system equivalent to this one, but there would be redondant data. Would someone have a better approach?

You can see the picture on the stackoverflow question.
Cheers,

Henri Colaux

[Henri Colaux]

As advised by Jonathan Guyer, I added an issue to GitHub:
https://github.com/usnistgov/fipy/issues/819

[Henri Colaux]

As a follow-up question: I could simulate my system using a triangular mesh this way. Is there a problem doing so ?

[Henri Colaux]

[Henri Colaux]

My apologies for the spam. I did not realize that it was not possible to include pictures.
I added the schematic in the stackoverflow question.
https://stackoverflow.com/questions/69196642/define-a-non-cartesian-mesh-in-fipy