GDSII import for Scheme

Author: stevengj

scheme/meep-ctl-swig.hpp

@@ -55,6 +55,10 @@ ctlio::cnumber_list make_casimir_g(double T, double dt, double sigma, meep::fiel
 
 ctlio::cnumber_list make_casimir_g_kz(double T, double dt, double sigma, meep::field_type ft);
 
+ctlio::integer_list GDSII_layers(const char *fname);
+meep::volume GDSII_vol(const char *fname, int layer = -1, double zmin = 0, double zmax = 0);
+ctlio::simple_prism_list get_GDSII_prisms(const char *GDSIIFile, int Layer = -1, double zmin = 0.0, double zmax = 0.0);
+
 // wrapper around constructor to fool SWIG
 meep::volume_list *make_volume_list(const meep::volume &v, int c, std::complex<double> weight,
                                     meep::volume_list *next);

scheme/meep.cpp

@@ -1,5 +1,12 @@
 #include "meep-ctl.hpp"
 
+/* hack to tell meepgeom.hpp, which includes ctlgeom.h,
+   to use the ctlgeom.h object definitions and not the
+   ctl-io.h definitions */
+#undef CXX_CTL_IO
+#undef CTL_IO_H
+#include "meepgeom.hpp"
+
 using namespace meep;
 using namespace std;
 
@@ -198,3 +205,36 @@ ctlio::number_list std_vector_double_to_scm(std::vector<double> *v) {
     res.items[i] = (*v)[i];
   return res;
 }
+/***************************************************************************/
+
+ctlio::integer_list GDSII_layers(const char *fname) {
+  std::vector<int> layers = meep_geom::get_GDSII_layers(fname);
+  ctlio::integer_list res;
+  res.num_items = int(layers.size());
+  res.items = new integer[res.num_items];
+  for (int i = 0; i < res.num_items; ++i) {
+    res.items[i] = layers[i];
+  }
+  return res;
+}
+
+volume GDSII_vol(const char *fname, int layer, double zmin, double zmax) {
+  return meep_geom::get_GDSII_volume(fname, layer, zmin, zmax);
+}
+
+ctlio::simple_prism_list get_GDSII_prisms(const char *GDSIIFile, int Layer, double zmin, double zmax) {
+  geometric_object_list go = meep_geom::get_GDSII_prisms(NULL, GDSIIFile, Layer, zmin, zmax);
+  ctlio::simple_prism_list res;
+  res.num_items = go.num_items;
+  res.items = new ctlio::simple_prism[res.num_items];
+  for (int i = 0; i < res.num_items; ++i) {
+    prism *p = go.items[i].subclass.prism_data;
+    res.items[i].vertices.num_items = p->vertices.num_items;
+    res.items[i].vertices.items = p->vertices.items;
+    res.items[i].height = p->height;
+    res.items[i].axis = p->axis;
+    free(p);
+  }
+  delete[] go.items;
+  return res;
+}

scheme/meep.scm.in

@@ -1105,6 +1105,18 @@
 (define (output-png c options) (output-png-rm? true c options))
 (define (output-png+h5 c options) (output-png-rm? false c options))
 
+; ****************************************************************
+; GDSII import functions
+
+(export-type (make-list-type 'integer))
+
+(define-class simple-prism no-parent
+  (define-property vertices no-default (make-list-type 'vector3))
+  (define-property height no-default 'number)
+  (define-property axis no-default 'vector3))
+(export-type (make-list-type 'simple-prism))
+
+
 ; ****************************************************************
 ; harminv functions for extracting bands, etcetera