1function plan = plate_rr_plan(params)
2% PLATE_RR_PLAN Sizes and estimates of a Rayleigh-Ritz run, before it runs.
3%
4% plan = plate_rr_plan(params) reports the problem size (Legendre
5% coefficients and quadrature points), the number of modes to expect below
6% 0.47 fs from the modal density of a thin plate, n(f) = (A/2) sqrt(rho h / D),
7% the pendulum frequency of the thread, the Hertz estimates of the contact
8% and the plate's coincidence frequency. Cheap: no matrix is assembled.
10P = plate_setup(params);
11plan = plate_plan_common(params, P);
12n = (params.rr.basis_order + 1)^2;
13plan.method = 'rr';
14plan.ndof = n;
15plan.nquad = params.rr.quadrature_order^2;
16plan.work = n^3;
17plan.size_note = sprintf('%d Legendre coefficients (order %d), %d quadrature points', n, params.rr.basis_order, plan.nquad);
18if params.rr.basis_order < 4
19 plan.notes{end + 1} = 'a basis order below 4 resolves almost no bending mode';
20end
21end