| 49 | |
| 50 | private: |
| 51 | const vec_vec_BaseDescMod DescMods_( // utility for Mod constructor |
| 52 | const BaseDescMod** const descmods_arg, // in: descriptor models |
| 53 | int ndescmods) // in: sanity check |
| 54 | { |
| 55 | if (stasm_NLANDMARKS != ndescmods / N_PYR_LEVS) // sanity check |
| 56 | { |
| 57 | Err("stasm_NLANDMARKS != ndescmods / N_PYR_LEVS\n\n" |
| 58 | " stasm_NLANDMARKS %d\n ndescmods / N_PYR_LEVS %d\n" |
| 59 | " ndescmods %d\n N_PYR_LEVS %d", |
| 60 | stasm_NLANDMARKS, ndescmods / N_PYR_LEVS, |
| 61 | ndescmods, N_PYR_LEVS); |
| 62 | } |
| 63 | vec_vec_BaseDescMod descmods(N_PYR_LEVS); |
| 64 | |
| 65 | for (int ilev = 0; ilev < N_PYR_LEVS; ilev++) |
| 66 | { |
| 67 | descmods[ilev].resize(stasm_NLANDMARKS); |
| 68 | for (int i = 0; i < stasm_NLANDMARKS; i++) |
| 69 | descmods[ilev][i] = |
| 70 | descmods_arg[ilev * stasm_NLANDMARKS + i]; |
| 71 | } |
| 72 | return descmods; |
| 73 | } |
| 74 | |
| 75 | public: |
| 76 | Mod(EYAW eyaw, // constructor |