diff --git a/src/meepgeom.cpp b/src/meepgeom.cpp index 2b5afdad8..1e90e728a 100644 --- a/src/meepgeom.cpp +++ b/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; diff --git a/src/susceptibility.cpp b/src/susceptibility.cpp index 7d927e590..28d54ec02 100644 --- a/src/susceptibility.cpp +++ b/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,