Wave Propagation in Linear Media | re-examined

PhiArg[a, b, c,d],1,opts] +
OscInt[RectTransNeg[a,-b,-c,d,w,k],
PhiArg[a,-b,-c,d],1,opts])
) RectConst[w,k]
];
GaussTransTemp[fpos_, fneg_, x_, t_, w_, k_, n_, opts___Rule] :=
Module[{a,b,c,d,o1,o2,
wp = WorkingPrecision/.{opts}/.Options[OscInt]},
a = -t;
b = N[Pi k/Sqrt[w],wp+2];
c = x;
d = N[-Pi k,wp+2];
o1 = QuadOffset[a, b+c,d];
o2 = QuadOffset[a,-b-c,d];
(OscInt[fpos[a, b, c,d,w,k,n],
QuadZero[a, b +c,d,o1],1,FunctionType->ZeroList,opts] +
OscInt[fneg[a,-b,-c,d,w,k,n],
QuadZero[a,-b -c,d,o2],1,FunctionType->ZeroList,opts]
) GaussConst[w,k,n]
];
GaussTransTruncTemp[fpos_, fneg_, x_, t_, w_, k_, n_, opts___Rule] :=
Module[{wp = WorkingPrecision/.{opts}/.Options[OscInt],
ag = AccuracyGoal/. {opts}/.Options[OscInt],
pg = PrecisionGoal/. {opts}/.Options[OscInt],
mi = MinRecursion/. {opts}/.Options[OscInt],
ma = MaxRecursion/. {opts}/.Options[OscInt],
limpos, limneg},
{limneg,limpos} = GaussTruncLims[w,k,n,wp];
(If[limpos DoubleExponential,
WorkingPrecision->wp,
AccuracyGoal->ag,
PrecisionGoal->pg,
MinRecursion->mi,
MaxRecursion->ma]] +
If[limneg >= -1, 0,
NIntegrate[fneg[-t,-Pi k/Sqrt[w],-x,-Pi k,w,k,n][xi],
{xi,1,-limneg},
Method->DoubleExponential,
WorkingPrecision->wp,
AccuracyGoal->ag,
PrecisionGoal->pg,
MinRecursion->mi,
MaxRecursion->ma]]
) GaussConst[w,k,n]
];
GaussEvanTruncTemp[f_, w_, k_, n_, opts___Rule] :=
Module[{limpos,limneg,
wp = WorkingPrecision/.{opts}/.Options[OscInt],
ag = AccuracyGoal/. {opts}/.Options[OscInt],
pg = PrecisionGoal/. {opts}/.Options[OscInt],
mi = MinRecursion/. {opts}/.Options[OscInt],
ma = MaxRecursion/. {opts}/.Options[OscInt]},
{limneg,limpos} = GaussTruncLims[w,k,n,wp];
limneg = Max[limneg,-1];
limpos = Min[limpos,1];
NIntegrate[f[xi],{xi,limneg,limpos},
Method->DoubleExponential,
WorkingPrecision->wp,
AccuracyGoal->ag,
PrecisionGoal->pg,
MinRecursion->mi,
MaxRecursion->ma] GaussConst[w,k,n]
];
(*-- Driver functions for the computation of the wave integrals --*)
PhiTrans[x_, t_, w_, k_, n_?Positive, opts___Rule] :=
Module[{sh = Shape/. {opts}/.Options[Phi],
pt = PPoints/.{opts}/.Options[Phi]},
Switch[sh,
Rect,
Switch[pt,
Approximate, RectTransTemp[RectTransPos,RectTransNeg,x,t,w,k,opts],
Zeros, RectTransExact[x,t,w,k,opts],
_, Message[Phi::invalidpart,pt]],
Tria,
Switch[pt,
226
A Mathematica packages