* apps/mesh-segm-skel/mesh-complex-max-curv-segm.cc: Here. --- milena/ChangeLog | 6 ++++++ .../mesh-segm-skel/mesh-complex-max-curv-segm.cc | 11 +++++++---- 2 files changed, 13 insertions(+), 4 deletions(-) diff --git a/milena/ChangeLog b/milena/ChangeLog index 00afa02..923f16b 100644 --- a/milena/ChangeLog +++ b/milena/ChangeLog @@ -1,5 +1,11 @@ 2010-06-14 Roland Levillain &lt;roland(a)lrde.epita.fr&gt; + Actually use the maximal curvature in mesh-complex-max-curv-segm. + + * apps/mesh-segm-skel/mesh-complex-max-curv-segm.cc: Here. + +2010-06-14 Roland Levillain &lt;roland(a)lrde.epita.fr&gt; + Improve the precision of some mesh-segm-skel apps. * apps/mesh-segm-skel/mesh-complex-segm.cc, diff --git a/milena/apps/mesh-segm-skel/mesh-complex-max-curv-segm.cc b/milena/apps/mesh-segm-skel/mesh-complex-max-curv-segm.cc index 1f43381..a7a741a 100644 --- a/milena/apps/mesh-segm-skel/mesh-complex-max-curv-segm.cc +++ b/milena/apps/mesh-segm-skel/mesh-complex-max-curv-segm.cc @@ -91,14 +91,17 @@ int main(int argc, char* argv[]) mln::p_n_faces_fwd_piter<D, G> v(float_ima.domain(), 0); for_all(v) { + // FIXME: The program should allow the user to choose the kind + // of measure (curvature). +#if 0 float h = (curv.first(v) + curv.second(v)) / 2; // Pseudo-inverse curvature. float h_inv = 1 / pi * (atan(-h) + pi / 2); float_ima(v) = h_inv; - // FIXME: The program should allow the user to choose the kind - // of measure. -// input(v) = mln::math::max(mln::math::sqr(curv.first(v)), -// mln::math::sqr(curv.second(v))); +#endif + // Max curvature. + float_ima(v) = mln::math::max(mln::math::sqr(curv.first(v)), + mln::math::sqr(curv.second(v))); } // Values on edges. -- 1.5.6.5