Eigenproblem Parameters: Difference between revisions
From WikiROMS
Jump to navigationJump to search
New page: <div class="title">Eigenproblem Parameters</div> {{#lst:Variables|NEV}} {{#lst:Variables|NCV}} '''Notice:''' At present, there is no a-priori analysis to guide the selection of <span cl... (change visibility) |
(No difference)
|
Latest revision as of 20:02, 16 April 2008
Eigenproblem Parameters
- NEV
- Number of eigenvalues to compute for the Lanczos/Arnoldi problem. Notice that the model memory requirement increases substantially as NEV increases. The GST requires NEV+1 copies of the model state vector. The memory requirements are decreased in distributed-memory applications.
- option =
- routine = mod_storage.F
- keyword = NEV
- input = roms.in
- NCV
- Number of eigenvectors to compute for the Lanczos/Arnoldi problem. NCV must be greater than NEV.
- option =
- routine = mod_storage.F
- keyword = NCV
- input = roms.in
Notice: At present, there is no a-priori analysis to guide the selection of NCV relative to NEV. The only formal requirement is that NCV > NEV. However in optimal perturbations, it is recommended to have NCV ≥ 2*NEV. In Finite Time Eigenmodes (FTE) and Adjoint Finite Time Eigenmodes (AFTE) the requirement is to have NCV ≥ 2*NEV+1. The efficiency of calculations depends critically on the combination of NEV and NCV. If NEV is large (greater than 10 say), you can use NCV=2*NEV+1 but for NEV small (less than 6) it will be inefficient to use NCV=2*NEV+1. In complicated applications, you can start with NEV=2 and NCV=10. Otherwise, it will iterate for very long time.
