search

HomerReid / scuff-em

A comprehensive and full-featured computational physics suite for boundary-element analysis of electromagnetic scattering, fluctuation-induced phenomena (Casimir forces and radiative heat transfer), nanophotonics, RF device engineering, electrostatics, and more. Includes a core library with C++ and python APIs as well as many command-line applications.
http://www.homerreid.com/scuff-em
GNU General Public License v2.0
128 stars 51 forks source link

Non-rectangular 2D lattices not yet supported ? #179

Closed CCherqui closed 6 years ago

CCherqui commented 6 years ago

So I've been working up to doing a triangular lattice of Ag nanodisks, when I built the lattice and tried it I got error: non-rectangular 2D lattices not yet supported (aborting).

Is there a timeline for this feature to be implemented?

HomerReid commented 6 years ago

No, at the moment there are several other things taking precedence over this, and you're the first person to ask for it, so I can't consider it a high priority---especially in the absence of any evidence that it will be particularly useful, or allow studies of systems that couldn't be done with other methods, or lead to high-impact publications, or otherwise justify the considerable expenditure of time it would require to implement.

But if you can furnish such evidence and would be willing to help with testing, let me know!

CCherqui commented 6 years ago

I can provide a lot of evidence, it's a very hot topic here at Northwestern and the entire reason I've been devoting so much time to understanding this wonderful code you've written. Also, I would be happy to help in any way I can. I've been using another BEM electrodynamics solver and what you've written is much faster.

Should I post such evidence in this issue thread or send a private message?

HomerReid commented 6 years ago

Sure, please describe what you have in mind. As a warmup for the case of non-rectangular lattices, is there any interesting calculation that can be done for rectangular lattices, or non-periodic geometries describing finite supercells of non-rectangular lattices? Note that, for something like Ag nanodiscs, it should be possible to do calculations involving thousands of unit cells, as in this example. Before discussing the possibility of adding new features it would be good to understand what can be done with the existing code and exactly why the new feature is required.

CCherqui commented 6 years ago

Yes, that's a good idea. Let me explore these large array examples first and see if this gets me to where I want to be.