in_subdomain.m

利用电磁场的源激发方法来计算光子晶体波导例如光子晶体光纤

function out = in_subdomain(mX, mY, sd)
if isempty(sd.insideCurve)
    out = ones(size(mX));
else
    out = in_shape(mX, mY, sd.insideCurve.xc, sd.insideCurve.yc, sd.insideCurve);
end
for iCurve = 1:length(sd.outsideCurveArray)
    curve = sd.outsideCurveArray(iCurve);
    out = out & ~in_shape(mX, mY, curve.xc, curve.yc, curve);
end