1function get_msOrIncidentField(M,tFinal,mu,solnType)
4if(strcmp(solnType,'ms'))
5 override = 0; % turn on '1' to avoid restrictions on chosen test function
7mu = linspace(40,50,M)';
8t0 = linspace(1,3,M)';
10% choose the initial time step
11dt0_vec = pi./(4*sqrt(mu*log(1/eps)));
12dt0 = min(dt0_vec);
14% ensure that the forcing function is bandlimited within the computational time domain
15if(any(2*t0>tFinal)&&(override == 0))
16 error('Forcing inappropriate for testing. For tFinal=%1.2f\nChoose other mu or tFinal',tFinal);
17end
19% Define the density functions to be Gaussians
20sig = @(t) exp(-mu.*((t - t0).^2));
22dataParam.sig = sig;
23dataParam.t0 = t0;
24dataParam.mu = mu;
26else
28t0 = -3;
30% choose the initial time step
31dt0 = pi/(4*sqrt(mu*log(1/eps)));
33dataParam.mu = mu;
34dataParam.t0 = t0;
35end
40end