2 logScale,printTimeStep,evalSol,order,numResolutions,solnType,tFinal,...
3 Nx,Nt_sol,M,maxNumNeighbors,src_dmn,ds,...
520bc24updated readme file and cleaned up codenalhassanieh 4 uniform_sgrid,betaMax,dataParam,dt0_min,uniform_beta,keepdt]...
10 case 1
11 % testing buttons
12 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
13 saveWorkspaceOpt = 0;
14 addErrResults = 0;
15 savePlot = 0;
16 logScale = 0; % '1' to use log scale in heat maps
17 printTimeStep = 0; % '1' to print out each time step
18 evalSol = 1;
20 % parameters
520bc24updated readme file and cleaned up codenalhassanieh 21 order = 3; % order of accuracy = interpolation
23 solnType = 'ms'; % type of solution 'ms' or 'true'
25 tFinal = 3*pi; % final time
27 %%% Choose number of spatial and temporal set for solution computation
520bc24updated readme file and cleaned up codenalhassanieh 28 Nx = 5; % size of spatial grid for testing and plotting
29 Nt_sol = 5; % size of the time grid for testing and plotting
ae869c2updated codenalhassanieh 30
31 %%% sources
ae869c2updated codenalhassanieh 37
38 %%% spring constants
ae869c2updated codenalhassanieh 40
41 %%% incident pulse parameter
ae869c2updated codenalhassanieh 43
ae869c2updated codenalhassanieh 45
48 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
49 saveWorkspaceOpt = 1;
50 addErrResults = 0;
51 savePlot = 0;
52 logScale = 0; % '1' to use log scale in heat maps
53 printTimeStep = 1; % '1' to print out each time step
54 evalSol = 1;
56 % parameters
57 order = 8; % order of accuracy = interpolation
58 numResolutions = 2; % num of grid resolutions
59 solnType = 'true'; % type of solution 'ms' or 'true'
61 tFinal = 10*pi; % final time
63 %%% Choose number of spatial and temporal set for solution computation
64 Nx = 400; % size of spatial grid for testing and plotting
65 Nt_sol = 400; % size of the time grid for testing and plotting
67 %%% sources
68 M = 150; % number of sources
69 maxNumNeighbors = 10; % set max number of neighbors
70 src_dmn = [-2,2];
71 ds = 1e-4; % set min distance between sources
72 uniform_sgrid = 0;
74 %%% spring constants
75 betaMax = 10; uniform_beta = 0;
77 %%% incident pulse parameter
81 dt0_min = 0.1;
ae869c2updated codenalhassanieh 84
85 % testing buttons
86 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
87 saveWorkspaceOpt = 1;
88 addErrResults = 0;
89 savePlot = 0;
90 logScale = 0; % '1' to use log scale in heat maps
91 printTimeStep = 1; % '1' to print out each time step
92 evalSol = 1;
94 % parameters
95 order = 8; % order of accuracy = interpolation
96 numResolutions = 2; % num of grid resolutions
97 solnType = 'true'; % type of solution 'ms' or 'true'
99 tFinal = 3*pi; % final time
101 %%% Choose number of spatial and temporal set for solution computation
102 Nx = 800; % size of spatial grid for testing and plotting
103 Nt_sol = 800; % size of the time grid for testing and plotting
105 %%% sources
106 M = 1000; % number of sources
107 maxNumNeighbors = 50; % set max number of neighbors
108 src_dmn = [-2,2];
109 ds = 1e-4; % set min distance between sources
110 uniform_sgrid = 0;
112 %%% spring constants
113 betaMax = 3; uniform_beta = 0;
115 %%% incident pulse parameter
520bc24updated readme file and cleaned up codenalhassanieh 116 dataParam.mu = M^2; dataParam.t0 = -3;
119 dt0_min = 0.1;
122 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
123 saveWorkspaceOpt = 1;
124 addErrResults = 0;
125 savePlot = 0;
126 logScale = 0; % '1' to use log scale in heat maps
520bc24updated readme file and cleaned up codenalhassanieh 127 printTimeStep = 1; % '1' to print out each time step
130 % parameters
131 order = 8; % order of accuracy = interpolation
132 numResolutions = 2; % num of grid resolutions
133 solnType = 'true'; % type of solution 'ms' or 'true'
135 tFinal = 3*pi; % final time
137 %%% Choose number of spatial and temporal set for solution computation
520bc24updated readme file and cleaned up codenalhassanieh 138 Nx = 400; % size of spatial grid for testing and plotting
141 %%% sources
142 M = 10000; % number of sources
520bc24updated readme file and cleaned up codenalhassanieh 143 maxNumNeighbors = 500; % set max number of neighbors
145 ds = 1e-4; % set min distance between sources
146 uniform_sgrid = 0;
148 %%% spring constants
149 betaMax = 3; uniform_beta = 0;
151 %%% incident pulse parameter
520bc24updated readme file and cleaned up codenalhassanieh 152 dataParam.mu = M^2; dataParam.t0 = -3;
153 dataParam.doubleTimeStep = 1;
156 dt0_min = 0.1;
157 case 5 % FP
158 % testing buttons
159 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
160 saveWorkspaceOpt = 1;
161 addErrResults = 0;
162 savePlot = 0;
163 logScale = 0; % '1' to use log scale in heat maps
164 printTimeStep = 0; % '1' to print out each time step
165 evalSol = 1;
167 % parameters
168 order = 8; % order of accuracy = interpolation
169 numResolutions = 2; % num of grid resolutions
170 solnType = 'true'; % type of solution 'ms' or 'true'
172 tFinal = 30*pi; % final time
174 %%% Choose number of spatial and temporal set for solution computation
175 Nx = 400; % size of spatial grid for testing and plotting
176 Nt_sol = 500; % size of the time grid for testing and plotting
178 %%% sources
179 M = 2; % number of sources
180 maxNumNeighbors = 2; % set max number of neighbors
181 src_dmn = [-0.5,0.5];
182 ds = 1e-4; % set min distance between sources
183 uniform_sgrid = 1;
185 %%% spring constants
186 betaMax = 100; uniform_beta = 1;
188 %%% incident pulse parameter
192 dt0_min = 0.02; % RBC constraint
193 dt0_min = min(0.01,dt0_min); % stability constraint
195 case 6 % FP M = 10
196 % testing buttons
197 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
200 savePlot = 0;
201 logScale = 0; % '1' to use log scale in heat maps
202 printTimeStep = 1; % '1' to print out each time step
203 evalSol = 1;
205 % parameters
206 order = 8; % order of accuracy = interpolation
520bc24updated readme file and cleaned up codenalhassanieh 207 numResolutions = 2; % num of grid resolutions
ae869c2updated codenalhassanieh 211
212 %%% Choose number of spatial and temporal set for solution computation
520bc24updated readme file and cleaned up codenalhassanieh 213 Nx = 800; % size of spatial grid for testing and plotting
214 Nt_sol = 800; % size of the time grid for testing and plotting
ae869c2updated codenalhassanieh 215
216 %%% sources
219 src_dmn = [-2,2];
220 ds = 1e-4; % set min distance between sources
221 uniform_sgrid = 1;
223 %%% spring constants
224 betaMax = 100; uniform_beta = 1;
226 %%% incident pulse parameter
230 dt0_min = 0.02; % RBC constraint
520bc24updated readme file and cleaned up codenalhassanieh 231 dt0_min = min(0.005,dt0_min); % stability constraint
233 % testing buttons
234 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
235 saveWorkspaceOpt = 1;
236 addErrResults = 0;
237 savePlot = 0;
238 logScale = 0; % '1' to use log scale in heat maps
239 printTimeStep = 0; % '1' to print out each time step
240 evalSol = 1;
242 % parameters
243 order = 8; % order of accuracy = interpolation
244 numResolutions = 2; % num of grid resolutions
245 solnType = 'true'; % type of solution 'ms' or 'true'
247 tFinal = 40*pi; % final time
249 %%% Choose number of spatial and temporal set for solution computation
250 Nx = 800; % size of spatial grid for testing and plotting
251 Nt_sol = 800; % size of the time grid for testing and plotting
253 %%% sources
254 M = 200; % number of sources
255 maxNumNeighbors = 20; % set max number of neighbors
256 src_dmn = [-2,2];
257 ds = 1e-4; % set min distance between sources
258 uniform_sgrid = 1;
260 %%% spring constants
261 betaMax = 1; uniform_beta = 1;
263 %%% incident pulse parameter
267 dt0_min = 0.02; % RBC constraint
520bc24updated readme file and cleaned up codenalhassanieh 268 dt0_min = min(0.005,dt0_min); % stability constraint
269 case 999
270 % testing buttons
271 plotOpt = 0; % '1' to plot sol at final time, '2' heat maps
272 saveWorkspaceOpt = 0;
273 addErrResults = 0;
274 savePlot = 0;
275 logScale = 0; % '1' to use log scale in heat maps
276 printTimeStep = 0; % '1' to print out each time step
277 evalSol = 1;
279 % parameters
280 order = 2; % order of accuracy = interpolation
281 numResolutions = 3; % num of grid resolutions
282 solnType = 'true'; % type of solution 'ms' or 'true'
284 tFinal = 3*pi; % final time
286 %%% Choose number of spatial and temporal set for solution computation
287 Nx = 5; % size of spatial grid for testing and plotting
288 Nt_sol = 5; % size of the time grid for testing and plotting
290 %%% sources
291 M = 10; % number of sources
292 maxNumNeighbors = NaN; % set max number of neighbors
293 src_dmn = [-2,2];
294 ds = 1e-4; % set min distance between sources
295 uniform_sgrid = 1;
297 %%% spring constants
298 betaMax = 3; uniform_beta = 1;
300 %%% incident pulse parameter
301 dataParam.mu = 30; dataParam.t0 = -3;
302 keepdt = 0;
304 dt0_min = 0.01; % RBC constraint
306 disp('experiment number inapplicable');
307end
309end