/opencascade-6.5.1/ros/src/Voxel/Voxel_ROctBoolDS.cxx

https://github.com/jehc/MondocosmOS · C++ · 717 lines · 554 code · 77 blank · 86 comment · 84 complexity · 33af78a4298adb5fa999bc947c5e1509 MD5 · raw file 

    

  1. // File:	Voxel_ROctBoolDS.cxx
    2. 

  2. // Created:	Mon Sep 01 11:00:08 2008
    3. 

  3. // Author:	Vladislav ROMASHKO
    4. 

  4. //		<vladislav.romashko@opencascade.com>
    5. 

  5. 

  6. #include <Voxel_ROctBoolDS.ixx>
    7. 

  7. #include <Voxel_SplitData.hxx>
    8. 

  8. 

  9. #include <stdlib.h>
    10. 

  10. 

  11. static Standard_Byte gbits[8] = {1, 2, 4, 8, 16, 32, 64, 128};
    12. 

  12. static Standard_Byte gnbits[8] = {255-1, 255-2, 255-4, 255-8, 255-16, 255-32, 255-64, 255-128};
    13. 

  13. 

  14. /* Data structure of the ROctBoolDS
    15. 

    16. 

  16. SplitData: 1 byte (8 values)
    17. 

  17.        (a) SplitData: 8 bytes (64 values)
    18. 

  18.                   (b) SplitData: 64 bytes (512 values)
    19. 

  19.                              (c) SplitData: ...
    20. 

  20.                                         (d) SplitData: ...
    21. 

  21. */
    22. 

  22. 

  23. // Empty constructor
    24. 

  24. Voxel_ROctBoolDS::Voxel_ROctBoolDS():Voxel_DS()
    25. 

  25. {
    26. 

  26. 

  27. }
    28. 

  28. 

  29. // Constructor with intialization.
    30. 

  30. Voxel_ROctBoolDS::Voxel_ROctBoolDS(const Standard_Real      x, const Standard_Real      y, const Standard_Real      z, 
    31. 

  31. 				   const Standard_Real   xlen, const Standard_Real   ylen, const Standard_Real   zlen,
    32. 

  32. 				   const Standard_Integer nbx, const Standard_Integer nby, const Standard_Integer nbz)
    33. 

  33. :Voxel_DS()
    34. 

  34. {
    35. 

  35.   Init(x, y, z, xlen, ylen, zlen, nbx, nby, nbz);
    36. 

  36. }
    37. 

  37. 

  38. // Initialization.
    39. 

  39. void Voxel_ROctBoolDS::Init(const Standard_Real      x, const Standard_Real      y, const Standard_Real      z, 
    40. 

  40. 			    const Standard_Real   xlen, const Standard_Real   ylen, const Standard_Real   zlen,
    41. 

  41. 			    const Standard_Integer nbx, const Standard_Integer nby, const Standard_Integer nbz)
    42. 

  42. {
    43. 

  43.   Destroy();
    44. 

  44. 

  45.   Voxel_DS::Init(x, y, z, xlen, ylen, zlen, nbx, nby, nbz);
    46. 

  46. 

  47.   if (!myNbX || !myNbY || !myNbZ)
    48. 

  48.     return;
    49. 

  49. 

  50.   Standard_Integer nb_slices = RealToInt(ceil(myNbXY * myNbZ / 8.0));
    51. 

  51.   myData = (Standard_Address) calloc(nb_slices, sizeof(Voxel_SplitData*));
    52. 

  52. }
    53. 

  53. 

  54. // Destructor
    55. 

  55. void Voxel_ROctBoolDS::Destroy()
    56. 

  56. {
    57. 

  57.   if (myData)
    58. 

  58.   {
    59. 

  59.     SetZero();
    60. 

  60.     free((Voxel_SplitData**)myData);
    61. 

  61.     myData = 0;
    62. 

  62.   }
    63. 

  63. }
    64. 

  64. 

  65. // A recursive method of deletion of data.
    66. 

  66. static void SetZeroSplitData(Voxel_SplitData* data)
    67. 

  67. {
    68. 

  68.   // Values:
    69. 

  69.   free((Standard_Byte*) data->GetValues());
    70. 

  70.   data->GetValues() = 0;
    71. 

  71.   if (data->GetSplitData())
    72. 

  72.   {
    73. 

  73.     SetZeroSplitData((Voxel_SplitData*) data->GetSplitData());
    74. 

  74.   }
    75. 

  75.   delete (data);
    76. 

  76.   data = 0;
    77. 

  77. }
    78. 

  78. 

  79. void Voxel_ROctBoolDS::SetZero()
    80. 

  80. {
    81. 

  81.   if (myData)
    82. 

  82.   {
    83. 

  83.     Standard_Integer ix = 0, nb_slices = RealToInt(ceil(myNbXY * myNbZ / 8.0));
    84. 

  84.     for (; ix < nb_slices; ix++)
    85. 

  85.     {
    86. 

  86.       if (((Voxel_SplitData**)myData)[ix])
    87. 

  87.       {
    88. 

  88. 	SetZeroSplitData((Voxel_SplitData*)((Voxel_SplitData**)myData)[ix]);
    89. 

  89.       }
    90. 

  90.     }
    91. 

  91.   }
    92. 

  92. }
    93. 

  93. 

  94. // Access to the boolean information attached to a particular voxel:
    95. 

  95. // Info: (ix >= 0 && ix < theNb_x), etc.
    96. 

  96. void Voxel_ROctBoolDS::Set(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz, 
    97. 

  97. 			   const Standard_Boolean data)
    98. 

  98. {
    99. 

  99.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    100. 

  100.   Standard_Integer islice = ibit >> 3;
    101. 

  101. 

  102.   if (!data && !((Voxel_SplitData**)myData)[islice])
    103. 

  103.     return; // don't allocate a slice of data for setting a 0 value
    104. 

  104. 

  105.   // Allocate the slice if it is not done yet.
    106. 

  106.   if (!((Voxel_SplitData**)myData)[islice])
    107. 

  107.   {
    108. 

  108.     ((Voxel_SplitData**)myData)[islice] = (Voxel_SplitData*) new Voxel_SplitData;
    109. 

  109.     // Values:
    110. 

  110.     ((Voxel_SplitData**)myData)[islice]->GetValues() = 
    111. 

  111.       (Standard_Byte*) calloc(1/*one byte: 8 1-bit values*/, sizeof(Standard_Byte));
    112. 

  112.     // Sub-voxels:
    113. 

  113.     ((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData() = 0;
    114. 

  114.   }
    115. 

  115. 

  116.   // Value
    117. 

  117.   Standard_Byte value = *((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues());
    118. 

  118. 

  119.   // Position of data in the 8 bit-"value".
    120. 

  120.   Standard_Integer shift = ibit - (islice << 3);
    121. 

  121. 

  122.   // Set data
    123. 

  123.   if (data != ((value & gbits[shift]) ? Standard_True : Standard_False))
    124. 

  124.   {
    125. 

  125.     if (data)
    126. 

  126.       value |= gbits[shift];
    127. 

  127.     else
    128. 

  128.       value &= gnbits[shift];
    129. 

  129.     *((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues()) = value;
    130. 

  130.   }
    131. 

  131. 

  132.   // Set the same value to sub-voxels.
    133. 

  133.   if (((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    134. 

  134.   {
    135. 

  135.     // Get sub-value
    136. 

  136.     Standard_Byte subvalue = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[shift];
    137. 

  137. 

  138.     // Set sub-value
    139. 

  139.     if (subvalue != (data ? 255 : 0))
    140. 

  140.     {
    141. 

  141.       subvalue = data ? 255 : 0;
    142. 

  142.       ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[shift] = subvalue;
    143. 

  143.     }
    144. 

  144. 

  145.     // Set the same value to sub-sub-voxels.
    146. 

  146.     if (((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())
    147. 

  147.     {
    148. 

  148.       // Start index of 64-bit value (index of byte of sub-sub-voxel).
    149. 

  149.       Standard_Integer ibyte2 = (shift << 3);
    150. 

  150.       for (Standard_Integer ioct2 = 0; ioct2 < 8; ioct2++)
    151. 

  151.       {
    152. 

  152. 	// Get sub-sub-value
    153. 

  153. 	Standard_Byte subsubvalue = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2 + ioct2];
    154. 

  154. 

  155. 	// Set sub-sub-value
    156. 

  156. 	if (subsubvalue != (data ? 255 : 0))
    157. 

  157. 	{
    158. 

  158. 	  subsubvalue = data ? 255 : 0;
    159. 

  159. 	  ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2 + ioct2] = subsubvalue;
    160. 

  160. 	}
    161. 

  161.       }
    162. 

  162.     }
    163. 

  163.   }
    164. 

  164. }
    165. 

  165. 

  166. void Voxel_ROctBoolDS::Set(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz, 
    167. 

  167. 			   const Standard_Integer ioct1, const Standard_Boolean data)
    168. 

  168. {
    169. 

  169.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    170. 

  170.   Standard_Integer islice = ibit >> 3;
    171. 

  171. 

  172.   if (!data && !((Voxel_SplitData**)myData)[islice])
    173. 

  173.     return; // don't allocate a slice of data for setting a 0 value
    174. 

  174. 

  175.   // Allocate the slice if it is not done yet.
    176. 

  176.   if (!((Voxel_SplitData**)myData)[islice])
    177. 

  177.   {
    178. 

  178.     ((Voxel_SplitData**)myData)[islice] = (Voxel_SplitData*) new Voxel_SplitData;
    179. 

  179.     // Values:
    180. 

  180.     ((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues() = 
    181. 

  181.       (Standard_Byte*) calloc(1/*one byte: 8 1-bit values*/, sizeof(Standard_Byte));
    182. 

  182.     // Sub-voxels:
    183. 

  183.     ((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData() = 0;
    184. 

  184.   }
    185. 

  185. 

  186.   // Check sub-voxels of the first level
    187. 

  187.   if (!((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    188. 

  188.   {
    189. 

  189.     // Sub-voxels:
    190. 

  190.     ((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData() = (Voxel_SplitData*) new Voxel_SplitData;
    191. 

  191.     // Value of sub-voxels:
    192. 

  192.     ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues() = 
    193. 

  193.       (Standard_Byte*) calloc(8/*eight bytes: 8 sub-voxels for each voxel*/, sizeof(Standard_Byte));
    194. 

  194. 

  195.     // Set parent value
    196. 

  196.     Standard_Byte parent_value = *((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues());
    197. 

  197.     if (parent_value)
    198. 

  198.     {
    199. 

  199.       for (Standard_Integer shift = 0; shift < 8; shift++)
    200. 

  200.       {
    201. 

  201. 	if ((parent_value & gbits[shift]) ? Standard_True : Standard_False)
    202. 

  202. 	{
    203. 

  203. 	  ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[shift] = 255;
    204. 

  204. 	}
    205. 

  205. 	else
    206. 

  206. 	{
    207. 

  207. 	  ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[shift] = 0;
    208. 

  208. 	}
    209. 

  209.       }
    210. 

  210.     }
    211. 

  211. 

  212.     // Sub-sub-voxels
    213. 

  213.     ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData() = 0;
    214. 

  214.   }
    215. 

  215. 

  216.   // Index of sub-voxel corresponding to ioct1: 8 voxels correspond to 64 sub-voxels.
    217. 

  217.   Standard_Integer ibyte = ibit - (islice << 3);
    218. 

  218. 

  219.   // Value
    220. 

  220.   Standard_Byte value = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[ibyte];
    221. 

  221. 

  222.   // Set data
    223. 

  223.   if (data != ((value & gbits[ioct1]) ? Standard_True : Standard_False))
    224. 

  224.   {
    225. 

  225.     if (data)
    226. 

  226.       value |= gbits[ioct1];
    227. 

  227.     else
    228. 

  228.       value &= gnbits[ioct1];
    229. 

  229.     ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[ibyte] = value;
    230. 

  230.   }
    231. 

  231. 

  232.   // Set the same value to sub-voxels.
    233. 

  233.   if (((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())
    234. 

  234.   {
    235. 

  235.     // Start index of 64-bit value (index of byte of sub-sub-voxel).
    236. 

  236.     Standard_Integer ibyte2 = (ibyte << 3) + ioct1;
    237. 

  237. 

  238.     // Get sub-sub-value
    239. 

  239.     Standard_Byte subsubvalue = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2];
    240. 

  240. 

  241.     // Set sub-sub-value
    242. 

  242.     if (subsubvalue != (data ? 255 : 0))
    243. 

  243.     {
    244. 

  244.       subsubvalue = data ? 255 : 0;
    245. 

  245.       ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2] = subsubvalue;
    246. 

  246.     }
    247. 

  247.   }
    248. 

  248. }
    249. 

  249. 

  250. void Voxel_ROctBoolDS::Set(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz, 
    251. 

  251. 			   const Standard_Integer ioct1, const Standard_Integer ioct2, const Standard_Boolean data)
    252. 

  252. {
    253. 

  253.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    254. 

  254.   Standard_Integer islice = ibit >> 3;
    255. 

  255. 

  256.   if (!data && !((Voxel_SplitData**)myData)[islice])
    257. 

  257.     return; // don't allocate a slice of data for setting a 0 value
    258. 

  258. 

  259.   // Allocate the slice if it is not done yet.
    260. 

  260.   if (!((Voxel_SplitData**)myData)[islice])
    261. 

  261.   {
    262. 

  262.     ((Voxel_SplitData**)myData)[islice] = (Voxel_SplitData*) new Voxel_SplitData;
    263. 

  263.     // Values:
    264. 

  264.     ((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues() = 
    265. 

  265.       (Standard_Byte*) calloc(1/*one byte: 8 1-bit values*/, sizeof(Standard_Byte));
    266. 

  266.     // Sub-voxels:
    267. 

  267.     ((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData() = 0;
    268. 

  268.   }
    269. 

  269. 

  270.   // Check sub-voxels of the first level
    271. 

  271.   if (!((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    272. 

  272.   {
    273. 

  273.     // Sub-voxels:
    274. 

  274.     ((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData() = (Voxel_SplitData*) new Voxel_SplitData;
    275. 

  275.     // Value of sub-voxels:
    276. 

  276.     ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues() = 
    277. 

  277.       (Standard_Byte*) calloc(8/*eight bytes: 8 sub-voxels for each voxel*/, sizeof(Standard_Byte));
    278. 

  278. 

  279.     // Set parent value
    280. 

  280.     Standard_Byte parent_value = *((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues());
    281. 

  281.     if (parent_value)
    282. 

  282.     {
    283. 

  283.       for (Standard_Integer shift = 0; shift < 8; shift++)
    284. 

  284.       {
    285. 

  285. 	if ((parent_value & gbits[shift]) ? Standard_True : Standard_False)
    286. 

  286. 	{
    287. 

  287. 	  ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[shift] = 255;
    288. 

  288. 	}
    289. 

  289. 	else
    290. 

  290. 	{
    291. 

  291. 	  ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[shift] = 0;
    292. 

  292. 	}
    293. 

  293.       }
    294. 

  294.     }
    295. 

  295. 

  296.      // Sub-sub-voxels
    297. 

  297.     ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData() = 0;
    298. 

  298.   }
    299. 

  299. 

  300.   // Check sub-voxels of the second level
    301. 

  301.   if (!((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())
    302. 

  302.   {
    303. 

  303.     // Sub-voxels 2:
    304. 

  304.     ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData() = 
    305. 

  305.       (Voxel_SplitData*) new Voxel_SplitData;
    306. 

  306.     // Value of sub-voxels 2:
    307. 

  307.     ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues() = 
    308. 

  308.       (Standard_Byte*) calloc(64/*sixty four bytes: 8 sub-voxels for each sub-voxel for each voxel*/, 
    309. 

  309. 			      sizeof(Standard_Byte));
    310. 

  310. 

  311.     // Set parent value
    312. 

  312.     for (Standard_Integer ibyte1 = 0; ibyte1 < 8; ibyte1++)
    313. 

  313.     {
    314. 

  314.       Standard_Byte parent_value = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[ibyte1];
    315. 

  315.       if (parent_value)
    316. 

  316.       {
    317. 

  317. 	Standard_Integer ibyte2 = (ibyte1 << 3);
    318. 

  318. 	for (Standard_Integer shift = 0; shift < 8; shift++)
    319. 

  319. 	{
    320. 

  320. 	  if ((parent_value & gbits[shift]) ? Standard_True : Standard_False)
    321. 

  321. 	  {
    322. 

  322. 	    ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2 + shift] = 255;
    323. 

  323. 	  }
    324. 

  324. 	  else
    325. 

  325. 	  {
    326. 

  326. 	    ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2 + shift] = 0;
    327. 

  327. 	  }
    328. 

  328. 	}
    329. 

  329.       }
    330. 

  330.     }
    331. 

  331. 

  332.     // Sub-sub-sub-voxels
    333. 

  333.     ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetSplitData() = 0;
    334. 

  334.   }
    335. 

  335. 

  336.   // Index of sub-voxel corresponding to ioct1: 8 voxels correspond to 64 sub-voxels.
    337. 

  337.   Standard_Integer ibyte1 = ibit - (islice << 3);  // imdex of byte of  8-byte value (sub-voxel 1).
    338. 

  338.   Standard_Integer ibyte2 = (ibyte1 << 3) + ioct1; // index of byte of 64-byte value (sub-voxel 2)
    339. 

  339. 

  340.   // Value
    341. 

  341.   Standard_Byte value = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2];
    342. 

  342. 

  343.   // Set data
    344. 

  344.   if (data != ((value & gbits[ioct2]) ? Standard_True : Standard_False))
    345. 

  345.   {
    346. 

  346.     if (data)
    347. 

  347.       value |= gbits[ioct2];
    348. 

  348.     else
    349. 

  349.       value &= gnbits[ioct2];
    350. 

  350.     ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2] = value;
    351. 

  351.   }
    352. 

  352. }
    353. 

  353. 

  354. Standard_Boolean Voxel_ROctBoolDS::Get(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz) const
    355. 

  355. {
    356. 

  356.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    357. 

  357.   Standard_Integer islice = ibit >> 3;
    358. 

  358. 

  359.   // If the slice of data is not allocated, it means that its values are 0.
    360. 

  360.   if (!((Voxel_SplitData**)myData)[islice])
    361. 

  361.     return Standard_False;
    362. 

  362. 

  363.   // Value (byte)
    364. 

  364.   Standard_Byte value = *((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues());
    365. 

  365. 

  366.   // Position of data in the 8 bit-"value".
    367. 

  367.   Standard_Integer shift = ibit - (islice << 3);
    368. 

  368.   return ((value & gbits[shift]) ? Standard_True : Standard_False);
    369. 

  369. }
    370. 

  370. 

  371. Standard_Boolean Voxel_ROctBoolDS::Get(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz,
    372. 

  372. 				      const Standard_Integer ioct1) const
    373. 

  373. {
    374. 

  374.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    375. 

  375.   Standard_Integer islice = ibit >> 3;
    376. 

  376. 

  377.   // If the slice of data is not allocated, it means that its values are 0.
    378. 

  378.   if (!((Voxel_SplitData**)myData)[islice])
    379. 

  379.     return Standard_False;
    380. 

  380. 

  381.   // If the voxel is not split, return the value of the voxel.
    382. 

  382.   if (!((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    383. 

  383.     return Get(ix, iy, iz);
    384. 

  384. 

  385.   // Index of sub-voxel corresponding to ioct1: 8 voxels correspond to 64 sub-voxels.
    386. 

  386.   Standard_Integer ibyte = ibit - (islice << 3);
    387. 

  387. 

  388.   // Value
    389. 

  389.   Standard_Byte value = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[ibyte];
    390. 

  390. 

  391.   return ((value & gbits[ioct1]) ? Standard_True : Standard_False);
    392. 

  392. }
    393. 

  393. 

  394. Standard_Boolean Voxel_ROctBoolDS::Get(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz,
    395. 

  395. 				       const Standard_Integer ioct1, const Standard_Integer ioct2) const
    396. 

  396. {
    397. 

  397.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    398. 

  398.   Standard_Integer islice = ibit >> 3;
    399. 

  399. 

  400.   // If the slice of data is not allocated, it means that its values are 0.
    401. 

  401.   if (!((Voxel_SplitData**)myData)[islice])
    402. 

  402.     return Standard_False;
    403. 

  403. 

  404.   // If the voxel is not split, return the value of the voxel.
    405. 

  405.   if (!((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    406. 

  406.     return Get(ix, iy, iz);
    407. 

  407. 

  408.   // If the split voxel (sub-voxel 1) is not split, return the value of the sub-voxel 1.
    409. 

  409.   if (!((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())
    410. 

  410.     return Get(ix, iy, iz, ioct1);
    411. 

  411. 

  412.   // Index of sub-voxel corresponding to ioct1: 8 voxels correspond to 64 sub-voxels.
    413. 

  413.   Standard_Integer ibyte1 = ibit - (islice << 3);  // index of byte of  8-byte value (sub-voxel 1).
    414. 

  414.   Standard_Integer ibyte2 = (ibyte1 << 3) + ioct1; // index of byte of 64-byte value (sub-voxel 2)
    415. 

  415. 

  416.   // Value
    417. 

  417.   Standard_Byte value = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2];
    418. 

  418. 

  419.   return ((value & gbits[ioct2]) ? Standard_True : Standard_False);
    420. 

  420. }
    421. 

  421. 

  422. Standard_Boolean Voxel_ROctBoolDS::IsSplit(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz) const
    423. 

  423. {
    424. 

  424.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    425. 

  425.   Standard_Integer islice = ibit >> 3;
    426. 

  426. 

  427.   // If the voxel has no value, it is not split.
    428. 

  428.   if (!((Voxel_SplitData**)myData)[islice])
    429. 

  429.     return Standard_False;
    430. 

  430. 

  431.   // Check existence of sub-voxels
    432. 

  432.   if (((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    433. 

  433.     return Standard_True;
    434. 

  434.   return Standard_False;
    435. 

  435. }
    436. 

  436. 

  437. Standard_Integer Voxel_ROctBoolDS::Deepness(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz) const
    438. 

  438. {
    439. 

  439.   Standard_Integer ibit = ix + myNbX * iy + myNbXY * iz;
    440. 

  440.   Standard_Integer islice = ibit >> 3;
    441. 

  441. 

  442.   // If the voxel has no value, it is not split.
    443. 

  443.   if (!((Voxel_SplitData**)myData)[islice])
    444. 

  444.     return 0;
    445. 

  445. 

  446.   // Test deepness.
    447. 

  447.   if (((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    448. 

  448.   {
    449. 

  449.     if (((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())
    450. 

  450.     {
    451. 

  451.       return 2;
    452. 

  452.     }
    453. 

  453.     else
    454. 

  454.     {
    455. 

  455.       return 1;
    456. 

  456.     }
    457. 

  457.   }
    458. 

  458.   return 0;
    459. 

  459. }
    460. 

  460. 

  461. void Voxel_ROctBoolDS::OptimizeMemory()
    462. 

  462. {
    463. 

  463.   // Iterate the array of voxels checking coincidence of values of sub-voxels.
    464. 

  464.   Standard_Integer islice = 0, nb_slices = RealToInt(ceil(myNbXY * myNbZ / 8.0));
    465. 

  465.   for (; islice < nb_slices; islice++)
    466. 

  466.   {
    467. 

  467.     if (!((Voxel_SplitData**)myData)[islice])
    468. 

  468.       continue;
    469. 

  469.     if (((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())
    470. 

  470.     {
    471. 

  471.       Standard_Boolean suppress = Standard_False;
    472. 

  472.       // Second level of sub-voxels
    473. 

  473.       if (((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())
    474. 

  474.       {
    475. 

  475. 	suppress = Standard_False;
    476. 

  476. 	Standard_Byte value1 = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[0];
    477. 

  477. 	if (value1 == 0 || value1 == 255)
    478. 

  478. 	{
    479. 

  479. 	  suppress = Standard_True;
    480. 

  480. 	  for (Standard_Integer ibyte2 = 1; ibyte2 < 64; ibyte2++)
    481. 

  481. 	  {
    482. 

  482. 	    Standard_Byte value2 = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData())->GetValues())[ibyte2];
    483. 

  483. 	    if (value2 != value1)
    484. 

  484. 	    {
    485. 

  485. 	      suppress = Standard_False;
    486. 

  486. 	      break;
    487. 

  487. 	    }
    488. 

  488. 	  }
    489. 

  489. 	}
    490. 

  490. 	if (suppress)
    491. 

  491. 	{
    492. 

  492. 	  SetZeroSplitData((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData());
    493. 

  493. 	  ((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetSplitData() = 0;
    494. 

  494. 	  // Set value to upper level
    495. 

  495. 	  for (Standard_Integer ibyte1 = 0; ibyte1 < 8; ibyte1++)
    496. 

  496. 	  {
    497. 

  497. 	    ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[ibyte1] = value1;
    498. 

  498. 	  }
    499. 

  499. 	}
    500. 

  500. 	else
    501. 

  501. 	{
    502. 

  502. 	  // If we don't suppress sub-sub-voxels, we don't touch sub-voxels.
    503. 

  503. 	  continue;
    504. 

  504. 	}
    505. 

  505.       }
    506. 

  506.       // First level of sub-voxels
    507. 

  507.       suppress = Standard_False;
    508. 

  508.       Standard_Byte value1 = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[0];
    509. 

  509.       if (value1 == 0 || value1 == 255)
    510. 

  510.       {
    511. 

  511. 	suppress = Standard_True;
    512. 

  512. 	for (Standard_Integer ibyte1 = 1; ibyte1 < 8; ibyte1++)
    513. 

  513. 	{
    514. 

  514. 	  Standard_Byte value2 = ((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData())->GetValues())[ibyte1];
    515. 

  515. 	  if (value2 != value1)
    516. 

  516. 	  {
    517. 

  517. 	    suppress = Standard_False;
    518. 

  518. 	    break;
    519. 

  519. 	  }
    520. 

  520. 	}
    521. 

  521.       }
    522. 

  522.       if (suppress)
    523. 

  523.       {
    524. 

  524. 	SetZeroSplitData((Voxel_SplitData*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData());
    525. 

  525. 	((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetSplitData() = 0;
    526. 

  526. 	// Set value to upper level
    527. 

  527. 	*((Standard_Byte*)((Voxel_SplitData*)((Voxel_SplitData**)myData)[islice])->GetValues()) = value1;
    528. 

  528.       }
    529. 

  529.     }
    530. 

  530.   }
    531. 

  531. }
    532. 

  532. 

  533. void Voxel_ROctBoolDS::GetCenter(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz, 
    534. 

  534. 				 const Standard_Integer i,
    535. 

  535. 				 Standard_Real& xc, Standard_Real& yc, Standard_Real& zc) const
    536. 

  536. {
    537. 

  537.   xc = myX + ix * myDX;
    538. 

  538.   yc = myY + iy * myDY;
    539. 

  539.   zc = myZ + iz * myDZ;
    540. 

  540. 

  541.   switch (i)
    542. 

  542.   {
    543. 

  543.     case 0:
    544. 

  544.     {
    545. 

  545.       xc += 0.5 * myHalfDX;
    546. 

  546.       yc += 0.5 * myHalfDY;
    547. 

  547.       zc += 0.5 * myHalfDZ;
    548. 

  548.       break;
    549. 

  549.     }
    550. 

  550.     case 1:
    551. 

  551.     {
    552. 

  552.       xc += 1.5 * myHalfDX;
    553. 

  553.       yc += 0.5 * myHalfDY;
    554. 

  554.       zc += 0.5 * myHalfDZ;
    555. 

  555.       break;
    556. 

  556.     }
    557. 

  557.     case 2:
    558. 

  558.     {
    559. 

  559.       xc += 0.5 * myHalfDX;
    560. 

  560.       yc += 1.5 * myHalfDY;
    561. 

  561.       zc += 0.5 * myHalfDZ;
    562. 

  562.       break;
    563. 

  563.     }
    564. 

  564.     case 3:
    565. 

  565.     {
    566. 

  566.       xc += 1.5 * myHalfDX;
    567. 

  567.       yc += 1.5 * myHalfDY;
    568. 

  568.       zc += 0.5 * myHalfDZ;
    569. 

  569.       break;
    570. 

  570.     }
    571. 

  571.     case 4:
    572. 

  572.     {
    573. 

  573.       xc += 0.5 * myHalfDX;
    574. 

  574.       yc += 0.5 * myHalfDY;
    575. 

  575.       zc += 1.5 * myHalfDZ;
    576. 

  576.       break;
    577. 

  577.     }
    578. 

  578.     case 5:
    579. 

  579.     {
    580. 

  580.       xc += 1.5 * myHalfDX;
    581. 

  581.       yc += 0.5 * myHalfDY;
    582. 

  582.       zc += 1.5 * myHalfDZ;
    583. 

  583.       break;
    584. 

  584.     }
    585. 

  585.     case 6:
    586. 

  586.     {
    587. 

  587.       xc += 0.5 * myHalfDX;
    588. 

  588.       yc += 1.5 * myHalfDY;
    589. 

  589.       zc += 1.5 * myHalfDZ;
    590. 

  590.       break;
    591. 

  591.     }
    592. 

  592.     case 7:
    593. 

  593.     {
    594. 

  594.       xc += 1.5 * myHalfDX;
    595. 

  595.       yc += 1.5 * myHalfDY;
    596. 

  596.       zc += 1.5 * myHalfDZ;
    597. 

  597.       break;
    598. 

  598.     }
    599. 

  599.   }
    600. 

  600. }
    601. 

  601. 

  602. void Voxel_ROctBoolDS::GetCenter(const Standard_Integer ix, const Standard_Integer iy, const Standard_Integer iz, 
    603. 

  603. 				 const Standard_Integer i, const Standard_Integer j,
    604. 

  604. 				 Standard_Real& xc, Standard_Real& yc, Standard_Real& zc) const
    605. 

  605. {
    606. 

  606.   xc = myX + ix * myDX;
    607. 

  607.   yc = myY + iy * myDY;
    608. 

  608.   zc = myZ + iz * myDZ;
    609. 

  609. 

  610.   switch (i)
    611. 

  611.   {
    612. 

  612.     case 0:
    613. 

  613.     {
    614. 

  614.       break;
    615. 

  615.     }
    616. 

  616.     case 1:
    617. 

  617.     {
    618. 

  618.       xc += myHalfDX;
    619. 

  619.       break;
    620. 

  620.     }
    621. 

  621.     case 2:
    622. 

  622.     {
    623. 

  623.       yc += myHalfDY;
    624. 

  624.       break;
    625. 

  625.     }
    626. 

  626.     case 3:
    627. 

  627.     {
    628. 

  628.       xc += myHalfDX;
    629. 

  629.       yc += myHalfDY;
    630. 

  630.       break;
    631. 

  631.     }
    632. 

  632.     case 4:
    633. 

  633.     {
    634. 

  634.       zc += myHalfDZ;
    635. 

  635.       break;
    636. 

  636.     }
    637. 

  637.     case 5:
    638. 

  638.     {
    639. 

  639.       xc += myHalfDX;
    640. 

  640.       zc += myHalfDZ;
    641. 

  641.       break;
    642. 

  642.     }
    643. 

  643.     case 6:
    644. 

  644.     {
    645. 

  645.       yc += myHalfDY;
    646. 

  646.       zc += myHalfDZ;
    647. 

  647.       break;
    648. 

  648.     }
    649. 

  649.     case 7:
    650. 

  650.     {
    651. 

  651.       xc += myHalfDX;
    652. 

  652.       yc += myHalfDY;
    653. 

  653.       zc += myHalfDZ;
    654. 

  654.       break;
    655. 

  655.     }
    656. 

  656.   }
    657. 

  657. 

  658.   switch (j)
    659. 

  659.   {
    660. 

  660.     case 0:
    661. 

  661.     {
    662. 

  662.       xc += 0.25 * myHalfDX;
    663. 

  663.       yc += 0.25 * myHalfDY;
    664. 

  664.       zc += 0.25 * myHalfDZ;
    665. 

  665.       break;
    666. 

  666.     }
    667. 

  667.     case 1:
    668. 

  668.     {
    669. 

  669.       xc += 0.75 * myHalfDX;
    670. 

  670.       yc += 0.25 * myHalfDY;
    671. 

  671.       zc += 0.25 * myHalfDZ;
    672. 

  672.       break;
    673. 

  673.     }
    674. 

  674.     case 2:
    675. 

  675.     {
    676. 

  676.       xc += 0.25 * myHalfDX;
    677. 

  677.       yc += 0.75 * myHalfDY;
    678. 

  678.       zc += 0.25 * myHalfDZ;
    679. 

  679.       break;
    680. 

  680.     }
    681. 

  681.     case 3:
    682. 

  682.     {
    683. 

  683.       xc += 0.75 * myHalfDX;
    684. 

  684.       yc += 0.75 * myHalfDY;
    685. 

  685.       zc += 0.25 * myHalfDZ;
    686. 

  686.       break;
    687. 

  687.     }
    688. 

  688.     case 4:
    689. 

  689.     {
    690. 

  690.       xc += 0.25 * myHalfDX;
    691. 

  691.       yc += 0.25 * myHalfDY;
    692. 

  692.       zc += 0.75 * myHalfDZ;
    693. 

  693.       break;
    694. 

  694.     }
    695. 

  695.     case 5:
    696. 

  696.     {
    697. 

  697.       xc += 0.75 * myHalfDX;
    698. 

  698.       yc += 0.25 * myHalfDY;
    699. 

  699.       zc += 0.75 * myHalfDZ;
    700. 

  700.       break;
    701. 

  701.     }
    702. 

  702.     case 6:
    703. 

  703.     {
    704. 

  704.       xc += 0.25 * myHalfDX;
    705. 

  705.       yc += 0.75 * myHalfDY;
    706. 

  706.       zc += 0.75 * myHalfDZ;
    707. 

  707.       break;
    708. 

  708.     }
    709. 

  709.     case 7:
    710. 

  710.     {
    711. 

  711.       xc += 0.75 * myHalfDX;
    712. 

  712.       yc += 0.75 * myHalfDY;
    713. 

  713.       zc += 0.75 * myHalfDZ;
    714. 

  714.       break;
    715. 

  715.     }
    716. 

  716.   }
    717. 

  717. }
    718.