HUMIM Tracker: /home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/core/containers/mesh/common/util/mesh

Go to the documentation of this file.
00001 #ifndef OPENTISSUE_CORE_CONTAINERS_MESH_COMMON_UTIL_MESH_ISOSURFACE_H
00002 #define OPENTISSUE_CORE_CONTAINERS_MESH_COMMON_UTIL_MESH_ISOSURFACE_H
00003 //
00004 // OpenTissue Template Library
00005 // - A generic toolbox for physics-based modeling and simulation.
00006 // Copyright (C) 2008 Department of Computer Science, University of Copenhagen.
00007 //
00008 // OTTL is licensed under zlib: http://opensource.org/licenses/zlib-license.php
00009 //
00010 #include <OpenTissue/configuration.h>
00011 
00012 #include <OpenTissue/core/containers/mesh/common/util/mesh_compute_mesh_center.h>
00013 #include <OpenTissue/core/containers/mesh/common/util/mesh_compute_mesh_minimum_coord.h>
00014 #include <OpenTissue/core/containers/mesh/common/util/mesh_compute_mesh_maximum_coord.h>
00015 #include <OpenTissue/core/containers/mesh/common/util/mesh_deformation_modifiers.h>
00016 
00017 #include <cassert>
00018 #include <cmath>
00019 
00020 namespace OpenTissue
00021 {
00022   namespace mesh
00023   {
00024 
00025 
00026     namespace detail
00027     {
00028 
00039       class IsoSurfaceGenerator
00040       {
00041       protected:
00042 
00043         typedef int identifier_type;
00044 
00045       protected:
00046 
00057         identifier_type compute_vertex_id(size_t i, size_t j, size_t k,size_t I, size_t J)
00058         {
00059           return 3*(k*J*I + j*I + i);
00060         }
00061 
00073         identifier_type get_edge_id(size_t i, size_t j, size_t k, identifier_type edge_number,size_t I,size_t J)
00074         {
00075           //switch (edge_number)
00076           //{
00077           //  case 0:        return compute_vertex_id(    i,     j,     k, I, J) + 1;
00078           //  case 1:        return compute_vertex_id(    i, j + 1,     k, I, J);
00079           //  case 2:        return compute_vertex_id(i + 1,     j,     k, I, J) + 1;
00080           //  case 3:        return compute_vertex_id(    i,     j,     k, I, J);
00081           //  case 4:        return compute_vertex_id(    i,     j, k + 1, I, J) + 1;
00082           //  case 5:        return compute_vertex_id(    i, j + 1, k + 1, I, J);
00083           //  case 6:        return compute_vertex_id(i + 1,     j, k + 1, I, J) + 1;
00084           //  case 7:        return compute_vertex_id(    i,     j, k + 1, I, J);
00085           //  case 8:        return compute_vertex_id(    i,     j,     k, I, J) + 2;
00086           //  case 9:        return compute_vertex_id(    i, j + 1,     k, I, J) + 2;
00087           //  case 10:       return compute_vertex_id(i + 1, j + 1,     k, I, J) + 2;
00088           //  case 11:       return compute_vertex_id(i + 1,     j,     k, I, J) + 2;
00089           //  default:       return -1;// Invalid edge no.
00090           //}
00091           switch (edge_number)
00092           {
00093           case 0:        return 3*(    k*J*I +     j*I + i  ) + 1;
00094           case 1:        return 3*(    k*J*I + (j+1)*I + i  );
00095           case 2:        return 3*(    k*J*I +     j*I + i+1) + 1;
00096           case 3:        return 3*(    k*J*I +     j*I + i  );
00097           case 4:        return 3*((k+1)*J*I +     j*I + i  ) + 1;
00098           case 5:        return 3*((k+1)*J*I + (j+1)*I + i  );
00099           case 6:        return 3*((k+1)*J*I +     j*I + i+1) + 1;
00100           case 7:        return 3*((k+1)*J*I +     j*I + i  );
00101           case 8:        return 3*(    k*J*I +     j*I + i  ) + 2;
00102           case 9:        return 3*(    k*J*I + (j+1)*I + i  ) + 2;
00103           case 10:       return 3*(    k*J*I + (j+1)*I + i+1) + 2;
00104           case 11:       return 3*(    k*J*I +     j*I + i+1) + 2;
00105           default:       return -1;// Invalid edge no.
00106           }
00107         }
00108 
00119         template<
00120           typename grid_type
00121           , typename size_t
00122         >
00123         typename grid_type::math_types::vector3_type
00124         calculate_intersection(
00125         grid_type const & phi
00126         , typename grid_type::value_type const & isovalue
00127         , size_t  i
00128         , size_t  j
00129         , size_t  k
00130         , identifier_type edge_number
00131         )
00132         {
00133           typedef typename grid_type::value_type       value_type;
00134           typedef typename grid_type::math_types       math_types;
00135           typedef typename math_types::vector3_type    vector3_type;
00136           typedef typename math_types::real_type       real_type;
00137 
00138           size_t i1 = i, j1 = j, k1 = k;
00139           size_t i2 = i, j2 = j, k2 = k;
00140           switch (edge_number)
00141           {
00142           case 0:        j2 += 1;        break;
00143           case 1:        j1 += 1;        i2 += 1;        j2 += 1;        break;
00144           case 2:        i1 += 1;        j1 += 1;        i2 += 1;        break;
00145           case 3:        i1 += 1;        break;
00146           case 4:        k1 += 1;        j2 += 1;        k2 += 1;        break;
00147           case 5:        j1 += 1;        k1 += 1;        i2 += 1;        j2 += 1;        k2 += 1;        break;
00148           case 6:        i1 += 1;        j1 += 1;        k1 += 1;        i2 += 1;        k2 += 1;        break;
00149           case 7:        i1 += 1;        k1 += 1;        k2 += 1;        break;
00150           case 8:        k2 += 1;        break;
00151           case 9:        j1 += 1;        j2 += 1;        k2 += 1;        break;
00152           case 10:       i1 += 1;        j1 += 1;        i2 += 1;        j2 += 1;        k2 += 1;        break;
00153           case 11:       i1 += 1;        i2 += 1;        k2 += 1;        break;
00154           }
00155           vector3_type p1 = vector3_type( i1*phi.dx(), j1*phi.dy(), k1*phi.dz()) + phi.min_coord();
00156           vector3_type p2 = vector3_type( i2*phi.dx(), j2*phi.dy(), k2*phi.dz()) + phi.min_coord();
00157           value_type  value1 = phi(i1,j1,k1);
00158           value_type  value2 = phi(i2,j2,k2);
00159           real_type mu = static_cast<real_type>(((isovalue - value1))/(value2 - value1));
00160           vector3_type p = p1 + mu*(p2-p1);
00161           return p;
00162         }
00163 
00164       protected:
00165 
00172         static identifier_type edge_table(identifier_type const & idx)
00173         {
00174           static identifier_type table[256] = {
00175             0x0  , 0x109, 0x203, 0x30a, 0x406, 0x50f, 0x605, 0x70c,
00176             0x80c, 0x905, 0xa0f, 0xb06, 0xc0a, 0xd03, 0xe09, 0xf00,
00177             0x190, 0x99 , 0x393, 0x29a, 0x596, 0x49f, 0x795, 0x69c,
00178             0x99c, 0x895, 0xb9f, 0xa96, 0xd9a, 0xc93, 0xf99, 0xe90,
00179             0x230, 0x339, 0x33 , 0x13a, 0x636, 0x73f, 0x435, 0x53c,
00180             0xa3c, 0xb35, 0x83f, 0x936, 0xe3a, 0xf33, 0xc39, 0xd30,
00181             0x3a0, 0x2a9, 0x1a3, 0xaa , 0x7a6, 0x6af, 0x5a5, 0x4ac,
00182             0xbac, 0xaa5, 0x9af, 0x8a6, 0xfaa, 0xea3, 0xda9, 0xca0,
00183             0x460, 0x569, 0x663, 0x76a, 0x66 , 0x16f, 0x265, 0x36c,
00184             0xc6c, 0xd65, 0xe6f, 0xf66, 0x86a, 0x963, 0xa69, 0xb60,
00185             0x5f0, 0x4f9, 0x7f3, 0x6fa, 0x1f6, 0xff , 0x3f5, 0x2fc,
00186             0xdfc, 0xcf5, 0xfff, 0xef6, 0x9fa, 0x8f3, 0xbf9, 0xaf0,
00187             0x650, 0x759, 0x453, 0x55a, 0x256, 0x35f, 0x55 , 0x15c,
00188             0xe5c, 0xf55, 0xc5f, 0xd56, 0xa5a, 0xb53, 0x859, 0x950,
00189             0x7c0, 0x6c9, 0x5c3, 0x4ca, 0x3c6, 0x2cf, 0x1c5, 0xcc ,
00190             0xfcc, 0xec5, 0xdcf, 0xcc6, 0xbca, 0xac3, 0x9c9, 0x8c0,
00191             0x8c0, 0x9c9, 0xac3, 0xbca, 0xcc6, 0xdcf, 0xec5, 0xfcc,
00192             0xcc , 0x1c5, 0x2cf, 0x3c6, 0x4ca, 0x5c3, 0x6c9, 0x7c0,
00193             0x950, 0x859, 0xb53, 0xa5a, 0xd56, 0xc5f, 0xf55, 0xe5c,
00194             0x15c, 0x55 , 0x35f, 0x256, 0x55a, 0x453, 0x759, 0x650,
00195             0xaf0, 0xbf9, 0x8f3, 0x9fa, 0xef6, 0xfff, 0xcf5, 0xdfc,
00196             0x2fc, 0x3f5, 0xff , 0x1f6, 0x6fa, 0x7f3, 0x4f9, 0x5f0,
00197             0xb60, 0xa69, 0x963, 0x86a, 0xf66, 0xe6f, 0xd65, 0xc6c,
00198             0x36c, 0x265, 0x16f, 0x66 , 0x76a, 0x663, 0x569, 0x460,
00199             0xca0, 0xda9, 0xea3, 0xfaa, 0x8a6, 0x9af, 0xaa5, 0xbac,
00200             0x4ac, 0x5a5, 0x6af, 0x7a6, 0xaa , 0x1a3, 0x2a9, 0x3a0,
00201             0xd30, 0xc39, 0xf33, 0xe3a, 0x936, 0x83f, 0xb35, 0xa3c,
00202             0x53c, 0x435, 0x73f, 0x636, 0x13a, 0x33 , 0x339, 0x230,
00203             0xe90, 0xf99, 0xc93, 0xd9a, 0xa96, 0xb9f, 0x895, 0x99c,
00204             0x69c, 0x795, 0x49f, 0x596, 0x29a, 0x393, 0x99 , 0x190,
00205             0xf00, 0xe09, 0xd03, 0xc0a, 0xb06, 0xa0f, 0x905, 0x80c,
00206             0x70c, 0x605, 0x50f, 0x406, 0x30a, 0x203, 0x109, 0x0
00207           };
00208           return table[idx];
00209         }
00210 
00211 
00218         static identifier_type triangle_table(identifier_type const & idx1,identifier_type const & idx2)
00219         {
00220           static identifier_type table[256][16] = {
00221             {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00222             {0, 8, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00223             {0, 1, 9, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00224             {1, 8, 3, 9, 8, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00225             {1, 2, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00226             {0, 8, 3, 1, 2, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00227             {9, 2, 10, 0, 2, 9, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00228             {2, 8, 3, 2, 10, 8, 10, 9, 8, -1, -1, -1, -1, -1, -1, -1},
00229             {3, 11, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00230             {0, 11, 2, 8, 11, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00231             {1, 9, 0, 2, 3, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00232             {1, 11, 2, 1, 9, 11, 9, 8, 11, -1, -1, -1, -1, -1, -1, -1},
00233             {3, 10, 1, 11, 10, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00234             {0, 10, 1, 0, 8, 10, 8, 11, 10, -1, -1, -1, -1, -1, -1, -1},
00235             {3, 9, 0, 3, 11, 9, 11, 10, 9, -1, -1, -1, -1, -1, -1, -1},
00236             {9, 8, 10, 10, 8, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00237             {4, 7, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00238             {4, 3, 0, 7, 3, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00239             {0, 1, 9, 8, 4, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00240             {4, 1, 9, 4, 7, 1, 7, 3, 1, -1, -1, -1, -1, -1, -1, -1},
00241             {1, 2, 10, 8, 4, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00242             {3, 4, 7, 3, 0, 4, 1, 2, 10, -1, -1, -1, -1, -1, -1, -1},
00243             {9, 2, 10, 9, 0, 2, 8, 4, 7, -1, -1, -1, -1, -1, -1, -1},
00244             {2, 10, 9, 2, 9, 7, 2, 7, 3, 7, 9, 4, -1, -1, -1, -1},
00245             {8, 4, 7, 3, 11, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00246             {11, 4, 7, 11, 2, 4, 2, 0, 4, -1, -1, -1, -1, -1, -1, -1},
00247             {9, 0, 1, 8, 4, 7, 2, 3, 11, -1, -1, -1, -1, -1, -1, -1},
00248             {4, 7, 11, 9, 4, 11, 9, 11, 2, 9, 2, 1, -1, -1, -1, -1},
00249             {3, 10, 1, 3, 11, 10, 7, 8, 4, -1, -1, -1, -1, -1, -1, -1},
00250             {1, 11, 10, 1, 4, 11, 1, 0, 4, 7, 11, 4, -1, -1, -1, -1},
00251             {4, 7, 8, 9, 0, 11, 9, 11, 10, 11, 0, 3, -1, -1, -1, -1},
00252             {4, 7, 11, 4, 11, 9, 9, 11, 10, -1, -1, -1, -1, -1, -1, -1},
00253             {9, 5, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00254             {9, 5, 4, 0, 8, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00255             {0, 5, 4, 1, 5, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00256             {8, 5, 4, 8, 3, 5, 3, 1, 5, -1, -1, -1, -1, -1, -1, -1},
00257             {1, 2, 10, 9, 5, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00258             {3, 0, 8, 1, 2, 10, 4, 9, 5, -1, -1, -1, -1, -1, -1, -1},
00259             {5, 2, 10, 5, 4, 2, 4, 0, 2, -1, -1, -1, -1, -1, -1, -1},
00260             {2, 10, 5, 3, 2, 5, 3, 5, 4, 3, 4, 8, -1, -1, -1, -1},
00261             {9, 5, 4, 2, 3, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00262             {0, 11, 2, 0, 8, 11, 4, 9, 5, -1, -1, -1, -1, -1, -1, -1},
00263             {0, 5, 4, 0, 1, 5, 2, 3, 11, -1, -1, -1, -1, -1, -1, -1},
00264             {2, 1, 5, 2, 5, 8, 2, 8, 11, 4, 8, 5, -1, -1, -1, -1},
00265             {10, 3, 11, 10, 1, 3, 9, 5, 4, -1, -1, -1, -1, -1, -1, -1},
00266             {4, 9, 5, 0, 8, 1, 8, 10, 1, 8, 11, 10, -1, -1, -1, -1},
00267             {5, 4, 0, 5, 0, 11, 5, 11, 10, 11, 0, 3, -1, -1, -1, -1},
00268             {5, 4, 8, 5, 8, 10, 10, 8, 11, -1, -1, -1, -1, -1, -1, -1},
00269             {9, 7, 8, 5, 7, 9, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00270             {9, 3, 0, 9, 5, 3, 5, 7, 3, -1, -1, -1, -1, -1, -1, -1},
00271             {0, 7, 8, 0, 1, 7, 1, 5, 7, -1, -1, -1, -1, -1, -1, -1},
00272             {1, 5, 3, 3, 5, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00273             {9, 7, 8, 9, 5, 7, 10, 1, 2, -1, -1, -1, -1, -1, -1, -1},
00274             {10, 1, 2, 9, 5, 0, 5, 3, 0, 5, 7, 3, -1, -1, -1, -1},
00275             {8, 0, 2, 8, 2, 5, 8, 5, 7, 10, 5, 2, -1, -1, -1, -1},
00276             {2, 10, 5, 2, 5, 3, 3, 5, 7, -1, -1, -1, -1, -1, -1, -1},
00277             {7, 9, 5, 7, 8, 9, 3, 11, 2, -1, -1, -1, -1, -1, -1, -1},
00278             {9, 5, 7, 9, 7, 2, 9, 2, 0, 2, 7, 11, -1, -1, -1, -1},
00279             {2, 3, 11, 0, 1, 8, 1, 7, 8, 1, 5, 7, -1, -1, -1, -1},
00280             {11, 2, 1, 11, 1, 7, 7, 1, 5, -1, -1, -1, -1, -1, -1, -1},
00281             {9, 5, 8, 8, 5, 7, 10, 1, 3, 10, 3, 11, -1, -1, -1, -1},
00282             {5, 7, 0, 5, 0, 9, 7, 11, 0, 1, 0, 10, 11, 10, 0, -1},
00283             {11, 10, 0, 11, 0, 3, 10, 5, 0, 8, 0, 7, 5, 7, 0, -1},
00284             {11, 10, 5, 7, 11, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00285             {10, 6, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00286             {0, 8, 3, 5, 10, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00287             {9, 0, 1, 5, 10, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00288             {1, 8, 3, 1, 9, 8, 5, 10, 6, -1, -1, -1, -1, -1, -1, -1},
00289             {1, 6, 5, 2, 6, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00290             {1, 6, 5, 1, 2, 6, 3, 0, 8, -1, -1, -1, -1, -1, -1, -1},
00291             {9, 6, 5, 9, 0, 6, 0, 2, 6, -1, -1, -1, -1, -1, -1, -1},
00292             {5, 9, 8, 5, 8, 2, 5, 2, 6, 3, 2, 8, -1, -1, -1, -1},
00293             {2, 3, 11, 10, 6, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00294             {11, 0, 8, 11, 2, 0, 10, 6, 5, -1, -1, -1, -1, -1, -1, -1},
00295             {0, 1, 9, 2, 3, 11, 5, 10, 6, -1, -1, -1, -1, -1, -1, -1},
00296             {5, 10, 6, 1, 9, 2, 9, 11, 2, 9, 8, 11, -1, -1, -1, -1},
00297             {6, 3, 11, 6, 5, 3, 5, 1, 3, -1, -1, -1, -1, -1, -1, -1},
00298             {0, 8, 11, 0, 11, 5, 0, 5, 1, 5, 11, 6, -1, -1, -1, -1},
00299             {3, 11, 6, 0, 3, 6, 0, 6, 5, 0, 5, 9, -1, -1, -1, -1},
00300             {6, 5, 9, 6, 9, 11, 11, 9, 8, -1, -1, -1, -1, -1, -1, -1},
00301             {5, 10, 6, 4, 7, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00302             {4, 3, 0, 4, 7, 3, 6, 5, 10, -1, -1, -1, -1, -1, -1, -1},
00303             {1, 9, 0, 5, 10, 6, 8, 4, 7, -1, -1, -1, -1, -1, -1, -1},
00304             {10, 6, 5, 1, 9, 7, 1, 7, 3, 7, 9, 4, -1, -1, -1, -1},
00305             {6, 1, 2, 6, 5, 1, 4, 7, 8, -1, -1, -1, -1, -1, -1, -1},
00306             {1, 2, 5, 5, 2, 6, 3, 0, 4, 3, 4, 7, -1, -1, -1, -1},
00307             {8, 4, 7, 9, 0, 5, 0, 6, 5, 0, 2, 6, -1, -1, -1, -1},
00308             {7, 3, 9, 7, 9, 4, 3, 2, 9, 5, 9, 6, 2, 6, 9, -1},
00309             {3, 11, 2, 7, 8, 4, 10, 6, 5, -1, -1, -1, -1, -1, -1, -1},
00310             {5, 10, 6, 4, 7, 2, 4, 2, 0, 2, 7, 11, -1, -1, -1, -1},
00311             {0, 1, 9, 4, 7, 8, 2, 3, 11, 5, 10, 6, -1, -1, -1, -1},
00312             {9, 2, 1, 9, 11, 2, 9, 4, 11, 7, 11, 4, 5, 10, 6, -1},
00313             {8, 4, 7, 3, 11, 5, 3, 5, 1, 5, 11, 6, -1, -1, -1, -1},
00314             {5, 1, 11, 5, 11, 6, 1, 0, 11, 7, 11, 4, 0, 4, 11, -1},
00315             {0, 5, 9, 0, 6, 5, 0, 3, 6, 11, 6, 3, 8, 4, 7, -1},
00316             {6, 5, 9, 6, 9, 11, 4, 7, 9, 7, 11, 9, -1, -1, -1, -1},
00317             {10, 4, 9, 6, 4, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00318             {4, 10, 6, 4, 9, 10, 0, 8, 3, -1, -1, -1, -1, -1, -1, -1},
00319             {10, 0, 1, 10, 6, 0, 6, 4, 0, -1, -1, -1, -1, -1, -1, -1},
00320             {8, 3, 1, 8, 1, 6, 8, 6, 4, 6, 1, 10, -1, -1, -1, -1},
00321             {1, 4, 9, 1, 2, 4, 2, 6, 4, -1, -1, -1, -1, -1, -1, -1},
00322             {3, 0, 8, 1, 2, 9, 2, 4, 9, 2, 6, 4, -1, -1, -1, -1},
00323             {0, 2, 4, 4, 2, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00324             {8, 3, 2, 8, 2, 4, 4, 2, 6, -1, -1, -1, -1, -1, -1, -1},
00325             {10, 4, 9, 10, 6, 4, 11, 2, 3, -1, -1, -1, -1, -1, -1, -1},
00326             {0, 8, 2, 2, 8, 11, 4, 9, 10, 4, 10, 6, -1, -1, -1, -1},
00327             {3, 11, 2, 0, 1, 6, 0, 6, 4, 6, 1, 10, -1, -1, -1, -1},
00328             {6, 4, 1, 6, 1, 10, 4, 8, 1, 2, 1, 11, 8, 11, 1, -1},
00329             {9, 6, 4, 9, 3, 6, 9, 1, 3, 11, 6, 3, -1, -1, -1, -1},
00330             {8, 11, 1, 8, 1, 0, 11, 6, 1, 9, 1, 4, 6, 4, 1, -1},
00331             {3, 11, 6, 3, 6, 0, 0, 6, 4, -1, -1, -1, -1, -1, -1, -1},
00332             {6, 4, 8, 11, 6, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00333             {7, 10, 6, 7, 8, 10, 8, 9, 10, -1, -1, -1, -1, -1, -1, -1},
00334             {0, 7, 3, 0, 10, 7, 0, 9, 10, 6, 7, 10, -1, -1, -1, -1},
00335             {10, 6, 7, 1, 10, 7, 1, 7, 8, 1, 8, 0, -1, -1, -1, -1},
00336             {10, 6, 7, 10, 7, 1, 1, 7, 3, -1, -1, -1, -1, -1, -1, -1},
00337             {1, 2, 6, 1, 6, 8, 1, 8, 9, 8, 6, 7, -1, -1, -1, -1},
00338             {2, 6, 9, 2, 9, 1, 6, 7, 9, 0, 9, 3, 7, 3, 9, -1},
00339             {7, 8, 0, 7, 0, 6, 6, 0, 2, -1, -1, -1, -1, -1, -1, -1},
00340             {7, 3, 2, 6, 7, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00341             {2, 3, 11, 10, 6, 8, 10, 8, 9, 8, 6, 7, -1, -1, -1, -1},
00342             {2, 0, 7, 2, 7, 11, 0, 9, 7, 6, 7, 10, 9, 10, 7, -1},
00343             {1, 8, 0, 1, 7, 8, 1, 10, 7, 6, 7, 10, 2, 3, 11, -1},
00344             {11, 2, 1, 11, 1, 7, 10, 6, 1, 6, 7, 1, -1, -1, -1, -1},
00345             {8, 9, 6, 8, 6, 7, 9, 1, 6, 11, 6, 3, 1, 3, 6, -1},
00346             {0, 9, 1, 11, 6, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00347             {7, 8, 0, 7, 0, 6, 3, 11, 0, 11, 6, 0, -1, -1, -1, -1},
00348             {7, 11, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00349             {7, 6, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00350             {3, 0, 8, 11, 7, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00351             {0, 1, 9, 11, 7, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00352             {8, 1, 9, 8, 3, 1, 11, 7, 6, -1, -1, -1, -1, -1, -1, -1},
00353             {10, 1, 2, 6, 11, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00354             {1, 2, 10, 3, 0, 8, 6, 11, 7, -1, -1, -1, -1, -1, -1, -1},
00355             {2, 9, 0, 2, 10, 9, 6, 11, 7, -1, -1, -1, -1, -1, -1, -1},
00356             {6, 11, 7, 2, 10, 3, 10, 8, 3, 10, 9, 8, -1, -1, -1, -1},
00357             {7, 2, 3, 6, 2, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00358             {7, 0, 8, 7, 6, 0, 6, 2, 0, -1, -1, -1, -1, -1, -1, -1},
00359             {2, 7, 6, 2, 3, 7, 0, 1, 9, -1, -1, -1, -1, -1, -1, -1},
00360             {1, 6, 2, 1, 8, 6, 1, 9, 8, 8, 7, 6, -1, -1, -1, -1},
00361             {10, 7, 6, 10, 1, 7, 1, 3, 7, -1, -1, -1, -1, -1, -1, -1},
00362             {10, 7, 6, 1, 7, 10, 1, 8, 7, 1, 0, 8, -1, -1, -1, -1},
00363             {0, 3, 7, 0, 7, 10, 0, 10, 9, 6, 10, 7, -1, -1, -1, -1},
00364             {7, 6, 10, 7, 10, 8, 8, 10, 9, -1, -1, -1, -1, -1, -1, -1},
00365             {6, 8, 4, 11, 8, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00366             {3, 6, 11, 3, 0, 6, 0, 4, 6, -1, -1, -1, -1, -1, -1, -1},
00367             {8, 6, 11, 8, 4, 6, 9, 0, 1, -1, -1, -1, -1, -1, -1, -1},
00368             {9, 4, 6, 9, 6, 3, 9, 3, 1, 11, 3, 6, -1, -1, -1, -1},
00369             {6, 8, 4, 6, 11, 8, 2, 10, 1, -1, -1, -1, -1, -1, -1, -1},
00370             {1, 2, 10, 3, 0, 11, 0, 6, 11, 0, 4, 6, -1, -1, -1, -1},
00371             {4, 11, 8, 4, 6, 11, 0, 2, 9, 2, 10, 9, -1, -1, -1, -1},
00372             {10, 9, 3, 10, 3, 2, 9, 4, 3, 11, 3, 6, 4, 6, 3, -1},
00373             {8, 2, 3, 8, 4, 2, 4, 6, 2, -1, -1, -1, -1, -1, -1, -1},
00374             {0, 4, 2, 4, 6, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00375             {1, 9, 0, 2, 3, 4, 2, 4, 6, 4, 3, 8, -1, -1, -1, -1},
00376             {1, 9, 4, 1, 4, 2, 2, 4, 6, -1, -1, -1, -1, -1, -1, -1},
00377             {8, 1, 3, 8, 6, 1, 8, 4, 6, 6, 10, 1, -1, -1, -1, -1},
00378             {10, 1, 0, 10, 0, 6, 6, 0, 4, -1, -1, -1, -1, -1, -1, -1},
00379             {4, 6, 3, 4, 3, 8, 6, 10, 3, 0, 3, 9, 10, 9, 3, -1},
00380             {10, 9, 4, 6, 10, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00381             {4, 9, 5, 7, 6, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00382             {0, 8, 3, 4, 9, 5, 11, 7, 6, -1, -1, -1, -1, -1, -1, -1},
00383             {5, 0, 1, 5, 4, 0, 7, 6, 11, -1, -1, -1, -1, -1, -1, -1},
00384             {11, 7, 6, 8, 3, 4, 3, 5, 4, 3, 1, 5, -1, -1, -1, -1},
00385             {9, 5, 4, 10, 1, 2, 7, 6, 11, -1, -1, -1, -1, -1, -1, -1},
00386             {6, 11, 7, 1, 2, 10, 0, 8, 3, 4, 9, 5, -1, -1, -1, -1},
00387             {7, 6, 11, 5, 4, 10, 4, 2, 10, 4, 0, 2, -1, -1, -1, -1},
00388             {3, 4, 8, 3, 5, 4, 3, 2, 5, 10, 5, 2, 11, 7, 6, -1},
00389             {7, 2, 3, 7, 6, 2, 5, 4, 9, -1, -1, -1, -1, -1, -1, -1},
00390             {9, 5, 4, 0, 8, 6, 0, 6, 2, 6, 8, 7, -1, -1, -1, -1},
00391             {3, 6, 2, 3, 7, 6, 1, 5, 0, 5, 4, 0, -1, -1, -1, -1},
00392             {6, 2, 8, 6, 8, 7, 2, 1, 8, 4, 8, 5, 1, 5, 8, -1},
00393             {9, 5, 4, 10, 1, 6, 1, 7, 6, 1, 3, 7, -1, -1, -1, -1},
00394             {1, 6, 10, 1, 7, 6, 1, 0, 7, 8, 7, 0, 9, 5, 4, -1},
00395             {4, 0, 10, 4, 10, 5, 0, 3, 10, 6, 10, 7, 3, 7, 10, -1},
00396             {7, 6, 10, 7, 10, 8, 5, 4, 10, 4, 8, 10, -1, -1, -1, -1},
00397             {6, 9, 5, 6, 11, 9, 11, 8, 9, -1, -1, -1, -1, -1, -1, -1},
00398             {3, 6, 11, 0, 6, 3, 0, 5, 6, 0, 9, 5, -1, -1, -1, -1},
00399             {0, 11, 8, 0, 5, 11, 0, 1, 5, 5, 6, 11, -1, -1, -1, -1},
00400             {6, 11, 3, 6, 3, 5, 5, 3, 1, -1, -1, -1, -1, -1, -1, -1},
00401             {1, 2, 10, 9, 5, 11, 9, 11, 8, 11, 5, 6, -1, -1, -1, -1},
00402             {0, 11, 3, 0, 6, 11, 0, 9, 6, 5, 6, 9, 1, 2, 10, -1},
00403             {11, 8, 5, 11, 5, 6, 8, 0, 5, 10, 5, 2, 0, 2, 5, -1},
00404             {6, 11, 3, 6, 3, 5, 2, 10, 3, 10, 5, 3, -1, -1, -1, -1},
00405             {5, 8, 9, 5, 2, 8, 5, 6, 2, 3, 8, 2, -1, -1, -1, -1},
00406             {9, 5, 6, 9, 6, 0, 0, 6, 2, -1, -1, -1, -1, -1, -1, -1},
00407             {1, 5, 8, 1, 8, 0, 5, 6, 8, 3, 8, 2, 6, 2, 8, -1},
00408             {1, 5, 6, 2, 1, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00409             {1, 3, 6, 1, 6, 10, 3, 8, 6, 5, 6, 9, 8, 9, 6, -1},
00410             {10, 1, 0, 10, 0, 6, 9, 5, 0, 5, 6, 0, -1, -1, -1, -1},
00411             {0, 3, 8, 5, 6, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00412             {10, 5, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00413             {11, 5, 10, 7, 5, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00414             {11, 5, 10, 11, 7, 5, 8, 3, 0, -1, -1, -1, -1, -1, -1, -1},
00415             {5, 11, 7, 5, 10, 11, 1, 9, 0, -1, -1, -1, -1, -1, -1, -1},
00416             {10, 7, 5, 10, 11, 7, 9, 8, 1, 8, 3, 1, -1, -1, -1, -1},
00417             {11, 1, 2, 11, 7, 1, 7, 5, 1, -1, -1, -1, -1, -1, -1, -1},
00418             {0, 8, 3, 1, 2, 7, 1, 7, 5, 7, 2, 11, -1, -1, -1, -1},
00419             {9, 7, 5, 9, 2, 7, 9, 0, 2, 2, 11, 7, -1, -1, -1, -1},
00420             {7, 5, 2, 7, 2, 11, 5, 9, 2, 3, 2, 8, 9, 8, 2, -1},
00421             {2, 5, 10, 2, 3, 5, 3, 7, 5, -1, -1, -1, -1, -1, -1, -1},
00422             {8, 2, 0, 8, 5, 2, 8, 7, 5, 10, 2, 5, -1, -1, -1, -1},
00423             {9, 0, 1, 5, 10, 3, 5, 3, 7, 3, 10, 2, -1, -1, -1, -1},
00424             {9, 8, 2, 9, 2, 1, 8, 7, 2, 10, 2, 5, 7, 5, 2, -1},
00425             {1, 3, 5, 3, 7, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00426             {0, 8, 7, 0, 7, 1, 1, 7, 5, -1, -1, -1, -1, -1, -1, -1},
00427             {9, 0, 3, 9, 3, 5, 5, 3, 7, -1, -1, -1, -1, -1, -1, -1},
00428             {9, 8, 7, 5, 9, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00429             {5, 8, 4, 5, 10, 8, 10, 11, 8, -1, -1, -1, -1, -1, -1, -1},
00430             {5, 0, 4, 5, 11, 0, 5, 10, 11, 11, 3, 0, -1, -1, -1, -1},
00431             {0, 1, 9, 8, 4, 10, 8, 10, 11, 10, 4, 5, -1, -1, -1, -1},
00432             {10, 11, 4, 10, 4, 5, 11, 3, 4, 9, 4, 1, 3, 1, 4, -1},
00433             {2, 5, 1, 2, 8, 5, 2, 11, 8, 4, 5, 8, -1, -1, -1, -1},
00434             {0, 4, 11, 0, 11, 3, 4, 5, 11, 2, 11, 1, 5, 1, 11, -1},
00435             {0, 2, 5, 0, 5, 9, 2, 11, 5, 4, 5, 8, 11, 8, 5, -1},
00436             {9, 4, 5, 2, 11, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00437             {2, 5, 10, 3, 5, 2, 3, 4, 5, 3, 8, 4, -1, -1, -1, -1},
00438             {5, 10, 2, 5, 2, 4, 4, 2, 0, -1, -1, -1, -1, -1, -1, -1},
00439             {3, 10, 2, 3, 5, 10, 3, 8, 5, 4, 5, 8, 0, 1, 9, -1},
00440             {5, 10, 2, 5, 2, 4, 1, 9, 2, 9, 4, 2, -1, -1, -1, -1},
00441             {8, 4, 5, 8, 5, 3, 3, 5, 1, -1, -1, -1, -1, -1, -1, -1},
00442             {0, 4, 5, 1, 0, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00443             {8, 4, 5, 8, 5, 3, 9, 0, 5, 0, 3, 5, -1, -1, -1, -1},
00444             {9, 4, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00445             {4, 11, 7, 4, 9, 11, 9, 10, 11, -1, -1, -1, -1, -1, -1, -1},
00446             {0, 8, 3, 4, 9, 7, 9, 11, 7, 9, 10, 11, -1, -1, -1, -1},
00447             {1, 10, 11, 1, 11, 4, 1, 4, 0, 7, 4, 11, -1, -1, -1, -1},
00448             {3, 1, 4, 3, 4, 8, 1, 10, 4, 7, 4, 11, 10, 11, 4, -1},
00449             {4, 11, 7, 9, 11, 4, 9, 2, 11, 9, 1, 2, -1, -1, -1, -1},
00450             {9, 7, 4, 9, 11, 7, 9, 1, 11, 2, 11, 1, 0, 8, 3, -1},
00451             {11, 7, 4, 11, 4, 2, 2, 4, 0, -1, -1, -1, -1, -1, -1, -1},
00452             {11, 7, 4, 11, 4, 2, 8, 3, 4, 3, 2, 4, -1, -1, -1, -1},
00453             {2, 9, 10, 2, 7, 9, 2, 3, 7, 7, 4, 9, -1, -1, -1, -1},
00454             {9, 10, 7, 9, 7, 4, 10, 2, 7, 8, 7, 0, 2, 0, 7, -1},
00455             {3, 7, 10, 3, 10, 2, 7, 4, 10, 1, 10, 0, 4, 0, 10, -1},
00456             {1, 10, 2, 8, 7, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00457             {4, 9, 1, 4, 1, 7, 7, 1, 3, -1, -1, -1, -1, -1, -1, -1},
00458             {4, 9, 1, 4, 1, 7, 0, 8, 1, 8, 7, 1, -1, -1, -1, -1},
00459             {4, 0, 3, 7, 4, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00460             {4, 8, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00461             {9, 10, 8, 10, 11, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00462             {3, 0, 9, 3, 9, 11, 11, 9, 10, -1, -1, -1, -1, -1, -1, -1},
00463             {0, 1, 10, 0, 10, 8, 8, 10, 11, -1, -1, -1, -1, -1, -1, -1},
00464             {3, 1, 10, 11, 3, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00465             {1, 2, 11, 1, 11, 9, 9, 11, 8, -1, -1, -1, -1, -1, -1, -1},
00466             {3, 0, 9, 3, 9, 11, 1, 2, 9, 2, 11, 9, -1, -1, -1, -1},
00467             {0, 2, 11, 8, 0, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00468             {3, 2, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00469             {2, 3, 8, 2, 8, 10, 10, 8, 9, -1, -1, -1, -1, -1, -1, -1},
00470             {9, 10, 2, 0, 9, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00471             {2, 3, 8, 2, 8, 10, 0, 1, 8, 1, 10, 8, -1, -1, -1, -1},
00472             {1, 10, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00473             {1, 3, 8, 9, 1, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00474             {0, 9, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00475             {0, 3, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00476             {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1}
00477           };
00478           return table[idx1][idx2];
00479         };
00480 
00481       public:
00482 
00489         template<typename grid_type,typename mesh_type>
00490         void operator()(grid_type const & phi, typename grid_type::value_type const & isolevel, mesh_type & mesh)
00491         {
00492           typedef typename mesh_type::face_handle     face_handle;
00493           typedef typename mesh_type::vertex_handle   vertex_handle;
00494           typedef typename mesh_type::math_types      math_types;
00495           typedef typename math_types::value_traits   value_traits;
00496           typedef typename math_types::vector3_type   vector3_type;
00497           typedef typename math_types::real_type      real_type;
00498 
00499           typedef typename grid_type::iterator         iterator;
00500           typedef typename grid_type::value_type       value_type;
00501 
00502           mesh.clear();
00503 
00504           size_t I = phi.I();
00505           size_t J = phi.J();
00506           size_t K = phi.K();
00507 
00508           //typedef std::map<identifier_type,vertex_handle> lut_type;
00509           //lut_type lut;
00510           typedef std::vector<vertex_handle> lut_type;
00511           lut_type lut(3*I*J*K);
00512 
00513           static value_type const unknown = phi.unused();// Convenience for better readability
00514 
00515           for (size_t k = 0; k < K-1; ++k)
00516           {
00517             for (size_t j = 0; j < J-1; ++j)
00518             {
00519               value_type const *  i0j0k0 = &(phi(0,j,k));
00520               value_type const *  i1j0k0 = &(phi(1,j,k));
00521               value_type const *  i0j1k0 = &(phi(0,j+1,k));
00522               value_type const *  i1j1k0 = &(phi(1,j+1,k));
00523               value_type const *  i0j0k1 = &(phi(0,j,k+1));
00524               value_type const *  i1j0k1 = &(phi(1,j,k+1));
00525               value_type const *  i0j1k1 = &(phi(0,j+1,k+1));
00526               value_type const *  i1j1k1 = &(phi(1,j+1,k+1));
00527               for (size_t i = 0; i < I-1; ++i)
00528               {
00529                 identifier_type table_index = static_cast<identifier_type>(0);
00530 
00531                 assert( phi(i  ,  j,  k)  == *i0j0k0);
00532                 assert( phi(i+1,  j,  k)  == *i1j0k0);
00533                 assert( phi(i  ,j+1,  k)  == *i0j1k0);
00534                 assert( phi(i+1,j+1,  k)  == *i1j1k0);
00535                 assert( phi(i  ,  j,k+1)  == *i0j0k1);
00536                 assert( phi(i+1,  j,k+1)  == *i1j0k1);
00537                 assert( phi(i  ,j+1,k+1)  == *i0j1k1);
00538                 assert( phi(i+1,j+1,k+1)  == *i1j1k1);
00539 
00540 
00541                 if( *i0j0k0 < isolevel ) table_index |= 1;
00542                 if( *i0j1k0 < isolevel ) table_index |= 2;
00543                 if( *i1j1k0 < isolevel ) table_index |= 4;
00544                 if( *i1j0k0 < isolevel ) table_index |= 8;
00545                 if( *i0j0k1 < isolevel ) table_index |= 16;
00546                 if( *i0j1k1 < isolevel ) table_index |= 32;
00547                 if( *i1j1k1 < isolevel ) table_index |= 64;
00548                 if( *i1j0k1 < isolevel ) table_index |= 128;
00549 
00550                 // Part of current cell was unknown so we just skip
00551                 // ahead to next cell
00552                 bool goto_next = ( ( *i0j0k0 == unknown )||( *i0j1k0 == unknown )||( *i1j1k0 == unknown )||( *i1j0k0 == unknown )||( *i0j0k1 == unknown )||( *i0j1k1 == unknown )||( *i1j1k1 == unknown )||( *i1j0k1 == unknown )  );
00553 
00554                 ++i0j0k0;
00555                 ++i1j0k0;
00556                 ++i0j1k0;
00557                 ++i1j1k0;
00558                 ++i0j0k1;
00559                 ++i1j0k1;
00560                 ++i0j1k1;
00561                 ++i1j1k1;
00562 
00563                 if(goto_next) continue;
00564 
00565                 assert(table_index >= 0  || !"table idx out of range");
00566                 assert(table_index < 256 || !"table idx out of range");
00567                 identifier_type  cell_case = edge_table(table_index);
00568 
00569                 // Now create a triangulation of the isosurface in this cell.
00570                 if (cell_case == 0)
00571                   continue;
00572 
00573                 //--- unrolling calculate_intersection may speed up things???
00574 
00575 #define test_edge(binary_edge_code,edge_number) \
00576   if (cell_case & binary_edge_code) \
00577                 {\
00578                 identifier_type  id = get_edge_id(i, j, k, edge_number, I, J);\
00579                 vector3_type p = calculate_intersection(phi,isolevel,i, j, k, edge_number);\
00580                 vertex_handle h = mesh.add_vertex( p ); \
00581                 assert(!h.is_null() || !"could not create vertex");\
00582                 lut[id] = h; \
00583                 }
00584 
00585                 test_edge(8,3);
00586                 test_edge(1,0);
00587                 test_edge(256,8);
00588                 if (i == I - 2)
00589                 {
00590                   test_edge(4,2);
00591                   test_edge(2048,11);
00592                 }
00593                 if (j == J - 2)
00594                 {
00595                   test_edge(2,1);
00596                   test_edge(512,9);
00597                 }
00598                 if (k == K - 2)
00599                 {
00600                   test_edge(16,4);
00601                   test_edge(128,7);
00602                 }
00603                 if ((i==I - 2) && (j==J - 2))
00604                 {
00605                   test_edge(1024,10);
00606                 }
00607                 if ((i==I - 2) && (k==K - 2))
00608                 {
00609                   test_edge(64,6);
00610                 }
00611                 if ((j==J - 2) && (k==K - 2))
00612                 {
00613                   test_edge(32,5);
00614                 }
00615 #undef test_edge \
00616 
00617               }
00618             }
00619           }
00620 
00621           for (size_t k = 0; k < (K-1); ++k)
00622           {
00623             for (size_t j = 0; j < (J-1); ++j)
00624             {
00625               value_type const *  i0j0k0 = &(phi(0,j,k));
00626               value_type const *  i1j0k0 = &(phi(1,j,k));
00627               value_type const *  i0j1k0 = &(phi(0,j+1,k));
00628               value_type const *  i1j1k0 = &(phi(1,j+1,k));
00629               value_type const *  i0j0k1 = &(phi(0,j,k+1));
00630               value_type const *  i1j0k1 = &(phi(1,j,k+1));
00631               value_type const *  i0j1k1 = &(phi(0,j+1,k+1));
00632               value_type const *  i1j1k1 = &(phi(1,j+1,k+1));
00633               for (size_t i = 0; i < (I-1); ++i)
00634               {
00635                 identifier_type table_index = static_cast<identifier_type>(0);
00636 
00637                 if( *i0j0k0 < isolevel ) table_index |= 1;
00638                 if( *i0j1k0 < isolevel ) table_index |= 2;
00639                 if( *i1j1k0 < isolevel ) table_index |= 4;
00640                 if( *i1j0k0 < isolevel ) table_index |= 8;
00641                 if( *i0j0k1 < isolevel ) table_index |= 16;
00642                 if( *i0j1k1 < isolevel ) table_index |= 32;
00643                 if( *i1j1k1 < isolevel ) table_index |= 64;
00644                 if( *i1j0k1 < isolevel ) table_index |= 128;
00645 
00646                 // Part of current cell was unknown so we just skip
00647                 // ahead to next cell
00648                 bool goto_next = ( ( *i0j0k0 == unknown )||( *i0j1k0 == unknown )||( *i1j1k0 == unknown )||( *i1j0k0 == unknown )||( *i0j0k1 == unknown )||( *i0j1k1 == unknown )||( *i1j1k1 == unknown )||( *i1j0k1 == unknown )  );
00649 
00650                 ++i0j0k0;
00651                 ++i1j0k0;
00652                 ++i0j1k0;
00653                 ++i1j1k0;
00654                 ++i0j0k1;
00655                 ++i1j0k1;
00656                 ++i0j1k1;
00657                 ++i1j1k1;
00658 
00659                 if(goto_next) continue;
00660 
00661                 assert(table_index >= 0  || !"table idx out of range");
00662                 assert(table_index < 256 || !"table idx out of range");
00663                 if (table_index==0 || table_index==255)  //--- if (edge_table(table_index) == 0)
00664                   continue;
00665                 //--- It may be a bit expensive to actually use a function call to look up in the triangle table??
00666                 for (identifier_type t = 0; triangle_table(table_index,t) != -1; t += 3)
00667                 {
00668                   identifier_type id0 = get_edge_id(i, j, k, triangle_table(table_index,t  ), I, J);
00669                   identifier_type id1 = get_edge_id(i, j, k, triangle_table(table_index,t+1), I, J);
00670                   identifier_type id2 = get_edge_id(i, j, k, triangle_table(table_index,t+2), I, J);
00671                   mesh.add_face(lut[id0],lut[id1],lut[id2]);
00672                 }
00673               }
00674             }
00675           }
00676 
00677         }
00678 
00679       };
00680 
00681     }// namespace detail
00682 
00683 
00691     template<typename grid_type,typename mesh_type>
00692     void isosurface(grid_type const & phi, typename grid_type::value_type const & isolevel, mesh_type & mesh)
00693     {
00694       detail::IsoSurfaceGenerator iso;
00695       iso(phi,isolevel,mesh);
00696 
00697       std::cout << "isosurface() : extracted " << mesh.size_vertices() << " vertices " << mesh.size_faces() << " faces" << std::endl;
00698     }
00699 
00700   } // namespace mesh
00701 } // namespace OpenTissue
00702 
00703 //OPENTISSUE_CORE_CONTAINERS_MESH_COMMON_UTIL_MESH_ISOSURFACE_H
00704 #endif
/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/core/containers/mesh/common/util/mesh_isosurface.h
