MarchingCubes.js 34 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210310410510610710810911011111211311411511611711811912012112212312412512612712812913013113213313413513613713813914014114214314414514614714814915015115215315415515615715815916016116216316416516616716816917017117217317417517617717817918018118218318418518618718818919019119219319419519619719819920020120220320420520620720820921021121221321421521621721821922022122222322422522622722822923023123223323423523623723823924024124224324424524624724824925025125225325425525625725825926026126226326426526626726826927027127227327427527627727827928028128228328428528628728828929029129229329429529629729829930030130230330430530630730830931031131231331431531631731831932032132232332432532632732832933033133233333433533633733833934034134234334434534634734834935035135235335435535635735835936036136236336436536636736836937037137237337437537637737837938038138238338438538638738838939039139239339439539639739839940040140240340440540640740840941041141241341441541641741841942042142242342442542642742842943043143243343443543643743843944044144244344444544644744844945045145245345445545645745845946046146246346446546646746846947047147247347447547647747847948048148248348448548648748848949049149249349449549649749849950050150250350450550650750850951051151251351451551651751851952052152252352452552652752852953053153253353453553653753853954054154254354454554654754854955055155255355455555655755855956056156256356456556656756856957057157257357457557657757857958058158258358458558658758858959059159259359459559659759859960060160260360460560660760860961061161261361461561661761861962062162262362462562662762862963063163263363463563663763863964064164264364464564664764864965065165265365465565665765865966066166266366466566666766866967067167267367467567667767867968068168268368468568668768868969069169269369469569669769869970070170270370470570670770870971071171271371471571671771871972072172272372472572672772872973073173273373473573673773873974074174274374474574674774874975075175275375475575675775875976076176276376476576676776876977077177277377477577677777877978078178278378478578678778878979079179279379479579679779879980080180280380480580680780880981081181281381481581681781881982082182282382482582682782882983083183283383483583683783883984084184284384484584684784884985085185285385485585685785885986086186286386486586686786886987087187287387487587687787887988088188288388488588688788888989089189289389489589689789889990090190290390490590690790890991091191291391491591691791891992092192292392492592692792892993093193293393493593693793893994094194294394494594694794894995095195295395495595695795895996096196296396496596696796896997097197297397497597697797897998098198298398498598698798898999099199299399499599699799899910001001100210031004100510061007100810091010101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010511052105310541055105610571058105910601061106210631064106510661067
  1. /**
  2. * @author alteredq / http://alteredqualia.com/
  3. * @author mrdoob / http://mrdoob.com
  4. * Port of http://webglsamples.org/blob/blob.html
  5. */
  6. THREE.MarchingCubes = function ( resolution, material, enableUvs, enableColors ) {
  7. THREE.ImmediateRenderObject.call( this, material );
  8. var scope = this;
  9. // temp buffers used in polygonize
  10. var vlist = new Float32Array( 12 * 3 );
  11. var nlist = new Float32Array( 12 * 3 );
  12. this.enableUvs = enableUvs !== undefined ? enableUvs : false;
  13. this.enableColors = enableColors !== undefined ? enableColors : false;
  14. // functions have to be object properties
  15. // prototype functions kill performance
  16. // (tested and it was 4x slower !!!)
  17. this.init = function ( resolution ) {
  18. this.resolution = resolution;
  19. // parameters
  20. this.isolation = 80.0;
  21. // size of field, 32 is pushing it in Javascript :)
  22. this.size = resolution;
  23. this.size2 = this.size * this.size;
  24. this.size3 = this.size2 * this.size;
  25. this.halfsize = this.size / 2.0;
  26. // deltas
  27. this.delta = 2.0 / this.size;
  28. this.yd = this.size;
  29. this.zd = this.size2;
  30. this.field = new Float32Array( this.size3 );
  31. this.normal_cache = new Float32Array( this.size3 * 3 );
  32. // immediate render mode simulator
  33. this.maxCount = 4096; // TODO: find the fastest size for this buffer
  34. this.count = 0;
  35. this.hasPositions = false;
  36. this.hasNormals = false;
  37. this.hasColors = false;
  38. this.hasUvs = false;
  39. this.positionArray = new Float32Array( this.maxCount * 3 );
  40. this.normalArray = new Float32Array( this.maxCount * 3 );
  41. if ( this.enableUvs ) {
  42. this.uvArray = new Float32Array( this.maxCount * 2 );
  43. }
  44. if ( this.enableColors ) {
  45. this.colorArray = new Float32Array( this.maxCount * 3 );
  46. }
  47. };
  48. ///////////////////////
  49. // Polygonization
  50. ///////////////////////
  51. function lerp( a, b, t ) {
  52. return a + ( b - a ) * t;
  53. }
  54. function VIntX( q, offset, isol, x, y, z, valp1, valp2 ) {
  55. var mu = ( isol - valp1 ) / ( valp2 - valp1 ),
  56. nc = scope.normal_cache;
  57. vlist[ offset + 0 ] = x + mu * scope.delta;
  58. vlist[ offset + 1 ] = y;
  59. vlist[ offset + 2 ] = z;
  60. nlist[ offset + 0 ] = lerp( nc[ q + 0 ], nc[ q + 3 ], mu );
  61. nlist[ offset + 1 ] = lerp( nc[ q + 1 ], nc[ q + 4 ], mu );
  62. nlist[ offset + 2 ] = lerp( nc[ q + 2 ], nc[ q + 5 ], mu );
  63. }
  64. function VIntY( q, offset, isol, x, y, z, valp1, valp2 ) {
  65. var mu = ( isol - valp1 ) / ( valp2 - valp1 ),
  66. nc = scope.normal_cache;
  67. vlist[ offset + 0 ] = x;
  68. vlist[ offset + 1 ] = y + mu * scope.delta;
  69. vlist[ offset + 2 ] = z;
  70. var q2 = q + scope.yd * 3;
  71. nlist[ offset + 0 ] = lerp( nc[ q + 0 ], nc[ q2 + 0 ], mu );
  72. nlist[ offset + 1 ] = lerp( nc[ q + 1 ], nc[ q2 + 1 ], mu );
  73. nlist[ offset + 2 ] = lerp( nc[ q + 2 ], nc[ q2 + 2 ], mu );
  74. }
  75. function VIntZ( q, offset, isol, x, y, z, valp1, valp2 ) {
  76. var mu = ( isol - valp1 ) / ( valp2 - valp1 ),
  77. nc = scope.normal_cache;
  78. vlist[ offset + 0 ] = x;
  79. vlist[ offset + 1 ] = y;
  80. vlist[ offset + 2 ] = z + mu * scope.delta;
  81. var q2 = q + scope.zd * 3;
  82. nlist[ offset + 0 ] = lerp( nc[ q + 0 ], nc[ q2 + 0 ], mu );
  83. nlist[ offset + 1 ] = lerp( nc[ q + 1 ], nc[ q2 + 1 ], mu );
  84. nlist[ offset + 2 ] = lerp( nc[ q + 2 ], nc[ q2 + 2 ], mu );
  85. }
  86. function compNorm( q ) {
  87. var q3 = q * 3;
  88. if ( scope.normal_cache[ q3 ] === 0.0 ) {
  89. scope.normal_cache[ q3 + 0 ] = scope.field[ q - 1 ] - scope.field[ q + 1 ];
  90. scope.normal_cache[ q3 + 1 ] = scope.field[ q - scope.yd ] - scope.field[ q + scope.yd ];
  91. scope.normal_cache[ q3 + 2 ] = scope.field[ q - scope.zd ] - scope.field[ q + scope.zd ];
  92. }
  93. }
  94. // Returns total number of triangles. Fills triangles.
  95. // (this is where most of time is spent - it's inner work of O(n3) loop )
  96. function polygonize( fx, fy, fz, q, isol, renderCallback ) {
  97. // cache indices
  98. var q1 = q + 1,
  99. qy = q + scope.yd,
  100. qz = q + scope.zd,
  101. q1y = q1 + scope.yd,
  102. q1z = q1 + scope.zd,
  103. qyz = q + scope.yd + scope.zd,
  104. q1yz = q1 + scope.yd + scope.zd;
  105. var cubeindex = 0,
  106. field0 = scope.field[ q ],
  107. field1 = scope.field[ q1 ],
  108. field2 = scope.field[ qy ],
  109. field3 = scope.field[ q1y ],
  110. field4 = scope.field[ qz ],
  111. field5 = scope.field[ q1z ],
  112. field6 = scope.field[ qyz ],
  113. field7 = scope.field[ q1yz ];
  114. if ( field0 < isol ) cubeindex |= 1;
  115. if ( field1 < isol ) cubeindex |= 2;
  116. if ( field2 < isol ) cubeindex |= 8;
  117. if ( field3 < isol ) cubeindex |= 4;
  118. if ( field4 < isol ) cubeindex |= 16;
  119. if ( field5 < isol ) cubeindex |= 32;
  120. if ( field6 < isol ) cubeindex |= 128;
  121. if ( field7 < isol ) cubeindex |= 64;
  122. // if cube is entirely in/out of the surface - bail, nothing to draw
  123. var bits = THREE.edgeTable[ cubeindex ];
  124. if ( bits === 0 ) return 0;
  125. var d = scope.delta,
  126. fx2 = fx + d,
  127. fy2 = fy + d,
  128. fz2 = fz + d;
  129. // top of the cube
  130. if ( bits & 1 ) {
  131. compNorm( q );
  132. compNorm( q1 );
  133. VIntX( q * 3, 0, isol, fx, fy, fz, field0, field1 );
  134. }
  135. if ( bits & 2 ) {
  136. compNorm( q1 );
  137. compNorm( q1y );
  138. VIntY( q1 * 3, 3, isol, fx2, fy, fz, field1, field3 );
  139. }
  140. if ( bits & 4 ) {
  141. compNorm( qy );
  142. compNorm( q1y );
  143. VIntX( qy * 3, 6, isol, fx, fy2, fz, field2, field3 );
  144. }
  145. if ( bits & 8 ) {
  146. compNorm( q );
  147. compNorm( qy );
  148. VIntY( q * 3, 9, isol, fx, fy, fz, field0, field2 );
  149. }
  150. // bottom of the cube
  151. if ( bits & 16 ) {
  152. compNorm( qz );
  153. compNorm( q1z );
  154. VIntX( qz * 3, 12, isol, fx, fy, fz2, field4, field5 );
  155. }
  156. if ( bits & 32 ) {
  157. compNorm( q1z );
  158. compNorm( q1yz );
  159. VIntY( q1z * 3, 15, isol, fx2, fy, fz2, field5, field7 );
  160. }
  161. if ( bits & 64 ) {
  162. compNorm( qyz );
  163. compNorm( q1yz );
  164. VIntX( qyz * 3, 18, isol, fx, fy2, fz2, field6, field7 );
  165. }
  166. if ( bits & 128 ) {
  167. compNorm( qz );
  168. compNorm( qyz );
  169. VIntY( qz * 3, 21, isol, fx, fy, fz2, field4, field6 );
  170. }
  171. // vertical lines of the cube
  172. if ( bits & 256 ) {
  173. compNorm( q );
  174. compNorm( qz );
  175. VIntZ( q * 3, 24, isol, fx, fy, fz, field0, field4 );
  176. }
  177. if ( bits & 512 ) {
  178. compNorm( q1 );
  179. compNorm( q1z );
  180. VIntZ( q1 * 3, 27, isol, fx2, fy, fz, field1, field5 );
  181. }
  182. if ( bits & 1024 ) {
  183. compNorm( q1y );
  184. compNorm( q1yz );
  185. VIntZ( q1y * 3, 30, isol, fx2, fy2, fz, field3, field7 );
  186. }
  187. if ( bits & 2048 ) {
  188. compNorm( qy );
  189. compNorm( qyz );
  190. VIntZ( qy * 3, 33, isol, fx, fy2, fz, field2, field6 );
  191. }
  192. cubeindex <<= 4; // re-purpose cubeindex into an offset into triTable
  193. var o1, o2, o3, numtris = 0, i = 0;
  194. // here is where triangles are created
  195. while ( THREE.triTable[ cubeindex + i ] != - 1 ) {
  196. o1 = cubeindex + i;
  197. o2 = o1 + 1;
  198. o3 = o1 + 2;
  199. posnormtriv( vlist, nlist,
  200. 3 * THREE.triTable[ o1 ],
  201. 3 * THREE.triTable[ o2 ],
  202. 3 * THREE.triTable[ o3 ],
  203. renderCallback );
  204. i += 3;
  205. numtris ++;
  206. }
  207. return numtris;
  208. }
  209. /////////////////////////////////////
  210. // Immediate render mode simulator
  211. /////////////////////////////////////
  212. function posnormtriv( pos, norm, o1, o2, o3, renderCallback ) {
  213. var c = scope.count * 3;
  214. // positions
  215. scope.positionArray[ c + 0 ] = pos[ o1 ];
  216. scope.positionArray[ c + 1 ] = pos[ o1 + 1 ];
  217. scope.positionArray[ c + 2 ] = pos[ o1 + 2 ];
  218. scope.positionArray[ c + 3 ] = pos[ o2 ];
  219. scope.positionArray[ c + 4 ] = pos[ o2 + 1 ];
  220. scope.positionArray[ c + 5 ] = pos[ o2 + 2 ];
  221. scope.positionArray[ c + 6 ] = pos[ o3 ];
  222. scope.positionArray[ c + 7 ] = pos[ o3 + 1 ];
  223. scope.positionArray[ c + 8 ] = pos[ o3 + 2 ];
  224. // normals
  225. if ( scope.material.flatShading === true ) {
  226. var nx = ( norm[ o1 + 0 ] + norm[ o2 + 0 ] + norm[ o3 + 0 ] ) / 3;
  227. var ny = ( norm[ o1 + 1 ] + norm[ o2 + 1 ] + norm[ o3 + 1 ] ) / 3;
  228. var nz = ( norm[ o1 + 2 ] + norm[ o2 + 2 ] + norm[ o3 + 2 ] ) / 3;
  229. scope.normalArray[ c + 0 ] = nx;
  230. scope.normalArray[ c + 1 ] = ny;
  231. scope.normalArray[ c + 2 ] = nz;
  232. scope.normalArray[ c + 3 ] = nx;
  233. scope.normalArray[ c + 4 ] = ny;
  234. scope.normalArray[ c + 5 ] = nz;
  235. scope.normalArray[ c + 6 ] = nx;
  236. scope.normalArray[ c + 7 ] = ny;
  237. scope.normalArray[ c + 8 ] = nz;
  238. } else {
  239. scope.normalArray[ c + 0 ] = norm[ o1 + 0 ];
  240. scope.normalArray[ c + 1 ] = norm[ o1 + 1 ];
  241. scope.normalArray[ c + 2 ] = norm[ o1 + 2 ];
  242. scope.normalArray[ c + 3 ] = norm[ o2 + 0 ];
  243. scope.normalArray[ c + 4 ] = norm[ o2 + 1 ];
  244. scope.normalArray[ c + 5 ] = norm[ o2 + 2 ];
  245. scope.normalArray[ c + 6 ] = norm[ o3 + 0 ];
  246. scope.normalArray[ c + 7 ] = norm[ o3 + 1 ];
  247. scope.normalArray[ c + 8 ] = norm[ o3 + 2 ];
  248. }
  249. // uvs
  250. if ( scope.enableUvs ) {
  251. var d = scope.count * 2;
  252. scope.uvArray[ d + 0 ] = pos[ o1 + 0 ];
  253. scope.uvArray[ d + 1 ] = pos[ o1 + 2 ];
  254. scope.uvArray[ d + 2 ] = pos[ o2 + 0 ];
  255. scope.uvArray[ d + 3 ] = pos[ o2 + 2 ];
  256. scope.uvArray[ d + 4 ] = pos[ o3 + 0 ];
  257. scope.uvArray[ d + 5 ] = pos[ o3 + 2 ];
  258. }
  259. // colors
  260. if ( scope.enableColors ) {
  261. scope.colorArray[ c + 0 ] = pos[ o1 + 0 ];
  262. scope.colorArray[ c + 1 ] = pos[ o1 + 1 ];
  263. scope.colorArray[ c + 2 ] = pos[ o1 + 2 ];
  264. scope.colorArray[ c + 3 ] = pos[ o2 + 0 ];
  265. scope.colorArray[ c + 4 ] = pos[ o2 + 1 ];
  266. scope.colorArray[ c + 5 ] = pos[ o2 + 2 ];
  267. scope.colorArray[ c + 6 ] = pos[ o3 + 0 ];
  268. scope.colorArray[ c + 7 ] = pos[ o3 + 1 ];
  269. scope.colorArray[ c + 8 ] = pos[ o3 + 2 ];
  270. }
  271. scope.count += 3;
  272. if ( scope.count >= scope.maxCount - 3 ) {
  273. scope.hasPositions = true;
  274. scope.hasNormals = true;
  275. if ( scope.enableUvs ) {
  276. scope.hasUvs = true;
  277. }
  278. if ( scope.enableColors ) {
  279. scope.hasColors = true;
  280. }
  281. renderCallback( scope );
  282. }
  283. }
  284. this.begin = function () {
  285. this.count = 0;
  286. this.hasPositions = false;
  287. this.hasNormals = false;
  288. this.hasUvs = false;
  289. this.hasColors = false;
  290. };
  291. this.end = function ( renderCallback ) {
  292. if ( this.count === 0 ) return;
  293. for ( var i = this.count * 3; i < this.positionArray.length; i ++ ) {
  294. this.positionArray[ i ] = 0.0;
  295. }
  296. this.hasPositions = true;
  297. this.hasNormals = true;
  298. if ( this.enableUvs && this.material.map ) {
  299. this.hasUvs = true;
  300. }
  301. if ( this.enableColors && this.material.vertexColors !== THREE.NoColors ) {
  302. this.hasColors = true;
  303. }
  304. renderCallback( this );
  305. };
  306. /////////////////////////////////////
  307. // Metaballs
  308. /////////////////////////////////////
  309. // Adds a reciprocal ball (nice and blobby) that, to be fast, fades to zero after
  310. // a fixed distance, determined by strength and subtract.
  311. this.addBall = function ( ballx, bally, ballz, strength, subtract ) {
  312. var sign = Math.sign( strength );
  313. strength = Math.abs( strength );
  314. // Let's solve the equation to find the radius:
  315. // 1.0 / (0.000001 + radius^2) * strength - subtract = 0
  316. // strength / (radius^2) = subtract
  317. // strength = subtract * radius^2
  318. // radius^2 = strength / subtract
  319. // radius = sqrt(strength / subtract)
  320. var radius = this.size * Math.sqrt( strength / subtract ),
  321. zs = ballz * this.size,
  322. ys = bally * this.size,
  323. xs = ballx * this.size;
  324. var min_z = Math.floor( zs - radius ); if ( min_z < 1 ) min_z = 1;
  325. var max_z = Math.floor( zs + radius ); if ( max_z > this.size - 1 ) max_z = this.size - 1;
  326. var min_y = Math.floor( ys - radius ); if ( min_y < 1 ) min_y = 1;
  327. var max_y = Math.floor( ys + radius ); if ( max_y > this.size - 1 ) max_y = this.size - 1;
  328. var min_x = Math.floor( xs - radius ); if ( min_x < 1 ) min_x = 1;
  329. var max_x = Math.floor( xs + radius ); if ( max_x > this.size - 1 ) max_x = this.size - 1;
  330. // Don't polygonize in the outer layer because normals aren't
  331. // well-defined there.
  332. var x, y, z, y_offset, z_offset, fx, fy, fz, fz2, fy2, val;
  333. for ( z = min_z; z < max_z; z ++ ) {
  334. z_offset = this.size2 * z;
  335. fz = z / this.size - ballz;
  336. fz2 = fz * fz;
  337. for ( y = min_y; y < max_y; y ++ ) {
  338. y_offset = z_offset + this.size * y;
  339. fy = y / this.size - bally;
  340. fy2 = fy * fy;
  341. for ( x = min_x; x < max_x; x ++ ) {
  342. fx = x / this.size - ballx;
  343. val = strength / ( 0.000001 + fx * fx + fy2 + fz2 ) - subtract;
  344. if ( val > 0.0 ) this.field[ y_offset + x ] += val * sign;
  345. }
  346. }
  347. }
  348. };
  349. this.addPlaneX = function ( strength, subtract ) {
  350. var x, y, z, xx, val, xdiv, cxy,
  351. // cache attribute lookups
  352. size = this.size,
  353. yd = this.yd,
  354. zd = this.zd,
  355. field = this.field,
  356. dist = size * Math.sqrt( strength / subtract );
  357. if ( dist > size ) dist = size;
  358. for ( x = 0; x < dist; x ++ ) {
  359. xdiv = x / size;
  360. xx = xdiv * xdiv;
  361. val = strength / ( 0.0001 + xx ) - subtract;
  362. if ( val > 0.0 ) {
  363. for ( y = 0; y < size; y ++ ) {
  364. cxy = x + y * yd;
  365. for ( z = 0; z < size; z ++ ) {
  366. field[ zd * z + cxy ] += val;
  367. }
  368. }
  369. }
  370. }
  371. };
  372. this.addPlaneY = function ( strength, subtract ) {
  373. var x, y, z, yy, val, ydiv, cy, cxy,
  374. // cache attribute lookups
  375. size = this.size,
  376. yd = this.yd,
  377. zd = this.zd,
  378. field = this.field,
  379. dist = size * Math.sqrt( strength / subtract );
  380. if ( dist > size ) dist = size;
  381. for ( y = 0; y < dist; y ++ ) {
  382. ydiv = y / size;
  383. yy = ydiv * ydiv;
  384. val = strength / ( 0.0001 + yy ) - subtract;
  385. if ( val > 0.0 ) {
  386. cy = y * yd;
  387. for ( x = 0; x < size; x ++ ) {
  388. cxy = cy + x;
  389. for ( z = 0; z < size; z ++ )
  390. field[ zd * z + cxy ] += val;
  391. }
  392. }
  393. }
  394. };
  395. this.addPlaneZ = function ( strength, subtract ) {
  396. var x, y, z, zz, val, zdiv, cz, cyz,
  397. // cache attribute lookups
  398. size = this.size,
  399. yd = this.yd,
  400. zd = this.zd,
  401. field = this.field,
  402. dist = size * Math.sqrt( strength / subtract );
  403. if ( dist > size ) dist = size;
  404. for ( z = 0; z < dist; z ++ ) {
  405. zdiv = z / size;
  406. zz = zdiv * zdiv;
  407. val = strength / ( 0.0001 + zz ) - subtract;
  408. if ( val > 0.0 ) {
  409. cz = zd * z;
  410. for ( y = 0; y < size; y ++ ) {
  411. cyz = cz + y * yd;
  412. for ( x = 0; x < size; x ++ )
  413. field[ cyz + x ] += val;
  414. }
  415. }
  416. }
  417. };
  418. /////////////////////////////////////
  419. // Updates
  420. /////////////////////////////////////
  421. this.reset = function () {
  422. var i;
  423. // wipe the normal cache
  424. for ( i = 0; i < this.size3; i ++ ) {
  425. this.normal_cache[ i * 3 ] = 0.0;
  426. this.field[ i ] = 0.0;
  427. }
  428. };
  429. this.render = function ( renderCallback ) {
  430. this.begin();
  431. // Triangulate. Yeah, this is slow.
  432. var smin2 = this.size - 2;
  433. for ( var z = 1; z < smin2; z ++ ) {
  434. var z_offset = this.size2 * z;
  435. var fz = ( z - this.halfsize ) / this.halfsize; //+ 1
  436. for ( var y = 1; y < smin2; y ++ ) {
  437. var y_offset = z_offset + this.size * y;
  438. var fy = ( y - this.halfsize ) / this.halfsize; //+ 1
  439. for ( var x = 1; x < smin2; x ++ ) {
  440. var fx = ( x - this.halfsize ) / this.halfsize; //+ 1
  441. var q = y_offset + x;
  442. polygonize( fx, fy, fz, q, this.isolation, renderCallback );
  443. }
  444. }
  445. }
  446. this.end( renderCallback );
  447. };
  448. this.generateGeometry = function () {
  449. console.warn( 'THREE.MarchingCubes: generateGeometry() now returns THREE.BufferGeometry' );
  450. return this.generateBufferGeometry();
  451. };
  452. function concatenate( a, b, length ) {
  453. var result = new Float32Array( a.length + length );
  454. result.set( a, 0 );
  455. result.set( b.slice( 0, length ), a.length );
  456. return result;
  457. }
  458. this.generateBufferGeometry = function () {
  459. var geo = new THREE.BufferGeometry();
  460. var posArray = new Float32Array();
  461. var normArray = new Float32Array();
  462. var colorArray = new Float32Array();
  463. var uvArray = new Float32Array();
  464. var scope = this;
  465. var geo_callback = function ( object ) {
  466. if ( scope.hasPositions ) posArray = concatenate( posArray, object.positionArray, object.count * 3 );
  467. if ( scope.hasNormals ) normArray = concatenate( normArray, object.normalArray, object.count * 3 );
  468. if ( scope.hasColors ) colorArray = concatenate( colorArray, object.colorArray, object.count * 3 );
  469. if ( scope.hasUvs ) uvArray = concatenate( uvArray, object.uvArray, object.count * 2 );
  470. object.count = 0;
  471. };
  472. this.render( geo_callback );
  473. if ( this.hasPositions ) geo.addAttribute( 'position', new THREE.BufferAttribute( posArray, 3 ) );
  474. if ( this.hasNormals ) geo.addAttribute( 'normal', new THREE.BufferAttribute( normArray, 3 ) );
  475. if ( this.hasColors ) geo.addAttribute( 'color', new THREE.BufferAttribute( colorArray, 3 ) );
  476. if ( this.hasUvs ) geo.addAttribute( 'uv', new THREE.BufferAttribute( uvArray, 2 ) );
  477. return geo;
  478. };
  479. this.init( resolution );
  480. };
  481. THREE.MarchingCubes.prototype = Object.create( THREE.ImmediateRenderObject.prototype );
  482. THREE.MarchingCubes.prototype.constructor = THREE.MarchingCubes;
  483. /////////////////////////////////////
  484. // Marching cubes lookup tables
  485. /////////////////////////////////////
  486. // These tables are straight from Paul Bourke's page:
  487. // http://local.wasp.uwa.edu.au/~pbourke/geometry/polygonise/
  488. // who in turn got them from Cory Gene Bloyd.
  489. THREE.edgeTable = new Int32Array( [
  490. 0x0, 0x109, 0x203, 0x30a, 0x406, 0x50f, 0x605, 0x70c,
  491. 0x80c, 0x905, 0xa0f, 0xb06, 0xc0a, 0xd03, 0xe09, 0xf00,
  492. 0x190, 0x99, 0x393, 0x29a, 0x596, 0x49f, 0x795, 0x69c,
  493. 0x99c, 0x895, 0xb9f, 0xa96, 0xd9a, 0xc93, 0xf99, 0xe90,
  494. 0x230, 0x339, 0x33, 0x13a, 0x636, 0x73f, 0x435, 0x53c,
  495. 0xa3c, 0xb35, 0x83f, 0x936, 0xe3a, 0xf33, 0xc39, 0xd30,
  496. 0x3a0, 0x2a9, 0x1a3, 0xaa, 0x7a6, 0x6af, 0x5a5, 0x4ac,
  497. 0xbac, 0xaa5, 0x9af, 0x8a6, 0xfaa, 0xea3, 0xda9, 0xca0,
  498. 0x460, 0x569, 0x663, 0x76a, 0x66, 0x16f, 0x265, 0x36c,
  499. 0xc6c, 0xd65, 0xe6f, 0xf66, 0x86a, 0x963, 0xa69, 0xb60,
  500. 0x5f0, 0x4f9, 0x7f3, 0x6fa, 0x1f6, 0xff, 0x3f5, 0x2fc,
  501. 0xdfc, 0xcf5, 0xfff, 0xef6, 0x9fa, 0x8f3, 0xbf9, 0xaf0,
  502. 0x650, 0x759, 0x453, 0x55a, 0x256, 0x35f, 0x55, 0x15c,
  503. 0xe5c, 0xf55, 0xc5f, 0xd56, 0xa5a, 0xb53, 0x859, 0x950,
  504. 0x7c0, 0x6c9, 0x5c3, 0x4ca, 0x3c6, 0x2cf, 0x1c5, 0xcc,
  505. 0xfcc, 0xec5, 0xdcf, 0xcc6, 0xbca, 0xac3, 0x9c9, 0x8c0,
  506. 0x8c0, 0x9c9, 0xac3, 0xbca, 0xcc6, 0xdcf, 0xec5, 0xfcc,
  507. 0xcc, 0x1c5, 0x2cf, 0x3c6, 0x4ca, 0x5c3, 0x6c9, 0x7c0,
  508. 0x950, 0x859, 0xb53, 0xa5a, 0xd56, 0xc5f, 0xf55, 0xe5c,
  509. 0x15c, 0x55, 0x35f, 0x256, 0x55a, 0x453, 0x759, 0x650,
  510. 0xaf0, 0xbf9, 0x8f3, 0x9fa, 0xef6, 0xfff, 0xcf5, 0xdfc,
  511. 0x2fc, 0x3f5, 0xff, 0x1f6, 0x6fa, 0x7f3, 0x4f9, 0x5f0,
  512. 0xb60, 0xa69, 0x963, 0x86a, 0xf66, 0xe6f, 0xd65, 0xc6c,
  513. 0x36c, 0x265, 0x16f, 0x66, 0x76a, 0x663, 0x569, 0x460,
  514. 0xca0, 0xda9, 0xea3, 0xfaa, 0x8a6, 0x9af, 0xaa5, 0xbac,
  515. 0x4ac, 0x5a5, 0x6af, 0x7a6, 0xaa, 0x1a3, 0x2a9, 0x3a0,
  516. 0xd30, 0xc39, 0xf33, 0xe3a, 0x936, 0x83f, 0xb35, 0xa3c,
  517. 0x53c, 0x435, 0x73f, 0x636, 0x13a, 0x33, 0x339, 0x230,
  518. 0xe90, 0xf99, 0xc93, 0xd9a, 0xa96, 0xb9f, 0x895, 0x99c,
  519. 0x69c, 0x795, 0x49f, 0x596, 0x29a, 0x393, 0x99, 0x190,
  520. 0xf00, 0xe09, 0xd03, 0xc0a, 0xb06, 0xa0f, 0x905, 0x80c,
  521. 0x70c, 0x605, 0x50f, 0x406, 0x30a, 0x203, 0x109, 0x0 ] );
  522. THREE.triTable = new Int32Array( [
  523. - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  524. 0, 8, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  525. 0, 1, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  526. 1, 8, 3, 9, 8, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  527. 1, 2, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  528. 0, 8, 3, 1, 2, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  529. 9, 2, 10, 0, 2, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  530. 2, 8, 3, 2, 10, 8, 10, 9, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  531. 3, 11, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  532. 0, 11, 2, 8, 11, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  533. 1, 9, 0, 2, 3, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  534. 1, 11, 2, 1, 9, 11, 9, 8, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  535. 3, 10, 1, 11, 10, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  536. 0, 10, 1, 0, 8, 10, 8, 11, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  537. 3, 9, 0, 3, 11, 9, 11, 10, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  538. 9, 8, 10, 10, 8, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  539. 4, 7, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  540. 4, 3, 0, 7, 3, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  541. 0, 1, 9, 8, 4, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  542. 4, 1, 9, 4, 7, 1, 7, 3, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  543. 1, 2, 10, 8, 4, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  544. 3, 4, 7, 3, 0, 4, 1, 2, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  545. 9, 2, 10, 9, 0, 2, 8, 4, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  546. 2, 10, 9, 2, 9, 7, 2, 7, 3, 7, 9, 4, - 1, - 1, - 1, - 1,
  547. 8, 4, 7, 3, 11, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  548. 11, 4, 7, 11, 2, 4, 2, 0, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  549. 9, 0, 1, 8, 4, 7, 2, 3, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  550. 4, 7, 11, 9, 4, 11, 9, 11, 2, 9, 2, 1, - 1, - 1, - 1, - 1,
  551. 3, 10, 1, 3, 11, 10, 7, 8, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  552. 1, 11, 10, 1, 4, 11, 1, 0, 4, 7, 11, 4, - 1, - 1, - 1, - 1,
  553. 4, 7, 8, 9, 0, 11, 9, 11, 10, 11, 0, 3, - 1, - 1, - 1, - 1,
  554. 4, 7, 11, 4, 11, 9, 9, 11, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  555. 9, 5, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  556. 9, 5, 4, 0, 8, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  557. 0, 5, 4, 1, 5, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  558. 8, 5, 4, 8, 3, 5, 3, 1, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  559. 1, 2, 10, 9, 5, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  560. 3, 0, 8, 1, 2, 10, 4, 9, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  561. 5, 2, 10, 5, 4, 2, 4, 0, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  562. 2, 10, 5, 3, 2, 5, 3, 5, 4, 3, 4, 8, - 1, - 1, - 1, - 1,
  563. 9, 5, 4, 2, 3, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  564. 0, 11, 2, 0, 8, 11, 4, 9, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  565. 0, 5, 4, 0, 1, 5, 2, 3, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  566. 2, 1, 5, 2, 5, 8, 2, 8, 11, 4, 8, 5, - 1, - 1, - 1, - 1,
  567. 10, 3, 11, 10, 1, 3, 9, 5, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  568. 4, 9, 5, 0, 8, 1, 8, 10, 1, 8, 11, 10, - 1, - 1, - 1, - 1,
  569. 5, 4, 0, 5, 0, 11, 5, 11, 10, 11, 0, 3, - 1, - 1, - 1, - 1,
  570. 5, 4, 8, 5, 8, 10, 10, 8, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  571. 9, 7, 8, 5, 7, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  572. 9, 3, 0, 9, 5, 3, 5, 7, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  573. 0, 7, 8, 0, 1, 7, 1, 5, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  574. 1, 5, 3, 3, 5, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  575. 9, 7, 8, 9, 5, 7, 10, 1, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  576. 10, 1, 2, 9, 5, 0, 5, 3, 0, 5, 7, 3, - 1, - 1, - 1, - 1,
  577. 8, 0, 2, 8, 2, 5, 8, 5, 7, 10, 5, 2, - 1, - 1, - 1, - 1,
  578. 2, 10, 5, 2, 5, 3, 3, 5, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  579. 7, 9, 5, 7, 8, 9, 3, 11, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  580. 9, 5, 7, 9, 7, 2, 9, 2, 0, 2, 7, 11, - 1, - 1, - 1, - 1,
  581. 2, 3, 11, 0, 1, 8, 1, 7, 8, 1, 5, 7, - 1, - 1, - 1, - 1,
  582. 11, 2, 1, 11, 1, 7, 7, 1, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  583. 9, 5, 8, 8, 5, 7, 10, 1, 3, 10, 3, 11, - 1, - 1, - 1, - 1,
  584. 5, 7, 0, 5, 0, 9, 7, 11, 0, 1, 0, 10, 11, 10, 0, - 1,
  585. 11, 10, 0, 11, 0, 3, 10, 5, 0, 8, 0, 7, 5, 7, 0, - 1,
  586. 11, 10, 5, 7, 11, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  587. 10, 6, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  588. 0, 8, 3, 5, 10, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  589. 9, 0, 1, 5, 10, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  590. 1, 8, 3, 1, 9, 8, 5, 10, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  591. 1, 6, 5, 2, 6, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  592. 1, 6, 5, 1, 2, 6, 3, 0, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  593. 9, 6, 5, 9, 0, 6, 0, 2, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  594. 5, 9, 8, 5, 8, 2, 5, 2, 6, 3, 2, 8, - 1, - 1, - 1, - 1,
  595. 2, 3, 11, 10, 6, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  596. 11, 0, 8, 11, 2, 0, 10, 6, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  597. 0, 1, 9, 2, 3, 11, 5, 10, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  598. 5, 10, 6, 1, 9, 2, 9, 11, 2, 9, 8, 11, - 1, - 1, - 1, - 1,
  599. 6, 3, 11, 6, 5, 3, 5, 1, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  600. 0, 8, 11, 0, 11, 5, 0, 5, 1, 5, 11, 6, - 1, - 1, - 1, - 1,
  601. 3, 11, 6, 0, 3, 6, 0, 6, 5, 0, 5, 9, - 1, - 1, - 1, - 1,
  602. 6, 5, 9, 6, 9, 11, 11, 9, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  603. 5, 10, 6, 4, 7, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  604. 4, 3, 0, 4, 7, 3, 6, 5, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  605. 1, 9, 0, 5, 10, 6, 8, 4, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  606. 10, 6, 5, 1, 9, 7, 1, 7, 3, 7, 9, 4, - 1, - 1, - 1, - 1,
  607. 6, 1, 2, 6, 5, 1, 4, 7, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  608. 1, 2, 5, 5, 2, 6, 3, 0, 4, 3, 4, 7, - 1, - 1, - 1, - 1,
  609. 8, 4, 7, 9, 0, 5, 0, 6, 5, 0, 2, 6, - 1, - 1, - 1, - 1,
  610. 7, 3, 9, 7, 9, 4, 3, 2, 9, 5, 9, 6, 2, 6, 9, - 1,
  611. 3, 11, 2, 7, 8, 4, 10, 6, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  612. 5, 10, 6, 4, 7, 2, 4, 2, 0, 2, 7, 11, - 1, - 1, - 1, - 1,
  613. 0, 1, 9, 4, 7, 8, 2, 3, 11, 5, 10, 6, - 1, - 1, - 1, - 1,
  614. 9, 2, 1, 9, 11, 2, 9, 4, 11, 7, 11, 4, 5, 10, 6, - 1,
  615. 8, 4, 7, 3, 11, 5, 3, 5, 1, 5, 11, 6, - 1, - 1, - 1, - 1,
  616. 5, 1, 11, 5, 11, 6, 1, 0, 11, 7, 11, 4, 0, 4, 11, - 1,
  617. 0, 5, 9, 0, 6, 5, 0, 3, 6, 11, 6, 3, 8, 4, 7, - 1,
  618. 6, 5, 9, 6, 9, 11, 4, 7, 9, 7, 11, 9, - 1, - 1, - 1, - 1,
  619. 10, 4, 9, 6, 4, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  620. 4, 10, 6, 4, 9, 10, 0, 8, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  621. 10, 0, 1, 10, 6, 0, 6, 4, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  622. 8, 3, 1, 8, 1, 6, 8, 6, 4, 6, 1, 10, - 1, - 1, - 1, - 1,
  623. 1, 4, 9, 1, 2, 4, 2, 6, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  624. 3, 0, 8, 1, 2, 9, 2, 4, 9, 2, 6, 4, - 1, - 1, - 1, - 1,
  625. 0, 2, 4, 4, 2, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  626. 8, 3, 2, 8, 2, 4, 4, 2, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  627. 10, 4, 9, 10, 6, 4, 11, 2, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  628. 0, 8, 2, 2, 8, 11, 4, 9, 10, 4, 10, 6, - 1, - 1, - 1, - 1,
  629. 3, 11, 2, 0, 1, 6, 0, 6, 4, 6, 1, 10, - 1, - 1, - 1, - 1,
  630. 6, 4, 1, 6, 1, 10, 4, 8, 1, 2, 1, 11, 8, 11, 1, - 1,
  631. 9, 6, 4, 9, 3, 6, 9, 1, 3, 11, 6, 3, - 1, - 1, - 1, - 1,
  632. 8, 11, 1, 8, 1, 0, 11, 6, 1, 9, 1, 4, 6, 4, 1, - 1,
  633. 3, 11, 6, 3, 6, 0, 0, 6, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  634. 6, 4, 8, 11, 6, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  635. 7, 10, 6, 7, 8, 10, 8, 9, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  636. 0, 7, 3, 0, 10, 7, 0, 9, 10, 6, 7, 10, - 1, - 1, - 1, - 1,
  637. 10, 6, 7, 1, 10, 7, 1, 7, 8, 1, 8, 0, - 1, - 1, - 1, - 1,
  638. 10, 6, 7, 10, 7, 1, 1, 7, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  639. 1, 2, 6, 1, 6, 8, 1, 8, 9, 8, 6, 7, - 1, - 1, - 1, - 1,
  640. 2, 6, 9, 2, 9, 1, 6, 7, 9, 0, 9, 3, 7, 3, 9, - 1,
  641. 7, 8, 0, 7, 0, 6, 6, 0, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  642. 7, 3, 2, 6, 7, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  643. 2, 3, 11, 10, 6, 8, 10, 8, 9, 8, 6, 7, - 1, - 1, - 1, - 1,
  644. 2, 0, 7, 2, 7, 11, 0, 9, 7, 6, 7, 10, 9, 10, 7, - 1,
  645. 1, 8, 0, 1, 7, 8, 1, 10, 7, 6, 7, 10, 2, 3, 11, - 1,
  646. 11, 2, 1, 11, 1, 7, 10, 6, 1, 6, 7, 1, - 1, - 1, - 1, - 1,
  647. 8, 9, 6, 8, 6, 7, 9, 1, 6, 11, 6, 3, 1, 3, 6, - 1,
  648. 0, 9, 1, 11, 6, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  649. 7, 8, 0, 7, 0, 6, 3, 11, 0, 11, 6, 0, - 1, - 1, - 1, - 1,
  650. 7, 11, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  651. 7, 6, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  652. 3, 0, 8, 11, 7, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  653. 0, 1, 9, 11, 7, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  654. 8, 1, 9, 8, 3, 1, 11, 7, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  655. 10, 1, 2, 6, 11, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  656. 1, 2, 10, 3, 0, 8, 6, 11, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  657. 2, 9, 0, 2, 10, 9, 6, 11, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  658. 6, 11, 7, 2, 10, 3, 10, 8, 3, 10, 9, 8, - 1, - 1, - 1, - 1,
  659. 7, 2, 3, 6, 2, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  660. 7, 0, 8, 7, 6, 0, 6, 2, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  661. 2, 7, 6, 2, 3, 7, 0, 1, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  662. 1, 6, 2, 1, 8, 6, 1, 9, 8, 8, 7, 6, - 1, - 1, - 1, - 1,
  663. 10, 7, 6, 10, 1, 7, 1, 3, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  664. 10, 7, 6, 1, 7, 10, 1, 8, 7, 1, 0, 8, - 1, - 1, - 1, - 1,
  665. 0, 3, 7, 0, 7, 10, 0, 10, 9, 6, 10, 7, - 1, - 1, - 1, - 1,
  666. 7, 6, 10, 7, 10, 8, 8, 10, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  667. 6, 8, 4, 11, 8, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  668. 3, 6, 11, 3, 0, 6, 0, 4, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  669. 8, 6, 11, 8, 4, 6, 9, 0, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  670. 9, 4, 6, 9, 6, 3, 9, 3, 1, 11, 3, 6, - 1, - 1, - 1, - 1,
  671. 6, 8, 4, 6, 11, 8, 2, 10, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  672. 1, 2, 10, 3, 0, 11, 0, 6, 11, 0, 4, 6, - 1, - 1, - 1, - 1,
  673. 4, 11, 8, 4, 6, 11, 0, 2, 9, 2, 10, 9, - 1, - 1, - 1, - 1,
  674. 10, 9, 3, 10, 3, 2, 9, 4, 3, 11, 3, 6, 4, 6, 3, - 1,
  675. 8, 2, 3, 8, 4, 2, 4, 6, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  676. 0, 4, 2, 4, 6, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  677. 1, 9, 0, 2, 3, 4, 2, 4, 6, 4, 3, 8, - 1, - 1, - 1, - 1,
  678. 1, 9, 4, 1, 4, 2, 2, 4, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  679. 8, 1, 3, 8, 6, 1, 8, 4, 6, 6, 10, 1, - 1, - 1, - 1, - 1,
  680. 10, 1, 0, 10, 0, 6, 6, 0, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  681. 4, 6, 3, 4, 3, 8, 6, 10, 3, 0, 3, 9, 10, 9, 3, - 1,
  682. 10, 9, 4, 6, 10, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  683. 4, 9, 5, 7, 6, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  684. 0, 8, 3, 4, 9, 5, 11, 7, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  685. 5, 0, 1, 5, 4, 0, 7, 6, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  686. 11, 7, 6, 8, 3, 4, 3, 5, 4, 3, 1, 5, - 1, - 1, - 1, - 1,
  687. 9, 5, 4, 10, 1, 2, 7, 6, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  688. 6, 11, 7, 1, 2, 10, 0, 8, 3, 4, 9, 5, - 1, - 1, - 1, - 1,
  689. 7, 6, 11, 5, 4, 10, 4, 2, 10, 4, 0, 2, - 1, - 1, - 1, - 1,
  690. 3, 4, 8, 3, 5, 4, 3, 2, 5, 10, 5, 2, 11, 7, 6, - 1,
  691. 7, 2, 3, 7, 6, 2, 5, 4, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  692. 9, 5, 4, 0, 8, 6, 0, 6, 2, 6, 8, 7, - 1, - 1, - 1, - 1,
  693. 3, 6, 2, 3, 7, 6, 1, 5, 0, 5, 4, 0, - 1, - 1, - 1, - 1,
  694. 6, 2, 8, 6, 8, 7, 2, 1, 8, 4, 8, 5, 1, 5, 8, - 1,
  695. 9, 5, 4, 10, 1, 6, 1, 7, 6, 1, 3, 7, - 1, - 1, - 1, - 1,
  696. 1, 6, 10, 1, 7, 6, 1, 0, 7, 8, 7, 0, 9, 5, 4, - 1,
  697. 4, 0, 10, 4, 10, 5, 0, 3, 10, 6, 10, 7, 3, 7, 10, - 1,
  698. 7, 6, 10, 7, 10, 8, 5, 4, 10, 4, 8, 10, - 1, - 1, - 1, - 1,
  699. 6, 9, 5, 6, 11, 9, 11, 8, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  700. 3, 6, 11, 0, 6, 3, 0, 5, 6, 0, 9, 5, - 1, - 1, - 1, - 1,
  701. 0, 11, 8, 0, 5, 11, 0, 1, 5, 5, 6, 11, - 1, - 1, - 1, - 1,
  702. 6, 11, 3, 6, 3, 5, 5, 3, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  703. 1, 2, 10, 9, 5, 11, 9, 11, 8, 11, 5, 6, - 1, - 1, - 1, - 1,
  704. 0, 11, 3, 0, 6, 11, 0, 9, 6, 5, 6, 9, 1, 2, 10, - 1,
  705. 11, 8, 5, 11, 5, 6, 8, 0, 5, 10, 5, 2, 0, 2, 5, - 1,
  706. 6, 11, 3, 6, 3, 5, 2, 10, 3, 10, 5, 3, - 1, - 1, - 1, - 1,
  707. 5, 8, 9, 5, 2, 8, 5, 6, 2, 3, 8, 2, - 1, - 1, - 1, - 1,
  708. 9, 5, 6, 9, 6, 0, 0, 6, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  709. 1, 5, 8, 1, 8, 0, 5, 6, 8, 3, 8, 2, 6, 2, 8, - 1,
  710. 1, 5, 6, 2, 1, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  711. 1, 3, 6, 1, 6, 10, 3, 8, 6, 5, 6, 9, 8, 9, 6, - 1,
  712. 10, 1, 0, 10, 0, 6, 9, 5, 0, 5, 6, 0, - 1, - 1, - 1, - 1,
  713. 0, 3, 8, 5, 6, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  714. 10, 5, 6, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  715. 11, 5, 10, 7, 5, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  716. 11, 5, 10, 11, 7, 5, 8, 3, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  717. 5, 11, 7, 5, 10, 11, 1, 9, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  718. 10, 7, 5, 10, 11, 7, 9, 8, 1, 8, 3, 1, - 1, - 1, - 1, - 1,
  719. 11, 1, 2, 11, 7, 1, 7, 5, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  720. 0, 8, 3, 1, 2, 7, 1, 7, 5, 7, 2, 11, - 1, - 1, - 1, - 1,
  721. 9, 7, 5, 9, 2, 7, 9, 0, 2, 2, 11, 7, - 1, - 1, - 1, - 1,
  722. 7, 5, 2, 7, 2, 11, 5, 9, 2, 3, 2, 8, 9, 8, 2, - 1,
  723. 2, 5, 10, 2, 3, 5, 3, 7, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  724. 8, 2, 0, 8, 5, 2, 8, 7, 5, 10, 2, 5, - 1, - 1, - 1, - 1,
  725. 9, 0, 1, 5, 10, 3, 5, 3, 7, 3, 10, 2, - 1, - 1, - 1, - 1,
  726. 9, 8, 2, 9, 2, 1, 8, 7, 2, 10, 2, 5, 7, 5, 2, - 1,
  727. 1, 3, 5, 3, 7, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  728. 0, 8, 7, 0, 7, 1, 1, 7, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  729. 9, 0, 3, 9, 3, 5, 5, 3, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  730. 9, 8, 7, 5, 9, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  731. 5, 8, 4, 5, 10, 8, 10, 11, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  732. 5, 0, 4, 5, 11, 0, 5, 10, 11, 11, 3, 0, - 1, - 1, - 1, - 1,
  733. 0, 1, 9, 8, 4, 10, 8, 10, 11, 10, 4, 5, - 1, - 1, - 1, - 1,
  734. 10, 11, 4, 10, 4, 5, 11, 3, 4, 9, 4, 1, 3, 1, 4, - 1,
  735. 2, 5, 1, 2, 8, 5, 2, 11, 8, 4, 5, 8, - 1, - 1, - 1, - 1,
  736. 0, 4, 11, 0, 11, 3, 4, 5, 11, 2, 11, 1, 5, 1, 11, - 1,
  737. 0, 2, 5, 0, 5, 9, 2, 11, 5, 4, 5, 8, 11, 8, 5, - 1,
  738. 9, 4, 5, 2, 11, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  739. 2, 5, 10, 3, 5, 2, 3, 4, 5, 3, 8, 4, - 1, - 1, - 1, - 1,
  740. 5, 10, 2, 5, 2, 4, 4, 2, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  741. 3, 10, 2, 3, 5, 10, 3, 8, 5, 4, 5, 8, 0, 1, 9, - 1,
  742. 5, 10, 2, 5, 2, 4, 1, 9, 2, 9, 4, 2, - 1, - 1, - 1, - 1,
  743. 8, 4, 5, 8, 5, 3, 3, 5, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  744. 0, 4, 5, 1, 0, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  745. 8, 4, 5, 8, 5, 3, 9, 0, 5, 0, 3, 5, - 1, - 1, - 1, - 1,
  746. 9, 4, 5, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  747. 4, 11, 7, 4, 9, 11, 9, 10, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  748. 0, 8, 3, 4, 9, 7, 9, 11, 7, 9, 10, 11, - 1, - 1, - 1, - 1,
  749. 1, 10, 11, 1, 11, 4, 1, 4, 0, 7, 4, 11, - 1, - 1, - 1, - 1,
  750. 3, 1, 4, 3, 4, 8, 1, 10, 4, 7, 4, 11, 10, 11, 4, - 1,
  751. 4, 11, 7, 9, 11, 4, 9, 2, 11, 9, 1, 2, - 1, - 1, - 1, - 1,
  752. 9, 7, 4, 9, 11, 7, 9, 1, 11, 2, 11, 1, 0, 8, 3, - 1,
  753. 11, 7, 4, 11, 4, 2, 2, 4, 0, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  754. 11, 7, 4, 11, 4, 2, 8, 3, 4, 3, 2, 4, - 1, - 1, - 1, - 1,
  755. 2, 9, 10, 2, 7, 9, 2, 3, 7, 7, 4, 9, - 1, - 1, - 1, - 1,
  756. 9, 10, 7, 9, 7, 4, 10, 2, 7, 8, 7, 0, 2, 0, 7, - 1,
  757. 3, 7, 10, 3, 10, 2, 7, 4, 10, 1, 10, 0, 4, 0, 10, - 1,
  758. 1, 10, 2, 8, 7, 4, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  759. 4, 9, 1, 4, 1, 7, 7, 1, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  760. 4, 9, 1, 4, 1, 7, 0, 8, 1, 8, 7, 1, - 1, - 1, - 1, - 1,
  761. 4, 0, 3, 7, 4, 3, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  762. 4, 8, 7, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  763. 9, 10, 8, 10, 11, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  764. 3, 0, 9, 3, 9, 11, 11, 9, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  765. 0, 1, 10, 0, 10, 8, 8, 10, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  766. 3, 1, 10, 11, 3, 10, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  767. 1, 2, 11, 1, 11, 9, 9, 11, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  768. 3, 0, 9, 3, 9, 11, 1, 2, 9, 2, 11, 9, - 1, - 1, - 1, - 1,
  769. 0, 2, 11, 8, 0, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  770. 3, 2, 11, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  771. 2, 3, 8, 2, 8, 10, 10, 8, 9, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  772. 9, 10, 2, 0, 9, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  773. 2, 3, 8, 2, 8, 10, 0, 1, 8, 1, 10, 8, - 1, - 1, - 1, - 1,
  774. 1, 10, 2, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  775. 1, 3, 8, 9, 1, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  776. 0, 9, 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  777. 0, 3, 8, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1,
  778. - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1, - 1 ] );