| 13 | // "b" is the name used for the eigenvector weights in Cootes' papers. |
| 14 | |
| 15 | static void LimitB( |
| 16 | VEC& b, // io: eigvec weights |
| 17 | const VEC& eigvals, // in |
| 18 | double bmax) // in |
| 19 | { |
| 20 | for (int i = 0; i < NSIZE(eigvals); i++) |
| 21 | { |
| 22 | const double limit = bmax * sqrt(eigvals(i)); |
| 23 | b(i) = Clamp(b(i), -limit, limit); |
| 24 | } |
| 25 | } |
| 26 | |
| 27 | // This implements Section 4.8 of CootesTaylor 2004 |
| 28 | // www.isbe.man.ac.uk/~bim/Mods/app_models.pdf. |
no test coverage detected