Question about Non-Cubic (Spherical) Mesh Support

[hz-xiaxz]

Question about Non-Cubic (Spherical) Mesh Support

Thanks for creating this excellent package!

Current Need

I'm working with a 3D Green's function that has a UV cutoff |k| < Λ, requiring a spherical mesh in k-space.

Question

1. Does MeshArray support non-cubic meshes (specifically spherical coordinates)? I've reviewed the documentation and source code but only found cubic-like mesh examples.

2. If spherical meshes aren't directly supported, what's the practical impact of using a cubic mesh with a spherical cutoff in physical calculations? Specifically:

   • Would using a cubic mesh be a reasonable approximation?
   • Are there known numerical issues or artifacts when using this approach?

Thanks for your time and help!

[kunyuan]

Thank you for your post.

The package is compatible with non-cubic meshes.

You can use the CompositeGrids.jl package to create three 1D meshes for ( k ), ( \theta ), and ( \phi ), respectively. Then, you can use these three meshes to create a Green's function object using the syntax: MeshArray(k, theta, phi, tau; dtype=Float64).

The Green's function object, "MeshArray," is essentially a high-dimensional array, with each dimension corresponding to a specific grid.

If you have any further questions, please let me know.

[hz-xiaxz]

Ah, yeah! You are definitely right. Thanks for your help!