add support for prisms (new type of geometric_object) to libctlgeom (#13)

* add support for prisms (new type of geometric_object) to libctlgeom

* updates

* updates to point_in_prism and other functions

* implemented intersect_line_with_object for prisms; added working unit test to test-prism.c

* updates

* new function intersect_line_segment_with_object

* finished implementing all geometry primitives for prisms; added passing unit test utils/test-prism.c

* finished implementing all geometry primitives for prisms; added passing unit test utils/test-prism.c

* inserted content of prism.c into geom.c, made updates to enable mpb and meep builds

* fixed compilation issues for test-prism

Author: Homer Reid

utils/Makefile.am

@@ -14,6 +14,13 @@ geomtst_SOURCES = geomtst.c
 geomtst_LDADD = libctlgeom.la
 geomtst_CPPFLAGS = -I$(top_srcdir)/src
 
+check_PROGRAMS = test-prism
+
+test_prism_SOURCES = test-prism.c
+test_prism_LDADD   = libctlgeom.la
+
+TESTS = test-prism
+
 dist_man_MANS = gen-ctl-io.1
 
 BUILT_SOURCES = gen-ctl-io geom-ctl-io.c ctlgeom-types.h nlopt-constants.scm

utils/ctlgeom.h

@@ -69,6 +69,9 @@ extern vector3 normal_to_object(vector3 p, GEOMETRIC_OBJECT o);
 extern vector3 normal_to_fixed_object(vector3 p, GEOMETRIC_OBJECT o);
 extern int intersect_line_with_object(vector3 p, vector3 d, GEOMETRIC_OBJECT o,
 				      double s[2]);
+extern double intersect_line_segment_with_object(vector3 p, vector3 d,
+                                                 GEOMETRIC_OBJECT o,
+                                                 double a, double b);
 extern MATERIAL_TYPE material_of_point_inobject(vector3 p, boolean *inobject);
 extern MATERIAL_TYPE material_of_point_inobject0(
      GEOMETRIC_OBJECT_LIST geometry, vector3 p, boolean *inobject);
@@ -144,7 +147,11 @@ GEOMETRIC_OBJECT make_block(MATERIAL_TYPE material, vector3 center,
 GEOMETRIC_OBJECT make_ellipsoid(MATERIAL_TYPE material, vector3 center,
 				vector3 e1, vector3 e2, vector3 e3,
 				vector3 size);
+GEOMETRIC_OBJECT make_prism(MATERIAL_TYPE material,
+			    const vector3 *vertices, int num_vertices,
+			    double height, vector3 axis);
 
+int vector3_nearly_equal(vector3 v1, vector3 v2, double tolerance);
 
 /**************************************************************************/