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