For Bug#1117

src/meepgeom.cpp

@@ -1190,6 +1190,7 @@ static bool susceptibility_equiv(const susceptibility *o0, const susceptibility
   if (!vector3_equal(o0->bias, o->bias)) return 0;
   if (o0->frequency != o->frequency) return 0;
   if (o0->gamma != o->gamma) return 0;
+  if (o0->alpha != o->alpha) return 0;
   if (o0->noise_amp != o->noise_amp) return 0;
   if (o0->drude != o->drude) return 0;
   if (o0->saturated_gyrotropy != o->saturated_gyrotropy) return 0;

src/susceptibility.cpp

@@ -603,8 +603,9 @@ void gyrotropic_susceptibility::subtract_P(field_type ft,
 
 int gyrotropic_susceptibility::num_cinternal_notowned_needed(component c,
                                                              void *P_internal_data) const {
-  gyrotropy_data *d = (gyrotropy_data *)P_internal_data;
-  return d->P[c][0][0] ? 3 : 0;
+    (void)c;
+    (void)P_internal_data;
+    return 0;
 }
 
 realnum *gyrotropic_susceptibility::cinternal_notowned_ptr(int inotowned, component c, int cmp,