1 #ifndef INTERPOLATION_HH
2 # define INTERPOLATION_HH
4 # include "definitions.hh"
5 # include "topology.hh"
20 return point2d((p.row() + 1) * 4, (p.col() + 1) * 4);
28 return point2d(p.row() / 4 - 1, p.col() / 4 - 1);
37 immerse_primary_2_faces(F_type& F,
47 if (f.domain().has(p_))
57 immerse_separating_2_faces(F_type& F, A a)
59 bool c25_vals[] = { 1, 1, 1, 1, 1,
64 neighb2d c25 = make::neighb2d(c25_vals);
68 if (is_2_face(p) && ! is_primary(p))
73 if (D.has(n) && is_primary(n))
82 immerse_0_and_1_faces(F_type& F, A a)
93 if (D.has(n) && is_2_face(n))
134 immerse(const image2d<value::int_u8>& f)
145 range<unsignedh> b = border_median(f);
146 immerse_primary_2_faces(F, f, b);
149 immerse_separating_2_faces(F, accu_max_range());
152 immerse_0_and_1_faces(F, accu_span_range());
180 #endif // ndef INTERPOLATION_HH