├── LICENSE
├── Main.cpp
├── Makefile
├── README.md
├── Simplify.h
├── a.out.js
├── a.out.wasm
├── index.html
├── vertex.h
└── worker.js


/LICENSE:
--------------------------------------------------------------------------------
 1 | MIT License
 2 | 
 3 | Copyright (c) 2018 MyMiniFactory
 4 | 
 5 | Permission is hereby granted, free of charge, to any person obtaining a copy
 6 | of this software and associated documentation files (the "Software"), to deal
 7 | in the Software without restriction, including without limitation the rights
 8 | to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 9 | copies of the Software, and to permit persons to whom the Software is
10 | furnished to do so, subject to the following conditions:
11 | 
12 | The above copyright notice and this permission notice shall be included in all
13 | copies or substantial portions of the Software.
14 | 
15 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
16 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
17 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
18 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
19 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
20 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
21 | SOFTWARE.
22 | 


--------------------------------------------------------------------------------
/Main.cpp:
--------------------------------------------------------------------------------
  1 | // Simple wrapper for Sven Forstmann's mesh simplification tool
  2 | //
  3 | // Loads a OBJ format mesh, decimates mesh, saves decimated mesh as OBJ format
  4 | // http://voxels.blogspot.com/2014/05/quadric-mesh-simplification-with-source.html
  5 | // https://github.com/sp4cerat/Fast-Quadric-Mesh-Simplification
  6 | //To compile for Linux/OSX (GCC/LLVM)
  7 | //  g++ Main.cpp -O3 -o simplify
  8 | //To compile for Windows (Visual Studio)
  9 | // vcvarsall amd64
 10 | // cl /EHsc Main.cpp /osimplify
 11 | //To execute
 12 | //  ./simplify wall.obj out.obj 0.04
 13 | //
 14 | // Pascal Version by Chris Roden:
 15 | // https://github.com/neurolabusc/Fast-Quadric-Mesh-Simplification-Pascal-
 16 | //
 17 | 
 18 | #include "Simplify.h"
 19 | #include <stdio.h>
 20 | #include <time.h>  // clock_t, clock, CLOCKS_PER_SEC
 21 | 
 22 | void showHelp(const char * argv[]) {
 23 |     const char *cstr = (argv[0]);
 24 |     printf("Usage: %s <input> <output> <ratio> <agressiveness)\n", cstr);
 25 |     printf(" Input: name of existing OBJ format mesh\n");
 26 |     printf(" Output: name for decimated OBJ format mesh\n");
 27 |     printf(" Ratio: (default = 0.5) for example 0.2 will decimate 80%% of triangles\n");
 28 |     printf(" Agressiveness: (default = 7.0) faster or better decimation\n");
 29 |     printf("Examples :\n");
 30 | #if defined(_WIN64) || defined(_WIN32)
 31 |     printf("  %s c:\\dir\\in.obj c:\\dir\\out.obj 0.2\n", cstr);
 32 | #else
 33 |     printf("  %s ~/dir/in.obj ~/dir/out.obj 0.2\n", cstr);
 34 | #endif
 35 | } //showHelp()
 36 | 
 37 | bool is_extension(const char* file_path, const char* extension) {
 38 |     char file_extension[3];
 39 | 
 40 |     file_extension[0] = file_path[strlen(file_path)-3];
 41 |     file_extension[1] = file_path[strlen(file_path)-2];
 42 |     file_extension[2] = file_path[strlen(file_path)-1];
 43 | 
 44 |     return (file_extension[0] == extension[0]
 45 |             and file_extension[1] == extension[1]
 46 |             and file_extension[2] == extension[2]);
 47 | }
 48 | 
 49 | bool is_obj(const char* file_path) {
 50 |     return is_extension(file_path, "obj");
 51 | }
 52 | 
 53 | bool is_stl(const char* file_path) {
 54 |     return is_extension(file_path, "stl");
 55 | }
 56 | 
 57 | int simplify(const char* file_path, const char* export_path, float reduceFraction, float agressiveness) {
 58 | 
 59 |     printf("Mesh Simplification (C)2014 by Sven Forstmann in 2014, MIT License (%zu-bit)\n", sizeof(size_t)*8);
 60 | 
 61 |     if (is_obj(file_path)) {
 62 |         Simplify::load_obj(file_path);
 63 |         printf("loading obj\n");
 64 |     }
 65 |     else if (is_stl(file_path)) {
 66 |         printf("loading stl\n");
 67 |         Simplify::load_stl(file_path);
 68 |     } else {
 69 |         printf("file is not obj or stl %s\n", file_path);
 70 |         return EXIT_FAILURE;
 71 |     }
 72 | 
 73 |     if ((Simplify::triangles.size() < 3) || (Simplify::vertices.size() < 3)) {
 74 |         printf("triangles size or vertices size less than 3\n");
 75 |         return EXIT_FAILURE;
 76 |     }
 77 | 
 78 |     int target_count =  Simplify::triangles.size() >> 1;
 79 | 
 80 |     if (reduceFraction > 1.0) reduceFraction = 1.0; //lossless only
 81 |     if (reduceFraction <= 0.0) {
 82 |         printf("Ratio must be BETWEEN zero and one.\n");
 83 |         return EXIT_FAILURE;
 84 |     }
 85 |     target_count = round((float)Simplify::triangles.size() * reduceFraction);
 86 | 
 87 |     if (target_count < 4) {
 88 |         printf("Object will not survive such extreme decimation\n");
 89 |         return EXIT_FAILURE;
 90 |     }
 91 |     clock_t start = clock();
 92 |     printf("Input: %zu vertices, %zu triangles (target %d)\n", Simplify::vertices.size(), Simplify::triangles.size(), target_count);
 93 |     int startSize = Simplify::triangles.size();
 94 |     Simplify::simplify_mesh(target_count, agressiveness, true);
 95 |     //Simplify::simplify_mesh_lossless( false);
 96 |     if ( Simplify::triangles.size() >= startSize) {
 97 |         printf("Unable to reduce mesh.\n");
 98 |         return EXIT_FAILURE;
 99 |     }
100 | 
101 |     if (is_obj(export_path)) {
102 |         printf("exporting obj\n");
103 |         Simplify::write_obj(export_path);
104 |     }
105 |     else if (is_stl(export_path)) {
106 |         printf("exporting stl\n");
107 |         Simplify::write_stl(export_path);
108 |     } else {
109 |         printf("export file is not obj or stl %s\n", export_path);
110 |         return EXIT_FAILURE;
111 |     }
112 | 
113 |     printf("Output: %zu vertices, %zu triangles (%f reduction; %.4f sec)\n",Simplify::vertices.size(), Simplify::triangles.size()
114 |         , (float)Simplify::triangles.size()/ (float) startSize  , ((float)(clock()-start))/CLOCKS_PER_SEC );
115 |     return EXIT_SUCCESS;
116 | }
117 | 
118 | #ifdef __EMSCRIPTEN__
119 | 
120 | extern "C" {
121 | int simplify(const char* file_path, float reduceFraction, const char* export_path) {
122 |     printf("Going to simplify %s\n", file_path);
123 |     return simplify(file_path, export_path, reduceFraction, 7.0);// aggressive
124 | }
125 | }
126 | 
127 | #else
128 | 
129 | int main(int argc, const char * argv[]) {
130 |     if (argc < 3) {
131 |         showHelp(argv);
132 |         return EXIT_SUCCESS;
133 |     }
134 | 
135 |     const char* file_path = argv[1];
136 |     const char* export_path = argv[2];
137 |     float reduceFraction = 0.5;
138 |     if (argc > 3) {
139 |         reduceFraction = atof(argv[3]);
140 |     }
141 | 
142 |     float agressiveness = 7.0;
143 |     if (argc > 4) {
144 |         agressiveness = atof(argv[4]);
145 |     }
146 |     return simplify(file_path, export_path, reduceFraction, agressiveness);
147 | }
148 | 
149 | #endif
150 | 


--------------------------------------------------------------------------------
/Makefile:
--------------------------------------------------------------------------------
 1 | CXX := Main.cpp
 2 | O3 := -O3
 3 | #WASM := -s DEMANGLE_SUPPORT=1 -s EXTRA_EXPORTED_RUNTIME_METHODS='["ccall", "cwrap", "FS_createDataFile", "FS_readFile", "FS_unlink"]' -s ALLOW_MEMORY_GROWTH=1 -s WASM=1 -s EXPORTED_FUNCTIONS='["_simplify"]'
 4 | WASM := -s DEMANGLE_SUPPORT=1 -s EXPORTED_RUNTIME_METHODS='["ccall", "cwrap", "FS_createDataFile", "FS_readFile", "FS_unlink"]' -s ALLOW_MEMORY_GROWTH=1 -s WASM=1 -s EXPORTED_FUNCTIONS='["_simplify"]'
 5 | 
 6 | 
 7 | all: release
 8 | 
 9 | release: CXX98 += ${O3}
10 | release: debug
11 | 
12 | debug:
13 | 	g++ ${CXX}
14 | 
15 | #make wasm
16 | wasm:
17 | 	em++ ${CXX} ${WASM} ${O3}
18 | wasm_debug:
19 | 	em++ ${CXX} ${WASM}
20 | 


--------------------------------------------------------------------------------
/README.md:
--------------------------------------------------------------------------------
1 | # Fast-Quadric-Mesh-Simplification
2 | WebAssembly version of Fast Quadric Mesh Simplification
3 | 
4 | Checkout https://myminifactory.github.io/Fast-Quadric-Mesh-Simplification/
5 | 


--------------------------------------------------------------------------------
/Simplify.h:
--------------------------------------------------------------------------------
   1 | /////////////////////////////////////////////
   2 | //
   3 | // Mesh Simplification Tutorial
   4 | //
   5 | // (C) by Sven Forstmann in 2014
   6 | //
   7 | // License : MIT
   8 | // http://opensource.org/licenses/MIT
   9 | //
  10 | //https://github.com/sp4cerat/Fast-Quadric-Mesh-Simplification
  11 | //
  12 | // 5/2016: Chris Rorden created minimal version for OSX/Linux/Windows compile
  13 | 
  14 | //#include <iostream>
  15 | //#include <stddef.h>
  16 | //#include <functional>
  17 | //#include <sys/stat.h>
  18 | //#include <stdbool.h>
  19 | #include <string.h>
  20 | #include <string>
  21 | #include <fstream>
  22 | #include <algorithm>
  23 | #include <iostream>
  24 | #include <cstring> // memcpy
  25 | //#include <ctype.h>
  26 | //#include <float.h>
  27 | #include <stdio.h>
  28 | #include <stdlib.h>
  29 | #include <vector>
  30 | #include <math.h>
  31 | #include <float.h> //FLT_EPSILON, DBL_EPSILON
  32 | #include <stdint.h>
  33 | 
  34 | #define loopi(start_l,end_l) for ( int i=start_l;i<end_l;++i )
  35 | #define loopi(start_l,end_l) for ( int i=start_l;i<end_l;++i )
  36 | #define loopj(start_l,end_l) for ( int j=start_l;j<end_l;++j )
  37 | #define loopk(start_l,end_l) for ( int k=start_l;k<end_l;++k )
  38 | 
  39 | /*
  40 |  *  Represents an optionally-indexed vertex in space
  41 |  */
  42 | struct VertexSTL
  43 | {
  44 |     VertexSTL() {}
  45 |     VertexSTL(float x, float y, float z) : x(x), y(y), z(z) {}
  46 | 
  47 |     float x, y, z;
  48 |     unsigned int i;
  49 | 
  50 |     bool operator!=(const VertexSTL& rhs) const
  51 |     {
  52 |         return x != rhs.x || y != rhs.y || z != rhs.z;
  53 |     }
  54 |     bool operator<(const VertexSTL& rhs) const
  55 |     {
  56 |         if      (x != rhs.x)    return x < rhs.x;
  57 |         else if (y != rhs.y)    return y < rhs.y;
  58 |         else if (z != rhs.z)    return z < rhs.z;
  59 |         else                    return false;
  60 |     }
  61 | };
  62 | inline std::string trim(std::string& str)
  63 | {
  64 |     str.erase(0, str.find_first_not_of(' '));       //prefixing spaces
  65 |     str.erase(str.find_last_not_of(' ')+1);         //surfixing spaces
  66 |     return str;
  67 | }
  68 | 
  69 | inline VertexSTL get_vector(std::string& str)
  70 | {
  71 |     // "vertex float float float"
  72 |     float x, y, z;
  73 | 
  74 |     // if(sscanf(str.c_str(),"%*[ \t]vertex%*[ \t]%f%*[ \t]%f%*[ \t]%f%*[ \t]",
  75 |         // &x,	&y,	&z) == 3)
  76 |         // printf("%f %f %f", x, y, z);
  77 |     if(sscanf(str.c_str(),"vertex %f %f %f",
  78 |         &x,	&y,	&z) != 3) {
  79 |         printf("weird format ascii stl exiting\n");
  80 |         exit(1);
  81 |     }
  82 |     VertexSTL v(x, y, z);
  83 | 
  84 |     // int space0 = str.find_first_of(' ');
  85 |     // str.erase(0, space0); // remove "vertex"
  86 |     // str.erase(0, str.find_first_not_of(' ')); //prefixing spaces
  87 |     // int space1 = str.find_first_of(' ');
  88 | 
  89 |     // float x = std::stof(str.substr(0, space1));
  90 | 
  91 |     // str.erase(0, space1+1); // remove x
  92 |     // str.erase(0, str.find_first_not_of(' ')); //prefixing spaces
  93 | 
  94 |     // int space2 = str.find_last_of(' ');
  95 | 
  96 |     // VertexSTL v(x,
  97 |         // std::stof(str.substr(0, space2)),
  98 |         // std::stof(str.substr(space2, str.size()))
  99 |     // );
 100 |     return v;
 101 | }
 102 | 
 103 | 
 104 | struct vector3
 105 | {
 106 | double x, y, z;
 107 | };
 108 | 
 109 | struct vec3f
 110 | {
 111 |     double x, y, z;
 112 | 
 113 |     inline vec3f( void ) {}
 114 | 
 115 |     //inline vec3f operator =( vector3 a )
 116 | 	// { vec3f b ; b.x = a.x; b.y = a.y; b.z = a.z; return b;}
 117 | 
 118 |     inline vec3f( vector3 a )
 119 | 	 { x = a.x; y = a.y; z = a.z; }
 120 | 
 121 |     inline vec3f( const double X, const double Y, const double Z )
 122 |     { x = X; y = Y; z = Z; }
 123 | 
 124 |     inline vec3f operator + ( const vec3f& a ) const
 125 |     { return vec3f( x + a.x, y + a.y, z + a.z ); }
 126 | 
 127 | 	inline vec3f operator += ( const vec3f& a ) const
 128 |     { return vec3f( x + a.x, y + a.y, z + a.z ); }
 129 | 
 130 |     inline vec3f operator * ( const double a ) const
 131 |     { return vec3f( x * a, y * a, z * a ); }
 132 | 
 133 |     inline vec3f operator * ( const vec3f a ) const
 134 |     { return vec3f( x * a.x, y * a.y, z * a.z ); }
 135 | 
 136 |     inline vec3f v3 () const
 137 |     { return vec3f( x , y, z ); }
 138 | 
 139 |     inline vec3f operator = ( const vector3 a )
 140 |     { x=a.x;y=a.y;z=a.z;return *this; }
 141 | 
 142 |     inline vec3f operator = ( const vec3f a )
 143 |     { x=a.x;y=a.y;z=a.z;return *this; }
 144 | 
 145 |     inline vec3f operator / ( const vec3f a ) const
 146 |     { return vec3f( x / a.x, y / a.y, z / a.z ); }
 147 | 
 148 |     inline vec3f operator - ( const vec3f& a ) const
 149 |     { return vec3f( x - a.x, y - a.y, z - a.z ); }
 150 | 
 151 |     inline vec3f operator / ( const double a ) const
 152 |     { return vec3f( x / a, y / a, z / a ); }
 153 | 
 154 |     inline double dot( const vec3f& a ) const
 155 |     { return a.x*x + a.y*y + a.z*z; }
 156 | 
 157 |     inline vec3f cross( const vec3f& a , const vec3f& b )
 158 |     {
 159 | 		x = a.y * b.z - a.z * b.y;
 160 | 		y = a.z * b.x - a.x * b.z;
 161 | 		z = a.x * b.y - a.y * b.x;
 162 | 		return *this;
 163 | 	}
 164 | 
 165 |     inline double angle( const vec3f& v )
 166 |     {
 167 | 		vec3f a = v , b = *this;
 168 | 		double dot = v.x*x + v.y*y + v.z*z;
 169 | 		double len = a.length() * b.length();
 170 | 		if(len==0)len=0.00001f;
 171 | 		double input = dot  / len;
 172 | 		if (input<-1) input=-1;
 173 | 		if (input>1) input=1;
 174 | 		return (double) acos ( input );
 175 | 	}
 176 | 
 177 |     inline double angle2( const vec3f& v , const vec3f& w )
 178 |     {
 179 | 		vec3f a = v , b= *this;
 180 | 		double dot = a.x*b.x + a.y*b.y + a.z*b.z;
 181 | 		double len = a.length() * b.length();
 182 | 		if(len==0)len=1;
 183 | 
 184 | 		vec3f plane; plane.cross( b,w );
 185 | 
 186 | 		if ( plane.x * a.x + plane.y * a.y + plane.z * a.z > 0 )
 187 | 			return (double) -acos ( dot  / len );
 188 | 
 189 | 		return (double) acos ( dot  / len );
 190 | 	}
 191 | 
 192 |     inline vec3f rot_x( double a )
 193 |     {
 194 | 		double yy = cos ( a ) * y + sin ( a ) * z;
 195 | 		double zz = cos ( a ) * z - sin ( a ) * y;
 196 | 		y = yy; z = zz;
 197 | 		return *this;
 198 | 	}
 199 |     inline vec3f rot_y( double a )
 200 |     {
 201 | 		double xx = cos ( -a ) * x + sin ( -a ) * z;
 202 | 		double zz = cos ( -a ) * z - sin ( -a ) * x;
 203 | 		x = xx; z = zz;
 204 | 		return *this;
 205 | 	}
 206 |     inline void clamp( double min, double max )
 207 |     {
 208 | 		if (x<min) x=min;
 209 | 		if (y<min) y=min;
 210 | 		if (z<min) z=min;
 211 | 		if (x>max) x=max;
 212 | 		if (y>max) y=max;
 213 | 		if (z>max) z=max;
 214 | 	}
 215 |     inline vec3f rot_z( double a )
 216 |     {
 217 | 		double yy = cos ( a ) * y + sin ( a ) * x;
 218 | 		double xx = cos ( a ) * x - sin ( a ) * y;
 219 | 		y = yy; x = xx;
 220 | 		return *this;
 221 | 	}
 222 |     inline vec3f invert()
 223 | 	{
 224 | 		x=-x;y=-y;z=-z;return *this;
 225 | 	}
 226 |     inline vec3f frac()
 227 | 	{
 228 | 		return vec3f(
 229 | 			x-double(int(x)),
 230 | 			y-double(int(y)),
 231 | 			z-double(int(z))
 232 | 			);
 233 | 	}
 234 | 
 235 |     inline vec3f integer()
 236 | 	{
 237 | 		return vec3f(
 238 | 			double(int(x)),
 239 | 			double(int(y)),
 240 | 			double(int(z))
 241 | 			);
 242 | 	}
 243 | 
 244 |     inline double length() const
 245 |     {
 246 | 		return (double)sqrt(x*x + y*y + z*z);
 247 | 	}
 248 | 
 249 |     inline vec3f normalize( double desired_length = 1 )
 250 |     {
 251 | 		double square = sqrt(x*x + y*y + z*z);
 252 | 		/*
 253 | 		if (square <= 0.00001f )
 254 | 		{
 255 | 			x=1;y=0;z=0;
 256 | 			return *this;
 257 | 		}*/
 258 | 		//double len = desired_length / square;
 259 | 		x/=square;y/=square;z/=square;
 260 | 
 261 | 		return *this;
 262 | 	}
 263 |     static vec3f normalize( vec3f a );
 264 | 
 265 | 	static void random_init();
 266 | 	static double random_double();
 267 | 	static vec3f random();
 268 | 
 269 | 	static int random_number;
 270 | 
 271 | 	double random_double_01(double a){
 272 | 		double rnf=a*14.434252+a*364.2343+a*4213.45352+a*2341.43255+a*254341.43535+a*223454341.3523534245+23453.423412;
 273 | 		int rni=((int)rnf)%100000;
 274 | 		return double(rni)/(100000.0f-1.0f);
 275 | 	}
 276 | 
 277 | 	vec3f random01_fxyz(){
 278 | 		x=(double)random_double_01(x);
 279 | 		y=(double)random_double_01(y);
 280 | 		z=(double)random_double_01(z);
 281 | 		return *this;
 282 | 	}
 283 | };
 284 | 
 285 | 
 286 | double min(double v1, double v2) {
 287 | 	return fmin(v1,v2);
 288 | }
 289 | 
 290 | 
 291 | class SymetricMatrix {
 292 | 
 293 | 	public:
 294 | 
 295 | 	// Constructor
 296 | 
 297 | 	SymetricMatrix(double c=0) { loopi(0,10) m[i] = c;  }
 298 | 
 299 | 	SymetricMatrix(	double m11, double m12, double m13, double m14,
 300 | 			            double m22, double m23, double m24,
 301 | 			                        double m33, double m34,
 302 | 			                                    double m44) {
 303 | 			 m[0] = m11;  m[1] = m12;  m[2] = m13;  m[3] = m14;
 304 | 			              m[4] = m22;  m[5] = m23;  m[6] = m24;
 305 | 			                           m[7] = m33;  m[8] = m34;
 306 | 			                                        m[9] = m44;
 307 | 	}
 308 | 
 309 | 	// Make plane
 310 | 
 311 | 	SymetricMatrix(double a,double b,double c,double d)
 312 | 	{
 313 | 		m[0] = a*a;  m[1] = a*b;  m[2] = a*c;  m[3] = a*d;
 314 | 		             m[4] = b*b;  m[5] = b*c;  m[6] = b*d;
 315 | 		                          m[7 ] =c*c; m[8 ] = c*d;
 316 | 		                                       m[9 ] = d*d;
 317 | 	}
 318 | 
 319 | 	double operator[](int c) const { return m[c]; }
 320 | 
 321 | 	// Determinant
 322 | 
 323 | 	double det(	int a11, int a12, int a13,
 324 | 				int a21, int a22, int a23,
 325 | 				int a31, int a32, int a33)
 326 | 	{
 327 | 		double det =  m[a11]*m[a22]*m[a33] + m[a13]*m[a21]*m[a32] + m[a12]*m[a23]*m[a31]
 328 | 					- m[a13]*m[a22]*m[a31] - m[a11]*m[a23]*m[a32]- m[a12]*m[a21]*m[a33];
 329 | 		return det;
 330 | 	}
 331 | 
 332 | 	const SymetricMatrix operator+(const SymetricMatrix& n) const
 333 | 	{
 334 | 		return SymetricMatrix( m[0]+n[0],   m[1]+n[1],   m[2]+n[2],   m[3]+n[3],
 335 | 						                    m[4]+n[4],   m[5]+n[5],   m[6]+n[6],
 336 | 						                                 m[ 7]+n[ 7], m[ 8]+n[8 ],
 337 | 						                                              m[ 9]+n[9 ]);
 338 | 	}
 339 | 
 340 | 	SymetricMatrix& operator+=(const SymetricMatrix& n)
 341 | 	{
 342 | 		 m[0]+=n[0];   m[1]+=n[1];   m[2]+=n[2];   m[3]+=n[3];
 343 | 		 m[4]+=n[4];   m[5]+=n[5];   m[6]+=n[6];   m[7]+=n[7];
 344 | 		 m[8]+=n[8];   m[9]+=n[9];
 345 | 		return *this;
 346 | 	}
 347 | 
 348 | 	double m[10];
 349 | };
 350 | ///////////////////////////////////////////
 351 | 
 352 | namespace Simplify
 353 | {
 354 | 	// Global Variables & Strctures
 355 | 
 356 | 	struct Triangle { int v[3];double err[4];int deleted,dirty;vec3f n; };
 357 | 	struct Vertex { vec3f p;int tstart,tcount;SymetricMatrix q;int border;};
 358 | 	struct Ref { int tid,tvertex; };
 359 | 	std::vector<Triangle> triangles;
 360 | 	std::vector<Vertex> vertices;
 361 | 	std::vector<Ref> refs;
 362 | 
 363 | 	// Helper functions
 364 | 
 365 | 	double vertex_error(SymetricMatrix q, double x, double y, double z);
 366 | 	double calculate_error(int id_v1, int id_v2, vec3f &p_result);
 367 | 	bool flipped(vec3f p,int i0,int i1,Vertex &v0,Vertex &v1,std::vector<int> &deleted);
 368 | 	void update_triangles(int i0,Vertex &v,std::vector<int> &deleted,int &deleted_triangles);
 369 | 	void update_mesh(int iteration);
 370 | 	void compact_mesh();
 371 | 	//
 372 | 	// Main simplification function
 373 | 	//
 374 | 	// target_count  : target nr. of triangles
 375 | 	// agressiveness : sharpness to increase the threashold.
 376 | 	//                 5..8 are good numbers
 377 | 	//                 more iterations yield higher quality
 378 | 	//
 379 | 
 380 | 	void simplify_mesh(int target_count, double agressiveness=7, bool verbose=false)
 381 | 	{
 382 | 		// init
 383 | 		loopi(0,triangles.size()) triangles[i].deleted=0;
 384 | 
 385 | 		// main iteration loop
 386 | 		int deleted_triangles=0;
 387 | 		std::vector<int> deleted0,deleted1;
 388 | 		int triangle_count=triangles.size();
 389 | 		//int iteration = 0;
 390 | 		//loop(iteration,0,100)
 391 | 		for (int iteration = 0; iteration < 100; iteration ++)
 392 | 		{
 393 | 			if(triangle_count-deleted_triangles<=target_count)break;
 394 | 
 395 | 			// update mesh once in a while
 396 | 			if(iteration%5==0)
 397 | 			{
 398 | 				update_mesh(iteration);
 399 | 			}
 400 | 
 401 | 			// clear dirty flag
 402 | 			loopi(0,triangles.size()) triangles[i].dirty=0;
 403 | 
 404 | 			//
 405 | 			// All triangles with edges below the threshold will be removed
 406 | 			//
 407 | 			// The following numbers works well for most models.
 408 | 			// If it does not, try to adjust the 3 parameters
 409 | 			//
 410 | 			double threshold = 0.000000001*pow(double(iteration+3),agressiveness);
 411 | 
 412 | 			// target number of triangles reached ? Then break
 413 | 			if ((verbose) && (iteration%5==0)) {
 414 | 				printf("iteration %d - triangles %d threshold %g\n",iteration,triangle_count-deleted_triangles, threshold);
 415 | 			}
 416 | 
 417 | 			// remove vertices & mark deleted triangles
 418 | 			loopi(0,triangles.size())
 419 | 			{
 420 | 				Triangle &t=triangles[i];
 421 | 				if(t.err[3]>threshold) continue;
 422 | 				if(t.deleted) continue;
 423 | 				if(t.dirty) continue;
 424 | 
 425 | 				loopj(0,3)if(t.err[j]<threshold)
 426 | 				{
 427 | 
 428 | 					int i0=t.v[ j     ]; Vertex &v0 = vertices[i0];
 429 | 					int i1=t.v[(j+1)%3]; Vertex &v1 = vertices[i1];
 430 | 					// Border check
 431 | 					if(v0.border != v1.border)  continue;
 432 | 
 433 | 					// Compute vertex to collapse to
 434 | 					vec3f p;
 435 | 					calculate_error(i0,i1,p);
 436 | 					deleted0.resize(v0.tcount); // normals temporarily
 437 | 					deleted1.resize(v1.tcount); // normals temporarily
 438 | 					// dont remove if flipped
 439 | 					if( flipped(p,i0,i1,v0,v1,deleted0) ) continue;
 440 | 
 441 | 					if( flipped(p,i1,i0,v1,v0,deleted1) ) continue;
 442 | 
 443 | 
 444 | 					// not flipped, so remove edge
 445 | 					v0.p=p;
 446 | 					v0.q=v1.q+v0.q;
 447 | 					int tstart=refs.size();
 448 | 
 449 | 					update_triangles(i0,v0,deleted0,deleted_triangles);
 450 | 					update_triangles(i0,v1,deleted1,deleted_triangles);
 451 | 
 452 | 					int tcount=refs.size()-tstart;
 453 | 
 454 | 					if(tcount<=v0.tcount)
 455 | 					{
 456 | 						// save ram
 457 | 						if(tcount)memcpy(&refs[v0.tstart],&refs[tstart],tcount*sizeof(Ref));
 458 | 					}
 459 | 					else
 460 | 						// append
 461 | 						v0.tstart=tstart;
 462 | 
 463 | 					v0.tcount=tcount;
 464 | 					break;
 465 | 				}
 466 | 				// done?
 467 | 				if(triangle_count-deleted_triangles<=target_count)break;
 468 | 			}
 469 | 		}
 470 | 		// clean up mesh
 471 | 		compact_mesh();
 472 | 	} //simplify_mesh()
 473 | 
 474 | 	void simplify_mesh_lossless(bool verbose=false)
 475 | 	{
 476 | 		// init
 477 | 		loopi(0,triangles.size()) triangles[i].deleted=0;
 478 | 
 479 | 		// main iteration loop
 480 | 		int deleted_triangles=0;
 481 | 		std::vector<int> deleted0,deleted1;
 482 | 		int triangle_count=triangles.size();
 483 | 		//int iteration = 0;
 484 | 		//loop(iteration,0,100)
 485 | 		for (int iteration = 0; iteration < 9999; iteration ++)
 486 | 		{
 487 | 			// update mesh constantly
 488 | 			update_mesh(iteration);
 489 | 			// clear dirty flag
 490 | 			loopi(0,triangles.size()) triangles[i].dirty=0;
 491 | 			//
 492 | 			// All triangles with edges below the threshold will be removed
 493 | 			//
 494 | 			// The following numbers works well for most models.
 495 | 			// If it does not, try to adjust the 3 parameters
 496 | 			//
 497 | 			double threshold = DBL_EPSILON; //1.0E-3 EPS;
 498 | 			if (verbose) {
 499 | 				printf("lossless iteration %d\n", iteration);
 500 | 			}
 501 | 
 502 | 			// remove vertices & mark deleted triangles
 503 | 			loopi(0,triangles.size())
 504 | 			{
 505 | 				Triangle &t=triangles[i];
 506 | 				if(t.err[3]>threshold) continue;
 507 | 				if(t.deleted) continue;
 508 | 				if(t.dirty) continue;
 509 | 
 510 | 				loopj(0,3)if(t.err[j]<threshold)
 511 | 				{
 512 | 					int i0=t.v[ j     ]; Vertex &v0 = vertices[i0];
 513 | 					int i1=t.v[(j+1)%3]; Vertex &v1 = vertices[i1];
 514 | 
 515 | 					// Border check
 516 | 					if(v0.border != v1.border)  continue;
 517 | 
 518 | 					// Compute vertex to collapse to
 519 | 					vec3f p;
 520 | 					calculate_error(i0,i1,p);
 521 | 
 522 | 					deleted0.resize(v0.tcount); // normals temporarily
 523 | 					deleted1.resize(v1.tcount); // normals temporarily
 524 | 
 525 | 					// dont remove if flipped
 526 | 					if( flipped(p,i0,i1,v0,v1,deleted0) ) continue;
 527 | 					if( flipped(p,i1,i0,v1,v0,deleted1) ) continue;
 528 | 
 529 | 					// not flipped, so remove edge
 530 | 					v0.p=p;
 531 | 					v0.q=v1.q+v0.q;
 532 | 					int tstart=refs.size();
 533 | 
 534 | 					update_triangles(i0,v0,deleted0,deleted_triangles);
 535 | 					update_triangles(i0,v1,deleted1,deleted_triangles);
 536 | 
 537 | 					int tcount=refs.size()-tstart;
 538 | 
 539 | 					if(tcount<=v0.tcount)
 540 | 					{
 541 | 						// save ram
 542 | 						if(tcount)memcpy(&refs[v0.tstart],&refs[tstart],tcount*sizeof(Ref));
 543 | 					}
 544 | 					else
 545 | 						// append
 546 | 						v0.tstart=tstart;
 547 | 
 548 | 					v0.tcount=tcount;
 549 | 					break;
 550 | 				}
 551 | 			}
 552 | 			if(deleted_triangles<=0)break;
 553 | 			deleted_triangles=0;
 554 | 		} //for each iteration
 555 | 		// clean up mesh
 556 | 		compact_mesh();
 557 | 	} //simplify_mesh_lossless()
 558 | 
 559 | 
 560 | 	// Check if a triangle flips when this edge is removed
 561 | 
 562 | 	bool flipped(vec3f p,int i0,int i1,Vertex &v0,Vertex &v1,std::vector<int> &deleted)
 563 | 	{
 564 | 
 565 | 		loopk(0,v0.tcount)
 566 | 		{
 567 | 			Triangle &t=triangles[refs[v0.tstart+k].tid];
 568 | 			if(t.deleted)continue;
 569 | 
 570 | 			int s=refs[v0.tstart+k].tvertex;
 571 | 			int id1=t.v[(s+1)%3];
 572 | 			int id2=t.v[(s+2)%3];
 573 | 
 574 | 			if(id1==i1 || id2==i1) // delete ?
 575 | 			{
 576 | 
 577 | 				deleted[k]=1;
 578 | 				continue;
 579 | 			}
 580 | 			vec3f d1 = vertices[id1].p-p; d1.normalize();
 581 | 			vec3f d2 = vertices[id2].p-p; d2.normalize();
 582 | 			if(fabs(d1.dot(d2))>0.999) return true;
 583 | 			vec3f n;
 584 | 			n.cross(d1,d2);
 585 | 			n.normalize();
 586 | 			deleted[k]=0;
 587 | 			if(n.dot(t.n)<0.2) return true;
 588 | 		}
 589 | 		return false;
 590 | 	}
 591 | 
 592 | 	// Update triangle connections and edge error after a edge is collapsed
 593 | 
 594 | 	void update_triangles(int i0,Vertex &v,std::vector<int> &deleted,int &deleted_triangles)
 595 | 	{
 596 | 		vec3f p;
 597 | 		loopk(0,v.tcount)
 598 | 		{
 599 | 			Ref &r=refs[v.tstart+k];
 600 | 			Triangle &t=triangles[r.tid];
 601 | 			if(t.deleted)continue;
 602 | 			if(deleted[k])
 603 | 			{
 604 | 				t.deleted=1;
 605 | 				deleted_triangles++;
 606 | 				continue;
 607 | 			}
 608 | 			t.v[r.tvertex]=i0;
 609 | 			t.dirty=1;
 610 | 			t.err[0]=calculate_error(t.v[0],t.v[1],p);
 611 | 			t.err[1]=calculate_error(t.v[1],t.v[2],p);
 612 | 			t.err[2]=calculate_error(t.v[2],t.v[0],p);
 613 | 			t.err[3]=min(t.err[0],min(t.err[1],t.err[2]));
 614 | 			refs.push_back(r);
 615 | 		}
 616 | 	}
 617 | 
 618 | 	// compact triangles, compute edge error and build reference list
 619 | 
 620 | 	void update_mesh(int iteration)
 621 | 	{
 622 | 		if(iteration>0) // compact triangles
 623 | 		{
 624 | 			int dst=0;
 625 | 			loopi(0,triangles.size())
 626 | 			if(!triangles[i].deleted)
 627 | 			{
 628 | 				triangles[dst++]=triangles[i];
 629 | 			}
 630 | 			triangles.resize(dst);
 631 | 		}
 632 | 		//
 633 | 		// Init Quadrics by Plane & Edge Errors
 634 | 		//
 635 | 		// required at the beginning ( iteration == 0 )
 636 | 		// recomputing during the simplification is not required,
 637 | 		// but mostly improves the result for closed meshes
 638 | 		//
 639 | 
 640 | 		// Init Reference ID list
 641 | 		loopi(0,vertices.size())
 642 | 		{
 643 | 			vertices[i].tstart=0;
 644 | 			vertices[i].tcount=0;
 645 | 		}
 646 | 		loopi(0,triangles.size())
 647 | 		{
 648 | 			Triangle &t=triangles[i];
 649 | 			loopj(0,3) vertices[t.v[j]].tcount++;
 650 | 		}
 651 | 		int tstart=0;
 652 | 		loopi(0,vertices.size())
 653 | 		{
 654 | 			Vertex &v=vertices[i];
 655 | 			v.tstart=tstart;
 656 | 			tstart+=v.tcount;
 657 | 			v.tcount=0;
 658 | 		}
 659 | 
 660 | 		// Write References
 661 | 		refs.resize(triangles.size()*3);
 662 | 		loopi(0,triangles.size())
 663 | 		{
 664 | 			Triangle &t=triangles[i];
 665 | 			loopj(0,3)
 666 | 			{
 667 | 				Vertex &v=vertices[t.v[j]];
 668 | 				refs[v.tstart+v.tcount].tid=i;
 669 | 				refs[v.tstart+v.tcount].tvertex=j;
 670 | 				v.tcount++;
 671 | 			}
 672 | 		}
 673 | 
 674 | 		if( iteration != 0 ) return;
 675 | 		//Initial setup: find border vertices and estimate cost function
 676 | 
 677 | 		// Identify boundary : vertices[].border=0,1
 678 | 
 679 | 		std::vector<int> vcount,vids;
 680 | 
 681 | 		loopi(0,vertices.size())
 682 | 			vertices[i].border=0;
 683 | 
 684 | 		loopi(0,vertices.size())
 685 | 		{
 686 | 			Vertex &v=vertices[i];
 687 | 			vcount.clear();
 688 | 			vids.clear();
 689 | 			loopj(0,v.tcount)
 690 | 			{
 691 | 				int k=refs[v.tstart+j].tid;
 692 | 				Triangle &t=triangles[k];
 693 | 				loopk(0,3)
 694 | 				{
 695 | 					int ofs=0,id=t.v[k];
 696 | 					while(ofs<vcount.size())
 697 | 					{
 698 | 						if(vids[ofs]==id)break;
 699 | 						ofs++;
 700 | 					}
 701 | 					if(ofs==vcount.size())
 702 | 					{
 703 | 						vcount.push_back(1);
 704 | 						vids.push_back(id);
 705 | 					}
 706 | 					else
 707 | 						vcount[ofs]++;
 708 | 				}
 709 | 			}
 710 | 			loopj(0,vcount.size()) if(vcount[j]==1)
 711 | 				vertices[vids[j]].border=1;
 712 | 		}
 713 | 
 714 | 		loopi(0,vertices.size())
 715 | 		vertices[i].q=SymetricMatrix(0.0);
 716 | 
 717 | 		loopi(0,triangles.size())
 718 | 		{
 719 | 			Triangle &t=triangles[i];
 720 | 			vec3f n,p[3];
 721 | 			loopj(0,3) p[j]=vertices[t.v[j]].p;
 722 | 			n.cross(p[1]-p[0],p[2]-p[0]);
 723 | 			n.normalize();
 724 | 			t.n=n;
 725 | 			loopj(0,3) vertices[t.v[j]].q =
 726 | 				vertices[t.v[j]].q+SymetricMatrix(n.x,n.y,n.z,-n.dot(p[0]));
 727 | 		}
 728 | 		loopi(0,triangles.size())
 729 | 		{
 730 | 			// Calc Edge Error
 731 | 			Triangle &t=triangles[i];vec3f p;
 732 | 			loopj(0,3) t.err[j]=calculate_error(t.v[j],t.v[(j+1)%3],p);
 733 | 			t.err[3]=min(t.err[0],min(t.err[1],t.err[2]));
 734 | 		}
 735 | 
 736 | 
 737 | 	}
 738 | 
 739 | 	// Finally compact mesh before exiting
 740 | 
 741 | 	void compact_mesh()
 742 | 	{
 743 | 		int dst=0;
 744 | 		loopi(0,vertices.size())
 745 | 		{
 746 | 			vertices[i].tcount=0;
 747 | 		}
 748 | 		loopi(0,triangles.size())
 749 | 		if(!triangles[i].deleted)
 750 | 		{
 751 | 			Triangle &t=triangles[i];
 752 | 			triangles[dst++]=t;
 753 | 			loopj(0,3)vertices[t.v[j]].tcount=1;
 754 | 		}
 755 | 		triangles.resize(dst);
 756 | 		dst=0;
 757 | 		loopi(0,vertices.size())
 758 | 		if(vertices[i].tcount)
 759 | 		{
 760 | 			vertices[i].tstart=dst;
 761 | 			vertices[dst].p=vertices[i].p;
 762 | 			dst++;
 763 | 		}
 764 | 		loopi(0,triangles.size())
 765 | 		{
 766 | 			Triangle &t=triangles[i];
 767 | 			loopj(0,3)t.v[j]=vertices[t.v[j]].tstart;
 768 | 		}
 769 | 		vertices.resize(dst);
 770 | 	}
 771 | 
 772 | 	// Error between vertex and Quadric
 773 | 
 774 | 	double vertex_error(SymetricMatrix q, double x, double y, double z)
 775 | 	{
 776 |  		return   q[0]*x*x + 2*q[1]*x*y + 2*q[2]*x*z + 2*q[3]*x + q[4]*y*y
 777 |  		     + 2*q[5]*y*z + 2*q[6]*y + q[7]*z*z + 2*q[8]*z + q[9];
 778 | 	}
 779 | 
 780 | 	// Error for one edge
 781 | 
 782 | 	double calculate_error(int id_v1, int id_v2, vec3f &p_result)
 783 | 	{
 784 | 		// compute interpolated vertex
 785 | 
 786 | 		SymetricMatrix q = vertices[id_v1].q + vertices[id_v2].q;
 787 | 		bool   border = vertices[id_v1].border & vertices[id_v2].border;
 788 | 		double error=0;
 789 | 		double det = q.det(0, 1, 2, 1, 4, 5, 2, 5, 7);
 790 | 		if ( det != 0 && !border )
 791 | 		{
 792 | 
 793 | 			// q_delta is invertible
 794 | 			p_result.x = -1/det*(q.det(1, 2, 3, 4, 5, 6, 5, 7 , 8));	// vx = A41/det(q_delta)
 795 | 			p_result.y =  1/det*(q.det(0, 2, 3, 1, 5, 6, 2, 7 , 8));	// vy = A42/det(q_delta)
 796 | 			p_result.z = -1/det*(q.det(0, 1, 3, 1, 4, 6, 2, 5,  8));	// vz = A43/det(q_delta)
 797 | 
 798 | 			error = vertex_error(q, p_result.x, p_result.y, p_result.z);
 799 | 		}
 800 | 		else
 801 | 		{
 802 | 			// det = 0 -> try to find best result
 803 | 			vec3f p1=vertices[id_v1].p;
 804 | 			vec3f p2=vertices[id_v2].p;
 805 | 			vec3f p3=(p1+p2)/2;
 806 | 			double error1 = vertex_error(q, p1.x,p1.y,p1.z);
 807 | 			double error2 = vertex_error(q, p2.x,p2.y,p2.z);
 808 | 			double error3 = vertex_error(q, p3.x,p3.y,p3.z);
 809 | 			error = min(error1, min(error2, error3));
 810 | 			if (error1 == error) p_result=p1;
 811 | 			if (error2 == error) p_result=p2;
 812 | 			if (error3 == error) p_result=p3;
 813 | 		}
 814 | 		return error;
 815 | 	}
 816 | 
 817 | 	//Option : Load OBJ
 818 | 	void load_obj(const char* filename){
 819 | 		vertices.clear();
 820 | 		triangles.clear();
 821 | 		//printf ( "Loading Objects %s ... \n",filename);
 822 | 		FILE* fn;
 823 | 		if(filename==NULL)		return ;
 824 | 		if((char)filename[0]==0)	return ;
 825 | 		if ((fn = fopen(filename, "rb")) == NULL)
 826 | 		{
 827 | 			printf ( "File %s not found!\n" ,filename );
 828 | 			return;
 829 | 		}
 830 | 		char line[1000];
 831 | 		memset ( line,0,1000 );
 832 | 		int vertex_cnt = 0;
 833 | 		while(fgets( line, 1000, fn ) != NULL)
 834 | 		{
 835 | 			Vertex v;
 836 | 			if ( line[0] == 'v' )
 837 | 			{
 838 | 				if ( line[1] == ' ' )
 839 | 				if(sscanf(line,"v %lf %lf %lf",
 840 | 					&v.p.x,	&v.p.y,	&v.p.z)==3)
 841 | 				{
 842 | 					vertices.push_back(v);
 843 | 				}
 844 | 			}
 845 | 			int integers[9];
 846 | 			if ( line[0] == 'f' )
 847 | 			{
 848 | 				Triangle t;
 849 | 				bool tri_ok = false;
 850 | 
 851 | 				if(sscanf(line,"f %d %d %d",
 852 | 					&integers[0],&integers[1],&integers[2])==3)
 853 | 				{
 854 | 					tri_ok = true;
 855 | 				}else
 856 | 				if(sscanf(line,"f %d// %d// %d//",
 857 | 					&integers[0],&integers[1],&integers[2])==3)
 858 | 				{
 859 | 					tri_ok = true;
 860 | 				}else
 861 | 				if(sscanf(line,"f %d//%d %d//%d %d//%d",
 862 | 					&integers[0],&integers[3],
 863 | 					&integers[1],&integers[4],
 864 | 					&integers[2],&integers[5])==6)
 865 | 				{
 866 | 					tri_ok = true;
 867 | 				}else
 868 | 				if(sscanf(line,"f %d/%d/%d %d/%d/%d %d/%d/%d",
 869 | 					&integers[0],&integers[6],&integers[3],
 870 | 					&integers[1],&integers[7],&integers[4],
 871 | 					&integers[2],&integers[8],&integers[5])==9)
 872 | 				{
 873 | 					tri_ok = true;
 874 | 				}
 875 | 				else
 876 | 				{
 877 | 					printf("unrecognized sequence exiting \n");
 878 | 					printf("%s\n",line);
 879 |                     exit(1);
 880 | 				}
 881 | 				if ( tri_ok )
 882 | 				{
 883 | 					t.v[0] = integers[0]-1-vertex_cnt;
 884 | 					t.v[1] = integers[1]-1-vertex_cnt;
 885 | 					t.v[2] = integers[2]-1-vertex_cnt;
 886 | 
 887 | 					//tri.material = material;
 888 | 					//geo.triangles.push_back ( tri );
 889 | 					triangles.push_back(t);
 890 | 					//state_before = state;
 891 | 					//state ='f';
 892 | 				}
 893 | 			}
 894 | 		}
 895 | 		fclose(fn);
 896 | 		//printf("load_obj: vertices = %lu, triangles = %lu\n", vertices.size(), triangles.size() );
 897 | 	} // load_obj()
 898 | 
 899 | 	// Optional : Store as OBJ
 900 | 
 901 | 	void write_obj(const char* filename)
 902 | 	{
 903 | 		FILE *file=fopen(filename, "w");
 904 | 		if (!file)
 905 | 		{
 906 | 			printf("write_obj: can't write data file \"%s\".\n", filename);
 907 | 			exit(0);
 908 | 		}
 909 | 		loopi(0,vertices.size())
 910 | 		{
 911 | 			//fprintf(file, "v %lf %lf %lf\n", vertices[i].p.x,vertices[i].p.y,vertices[i].p.z);
 912 | 			fprintf(file, "v %g %g %g\n", vertices[i].p.x,vertices[i].p.y,vertices[i].p.z); //more compact: remove trailing zeros
 913 | 		}
 914 | 		loopi(0,triangles.size()) if(!triangles[i].deleted)
 915 | 		{
 916 | 			fprintf(file, "f %d %d %d\n", triangles[i].v[0]+1, triangles[i].v[1]+1, triangles[i].v[2]+1);
 917 | 			//fprintf(file, "f %d// %d// %d//\n", triangles[i].v[0]+1, triangles[i].v[1]+1, triangles[i].v[2]+1); //more compact: remove trailing zeros
 918 | 		}
 919 | 		fclose(file);
 920 | 	}
 921 | 
 922 |     // std::vector<VertexSTL> load_binary(uint8_t* buf) {
 923 |         // std::uint32_t num_faces;
 924 |         // std::memcpy(&num_faces, &buf[80], 4);
 925 |         // printf("num faces %d\n", num_faces);
 926 | 
 927 |         // const unsigned int num_indices = num_faces*3;
 928 |         // std::vector<VertexSTL> all_vertices(num_indices);
 929 | 
 930 |         // for (int i=0;i<num_faces;i+=1) {
 931 |             // for (int j=0;j<3;j++) {
 932 |                 // const int index = i*3+j;
 933 |                 // const int position = 84 + 12 + i*50 + j*12;
 934 |                 // std::memcpy(&all_vertices[index], &buf[position], 12);
 935 |             // }
 936 |         // }
 937 |         // return all_vertices;
 938 |     // }
 939 | 
 940 |     std::vector<VertexSTL> load_binary(const char* filename) {
 941 |         printf("loading binary stl\n");
 942 |         std::fstream fbin;
 943 |         fbin.open(filename, std::ios::in | std::ios::binary);
 944 |         fbin.seekg(80);
 945 | 
 946 |         int num_faces;
 947 |         fbin.read(reinterpret_cast<char *>(&num_faces), 4);
 948 | 
 949 |         const unsigned int num_indices = num_faces*3;
 950 | 
 951 |         size_t len = num_faces*50;
 952 |         char *ret = new char[len];
 953 |         fbin.read(ret, len);
 954 |         std::vector<VertexSTL> all_vertices(num_indices);
 955 | 
 956 |         for (int i=0;i<num_faces;i+=1) {
 957 |             for (int j=0;j<3;j++) {
 958 |                 const int index = i*3+j;
 959 |                 std::memcpy(&all_vertices[index], &ret[12 + i*50 + j*12], 12);
 960 |             }
 961 |         }
 962 | 
 963 |         fbin.close();
 964 | 
 965 |         return all_vertices;
 966 |     }
 967 | 
 968 |     std::vector<VertexSTL> load_ascii(const char* filename) {
 969 |         printf("loading ascii\n");
 970 | 
 971 |         std::vector<VertexSTL> all_vertices;
 972 | 
 973 |         std::ifstream file;
 974 |         file.open(filename);
 975 | 
 976 |         std::string line;
 977 |         while (!file.eof()) {
 978 |             std::getline(file, line);
 979 |             line = trim(line);
 980 |             if (line.rfind("vertex", 0) == 0) {
 981 |                 all_vertices.push_back(get_vector(line));
 982 |             }
 983 |         }
 984 |         file.close();
 985 | 
 986 |         return all_vertices;
 987 |     }
 988 | 
 989 |     // thanks to https://github.com/mkeeter/fstl/blob/master/src/loader.cpp
 990 |     std::vector<VertexSTL> load_stl_vertices(const char* filename) {
 991 |         std::ifstream file(filename);
 992 |         std::string line;
 993 |         std::getline(file, line);
 994 |         if (line.rfind("solid ", 0) == 0) {
 995 |             std::getline(file, line);
 996 |             line = trim(line);
 997 |             if (line.rfind("facet", 0) == 0)
 998 |             {
 999 |                 file.close();
1000 |                 return load_ascii(filename);
1001 |             }
1002 |         }
1003 |         file.close();
1004 |         return load_binary(filename);
1005 |     }
1006 | 
1007 | 
1008 |     void load_stl(const char* filename) {
1009 |         vertices.clear();
1010 |         triangles.clear();
1011 | 
1012 |         std::vector<VertexSTL> all_vertices = load_stl_vertices(filename);
1013 |         const int32_t num_indices = all_vertices.size();
1014 |         for (int c=0;c<all_vertices.size();c++)
1015 |             all_vertices[c].i = c;
1016 | 
1017 |         int32_t *indices;
1018 |         indices = (int32_t *) malloc(num_indices * sizeof(int32_t));
1019 | 
1020 |         std::sort(all_vertices.begin(), all_vertices.end());
1021 | 
1022 |         float minx =  999999;
1023 |         float miny =  999999;
1024 |         float minz =  999999;
1025 |         float maxx = -999999;
1026 |         float maxy = -999999;
1027 |         float maxz = -999999;
1028 | 
1029 |         unsigned int num_vertices = 0;
1030 |         for (int i=0;i<all_vertices.size();i++)
1031 |         {
1032 |             VertexSTL v = all_vertices[i];
1033 |             if (!num_vertices || v != all_vertices[num_vertices-1])
1034 |             {
1035 |                 all_vertices[num_vertices++] = v;
1036 |                 if (v.x < minx) minx = v.x;
1037 |                 if (v.x > maxx) maxx = v.x;
1038 |                 if (v.y < miny) miny = v.y;
1039 |                 if (v.y > maxy) maxy = v.y;
1040 |                 if (v.z < minz) minz = v.z;
1041 |                 if (v.z > maxz) maxz = v.z;
1042 |             }
1043 |             indices[v.i] = num_vertices - 1;
1044 |         }
1045 |         all_vertices.resize(num_vertices);
1046 | 
1047 |         for (int i=0;i<all_vertices.size();i++)
1048 |         {
1049 |             VertexSTL _v = all_vertices[i];
1050 |             Vertex v;
1051 |             v.p.x = _v.x; v.p.y = _v.y; v.p.z = _v.z;
1052 |             vertices.push_back(v);
1053 |         }
1054 | 
1055 |         for (int i=0;i<num_indices;i+=3) {
1056 |             Triangle t;
1057 |             t.v[0] = indices[i];
1058 |             t.v[1] = indices[i+1];
1059 |             t.v[2] = indices[i+2];
1060 |             triangles.push_back(t);
1061 |         }
1062 | 
1063 |         printf("STL num of vertices %zu num of triangles %zu\n",
1064 |                 vertices.size(), triangles.size());
1065 |     }
1066 | 
1067 |     inline void write_float(float f, FILE* file) {
1068 |         fwrite(&f,sizeof(f),1,file);
1069 |     }
1070 | 
1071 |     inline void write_vertex_stl(vec3f v, FILE* file) {
1072 |         write_float(v.x, file);
1073 |         write_float(v.y, file);
1074 |         write_float(v.z, file);
1075 |     }
1076 | 
1077 |     void write_stl(const char* filename)
1078 |     {
1079 |         FILE *file=fopen(filename, "wb");
1080 |         if (!file)
1081 |         {
1082 |             printf("write_obj: can't write data file \"%s\".\n", filename);
1083 |             exit(0);
1084 |         }
1085 |         unsigned char buffer[80] = {'T', 'I', 'G', 'E', 'R'}; // write 80 empty
1086 |         fwrite(buffer,sizeof(buffer),1,file);
1087 |         unsigned char spacer[2] = {'\0', '\0'}; // null char
1088 | 
1089 |         // loopi(0,vertices.size())
1090 |         // {
1091 |             // fprintf(file, "v %g %g %g\n", vertices[i].p.x,vertices[i].p.y,vertices[i].p.z); //more compact: remove trailing zeros
1092 |         // }
1093 |         unsigned int number_triangles = 0;
1094 |         loopi(0,triangles.size()) if(!triangles[i].deleted)
1095 |             number_triangles += 1;
1096 |         fwrite(&number_triangles,sizeof(number_triangles),1,file);
1097 | 
1098 |         loopi(0,triangles.size()) if(!triangles[i].deleted)
1099 |         {
1100 |             vec3f v0 = vertices[triangles[i].v[0]].p;
1101 |             vec3f v1 = vertices[triangles[i].v[1]].p;
1102 |             vec3f v2 = vertices[triangles[i].v[2]].p;
1103 | 
1104 |             vec3f n;
1105 |             n.cross(v1-v0,v2-v0);
1106 |             n.normalize();
1107 | 
1108 |             write_vertex_stl(n, file);
1109 |             write_vertex_stl(v0, file);
1110 |             write_vertex_stl(v1, file);
1111 |             write_vertex_stl(v2, file);
1112 | 
1113 |             fwrite(spacer,sizeof(spacer),1,file);
1114 |         }
1115 |         fclose(file);
1116 |     }
1117 | 
1118 | };
1119 | ///////////////////////////////////////////
1120 | 


--------------------------------------------------------------------------------
/a.out.js:
--------------------------------------------------------------------------------
1 | var Module=typeof Module!=="undefined"?Module:{};var objAssign=Object.assign;var moduleOverrides=objAssign({},Module);var arguments_=[];var thisProgram="./this.program";var quit_=(status,toThrow)=>{throw toThrow};var ENVIRONMENT_IS_WEB=typeof window==="object";var ENVIRONMENT_IS_WORKER=typeof importScripts==="function";var ENVIRONMENT_IS_NODE=typeof process==="object"&&typeof process.versions==="object"&&typeof process.versions.node==="string";var scriptDirectory="";function locateFile(path){if(Module["locateFile"]){return Module["locateFile"](path,scriptDirectory)}return scriptDirectory+path}var read_,readAsync,readBinary,setWindowTitle;function logExceptionOnExit(e){if(e instanceof ExitStatus)return;let toLog=e;err("exiting due to exception: "+toLog)}var fs;var nodePath;var requireNodeFS;if(ENVIRONMENT_IS_NODE){if(ENVIRONMENT_IS_WORKER){scriptDirectory=require("path").dirname(scriptDirectory)+"/"}else{scriptDirectory=__dirname+"/"}requireNodeFS=(()=>{if(!nodePath){fs=require("fs");nodePath=require("path")}});read_=function shell_read(filename,binary){requireNodeFS();filename=nodePath["normalize"](filename);return fs.readFileSync(filename,binary?null:"utf8")};readBinary=(filename=>{var ret=read_(filename,true);if(!ret.buffer){ret=new Uint8Array(ret)}return ret});readAsync=((filename,onload,onerror)=>{requireNodeFS();filename=nodePath["normalize"](filename);fs.readFile(filename,function(err,data){if(err)onerror(err);else onload(data.buffer)})});if(process["argv"].length>1){thisProgram=process["argv"][1].replace(/\\/g,"/")}arguments_=process["argv"].slice(2);if(typeof module!=="undefined"){module["exports"]=Module}process["on"]("uncaughtException",function(ex){if(!(ex instanceof ExitStatus)){throw ex}});process["on"]("unhandledRejection",function(reason){throw reason});quit_=((status,toThrow)=>{if(keepRuntimeAlive()){process["exitCode"]=status;throw toThrow}logExceptionOnExit(toThrow);process["exit"](status)});Module["inspect"]=function(){return"[Emscripten Module object]"}}else if(ENVIRONMENT_IS_WEB||ENVIRONMENT_IS_WORKER){if(ENVIRONMENT_IS_WORKER){scriptDirectory=self.location.href}else if(typeof document!=="undefined"&&document.currentScript){scriptDirectory=document.currentScript.src}if(scriptDirectory.indexOf("blob:")!==0){scriptDirectory=scriptDirectory.substr(0,scriptDirectory.replace(/[?#].*/,"").lastIndexOf("/")+1)}else{scriptDirectory=""}{read_=(url=>{var xhr=new XMLHttpRequest;xhr.open("GET",url,false);xhr.send(null);return xhr.responseText});if(ENVIRONMENT_IS_WORKER){readBinary=(url=>{var xhr=new XMLHttpRequest;xhr.open("GET",url,false);xhr.responseType="arraybuffer";xhr.send(null);return new Uint8Array(xhr.response)})}readAsync=((url,onload,onerror)=>{var xhr=new XMLHttpRequest;xhr.open("GET",url,true);xhr.responseType="arraybuffer";xhr.onload=(()=>{if(xhr.status==200||xhr.status==0&&xhr.response){onload(xhr.response);return}onerror()});xhr.onerror=onerror;xhr.send(null)})}setWindowTitle=(title=>document.title=title)}else{}var out=Module["print"]||console.log.bind(console);var err=Module["printErr"]||console.warn.bind(console);objAssign(Module,moduleOverrides);moduleOverrides=null;if(Module["arguments"])arguments_=Module["arguments"];if(Module["thisProgram"])thisProgram=Module["thisProgram"];if(Module["quit"])quit_=Module["quit"];var wasmBinary;if(Module["wasmBinary"])wasmBinary=Module["wasmBinary"];var noExitRuntime=Module["noExitRuntime"]||true;if(typeof WebAssembly!=="object"){abort("no native wasm support detected")}var wasmMemory;var ABORT=false;var EXITSTATUS;function assert(condition,text){if(!condition){abort(text)}}function getCFunc(ident){var func=Module["_"+ident];return func}function ccall(ident,returnType,argTypes,args,opts){var toC={"string":function(str){var ret=0;if(str!==null&&str!==undefined&&str!==0){var len=(str.length<<2)+1;ret=stackAlloc(len);stringToUTF8(str,ret,len)}return ret},"array":function(arr){var ret=stackAlloc(arr.length);writeArrayToMemory(arr,ret);return ret}};function convertReturnValue(ret){if(returnType==="string")return UTF8ToString(ret);if(returnType==="boolean")return Boolean(ret);return ret}var func=getCFunc(ident);var cArgs=[];var stack=0;if(args){for(var i=0;i<args.length;i++){var converter=toC[argTypes[i]];if(converter){if(stack===0)stack=stackSave();cArgs[i]=converter(args[i])}else{cArgs[i]=args[i]}}}var ret=func.apply(null,cArgs);function onDone(ret){if(stack!==0)stackRestore(stack);return convertReturnValue(ret)}ret=onDone(ret);return ret}function cwrap(ident,returnType,argTypes,opts){argTypes=argTypes||[];var numericArgs=argTypes.every(function(type){return type==="number"});var numericRet=returnType!=="string";if(numericRet&&numericArgs&&!opts){return getCFunc(ident)}return function(){return ccall(ident,returnType,argTypes,arguments,opts)}}var UTF8Decoder=typeof TextDecoder!=="undefined"?new TextDecoder("utf8"):undefined;function UTF8ArrayToString(heap,idx,maxBytesToRead){var endIdx=idx+maxBytesToRead;var endPtr=idx;while(heap[endPtr]&&!(endPtr>=endIdx))++endPtr;if(endPtr-idx>16&&heap.subarray&&UTF8Decoder){return UTF8Decoder.decode(heap.subarray(idx,endPtr))}else{var str="";while(idx<endPtr){var u0=heap[idx++];if(!(u0&128)){str+=String.fromCharCode(u0);continue}var u1=heap[idx++]&63;if((u0&224)==192){str+=String.fromCharCode((u0&31)<<6|u1);continue}var u2=heap[idx++]&63;if((u0&240)==224){u0=(u0&15)<<12|u1<<6|u2}else{u0=(u0&7)<<18|u1<<12|u2<<6|heap[idx++]&63}if(u0<65536){str+=String.fromCharCode(u0)}else{var ch=u0-65536;str+=String.fromCharCode(55296|ch>>10,56320|ch&1023)}}}return str}function UTF8ToString(ptr,maxBytesToRead){return ptr?UTF8ArrayToString(HEAPU8,ptr,maxBytesToRead):""}function stringToUTF8Array(str,heap,outIdx,maxBytesToWrite){if(!(maxBytesToWrite>0))return 0;var startIdx=outIdx;var endIdx=outIdx+maxBytesToWrite-1;for(var i=0;i<str.length;++i){var u=str.charCodeAt(i);if(u>=55296&&u<=57343){var u1=str.charCodeAt(++i);u=65536+((u&1023)<<10)|u1&1023}if(u<=127){if(outIdx>=endIdx)break;heap[outIdx++]=u}else if(u<=2047){if(outIdx+1>=endIdx)break;heap[outIdx++]=192|u>>6;heap[outIdx++]=128|u&63}else if(u<=65535){if(outIdx+2>=endIdx)break;heap[outIdx++]=224|u>>12;heap[outIdx++]=128|u>>6&63;heap[outIdx++]=128|u&63}else{if(outIdx+3>=endIdx)break;heap[outIdx++]=240|u>>18;heap[outIdx++]=128|u>>12&63;heap[outIdx++]=128|u>>6&63;heap[outIdx++]=128|u&63}}heap[outIdx]=0;return outIdx-startIdx}function stringToUTF8(str,outPtr,maxBytesToWrite){return stringToUTF8Array(str,HEAPU8,outPtr,maxBytesToWrite)}function lengthBytesUTF8(str){var len=0;for(var i=0;i<str.length;++i){var u=str.charCodeAt(i);if(u>=55296&&u<=57343)u=65536+((u&1023)<<10)|str.charCodeAt(++i)&1023;if(u<=127)++len;else if(u<=2047)len+=2;else if(u<=65535)len+=3;else len+=4}return len}function writeArrayToMemory(array,buffer){HEAP8.set(array,buffer)}function writeAsciiToMemory(str,buffer,dontAddNull){for(var i=0;i<str.length;++i){HEAP8[buffer++>>0]=str.charCodeAt(i)}if(!dontAddNull)HEAP8[buffer>>0]=0}function alignUp(x,multiple){if(x%multiple>0){x+=multiple-x%multiple}return x}var buffer,HEAP8,HEAPU8,HEAP16,HEAPU16,HEAP32,HEAPU32,HEAPF32,HEAPF64;function updateGlobalBufferAndViews(buf){buffer=buf;Module["HEAP8"]=HEAP8=new Int8Array(buf);Module["HEAP16"]=HEAP16=new Int16Array(buf);Module["HEAP32"]=HEAP32=new Int32Array(buf);Module["HEAPU8"]=HEAPU8=new Uint8Array(buf);Module["HEAPU16"]=HEAPU16=new Uint16Array(buf);Module["HEAPU32"]=HEAPU32=new Uint32Array(buf);Module["HEAPF32"]=HEAPF32=new Float32Array(buf);Module["HEAPF64"]=HEAPF64=new Float64Array(buf)}var INITIAL_MEMORY=Module["INITIAL_MEMORY"]||16777216;var wasmTable;var __ATPRERUN__=[];var __ATINIT__=[];var __ATPOSTRUN__=[];var runtimeInitialized=false;var runtimeExited=false;var runtimeKeepaliveCounter=0;function keepRuntimeAlive(){return noExitRuntime||runtimeKeepaliveCounter>0}function preRun(){if(Module["preRun"]){if(typeof Module["preRun"]=="function")Module["preRun"]=[Module["preRun"]];while(Module["preRun"].length){addOnPreRun(Module["preRun"].shift())}}callRuntimeCallbacks(__ATPRERUN__)}function initRuntime(){runtimeInitialized=true;if(!Module["noFSInit"]&&!FS.init.initialized)FS.init();FS.ignorePermissions=false;TTY.init();callRuntimeCallbacks(__ATINIT__)}function exitRuntime(){runtimeExited=true}function postRun(){if(Module["postRun"]){if(typeof Module["postRun"]=="function")Module["postRun"]=[Module["postRun"]];while(Module["postRun"].length){addOnPostRun(Module["postRun"].shift())}}callRuntimeCallbacks(__ATPOSTRUN__)}function addOnPreRun(cb){__ATPRERUN__.unshift(cb)}function addOnInit(cb){__ATINIT__.unshift(cb)}function addOnPostRun(cb){__ATPOSTRUN__.unshift(cb)}var runDependencies=0;var runDependencyWatcher=null;var dependenciesFulfilled=null;function getUniqueRunDependency(id){return id}function addRunDependency(id){runDependencies++;if(Module["monitorRunDependencies"]){Module["monitorRunDependencies"](runDependencies)}}function removeRunDependency(id){runDependencies--;if(Module["monitorRunDependencies"]){Module["monitorRunDependencies"](runDependencies)}if(runDependencies==0){if(runDependencyWatcher!==null){clearInterval(runDependencyWatcher);runDependencyWatcher=null}if(dependenciesFulfilled){var callback=dependenciesFulfilled;dependenciesFulfilled=null;callback()}}}Module["preloadedImages"]={};Module["preloadedAudios"]={};function abort(what){{if(Module["onAbort"]){Module["onAbort"](what)}}what="Aborted("+what+")";err(what);ABORT=true;EXITSTATUS=1;what+=". Build with -s ASSERTIONS=1 for more info.";var e=new WebAssembly.RuntimeError(what);throw e}var dataURIPrefix="data:application/octet-stream;base64,";function isDataURI(filename){return filename.startsWith(dataURIPrefix)}function isFileURI(filename){return filename.startsWith("file://")}var wasmBinaryFile;wasmBinaryFile="a.out.wasm";if(!isDataURI(wasmBinaryFile)){wasmBinaryFile=locateFile(wasmBinaryFile)}function getBinary(file){try{if(file==wasmBinaryFile&&wasmBinary){return new Uint8Array(wasmBinary)}if(readBinary){return readBinary(file)}else{throw"both async and sync fetching of the wasm failed"}}catch(err){abort(err)}}function getBinaryPromise(){if(!wasmBinary&&(ENVIRONMENT_IS_WEB||ENVIRONMENT_IS_WORKER)){if(typeof fetch==="function"&&!isFileURI(wasmBinaryFile)){return fetch(wasmBinaryFile,{credentials:"same-origin"}).then(function(response){if(!response["ok"]){throw"failed to load wasm binary file at '"+wasmBinaryFile+"'"}return response["arrayBuffer"]()}).catch(function(){return getBinary(wasmBinaryFile)})}else{if(readAsync){return new Promise(function(resolve,reject){readAsync(wasmBinaryFile,function(response){resolve(new Uint8Array(response))},reject)})}}}return Promise.resolve().then(function(){return getBinary(wasmBinaryFile)})}function createWasm(){var info={"a":asmLibraryArg};function receiveInstance(instance,module){var exports=instance.exports;Module["asm"]=exports;wasmMemory=Module["asm"]["r"];updateGlobalBufferAndViews(wasmMemory.buffer);wasmTable=Module["asm"]["v"];addOnInit(Module["asm"]["s"]);removeRunDependency("wasm-instantiate")}addRunDependency("wasm-instantiate");function receiveInstantiationResult(result){receiveInstance(result["instance"])}function instantiateArrayBuffer(receiver){return getBinaryPromise().then(function(binary){return WebAssembly.instantiate(binary,info)}).then(function(instance){return instance}).then(receiver,function(reason){err("failed to asynchronously prepare wasm: "+reason);abort(reason)})}function instantiateAsync(){if(!wasmBinary&&typeof WebAssembly.instantiateStreaming==="function"&&!isDataURI(wasmBinaryFile)&&!isFileURI(wasmBinaryFile)&&typeof fetch==="function"){return fetch(wasmBinaryFile,{credentials:"same-origin"}).then(function(response){var result=WebAssembly.instantiateStreaming(response,info);return result.then(receiveInstantiationResult,function(reason){err("wasm streaming compile failed: "+reason);err("falling back to ArrayBuffer instantiation");return instantiateArrayBuffer(receiveInstantiationResult)})})}else{return instantiateArrayBuffer(receiveInstantiationResult)}}if(Module["instantiateWasm"]){try{var exports=Module["instantiateWasm"](info,receiveInstance);return exports}catch(e){err("Module.instantiateWasm callback failed with error: "+e);return false}}instantiateAsync();return{}}var tempDouble;var tempI64;function callRuntimeCallbacks(callbacks){while(callbacks.length>0){var callback=callbacks.shift();if(typeof callback=="function"){callback(Module);continue}var func=callback.func;if(typeof func==="number"){if(callback.arg===undefined){getWasmTableEntry(func)()}else{getWasmTableEntry(func)(callback.arg)}}else{func(callback.arg===undefined?null:callback.arg)}}}var wasmTableMirror=[];function getWasmTableEntry(funcPtr){var func=wasmTableMirror[funcPtr];if(!func){if(funcPtr>=wasmTableMirror.length)wasmTableMirror.length=funcPtr+1;wasmTableMirror[funcPtr]=func=wasmTable.get(funcPtr)}return func}function ___cxa_allocate_exception(size){return _malloc(size+16)+16}function ExceptionInfo(excPtr){this.excPtr=excPtr;this.ptr=excPtr-16;this.set_type=function(type){HEAP32[this.ptr+4>>2]=type};this.get_type=function(){return HEAP32[this.ptr+4>>2]};this.set_destructor=function(destructor){HEAP32[this.ptr+8>>2]=destructor};this.get_destructor=function(){return HEAP32[this.ptr+8>>2]};this.set_refcount=function(refcount){HEAP32[this.ptr>>2]=refcount};this.set_caught=function(caught){caught=caught?1:0;HEAP8[this.ptr+12>>0]=caught};this.get_caught=function(){return HEAP8[this.ptr+12>>0]!=0};this.set_rethrown=function(rethrown){rethrown=rethrown?1:0;HEAP8[this.ptr+13>>0]=rethrown};this.get_rethrown=function(){return HEAP8[this.ptr+13>>0]!=0};this.init=function(type,destructor){this.set_type(type);this.set_destructor(destructor);this.set_refcount(0);this.set_caught(false);this.set_rethrown(false)};this.add_ref=function(){var value=HEAP32[this.ptr>>2];HEAP32[this.ptr>>2]=value+1};this.release_ref=function(){var prev=HEAP32[this.ptr>>2];HEAP32[this.ptr>>2]=prev-1;return prev===1}}var exceptionLast=0;var uncaughtExceptionCount=0;function ___cxa_throw(ptr,type,destructor){var info=new ExceptionInfo(ptr);info.init(type,destructor);exceptionLast=ptr;uncaughtExceptionCount++;throw ptr}function setErrNo(value){HEAP32[___errno_location()>>2]=value;return value}var PATH={splitPath:function(filename){var splitPathRe=/^(\/?|)([\s\S]*?)((?:\.{1,2}|[^\/]+?|)(\.[^.\/]*|))(?:[\/]*)$/;return splitPathRe.exec(filename).slice(1)},normalizeArray:function(parts,allowAboveRoot){var up=0;for(var i=parts.length-1;i>=0;i--){var last=parts[i];if(last==="."){parts.splice(i,1)}else if(last===".."){parts.splice(i,1);up++}else if(up){parts.splice(i,1);up--}}if(allowAboveRoot){for(;up;up--){parts.unshift("..")}}return parts},normalize:function(path){var isAbsolute=path.charAt(0)==="/",trailingSlash=path.substr(-1)==="/";path=PATH.normalizeArray(path.split("/").filter(function(p){return!!p}),!isAbsolute).join("/");if(!path&&!isAbsolute){path="."}if(path&&trailingSlash){path+="/"}return(isAbsolute?"/":"")+path},dirname:function(path){var result=PATH.splitPath(path),root=result[0],dir=result[1];if(!root&&!dir){return"."}if(dir){dir=dir.substr(0,dir.length-1)}return root+dir},basename:function(path){if(path==="/")return"/";path=PATH.normalize(path);path=path.replace(/\/$/,"");var lastSlash=path.lastIndexOf("/");if(lastSlash===-1)return path;return path.substr(lastSlash+1)},extname:function(path){return PATH.splitPath(path)[3]},join:function(){var paths=Array.prototype.slice.call(arguments,0);return PATH.normalize(paths.join("/"))},join2:function(l,r){return PATH.normalize(l+"/"+r)}};function getRandomDevice(){if(typeof crypto==="object"&&typeof crypto["getRandomValues"]==="function"){var randomBuffer=new Uint8Array(1);return function(){crypto.getRandomValues(randomBuffer);return randomBuffer[0]}}else if(ENVIRONMENT_IS_NODE){try{var crypto_module=require("crypto");return function(){return crypto_module["randomBytes"](1)[0]}}catch(e){}}return function(){abort("randomDevice")}}var PATH_FS={resolve:function(){var resolvedPath="",resolvedAbsolute=false;for(var i=arguments.length-1;i>=-1&&!resolvedAbsolute;i--){var path=i>=0?arguments[i]:FS.cwd();if(typeof path!=="string"){throw new TypeError("Arguments to path.resolve must be strings")}else if(!path){return""}resolvedPath=path+"/"+resolvedPath;resolvedAbsolute=path.charAt(0)==="/"}resolvedPath=PATH.normalizeArray(resolvedPath.split("/").filter(function(p){return!!p}),!resolvedAbsolute).join("/");return(resolvedAbsolute?"/":"")+resolvedPath||"."},relative:function(from,to){from=PATH_FS.resolve(from).substr(1);to=PATH_FS.resolve(to).substr(1);function trim(arr){var start=0;for(;start<arr.length;start++){if(arr[start]!=="")break}var end=arr.length-1;for(;end>=0;end--){if(arr[end]!=="")break}if(start>end)return[];return arr.slice(start,end-start+1)}var fromParts=trim(from.split("/"));var toParts=trim(to.split("/"));var length=Math.min(fromParts.length,toParts.length);var samePartsLength=length;for(var i=0;i<length;i++){if(fromParts[i]!==toParts[i]){samePartsLength=i;break}}var outputParts=[];for(var i=samePartsLength;i<fromParts.length;i++){outputParts.push("..")}outputParts=outputParts.concat(toParts.slice(samePartsLength));return outputParts.join("/")}};var TTY={ttys:[],init:function(){},shutdown:function(){},register:function(dev,ops){TTY.ttys[dev]={input:[],output:[],ops:ops};FS.registerDevice(dev,TTY.stream_ops)},stream_ops:{open:function(stream){var tty=TTY.ttys[stream.node.rdev];if(!tty){throw new FS.ErrnoError(43)}stream.tty=tty;stream.seekable=false},close:function(stream){stream.tty.ops.flush(stream.tty)},flush:function(stream){stream.tty.ops.flush(stream.tty)},read:function(stream,buffer,offset,length,pos){if(!stream.tty||!stream.tty.ops.get_char){throw new FS.ErrnoError(60)}var bytesRead=0;for(var i=0;i<length;i++){var result;try{result=stream.tty.ops.get_char(stream.tty)}catch(e){throw new FS.ErrnoError(29)}if(result===undefined&&bytesRead===0){throw new FS.ErrnoError(6)}if(result===null||result===undefined)break;bytesRead++;buffer[offset+i]=result}if(bytesRead){stream.node.timestamp=Date.now()}return bytesRead},write:function(stream,buffer,offset,length,pos){if(!stream.tty||!stream.tty.ops.put_char){throw new FS.ErrnoError(60)}try{for(var i=0;i<length;i++){stream.tty.ops.put_char(stream.tty,buffer[offset+i])}}catch(e){throw new FS.ErrnoError(29)}if(length){stream.node.timestamp=Date.now()}return i}},default_tty_ops:{get_char:function(tty){if(!tty.input.length){var result=null;if(ENVIRONMENT_IS_NODE){var BUFSIZE=256;var buf=Buffer.alloc(BUFSIZE);var bytesRead=0;try{bytesRead=fs.readSync(process.stdin.fd,buf,0,BUFSIZE,null)}catch(e){if(e.toString().includes("EOF"))bytesRead=0;else throw e}if(bytesRead>0){result=buf.slice(0,bytesRead).toString("utf-8")}else{result=null}}else if(typeof window!="undefined"&&typeof window.prompt=="function"){result=window.prompt("Input: ");if(result!==null){result+="\n"}}else if(typeof readline=="function"){result=readline();if(result!==null){result+="\n"}}if(!result){return null}tty.input=intArrayFromString(result,true)}return tty.input.shift()},put_char:function(tty,val){if(val===null||val===10){out(UTF8ArrayToString(tty.output,0));tty.output=[]}else{if(val!=0)tty.output.push(val)}},flush:function(tty){if(tty.output&&tty.output.length>0){out(UTF8ArrayToString(tty.output,0));tty.output=[]}}},default_tty1_ops:{put_char:function(tty,val){if(val===null||val===10){err(UTF8ArrayToString(tty.output,0));tty.output=[]}else{if(val!=0)tty.output.push(val)}},flush:function(tty){if(tty.output&&tty.output.length>0){err(UTF8ArrayToString(tty.output,0));tty.output=[]}}}};function mmapAlloc(size){abort()}var MEMFS={ops_table:null,mount:function(mount){return MEMFS.createNode(null,"/",16384|511,0)},createNode:function(parent,name,mode,dev){if(FS.isBlkdev(mode)||FS.isFIFO(mode)){throw new FS.ErrnoError(63)}if(!MEMFS.ops_table){MEMFS.ops_table={dir:{node:{getattr:MEMFS.node_ops.getattr,setattr:MEMFS.node_ops.setattr,lookup:MEMFS.node_ops.lookup,mknod:MEMFS.node_ops.mknod,rename:MEMFS.node_ops.rename,unlink:MEMFS.node_ops.unlink,rmdir:MEMFS.node_ops.rmdir,readdir:MEMFS.node_ops.readdir,symlink:MEMFS.node_ops.symlink},stream:{llseek:MEMFS.stream_ops.llseek}},file:{node:{getattr:MEMFS.node_ops.getattr,setattr:MEMFS.node_ops.setattr},stream:{llseek:MEMFS.stream_ops.llseek,read:MEMFS.stream_ops.read,write:MEMFS.stream_ops.write,allocate:MEMFS.stream_ops.allocate,mmap:MEMFS.stream_ops.mmap,msync:MEMFS.stream_ops.msync}},link:{node:{getattr:MEMFS.node_ops.getattr,setattr:MEMFS.node_ops.setattr,readlink:MEMFS.node_ops.readlink},stream:{}},chrdev:{node:{getattr:MEMFS.node_ops.getattr,setattr:MEMFS.node_ops.setattr},stream:FS.chrdev_stream_ops}}}var node=FS.createNode(parent,name,mode,dev);if(FS.isDir(node.mode)){node.node_ops=MEMFS.ops_table.dir.node;node.stream_ops=MEMFS.ops_table.dir.stream;node.contents={}}else if(FS.isFile(node.mode)){node.node_ops=MEMFS.ops_table.file.node;node.stream_ops=MEMFS.ops_table.file.stream;node.usedBytes=0;node.contents=null}else if(FS.isLink(node.mode)){node.node_ops=MEMFS.ops_table.link.node;node.stream_ops=MEMFS.ops_table.link.stream}else if(FS.isChrdev(node.mode)){node.node_ops=MEMFS.ops_table.chrdev.node;node.stream_ops=MEMFS.ops_table.chrdev.stream}node.timestamp=Date.now();if(parent){parent.contents[name]=node;parent.timestamp=node.timestamp}return node},getFileDataAsTypedArray:function(node){if(!node.contents)return new Uint8Array(0);if(node.contents.subarray)return node.contents.subarray(0,node.usedBytes);return new Uint8Array(node.contents)},expandFileStorage:function(node,newCapacity){var prevCapacity=node.contents?node.contents.length:0;if(prevCapacity>=newCapacity)return;var CAPACITY_DOUBLING_MAX=1024*1024;newCapacity=Math.max(newCapacity,prevCapacity*(prevCapacity<CAPACITY_DOUBLING_MAX?2:1.125)>>>0);if(prevCapacity!=0)newCapacity=Math.max(newCapacity,256);var oldContents=node.contents;node.contents=new Uint8Array(newCapacity);if(node.usedBytes>0)node.contents.set(oldContents.subarray(0,node.usedBytes),0)},resizeFileStorage:function(node,newSize){if(node.usedBytes==newSize)return;if(newSize==0){node.contents=null;node.usedBytes=0}else{var oldContents=node.contents;node.contents=new Uint8Array(newSize);if(oldContents){node.contents.set(oldContents.subarray(0,Math.min(newSize,node.usedBytes)))}node.usedBytes=newSize}},node_ops:{getattr:function(node){var attr={};attr.dev=FS.isChrdev(node.mode)?node.id:1;attr.ino=node.id;attr.mode=node.mode;attr.nlink=1;attr.uid=0;attr.gid=0;attr.rdev=node.rdev;if(FS.isDir(node.mode)){attr.size=4096}else if(FS.isFile(node.mode)){attr.size=node.usedBytes}else if(FS.isLink(node.mode)){attr.size=node.link.length}else{attr.size=0}attr.atime=new Date(node.timestamp);attr.mtime=new Date(node.timestamp);attr.ctime=new Date(node.timestamp);attr.blksize=4096;attr.blocks=Math.ceil(attr.size/attr.blksize);return attr},setattr:function(node,attr){if(attr.mode!==undefined){node.mode=attr.mode}if(attr.timestamp!==undefined){node.timestamp=attr.timestamp}if(attr.size!==undefined){MEMFS.resizeFileStorage(node,attr.size)}},lookup:function(parent,name){throw FS.genericErrors[44]},mknod:function(parent,name,mode,dev){return MEMFS.createNode(parent,name,mode,dev)},rename:function(old_node,new_dir,new_name){if(FS.isDir(old_node.mode)){var new_node;try{new_node=FS.lookupNode(new_dir,new_name)}catch(e){}if(new_node){for(var i in new_node.contents){throw new FS.ErrnoError(55)}}}delete old_node.parent.contents[old_node.name];old_node.parent.timestamp=Date.now();old_node.name=new_name;new_dir.contents[new_name]=old_node;new_dir.timestamp=old_node.parent.timestamp;old_node.parent=new_dir},unlink:function(parent,name){delete parent.contents[name];parent.timestamp=Date.now()},rmdir:function(parent,name){var node=FS.lookupNode(parent,name);for(var i in node.contents){throw new FS.ErrnoError(55)}delete parent.contents[name];parent.timestamp=Date.now()},readdir:function(node){var entries=[".",".."];for(var key in node.contents){if(!node.contents.hasOwnProperty(key)){continue}entries.push(key)}return entries},symlink:function(parent,newname,oldpath){var node=MEMFS.createNode(parent,newname,511|40960,0);node.link=oldpath;return node},readlink:function(node){if(!FS.isLink(node.mode)){throw new FS.ErrnoError(28)}return node.link}},stream_ops:{read:function(stream,buffer,offset,length,position){var contents=stream.node.contents;if(position>=stream.node.usedBytes)return 0;var size=Math.min(stream.node.usedBytes-position,length);if(size>8&&contents.subarray){buffer.set(contents.subarray(position,position+size),offset)}else{for(var i=0;i<size;i++)buffer[offset+i]=contents[position+i]}return size},write:function(stream,buffer,offset,length,position,canOwn){if(buffer.buffer===HEAP8.buffer){canOwn=false}if(!length)return 0;var node=stream.node;node.timestamp=Date.now();if(buffer.subarray&&(!node.contents||node.contents.subarray)){if(canOwn){node.contents=buffer.subarray(offset,offset+length);node.usedBytes=length;return length}else if(node.usedBytes===0&&position===0){node.contents=buffer.slice(offset,offset+length);node.usedBytes=length;return length}else if(position+length<=node.usedBytes){node.contents.set(buffer.subarray(offset,offset+length),position);return length}}MEMFS.expandFileStorage(node,position+length);if(node.contents.subarray&&buffer.subarray){node.contents.set(buffer.subarray(offset,offset+length),position)}else{for(var i=0;i<length;i++){node.contents[position+i]=buffer[offset+i]}}node.usedBytes=Math.max(node.usedBytes,position+length);return length},llseek:function(stream,offset,whence){var position=offset;if(whence===1){position+=stream.position}else if(whence===2){if(FS.isFile(stream.node.mode)){position+=stream.node.usedBytes}}if(position<0){throw new FS.ErrnoError(28)}return position},allocate:function(stream,offset,length){MEMFS.expandFileStorage(stream.node,offset+length);stream.node.usedBytes=Math.max(stream.node.usedBytes,offset+length)},mmap:function(stream,address,length,position,prot,flags){if(address!==0){throw new FS.ErrnoError(28)}if(!FS.isFile(stream.node.mode)){throw new FS.ErrnoError(43)}var ptr;var allocated;var contents=stream.node.contents;if(!(flags&2)&&contents.buffer===buffer){allocated=false;ptr=contents.byteOffset}else{if(position>0||position+length<contents.length){if(contents.subarray){contents=contents.subarray(position,position+length)}else{contents=Array.prototype.slice.call(contents,position,position+length)}}allocated=true;ptr=mmapAlloc(length);if(!ptr){throw new FS.ErrnoError(48)}HEAP8.set(contents,ptr)}return{ptr:ptr,allocated:allocated}},msync:function(stream,buffer,offset,length,mmapFlags){if(!FS.isFile(stream.node.mode)){throw new FS.ErrnoError(43)}if(mmapFlags&2){return 0}var bytesWritten=MEMFS.stream_ops.write(stream,buffer,0,length,offset,false);return 0}}};function asyncLoad(url,onload,onerror,noRunDep){var dep=!noRunDep?getUniqueRunDependency("al "+url):"";readAsync(url,function(arrayBuffer){assert(arrayBuffer,'Loading data file "'+url+'" failed (no arrayBuffer).');onload(new Uint8Array(arrayBuffer));if(dep)removeRunDependency(dep)},function(event){if(onerror){onerror()}else{throw'Loading data file "'+url+'" failed.'}});if(dep)addRunDependency(dep)}var FS={root:null,mounts:[],devices:{},streams:[],nextInode:1,nameTable:null,currentPath:"/",initialized:false,ignorePermissions:true,ErrnoError:null,genericErrors:{},filesystems:null,syncFSRequests:0,lookupPath:(path,opts={})=>{path=PATH_FS.resolve(FS.cwd(),path);if(!path)return{path:"",node:null};var defaults={follow_mount:true,recurse_count:0};for(var key in defaults){if(opts[key]===undefined){opts[key]=defaults[key]}}if(opts.recurse_count>8){throw new FS.ErrnoError(32)}var parts=PATH.normalizeArray(path.split("/").filter(p=>!!p),false);var current=FS.root;var current_path="/";for(var i=0;i<parts.length;i++){var islast=i===parts.length-1;if(islast&&opts.parent){break}current=FS.lookupNode(current,parts[i]);current_path=PATH.join2(current_path,parts[i]);if(FS.isMountpoint(current)){if(!islast||islast&&opts.follow_mount){current=current.mounted.root}}if(!islast||opts.follow){var count=0;while(FS.isLink(current.mode)){var link=FS.readlink(current_path);current_path=PATH_FS.resolve(PATH.dirname(current_path),link);var lookup=FS.lookupPath(current_path,{recurse_count:opts.recurse_count});current=lookup.node;if(count++>40){throw new FS.ErrnoError(32)}}}}return{path:current_path,node:current}},getPath:node=>{var path;while(true){if(FS.isRoot(node)){var mount=node.mount.mountpoint;if(!path)return mount;return mount[mount.length-1]!=="/"?mount+"/"+path:mount+path}path=path?node.name+"/"+path:node.name;node=node.parent}},hashName:(parentid,name)=>{var hash=0;for(var i=0;i<name.length;i++){hash=(hash<<5)-hash+name.charCodeAt(i)|0}return(parentid+hash>>>0)%FS.nameTable.length},hashAddNode:node=>{var hash=FS.hashName(node.parent.id,node.name);node.name_next=FS.nameTable[hash];FS.nameTable[hash]=node},hashRemoveNode:node=>{var hash=FS.hashName(node.parent.id,node.name);if(FS.nameTable[hash]===node){FS.nameTable[hash]=node.name_next}else{var current=FS.nameTable[hash];while(current){if(current.name_next===node){current.name_next=node.name_next;break}current=current.name_next}}},lookupNode:(parent,name)=>{var errCode=FS.mayLookup(parent);if(errCode){throw new FS.ErrnoError(errCode,parent)}var hash=FS.hashName(parent.id,name);for(var node=FS.nameTable[hash];node;node=node.name_next){var nodeName=node.name;if(node.parent.id===parent.id&&nodeName===name){return node}}return FS.lookup(parent,name)},createNode:(parent,name,mode,rdev)=>{var node=new FS.FSNode(parent,name,mode,rdev);FS.hashAddNode(node);return node},destroyNode:node=>{FS.hashRemoveNode(node)},isRoot:node=>{return node===node.parent},isMountpoint:node=>{return!!node.mounted},isFile:mode=>{return(mode&61440)===32768},isDir:mode=>{return(mode&61440)===16384},isLink:mode=>{return(mode&61440)===40960},isChrdev:mode=>{return(mode&61440)===8192},isBlkdev:mode=>{return(mode&61440)===24576},isFIFO:mode=>{return(mode&61440)===4096},isSocket:mode=>{return(mode&49152)===49152},flagModes:{"r":0,"r+":2,"w":577,"w+":578,"a":1089,"a+":1090},modeStringToFlags:str=>{var flags=FS.flagModes[str];if(typeof flags==="undefined"){throw new Error("Unknown file open mode: "+str)}return flags},flagsToPermissionString:flag=>{var perms=["r","w","rw"][flag&3];if(flag&512){perms+="w"}return perms},nodePermissions:(node,perms)=>{if(FS.ignorePermissions){return 0}if(perms.includes("r")&&!(node.mode&292)){return 2}else if(perms.includes("w")&&!(node.mode&146)){return 2}else if(perms.includes("x")&&!(node.mode&73)){return 2}return 0},mayLookup:dir=>{var errCode=FS.nodePermissions(dir,"x");if(errCode)return errCode;if(!dir.node_ops.lookup)return 2;return 0},mayCreate:(dir,name)=>{try{var node=FS.lookupNode(dir,name);return 20}catch(e){}return FS.nodePermissions(dir,"wx")},mayDelete:(dir,name,isdir)=>{var node;try{node=FS.lookupNode(dir,name)}catch(e){return e.errno}var errCode=FS.nodePermissions(dir,"wx");if(errCode){return errCode}if(isdir){if(!FS.isDir(node.mode)){return 54}if(FS.isRoot(node)||FS.getPath(node)===FS.cwd()){return 10}}else{if(FS.isDir(node.mode)){return 31}}return 0},mayOpen:(node,flags)=>{if(!node){return 44}if(FS.isLink(node.mode)){return 32}else if(FS.isDir(node.mode)){if(FS.flagsToPermissionString(flags)!=="r"||flags&512){return 31}}return FS.nodePermissions(node,FS.flagsToPermissionString(flags))},MAX_OPEN_FDS:4096,nextfd:(fd_start=0,fd_end=FS.MAX_OPEN_FDS)=>{for(var fd=fd_start;fd<=fd_end;fd++){if(!FS.streams[fd]){return fd}}throw new FS.ErrnoError(33)},getStream:fd=>FS.streams[fd],createStream:(stream,fd_start,fd_end)=>{if(!FS.FSStream){FS.FSStream=function(){};FS.FSStream.prototype={object:{get:function(){return this.node},set:function(val){this.node=val}},isRead:{get:function(){return(this.flags&2097155)!==1}},isWrite:{get:function(){return(this.flags&2097155)!==0}},isAppend:{get:function(){return this.flags&1024}}}}var newStream=new FS.FSStream;for(var p in stream){newStream[p]=stream[p]}stream=newStream;var fd=FS.nextfd(fd_start,fd_end);stream.fd=fd;FS.streams[fd]=stream;return stream},closeStream:fd=>{FS.streams[fd]=null},chrdev_stream_ops:{open:stream=>{var device=FS.getDevice(stream.node.rdev);stream.stream_ops=device.stream_ops;if(stream.stream_ops.open){stream.stream_ops.open(stream)}},llseek:()=>{throw new FS.ErrnoError(70)}},major:dev=>dev>>8,minor:dev=>dev&255,makedev:(ma,mi)=>ma<<8|mi,registerDevice:(dev,ops)=>{FS.devices[dev]={stream_ops:ops}},getDevice:dev=>FS.devices[dev],getMounts:mount=>{var mounts=[];var check=[mount];while(check.length){var m=check.pop();mounts.push(m);check.push.apply(check,m.mounts)}return mounts},syncfs:(populate,callback)=>{if(typeof populate==="function"){callback=populate;populate=false}FS.syncFSRequests++;if(FS.syncFSRequests>1){err("warning: "+FS.syncFSRequests+" FS.syncfs operations in flight at once, probably just doing extra work")}var mounts=FS.getMounts(FS.root.mount);var completed=0;function doCallback(errCode){FS.syncFSRequests--;return callback(errCode)}function done(errCode){if(errCode){if(!done.errored){done.errored=true;return doCallback(errCode)}return}if(++completed>=mounts.length){doCallback(null)}}mounts.forEach(mount=>{if(!mount.type.syncfs){return done(null)}mount.type.syncfs(mount,populate,done)})},mount:(type,opts,mountpoint)=>{var root=mountpoint==="/";var pseudo=!mountpoint;var node;if(root&&FS.root){throw new FS.ErrnoError(10)}else if(!root&&!pseudo){var lookup=FS.lookupPath(mountpoint,{follow_mount:false});mountpoint=lookup.path;node=lookup.node;if(FS.isMountpoint(node)){throw new FS.ErrnoError(10)}if(!FS.isDir(node.mode)){throw new FS.ErrnoError(54)}}var mount={type:type,opts:opts,mountpoint:mountpoint,mounts:[]};var mountRoot=type.mount(mount);mountRoot.mount=mount;mount.root=mountRoot;if(root){FS.root=mountRoot}else if(node){node.mounted=mount;if(node.mount){node.mount.mounts.push(mount)}}return mountRoot},unmount:mountpoint=>{var lookup=FS.lookupPath(mountpoint,{follow_mount:false});if(!FS.isMountpoint(lookup.node)){throw new FS.ErrnoError(28)}var node=lookup.node;var mount=node.mounted;var mounts=FS.getMounts(mount);Object.keys(FS.nameTable).forEach(hash=>{var current=FS.nameTable[hash];while(current){var next=current.name_next;if(mounts.includes(current.mount)){FS.destroyNode(current)}current=next}});node.mounted=null;var idx=node.mount.mounts.indexOf(mount);node.mount.mounts.splice(idx,1)},lookup:(parent,name)=>{return parent.node_ops.lookup(parent,name)},mknod:(path,mode,dev)=>{var lookup=FS.lookupPath(path,{parent:true});var parent=lookup.node;var name=PATH.basename(path);if(!name||name==="."||name===".."){throw new FS.ErrnoError(28)}var errCode=FS.mayCreate(parent,name);if(errCode){throw new FS.ErrnoError(errCode)}if(!parent.node_ops.mknod){throw new FS.ErrnoError(63)}return parent.node_ops.mknod(parent,name,mode,dev)},create:(path,mode)=>{mode=mode!==undefined?mode:438;mode&=4095;mode|=32768;return FS.mknod(path,mode,0)},mkdir:(path,mode)=>{mode=mode!==undefined?mode:511;mode&=511|512;mode|=16384;return FS.mknod(path,mode,0)},mkdirTree:(path,mode)=>{var dirs=path.split("/");var d="";for(var i=0;i<dirs.length;++i){if(!dirs[i])continue;d+="/"+dirs[i];try{FS.mkdir(d,mode)}catch(e){if(e.errno!=20)throw e}}},mkdev:(path,mode,dev)=>{if(typeof dev==="undefined"){dev=mode;mode=438}mode|=8192;return FS.mknod(path,mode,dev)},symlink:(oldpath,newpath)=>{if(!PATH_FS.resolve(oldpath)){throw new FS.ErrnoError(44)}var lookup=FS.lookupPath(newpath,{parent:true});var parent=lookup.node;if(!parent){throw new FS.ErrnoError(44)}var newname=PATH.basename(newpath);var errCode=FS.mayCreate(parent,newname);if(errCode){throw new FS.ErrnoError(errCode)}if(!parent.node_ops.symlink){throw new FS.ErrnoError(63)}return parent.node_ops.symlink(parent,newname,oldpath)},rename:(old_path,new_path)=>{var old_dirname=PATH.dirname(old_path);var new_dirname=PATH.dirname(new_path);var old_name=PATH.basename(old_path);var new_name=PATH.basename(new_path);var lookup,old_dir,new_dir;lookup=FS.lookupPath(old_path,{parent:true});old_dir=lookup.node;lookup=FS.lookupPath(new_path,{parent:true});new_dir=lookup.node;if(!old_dir||!new_dir)throw new FS.ErrnoError(44);if(old_dir.mount!==new_dir.mount){throw new FS.ErrnoError(75)}var old_node=FS.lookupNode(old_dir,old_name);var relative=PATH_FS.relative(old_path,new_dirname);if(relative.charAt(0)!=="."){throw new FS.ErrnoError(28)}relative=PATH_FS.relative(new_path,old_dirname);if(relative.charAt(0)!=="."){throw new FS.ErrnoError(55)}var new_node;try{new_node=FS.lookupNode(new_dir,new_name)}catch(e){}if(old_node===new_node){return}var isdir=FS.isDir(old_node.mode);var errCode=FS.mayDelete(old_dir,old_name,isdir);if(errCode){throw new FS.ErrnoError(errCode)}errCode=new_node?FS.mayDelete(new_dir,new_name,isdir):FS.mayCreate(new_dir,new_name);if(errCode){throw new FS.ErrnoError(errCode)}if(!old_dir.node_ops.rename){throw new FS.ErrnoError(63)}if(FS.isMountpoint(old_node)||new_node&&FS.isMountpoint(new_node)){throw new FS.ErrnoError(10)}if(new_dir!==old_dir){errCode=FS.nodePermissions(old_dir,"w");if(errCode){throw new FS.ErrnoError(errCode)}}FS.hashRemoveNode(old_node);try{old_dir.node_ops.rename(old_node,new_dir,new_name)}catch(e){throw e}finally{FS.hashAddNode(old_node)}},rmdir:path=>{var lookup=FS.lookupPath(path,{parent:true});var parent=lookup.node;var name=PATH.basename(path);var node=FS.lookupNode(parent,name);var errCode=FS.mayDelete(parent,name,true);if(errCode){throw new FS.ErrnoError(errCode)}if(!parent.node_ops.rmdir){throw new FS.ErrnoError(63)}if(FS.isMountpoint(node)){throw new FS.ErrnoError(10)}parent.node_ops.rmdir(parent,name);FS.destroyNode(node)},readdir:path=>{var lookup=FS.lookupPath(path,{follow:true});var node=lookup.node;if(!node.node_ops.readdir){throw new FS.ErrnoError(54)}return node.node_ops.readdir(node)},unlink:path=>{var lookup=FS.lookupPath(path,{parent:true});var parent=lookup.node;if(!parent){throw new FS.ErrnoError(44)}var name=PATH.basename(path);var node=FS.lookupNode(parent,name);var errCode=FS.mayDelete(parent,name,false);if(errCode){throw new FS.ErrnoError(errCode)}if(!parent.node_ops.unlink){throw new FS.ErrnoError(63)}if(FS.isMountpoint(node)){throw new FS.ErrnoError(10)}parent.node_ops.unlink(parent,name);FS.destroyNode(node)},readlink:path=>{var lookup=FS.lookupPath(path);var link=lookup.node;if(!link){throw new FS.ErrnoError(44)}if(!link.node_ops.readlink){throw new FS.ErrnoError(28)}return PATH_FS.resolve(FS.getPath(link.parent),link.node_ops.readlink(link))},stat:(path,dontFollow)=>{var lookup=FS.lookupPath(path,{follow:!dontFollow});var node=lookup.node;if(!node){throw new FS.ErrnoError(44)}if(!node.node_ops.getattr){throw new FS.ErrnoError(63)}return node.node_ops.getattr(node)},lstat:path=>{return FS.stat(path,true)},chmod:(path,mode,dontFollow)=>{var node;if(typeof path==="string"){var lookup=FS.lookupPath(path,{follow:!dontFollow});node=lookup.node}else{node=path}if(!node.node_ops.setattr){throw new FS.ErrnoError(63)}node.node_ops.setattr(node,{mode:mode&4095|node.mode&~4095,timestamp:Date.now()})},lchmod:(path,mode)=>{FS.chmod(path,mode,true)},fchmod:(fd,mode)=>{var stream=FS.getStream(fd);if(!stream){throw new FS.ErrnoError(8)}FS.chmod(stream.node,mode)},chown:(path,uid,gid,dontFollow)=>{var node;if(typeof path==="string"){var lookup=FS.lookupPath(path,{follow:!dontFollow});node=lookup.node}else{node=path}if(!node.node_ops.setattr){throw new FS.ErrnoError(63)}node.node_ops.setattr(node,{timestamp:Date.now()})},lchown:(path,uid,gid)=>{FS.chown(path,uid,gid,true)},fchown:(fd,uid,gid)=>{var stream=FS.getStream(fd);if(!stream){throw new FS.ErrnoError(8)}FS.chown(stream.node,uid,gid)},truncate:(path,len)=>{if(len<0){throw new FS.ErrnoError(28)}var node;if(typeof path==="string"){var lookup=FS.lookupPath(path,{follow:true});node=lookup.node}else{node=path}if(!node.node_ops.setattr){throw new FS.ErrnoError(63)}if(FS.isDir(node.mode)){throw new FS.ErrnoError(31)}if(!FS.isFile(node.mode)){throw new FS.ErrnoError(28)}var errCode=FS.nodePermissions(node,"w");if(errCode){throw new FS.ErrnoError(errCode)}node.node_ops.setattr(node,{size:len,timestamp:Date.now()})},ftruncate:(fd,len)=>{var stream=FS.getStream(fd);if(!stream){throw new FS.ErrnoError(8)}if((stream.flags&2097155)===0){throw new FS.ErrnoError(28)}FS.truncate(stream.node,len)},utime:(path,atime,mtime)=>{var lookup=FS.lookupPath(path,{follow:true});var node=lookup.node;node.node_ops.setattr(node,{timestamp:Math.max(atime,mtime)})},open:(path,flags,mode,fd_start,fd_end)=>{if(path===""){throw new FS.ErrnoError(44)}flags=typeof flags==="string"?FS.modeStringToFlags(flags):flags;mode=typeof mode==="undefined"?438:mode;if(flags&64){mode=mode&4095|32768}else{mode=0}var node;if(typeof path==="object"){node=path}else{path=PATH.normalize(path);try{var lookup=FS.lookupPath(path,{follow:!(flags&131072)});node=lookup.node}catch(e){}}var created=false;if(flags&64){if(node){if(flags&128){throw new FS.ErrnoError(20)}}else{node=FS.mknod(path,mode,0);created=true}}if(!node){throw new FS.ErrnoError(44)}if(FS.isChrdev(node.mode)){flags&=~512}if(flags&65536&&!FS.isDir(node.mode)){throw new FS.ErrnoError(54)}if(!created){var errCode=FS.mayOpen(node,flags);if(errCode){throw new FS.ErrnoError(errCode)}}if(flags&512){FS.truncate(node,0)}flags&=~(128|512|131072);var stream=FS.createStream({node:node,path:FS.getPath(node),flags:flags,seekable:true,position:0,stream_ops:node.stream_ops,ungotten:[],error:false},fd_start,fd_end);if(stream.stream_ops.open){stream.stream_ops.open(stream)}if(Module["logReadFiles"]&&!(flags&1)){if(!FS.readFiles)FS.readFiles={};if(!(path in FS.readFiles)){FS.readFiles[path]=1}}return stream},close:stream=>{if(FS.isClosed(stream)){throw new FS.ErrnoError(8)}if(stream.getdents)stream.getdents=null;try{if(stream.stream_ops.close){stream.stream_ops.close(stream)}}catch(e){throw e}finally{FS.closeStream(stream.fd)}stream.fd=null},isClosed:stream=>{return stream.fd===null},llseek:(stream,offset,whence)=>{if(FS.isClosed(stream)){throw new FS.ErrnoError(8)}if(!stream.seekable||!stream.stream_ops.llseek){throw new FS.ErrnoError(70)}if(whence!=0&&whence!=1&&whence!=2){throw new FS.ErrnoError(28)}stream.position=stream.stream_ops.llseek(stream,offset,whence);stream.ungotten=[];return stream.position},read:(stream,buffer,offset,length,position)=>{if(length<0||position<0){throw new FS.ErrnoError(28)}if(FS.isClosed(stream)){throw new FS.ErrnoError(8)}if((stream.flags&2097155)===1){throw new FS.ErrnoError(8)}if(FS.isDir(stream.node.mode)){throw new FS.ErrnoError(31)}if(!stream.stream_ops.read){throw new FS.ErrnoError(28)}var seeking=typeof position!=="undefined";if(!seeking){position=stream.position}else if(!stream.seekable){throw new FS.ErrnoError(70)}var bytesRead=stream.stream_ops.read(stream,buffer,offset,length,position);if(!seeking)stream.position+=bytesRead;return bytesRead},write:(stream,buffer,offset,length,position,canOwn)=>{if(length<0||position<0){throw new FS.ErrnoError(28)}if(FS.isClosed(stream)){throw new FS.ErrnoError(8)}if((stream.flags&2097155)===0){throw new FS.ErrnoError(8)}if(FS.isDir(stream.node.mode)){throw new FS.ErrnoError(31)}if(!stream.stream_ops.write){throw new FS.ErrnoError(28)}if(stream.seekable&&stream.flags&1024){FS.llseek(stream,0,2)}var seeking=typeof position!=="undefined";if(!seeking){position=stream.position}else if(!stream.seekable){throw new FS.ErrnoError(70)}var bytesWritten=stream.stream_ops.write(stream,buffer,offset,length,position,canOwn);if(!seeking)stream.position+=bytesWritten;return bytesWritten},allocate:(stream,offset,length)=>{if(FS.isClosed(stream)){throw new FS.ErrnoError(8)}if(offset<0||length<=0){throw new FS.ErrnoError(28)}if((stream.flags&2097155)===0){throw new FS.ErrnoError(8)}if(!FS.isFile(stream.node.mode)&&!FS.isDir(stream.node.mode)){throw new FS.ErrnoError(43)}if(!stream.stream_ops.allocate){throw new FS.ErrnoError(138)}stream.stream_ops.allocate(stream,offset,length)},mmap:(stream,address,length,position,prot,flags)=>{if((prot&2)!==0&&(flags&2)===0&&(stream.flags&2097155)!==2){throw new FS.ErrnoError(2)}if((stream.flags&2097155)===1){throw new FS.ErrnoError(2)}if(!stream.stream_ops.mmap){throw new FS.ErrnoError(43)}return stream.stream_ops.mmap(stream,address,length,position,prot,flags)},msync:(stream,buffer,offset,length,mmapFlags)=>{if(!stream||!stream.stream_ops.msync){return 0}return stream.stream_ops.msync(stream,buffer,offset,length,mmapFlags)},munmap:stream=>0,ioctl:(stream,cmd,arg)=>{if(!stream.stream_ops.ioctl){throw new FS.ErrnoError(59)}return stream.stream_ops.ioctl(stream,cmd,arg)},readFile:(path,opts={})=>{opts.flags=opts.flags||0;opts.encoding=opts.encoding||"binary";if(opts.encoding!=="utf8"&&opts.encoding!=="binary"){throw new Error('Invalid encoding type "'+opts.encoding+'"')}var ret;var stream=FS.open(path,opts.flags);var stat=FS.stat(path);var length=stat.size;var buf=new Uint8Array(length);FS.read(stream,buf,0,length,0);if(opts.encoding==="utf8"){ret=UTF8ArrayToString(buf,0)}else if(opts.encoding==="binary"){ret=buf}FS.close(stream);return ret},writeFile:(path,data,opts={})=>{opts.flags=opts.flags||577;var stream=FS.open(path,opts.flags,opts.mode);if(typeof data==="string"){var buf=new Uint8Array(lengthBytesUTF8(data)+1);var actualNumBytes=stringToUTF8Array(data,buf,0,buf.length);FS.write(stream,buf,0,actualNumBytes,undefined,opts.canOwn)}else if(ArrayBuffer.isView(data)){FS.write(stream,data,0,data.byteLength,undefined,opts.canOwn)}else{throw new Error("Unsupported data type")}FS.close(stream)},cwd:()=>FS.currentPath,chdir:path=>{var lookup=FS.lookupPath(path,{follow:true});if(lookup.node===null){throw new FS.ErrnoError(44)}if(!FS.isDir(lookup.node.mode)){throw new FS.ErrnoError(54)}var errCode=FS.nodePermissions(lookup.node,"x");if(errCode){throw new FS.ErrnoError(errCode)}FS.currentPath=lookup.path},createDefaultDirectories:()=>{FS.mkdir("/tmp");FS.mkdir("/home");FS.mkdir("/home/web_user")},createDefaultDevices:()=>{FS.mkdir("/dev");FS.registerDevice(FS.makedev(1,3),{read:()=>0,write:(stream,buffer,offset,length,pos)=>length});FS.mkdev("/dev/null",FS.makedev(1,3));TTY.register(FS.makedev(5,0),TTY.default_tty_ops);TTY.register(FS.makedev(6,0),TTY.default_tty1_ops);FS.mkdev("/dev/tty",FS.makedev(5,0));FS.mkdev("/dev/tty1",FS.makedev(6,0));var random_device=getRandomDevice();FS.createDevice("/dev","random",random_device);FS.createDevice("/dev","urandom",random_device);FS.mkdir("/dev/shm");FS.mkdir("/dev/shm/tmp")},createSpecialDirectories:()=>{FS.mkdir("/proc");var proc_self=FS.mkdir("/proc/self");FS.mkdir("/proc/self/fd");FS.mount({mount:()=>{var node=FS.createNode(proc_self,"fd",16384|511,73);node.node_ops={lookup:(parent,name)=>{var fd=+name;var stream=FS.getStream(fd);if(!stream)throw new FS.ErrnoError(8);var ret={parent:null,mount:{mountpoint:"fake"},node_ops:{readlink:()=>stream.path}};ret.parent=ret;return ret}};return node}},{},"/proc/self/fd")},createStandardStreams:()=>{if(Module["stdin"]){FS.createDevice("/dev","stdin",Module["stdin"])}else{FS.symlink("/dev/tty","/dev/stdin")}if(Module["stdout"]){FS.createDevice("/dev","stdout",null,Module["stdout"])}else{FS.symlink("/dev/tty","/dev/stdout")}if(Module["stderr"]){FS.createDevice("/dev","stderr",null,Module["stderr"])}else{FS.symlink("/dev/tty1","/dev/stderr")}var stdin=FS.open("/dev/stdin",0);var stdout=FS.open("/dev/stdout",1);var stderr=FS.open("/dev/stderr",1)},ensureErrnoError:()=>{if(FS.ErrnoError)return;FS.ErrnoError=function ErrnoError(errno,node){this.node=node;this.setErrno=function(errno){this.errno=errno};this.setErrno(errno);this.message="FS error"};FS.ErrnoError.prototype=new Error;FS.ErrnoError.prototype.constructor=FS.ErrnoError;[44].forEach(code=>{FS.genericErrors[code]=new FS.ErrnoError(code);FS.genericErrors[code].stack="<generic error, no stack>"})},staticInit:()=>{FS.ensureErrnoError();FS.nameTable=new Array(4096);FS.mount(MEMFS,{},"/");FS.createDefaultDirectories();FS.createDefaultDevices();FS.createSpecialDirectories();FS.filesystems={"MEMFS":MEMFS}},init:(input,output,error)=>{FS.init.initialized=true;FS.ensureErrnoError();Module["stdin"]=input||Module["stdin"];Module["stdout"]=output||Module["stdout"];Module["stderr"]=error||Module["stderr"];FS.createStandardStreams()},quit:()=>{FS.init.initialized=false;for(var i=0;i<FS.streams.length;i++){var stream=FS.streams[i];if(!stream){continue}FS.close(stream)}},getMode:(canRead,canWrite)=>{var mode=0;if(canRead)mode|=292|73;if(canWrite)mode|=146;return mode},findObject:(path,dontResolveLastLink)=>{var ret=FS.analyzePath(path,dontResolveLastLink);if(ret.exists){return ret.object}else{return null}},analyzePath:(path,dontResolveLastLink)=>{try{var lookup=FS.lookupPath(path,{follow:!dontResolveLastLink});path=lookup.path}catch(e){}var ret={isRoot:false,exists:false,error:0,name:null,path:null,object:null,parentExists:false,parentPath:null,parentObject:null};try{var lookup=FS.lookupPath(path,{parent:true});ret.parentExists=true;ret.parentPath=lookup.path;ret.parentObject=lookup.node;ret.name=PATH.basename(path);lookup=FS.lookupPath(path,{follow:!dontResolveLastLink});ret.exists=true;ret.path=lookup.path;ret.object=lookup.node;ret.name=lookup.node.name;ret.isRoot=lookup.path==="/"}catch(e){ret.error=e.errno}return ret},createPath:(parent,path,canRead,canWrite)=>{parent=typeof parent==="string"?parent:FS.getPath(parent);var parts=path.split("/").reverse();while(parts.length){var part=parts.pop();if(!part)continue;var current=PATH.join2(parent,part);try{FS.mkdir(current)}catch(e){}parent=current}return current},createFile:(parent,name,properties,canRead,canWrite)=>{var path=PATH.join2(typeof parent==="string"?parent:FS.getPath(parent),name);var mode=FS.getMode(canRead,canWrite);return FS.create(path,mode)},createDataFile:(parent,name,data,canRead,canWrite,canOwn)=>{var path=name;if(parent){parent=typeof parent==="string"?parent:FS.getPath(parent);path=name?PATH.join2(parent,name):parent}var mode=FS.getMode(canRead,canWrite);var node=FS.create(path,mode);if(data){if(typeof data==="string"){var arr=new Array(data.length);for(var i=0,len=data.length;i<len;++i)arr[i]=data.charCodeAt(i);data=arr}FS.chmod(node,mode|146);var stream=FS.open(node,577);FS.write(stream,data,0,data.length,0,canOwn);FS.close(stream);FS.chmod(node,mode)}return node},createDevice:(parent,name,input,output)=>{var path=PATH.join2(typeof parent==="string"?parent:FS.getPath(parent),name);var mode=FS.getMode(!!input,!!output);if(!FS.createDevice.major)FS.createDevice.major=64;var dev=FS.makedev(FS.createDevice.major++,0);FS.registerDevice(dev,{open:stream=>{stream.seekable=false},close:stream=>{if(output&&output.buffer&&output.buffer.length){output(10)}},read:(stream,buffer,offset,length,pos)=>{var bytesRead=0;for(var i=0;i<length;i++){var result;try{result=input()}catch(e){throw new FS.ErrnoError(29)}if(result===undefined&&bytesRead===0){throw new FS.ErrnoError(6)}if(result===null||result===undefined)break;bytesRead++;buffer[offset+i]=result}if(bytesRead){stream.node.timestamp=Date.now()}return bytesRead},write:(stream,buffer,offset,length,pos)=>{for(var i=0;i<length;i++){try{output(buffer[offset+i])}catch(e){throw new FS.ErrnoError(29)}}if(length){stream.node.timestamp=Date.now()}return i}});return FS.mkdev(path,mode,dev)},forceLoadFile:obj=>{if(obj.isDevice||obj.isFolder||obj.link||obj.contents)return true;if(typeof XMLHttpRequest!=="undefined"){throw new Error("Lazy loading should have been performed (contents set) in createLazyFile, but it was not. Lazy loading only works in web workers. Use --embed-file or --preload-file in emcc on the main thread.")}else if(read_){try{obj.contents=intArrayFromString(read_(obj.url),true);obj.usedBytes=obj.contents.length}catch(e){throw new FS.ErrnoError(29)}}else{throw new Error("Cannot load without read() or XMLHttpRequest.")}},createLazyFile:(parent,name,url,canRead,canWrite)=>{function LazyUint8Array(){this.lengthKnown=false;this.chunks=[]}LazyUint8Array.prototype.get=function LazyUint8Array_get(idx){if(idx>this.length-1||idx<0){return undefined}var chunkOffset=idx%this.chunkSize;var chunkNum=idx/this.chunkSize|0;return this.getter(chunkNum)[chunkOffset]};LazyUint8Array.prototype.setDataGetter=function LazyUint8Array_setDataGetter(getter){this.getter=getter};LazyUint8Array.prototype.cacheLength=function LazyUint8Array_cacheLength(){var xhr=new XMLHttpRequest;xhr.open("HEAD",url,false);xhr.send(null);if(!(xhr.status>=200&&xhr.status<300||xhr.status===304))throw new Error("Couldn't load "+url+". Status: "+xhr.status);var datalength=Number(xhr.getResponseHeader("Content-length"));var header;var hasByteServing=(header=xhr.getResponseHeader("Accept-Ranges"))&&header==="bytes";var usesGzip=(header=xhr.getResponseHeader("Content-Encoding"))&&header==="gzip";var chunkSize=1024*1024;if(!hasByteServing)chunkSize=datalength;var doXHR=(from,to)=>{if(from>to)throw new Error("invalid range ("+from+", "+to+") or no bytes requested!");if(to>datalength-1)throw new Error("only "+datalength+" bytes available! programmer error!");var xhr=new XMLHttpRequest;xhr.open("GET",url,false);if(datalength!==chunkSize)xhr.setRequestHeader("Range","bytes="+from+"-"+to);if(typeof Uint8Array!="undefined")xhr.responseType="arraybuffer";if(xhr.overrideMimeType){xhr.overrideMimeType("text/plain; charset=x-user-defined")}xhr.send(null);if(!(xhr.status>=200&&xhr.status<300||xhr.status===304))throw new Error("Couldn't load "+url+". Status: "+xhr.status);if(xhr.response!==undefined){return new Uint8Array(xhr.response||[])}else{return intArrayFromString(xhr.responseText||"",true)}};var lazyArray=this;lazyArray.setDataGetter(chunkNum=>{var start=chunkNum*chunkSize;var end=(chunkNum+1)*chunkSize-1;end=Math.min(end,datalength-1);if(typeof lazyArray.chunks[chunkNum]==="undefined"){lazyArray.chunks[chunkNum]=doXHR(start,end)}if(typeof lazyArray.chunks[chunkNum]==="undefined")throw new Error("doXHR failed!");return lazyArray.chunks[chunkNum]});if(usesGzip||!datalength){chunkSize=datalength=1;datalength=this.getter(0).length;chunkSize=datalength;out("LazyFiles on gzip forces download of the whole file when length is accessed")}this._length=datalength;this._chunkSize=chunkSize;this.lengthKnown=true};if(typeof XMLHttpRequest!=="undefined"){if(!ENVIRONMENT_IS_WORKER)throw"Cannot do synchronous binary XHRs outside webworkers in modern browsers. Use --embed-file or --preload-file in emcc";var lazyArray=new LazyUint8Array;Object.defineProperties(lazyArray,{length:{get:function(){if(!this.lengthKnown){this.cacheLength()}return this._length}},chunkSize:{get:function(){if(!this.lengthKnown){this.cacheLength()}return this._chunkSize}}});var properties={isDevice:false,contents:lazyArray}}else{var properties={isDevice:false,url:url}}var node=FS.createFile(parent,name,properties,canRead,canWrite);if(properties.contents){node.contents=properties.contents}else if(properties.url){node.contents=null;node.url=properties.url}Object.defineProperties(node,{usedBytes:{get:function(){return this.contents.length}}});var stream_ops={};var keys=Object.keys(node.stream_ops);keys.forEach(key=>{var fn=node.stream_ops[key];stream_ops[key]=function forceLoadLazyFile(){FS.forceLoadFile(node);return fn.apply(null,arguments)}});stream_ops.read=((stream,buffer,offset,length,position)=>{FS.forceLoadFile(node);var contents=stream.node.contents;if(position>=contents.length)return 0;var size=Math.min(contents.length-position,length);if(contents.slice){for(var i=0;i<size;i++){buffer[offset+i]=contents[position+i]}}else{for(var i=0;i<size;i++){buffer[offset+i]=contents.get(position+i)}}return size});node.stream_ops=stream_ops;return node},createPreloadedFile:(parent,name,url,canRead,canWrite,onload,onerror,dontCreateFile,canOwn,preFinish)=>{var fullname=name?PATH_FS.resolve(PATH.join2(parent,name)):parent;var dep=getUniqueRunDependency("cp "+fullname);function processData(byteArray){function finish(byteArray){if(preFinish)preFinish();if(!dontCreateFile){FS.createDataFile(parent,name,byteArray,canRead,canWrite,canOwn)}if(onload)onload();removeRunDependency(dep)}if(Browser.handledByPreloadPlugin(byteArray,fullname,finish,()=>{if(onerror)onerror();removeRunDependency(dep)})){return}finish(byteArray)}addRunDependency(dep);if(typeof url=="string"){asyncLoad(url,byteArray=>processData(byteArray),onerror)}else{processData(url)}},indexedDB:()=>{return window.indexedDB||window.mozIndexedDB||window.webkitIndexedDB||window.msIndexedDB},DB_NAME:()=>{return"EM_FS_"+window.location.pathname},DB_VERSION:20,DB_STORE_NAME:"FILE_DATA",saveFilesToDB:(paths,onload,onerror)=>{onload=onload||(()=>{});onerror=onerror||(()=>{});var indexedDB=FS.indexedDB();try{var openRequest=indexedDB.open(FS.DB_NAME(),FS.DB_VERSION)}catch(e){return onerror(e)}openRequest.onupgradeneeded=(()=>{out("creating db");var db=openRequest.result;db.createObjectStore(FS.DB_STORE_NAME)});openRequest.onsuccess=(()=>{var db=openRequest.result;var transaction=db.transaction([FS.DB_STORE_NAME],"readwrite");var files=transaction.objectStore(FS.DB_STORE_NAME);var ok=0,fail=0,total=paths.length;function finish(){if(fail==0)onload();else onerror()}paths.forEach(path=>{var putRequest=files.put(FS.analyzePath(path).object.contents,path);putRequest.onsuccess=(()=>{ok++;if(ok+fail==total)finish()});putRequest.onerror=(()=>{fail++;if(ok+fail==total)finish()})});transaction.onerror=onerror});openRequest.onerror=onerror},loadFilesFromDB:(paths,onload,onerror)=>{onload=onload||(()=>{});onerror=onerror||(()=>{});var indexedDB=FS.indexedDB();try{var openRequest=indexedDB.open(FS.DB_NAME(),FS.DB_VERSION)}catch(e){return onerror(e)}openRequest.onupgradeneeded=onerror;openRequest.onsuccess=(()=>{var db=openRequest.result;try{var transaction=db.transaction([FS.DB_STORE_NAME],"readonly")}catch(e){onerror(e);return}var files=transaction.objectStore(FS.DB_STORE_NAME);var ok=0,fail=0,total=paths.length;function finish(){if(fail==0)onload();else onerror()}paths.forEach(path=>{var getRequest=files.get(path);getRequest.onsuccess=(()=>{if(FS.analyzePath(path).exists){FS.unlink(path)}FS.createDataFile(PATH.dirname(path),PATH.basename(path),getRequest.result,true,true,true);ok++;if(ok+fail==total)finish()});getRequest.onerror=(()=>{fail++;if(ok+fail==total)finish()})});transaction.onerror=onerror});openRequest.onerror=onerror}};var SYSCALLS={mappings:{},DEFAULT_POLLMASK:5,calculateAt:function(dirfd,path,allowEmpty){if(path[0]==="/"){return path}var dir;if(dirfd===-100){dir=FS.cwd()}else{var dirstream=FS.getStream(dirfd);if(!dirstream)throw new FS.ErrnoError(8);dir=dirstream.path}if(path.length==0){if(!allowEmpty){throw new FS.ErrnoError(44)}return dir}return PATH.join2(dir,path)},doStat:function(func,path,buf){try{var stat=func(path)}catch(e){if(e&&e.node&&PATH.normalize(path)!==PATH.normalize(FS.getPath(e.node))){return-54}throw e}HEAP32[buf>>2]=stat.dev;HEAP32[buf+4>>2]=0;HEAP32[buf+8>>2]=stat.ino;HEAP32[buf+12>>2]=stat.mode;HEAP32[buf+16>>2]=stat.nlink;HEAP32[buf+20>>2]=stat.uid;HEAP32[buf+24>>2]=stat.gid;HEAP32[buf+28>>2]=stat.rdev;HEAP32[buf+32>>2]=0;tempI64=[stat.size>>>0,(tempDouble=stat.size,+Math.abs(tempDouble)>=1?tempDouble>0?(Math.min(+Math.floor(tempDouble/4294967296),4294967295)|0)>>>0:~~+Math.ceil((tempDouble-+(~~tempDouble>>>0))/4294967296)>>>0:0)],HEAP32[buf+40>>2]=tempI64[0],HEAP32[buf+44>>2]=tempI64[1];HEAP32[buf+48>>2]=4096;HEAP32[buf+52>>2]=stat.blocks;HEAP32[buf+56>>2]=stat.atime.getTime()/1e3|0;HEAP32[buf+60>>2]=0;HEAP32[buf+64>>2]=stat.mtime.getTime()/1e3|0;HEAP32[buf+68>>2]=0;HEAP32[buf+72>>2]=stat.ctime.getTime()/1e3|0;HEAP32[buf+76>>2]=0;tempI64=[stat.ino>>>0,(tempDouble=stat.ino,+Math.abs(tempDouble)>=1?tempDouble>0?(Math.min(+Math.floor(tempDouble/4294967296),4294967295)|0)>>>0:~~+Math.ceil((tempDouble-+(~~tempDouble>>>0))/4294967296)>>>0:0)],HEAP32[buf+80>>2]=tempI64[0],HEAP32[buf+84>>2]=tempI64[1];return 0},doMsync:function(addr,stream,len,flags,offset){var buffer=HEAPU8.slice(addr,addr+len);FS.msync(stream,buffer,offset,len,flags)},doMkdir:function(path,mode){path=PATH.normalize(path);if(path[path.length-1]==="/")path=path.substr(0,path.length-1);FS.mkdir(path,mode,0);return 0},doMknod:function(path,mode,dev){switch(mode&61440){case 32768:case 8192:case 24576:case 4096:case 49152:break;default:return-28}FS.mknod(path,mode,dev);return 0},doReadlink:function(path,buf,bufsize){if(bufsize<=0)return-28;var ret=FS.readlink(path);var len=Math.min(bufsize,lengthBytesUTF8(ret));var endChar=HEAP8[buf+len];stringToUTF8(ret,buf,bufsize+1);HEAP8[buf+len]=endChar;return len},doAccess:function(path,amode){if(amode&~7){return-28}var lookup=FS.lookupPath(path,{follow:true});var node=lookup.node;if(!node){return-44}var perms="";if(amode&4)perms+="r";if(amode&2)perms+="w";if(amode&1)perms+="x";if(perms&&FS.nodePermissions(node,perms)){return-2}return 0},doDup:function(path,flags,suggestFD){var suggest=FS.getStream(suggestFD);if(suggest)FS.close(suggest);return FS.open(path,flags,0,suggestFD,suggestFD).fd},doReadv:function(stream,iov,iovcnt,offset){var ret=0;for(var i=0;i<iovcnt;i++){var ptr=HEAP32[iov+i*8>>2];var len=HEAP32[iov+(i*8+4)>>2];var curr=FS.read(stream,HEAP8,ptr,len,offset);if(curr<0)return-1;ret+=curr;if(curr<len)break}return ret},doWritev:function(stream,iov,iovcnt,offset){var ret=0;for(var i=0;i<iovcnt;i++){var ptr=HEAP32[iov+i*8>>2];var len=HEAP32[iov+(i*8+4)>>2];var curr=FS.write(stream,HEAP8,ptr,len,offset);if(curr<0)return-1;ret+=curr}return ret},varargs:undefined,get:function(){SYSCALLS.varargs+=4;var ret=HEAP32[SYSCALLS.varargs-4>>2];return ret},getStr:function(ptr){var ret=UTF8ToString(ptr);return ret},getStreamFromFD:function(fd){var stream=FS.getStream(fd);if(!stream)throw new FS.ErrnoError(8);return stream},get64:function(low,high){return low}};function ___syscall_fcntl64(fd,cmd,varargs){SYSCALLS.varargs=varargs;try{var stream=SYSCALLS.getStreamFromFD(fd);switch(cmd){case 0:{var arg=SYSCALLS.get();if(arg<0){return-28}var newStream;newStream=FS.open(stream.path,stream.flags,0,arg);return newStream.fd}case 1:case 2:return 0;case 3:return stream.flags;case 4:{var arg=SYSCALLS.get();stream.flags|=arg;return 0}case 5:{var arg=SYSCALLS.get();var offset=0;HEAP16[arg+offset>>1]=2;return 0}case 6:case 7:return 0;case 16:case 8:return-28;case 9:setErrNo(28);return-1;default:{return-28}}}catch(e){if(typeof FS==="undefined"||!(e instanceof FS.ErrnoError))throw e;return-e.errno}}function ___syscall_ioctl(fd,op,varargs){SYSCALLS.varargs=varargs;try{var stream=SYSCALLS.getStreamFromFD(fd);switch(op){case 21509:case 21505:{if(!stream.tty)return-59;return 0}case 21510:case 21511:case 21512:case 21506:case 21507:case 21508:{if(!stream.tty)return-59;return 0}case 21519:{if(!stream.tty)return-59;var argp=SYSCALLS.get();HEAP32[argp>>2]=0;return 0}case 21520:{if(!stream.tty)return-59;return-28}case 21531:{var argp=SYSCALLS.get();return FS.ioctl(stream,op,argp)}case 21523:{if(!stream.tty)return-59;return 0}case 21524:{if(!stream.tty)return-59;return 0}default:abort("bad ioctl syscall "+op)}}catch(e){if(typeof FS==="undefined"||!(e instanceof FS.ErrnoError))throw e;return-e.errno}}function ___syscall_open(path,flags,varargs){SYSCALLS.varargs=varargs;try{var pathname=SYSCALLS.getStr(path);var mode=varargs?SYSCALLS.get():0;var stream=FS.open(pathname,flags,mode);return stream.fd}catch(e){if(typeof FS==="undefined"||!(e instanceof FS.ErrnoError))throw e;return-e.errno}}function _abort(){abort("")}function _clock(){if(_clock.start===undefined)_clock.start=Date.now();return(Date.now()-_clock.start)*(1e6/1e3)|0}function _emscripten_memcpy_big(dest,src,num){HEAPU8.copyWithin(dest,src,src+num)}function emscripten_realloc_buffer(size){try{wasmMemory.grow(size-buffer.byteLength+65535>>>16);updateGlobalBufferAndViews(wasmMemory.buffer);return 1}catch(e){}}function _emscripten_resize_heap(requestedSize){var oldSize=HEAPU8.length;requestedSize=requestedSize>>>0;var maxHeapSize=2147483648;if(requestedSize>maxHeapSize){return false}for(var cutDown=1;cutDown<=4;cutDown*=2){var overGrownHeapSize=oldSize*(1+.2/cutDown);overGrownHeapSize=Math.min(overGrownHeapSize,requestedSize+100663296);var newSize=Math.min(maxHeapSize,alignUp(Math.max(requestedSize,overGrownHeapSize),65536));var replacement=emscripten_realloc_buffer(newSize);if(replacement){return true}}return false}var ENV={};function getExecutableName(){return thisProgram||"./this.program"}function getEnvStrings(){if(!getEnvStrings.strings){var lang=(typeof navigator==="object"&&navigator.languages&&navigator.languages[0]||"C").replace("-","_")+".UTF-8";var env={"USER":"web_user","LOGNAME":"web_user","PATH":"/","PWD":"/","HOME":"/home/web_user","LANG":lang,"_":getExecutableName()};for(var x in ENV){if(ENV[x]===undefined)delete env[x];else env[x]=ENV[x]}var strings=[];for(var x in env){strings.push(x+"="+env[x])}getEnvStrings.strings=strings}return getEnvStrings.strings}function _environ_get(__environ,environ_buf){var bufSize=0;getEnvStrings().forEach(function(string,i){var ptr=environ_buf+bufSize;HEAP32[__environ+i*4>>2]=ptr;writeAsciiToMemory(string,ptr);bufSize+=string.length+1});return 0}function _environ_sizes_get(penviron_count,penviron_buf_size){var strings=getEnvStrings();HEAP32[penviron_count>>2]=strings.length;var bufSize=0;strings.forEach(function(string){bufSize+=string.length+1});HEAP32[penviron_buf_size>>2]=bufSize;return 0}function _exit(status){exit(status)}function _fd_close(fd){try{var stream=SYSCALLS.getStreamFromFD(fd);FS.close(stream);return 0}catch(e){if(typeof FS==="undefined"||!(e instanceof FS.ErrnoError))throw e;return e.errno}}function _fd_read(fd,iov,iovcnt,pnum){try{var stream=SYSCALLS.getStreamFromFD(fd);var num=SYSCALLS.doReadv(stream,iov,iovcnt);HEAP32[pnum>>2]=num;return 0}catch(e){if(typeof FS==="undefined"||!(e instanceof FS.ErrnoError))throw e;return e.errno}}function _fd_seek(fd,offset_low,offset_high,whence,newOffset){try{var stream=SYSCALLS.getStreamFromFD(fd);var HIGH_OFFSET=4294967296;var offset=offset_high*HIGH_OFFSET+(offset_low>>>0);var DOUBLE_LIMIT=9007199254740992;if(offset<=-DOUBLE_LIMIT||offset>=DOUBLE_LIMIT){return-61}FS.llseek(stream,offset,whence);tempI64=[stream.position>>>0,(tempDouble=stream.position,+Math.abs(tempDouble)>=1?tempDouble>0?(Math.min(+Math.floor(tempDouble/4294967296),4294967295)|0)>>>0:~~+Math.ceil((tempDouble-+(~~tempDouble>>>0))/4294967296)>>>0:0)],HEAP32[newOffset>>2]=tempI64[0],HEAP32[newOffset+4>>2]=tempI64[1];if(stream.getdents&&offset===0&&whence===0)stream.getdents=null;return 0}catch(e){if(typeof FS==="undefined"||!(e instanceof FS.ErrnoError))throw e;return e.errno}}function _fd_write(fd,iov,iovcnt,pnum){try{var stream=SYSCALLS.getStreamFromFD(fd);var num=SYSCALLS.doWritev(stream,iov,iovcnt);HEAP32[pnum>>2]=num;return 0}catch(e){if(typeof FS==="undefined"||!(e instanceof FS.ErrnoError))throw e;return e.errno}}function __isLeapYear(year){return year%4===0&&(year%100!==0||year%400===0)}function __arraySum(array,index){var sum=0;for(var i=0;i<=index;sum+=array[i++]){}return sum}var __MONTH_DAYS_LEAP=[31,29,31,30,31,30,31,31,30,31,30,31];var __MONTH_DAYS_REGULAR=[31,28,31,30,31,30,31,31,30,31,30,31];function __addDays(date,days){var newDate=new Date(date.getTime());while(days>0){var leap=__isLeapYear(newDate.getFullYear());var currentMonth=newDate.getMonth();var daysInCurrentMonth=(leap?__MONTH_DAYS_LEAP:__MONTH_DAYS_REGULAR)[currentMonth];if(days>daysInCurrentMonth-newDate.getDate()){days-=daysInCurrentMonth-newDate.getDate()+1;newDate.setDate(1);if(currentMonth<11){newDate.setMonth(currentMonth+1)}else{newDate.setMonth(0);newDate.setFullYear(newDate.getFullYear()+1)}}else{newDate.setDate(newDate.getDate()+days);return newDate}}return newDate}function _strftime(s,maxsize,format,tm){var tm_zone=HEAP32[tm+40>>2];var date={tm_sec:HEAP32[tm>>2],tm_min:HEAP32[tm+4>>2],tm_hour:HEAP32[tm+8>>2],tm_mday:HEAP32[tm+12>>2],tm_mon:HEAP32[tm+16>>2],tm_year:HEAP32[tm+20>>2],tm_wday:HEAP32[tm+24>>2],tm_yday:HEAP32[tm+28>>2],tm_isdst:HEAP32[tm+32>>2],tm_gmtoff:HEAP32[tm+36>>2],tm_zone:tm_zone?UTF8ToString(tm_zone):""};var pattern=UTF8ToString(format);var EXPANSION_RULES_1={"%c":"%a %b %d %H:%M:%S %Y","%D":"%m/%d/%y","%F":"%Y-%m-%d","%h":"%b","%r":"%I:%M:%S %p","%R":"%H:%M","%T":"%H:%M:%S","%x":"%m/%d/%y","%X":"%H:%M:%S","%Ec":"%c","%EC":"%C","%Ex":"%m/%d/%y","%EX":"%H:%M:%S","%Ey":"%y","%EY":"%Y","%Od":"%d","%Oe":"%e","%OH":"%H","%OI":"%I","%Om":"%m","%OM":"%M","%OS":"%S","%Ou":"%u","%OU":"%U","%OV":"%V","%Ow":"%w","%OW":"%W","%Oy":"%y"};for(var rule in EXPANSION_RULES_1){pattern=pattern.replace(new RegExp(rule,"g"),EXPANSION_RULES_1[rule])}var WEEKDAYS=["Sunday","Monday","Tuesday","Wednesday","Thursday","Friday","Saturday"];var MONTHS=["January","February","March","April","May","June","July","August","September","October","November","December"];function leadingSomething(value,digits,character){var str=typeof value==="number"?value.toString():value||"";while(str.length<digits){str=character[0]+str}return str}function leadingNulls(value,digits){return leadingSomething(value,digits,"0")}function compareByDay(date1,date2){function sgn(value){return value<0?-1:value>0?1:0}var compare;if((compare=sgn(date1.getFullYear()-date2.getFullYear()))===0){if((compare=sgn(date1.getMonth()-date2.getMonth()))===0){compare=sgn(date1.getDate()-date2.getDate())}}return compare}function getFirstWeekStartDate(janFourth){switch(janFourth.getDay()){case 0:return new Date(janFourth.getFullYear()-1,11,29);case 1:return janFourth;case 2:return new Date(janFourth.getFullYear(),0,3);case 3:return new Date(janFourth.getFullYear(),0,2);case 4:return new Date(janFourth.getFullYear(),0,1);case 5:return new Date(janFourth.getFullYear()-1,11,31);case 6:return new Date(janFourth.getFullYear()-1,11,30)}}function getWeekBasedYear(date){var thisDate=__addDays(new Date(date.tm_year+1900,0,1),date.tm_yday);var janFourthThisYear=new Date(thisDate.getFullYear(),0,4);var janFourthNextYear=new Date(thisDate.getFullYear()+1,0,4);var firstWeekStartThisYear=getFirstWeekStartDate(janFourthThisYear);var firstWeekStartNextYear=getFirstWeekStartDate(janFourthNextYear);if(compareByDay(firstWeekStartThisYear,thisDate)<=0){if(compareByDay(firstWeekStartNextYear,thisDate)<=0){return thisDate.getFullYear()+1}else{return thisDate.getFullYear()}}else{return thisDate.getFullYear()-1}}var EXPANSION_RULES_2={"%a":function(date){return WEEKDAYS[date.tm_wday].substring(0,3)},"%A":function(date){return WEEKDAYS[date.tm_wday]},"%b":function(date){return MONTHS[date.tm_mon].substring(0,3)},"%B":function(date){return MONTHS[date.tm_mon]},"%C":function(date){var year=date.tm_year+1900;return leadingNulls(year/100|0,2)},"%d":function(date){return leadingNulls(date.tm_mday,2)},"%e":function(date){return leadingSomething(date.tm_mday,2," ")},"%g":function(date){return getWeekBasedYear(date).toString().substring(2)},"%G":function(date){return getWeekBasedYear(date)},"%H":function(date){return leadingNulls(date.tm_hour,2)},"%I":function(date){var twelveHour=date.tm_hour;if(twelveHour==0)twelveHour=12;else if(twelveHour>12)twelveHour-=12;return leadingNulls(twelveHour,2)},"%j":function(date){return leadingNulls(date.tm_mday+__arraySum(__isLeapYear(date.tm_year+1900)?__MONTH_DAYS_LEAP:__MONTH_DAYS_REGULAR,date.tm_mon-1),3)},"%m":function(date){return leadingNulls(date.tm_mon+1,2)},"%M":function(date){return leadingNulls(date.tm_min,2)},"%n":function(){return"\n"},"%p":function(date){if(date.tm_hour>=0&&date.tm_hour<12){return"AM"}else{return"PM"}},"%S":function(date){return leadingNulls(date.tm_sec,2)},"%t":function(){return"\t"},"%u":function(date){return date.tm_wday||7},"%U":function(date){var janFirst=new Date(date.tm_year+1900,0,1);var firstSunday=janFirst.getDay()===0?janFirst:__addDays(janFirst,7-janFirst.getDay());var endDate=new Date(date.tm_year+1900,date.tm_mon,date.tm_mday);if(compareByDay(firstSunday,endDate)<0){var februaryFirstUntilEndMonth=__arraySum(__isLeapYear(endDate.getFullYear())?__MONTH_DAYS_LEAP:__MONTH_DAYS_REGULAR,endDate.getMonth()-1)-31;var firstSundayUntilEndJanuary=31-firstSunday.getDate();var days=firstSundayUntilEndJanuary+februaryFirstUntilEndMonth+endDate.getDate();return leadingNulls(Math.ceil(days/7),2)}return compareByDay(firstSunday,janFirst)===0?"01":"00"},"%V":function(date){var janFourthThisYear=new Date(date.tm_year+1900,0,4);var janFourthNextYear=new Date(date.tm_year+1901,0,4);var firstWeekStartThisYear=getFirstWeekStartDate(janFourthThisYear);var firstWeekStartNextYear=getFirstWeekStartDate(janFourthNextYear);var endDate=__addDays(new Date(date.tm_year+1900,0,1),date.tm_yday);if(compareByDay(endDate,firstWeekStartThisYear)<0){return"53"}if(compareByDay(firstWeekStartNextYear,endDate)<=0){return"01"}var daysDifference;if(firstWeekStartThisYear.getFullYear()<date.tm_year+1900){daysDifference=date.tm_yday+32-firstWeekStartThisYear.getDate()}else{daysDifference=date.tm_yday+1-firstWeekStartThisYear.getDate()}return leadingNulls(Math.ceil(daysDifference/7),2)},"%w":function(date){return date.tm_wday},"%W":function(date){var janFirst=new Date(date.tm_year,0,1);var firstMonday=janFirst.getDay()===1?janFirst:__addDays(janFirst,janFirst.getDay()===0?1:7-janFirst.getDay()+1);var endDate=new Date(date.tm_year+1900,date.tm_mon,date.tm_mday);if(compareByDay(firstMonday,endDate)<0){var februaryFirstUntilEndMonth=__arraySum(__isLeapYear(endDate.getFullYear())?__MONTH_DAYS_LEAP:__MONTH_DAYS_REGULAR,endDate.getMonth()-1)-31;var firstMondayUntilEndJanuary=31-firstMonday.getDate();var days=firstMondayUntilEndJanuary+februaryFirstUntilEndMonth+endDate.getDate();return leadingNulls(Math.ceil(days/7),2)}return compareByDay(firstMonday,janFirst)===0?"01":"00"},"%y":function(date){return(date.tm_year+1900).toString().substring(2)},"%Y":function(date){return date.tm_year+1900},"%z":function(date){var off=date.tm_gmtoff;var ahead=off>=0;off=Math.abs(off)/60;off=off/60*100+off%60;return(ahead?"+":"-")+String("0000"+off).slice(-4)},"%Z":function(date){return date.tm_zone},"%%":function(){return"%"}};for(var rule in EXPANSION_RULES_2){if(pattern.includes(rule)){pattern=pattern.replace(new RegExp(rule,"g"),EXPANSION_RULES_2[rule](date))}}var bytes=intArrayFromString(pattern,false);if(bytes.length>maxsize){return 0}writeArrayToMemory(bytes,s);return bytes.length-1}function _strftime_l(s,maxsize,format,tm){return _strftime(s,maxsize,format,tm)}var FSNode=function(parent,name,mode,rdev){if(!parent){parent=this}this.parent=parent;this.mount=parent.mount;this.mounted=null;this.id=FS.nextInode++;this.name=name;this.mode=mode;this.node_ops={};this.stream_ops={};this.rdev=rdev};var readMode=292|73;var writeMode=146;Object.defineProperties(FSNode.prototype,{read:{get:function(){return(this.mode&readMode)===readMode},set:function(val){val?this.mode|=readMode:this.mode&=~readMode}},write:{get:function(){return(this.mode&writeMode)===writeMode},set:function(val){val?this.mode|=writeMode:this.mode&=~writeMode}},isFolder:{get:function(){return FS.isDir(this.mode)}},isDevice:{get:function(){return FS.isChrdev(this.mode)}}});FS.FSNode=FSNode;FS.staticInit();Module["FS_createDataFile"]=FS.createDataFile;Module["FS_readFile"]=FS.readFile;Module["FS_unlink"]=FS.unlink;function intArrayFromString(stringy,dontAddNull,length){var len=length>0?length:lengthBytesUTF8(stringy)+1;var u8array=new Array(len);var numBytesWritten=stringToUTF8Array(stringy,u8array,0,u8array.length);if(dontAddNull)u8array.length=numBytesWritten;return u8array}var asmLibraryArg={"l":___cxa_allocate_exception,"k":___cxa_throw,"c":___syscall_fcntl64,"h":___syscall_ioctl,"i":___syscall_open,"d":_abort,"f":_clock,"j":_emscripten_memcpy_big,"q":_emscripten_resize_heap,"o":_environ_get,"p":_environ_sizes_get,"a":_exit,"e":_fd_close,"g":_fd_read,"m":_fd_seek,"b":_fd_write,"n":_strftime_l};var asm=createWasm();var ___wasm_call_ctors=Module["___wasm_call_ctors"]=function(){return(___wasm_call_ctors=Module["___wasm_call_ctors"]=Module["asm"]["s"]).apply(null,arguments)};var _malloc=Module["_malloc"]=function(){return(_malloc=Module["_malloc"]=Module["asm"]["t"]).apply(null,arguments)};var _simplify=Module["_simplify"]=function(){return(_simplify=Module["_simplify"]=Module["asm"]["u"]).apply(null,arguments)};var ___errno_location=Module["___errno_location"]=function(){return(___errno_location=Module["___errno_location"]=Module["asm"]["w"]).apply(null,arguments)};var stackSave=Module["stackSave"]=function(){return(stackSave=Module["stackSave"]=Module["asm"]["x"]).apply(null,arguments)};var stackRestore=Module["stackRestore"]=function(){return(stackRestore=Module["stackRestore"]=Module["asm"]["y"]).apply(null,arguments)};var stackAlloc=Module["stackAlloc"]=function(){return(stackAlloc=Module["stackAlloc"]=Module["asm"]["z"]).apply(null,arguments)};Module["ccall"]=ccall;Module["cwrap"]=cwrap;Module["FS_createDataFile"]=FS.createDataFile;Module["FS_unlink"]=FS.unlink;var calledRun;function ExitStatus(status){this.name="ExitStatus";this.message="Program terminated with exit("+status+")";this.status=status}dependenciesFulfilled=function runCaller(){if(!calledRun)run();if(!calledRun)dependenciesFulfilled=runCaller};function run(args){args=args||arguments_;if(runDependencies>0){return}preRun();if(runDependencies>0){return}function doRun(){if(calledRun)return;calledRun=true;Module["calledRun"]=true;if(ABORT)return;initRuntime();if(Module["onRuntimeInitialized"])Module["onRuntimeInitialized"]();postRun()}if(Module["setStatus"]){Module["setStatus"]("Running...");setTimeout(function(){setTimeout(function(){Module["setStatus"]("")},1);doRun()},1)}else{doRun()}}Module["run"]=run;function exit(status,implicit){EXITSTATUS=status;if(keepRuntimeAlive()){}else{exitRuntime()}procExit(status)}function procExit(code){EXITSTATUS=code;if(!keepRuntimeAlive()){if(Module["onExit"])Module["onExit"](code);ABORT=true}quit_(code,new ExitStatus(code))}if(Module["preInit"]){if(typeof Module["preInit"]=="function")Module["preInit"]=[Module["preInit"]];while(Module["preInit"].length>0){Module["preInit"].pop()()}}run();
2 | 


--------------------------------------------------------------------------------
/a.out.wasm:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/MyMiniFactory/Fast-Quadric-Mesh-Simplification/9bd297a5f840aa0178235469cc29114b0fe191a6/a.out.wasm


--------------------------------------------------------------------------------
/index.html:
--------------------------------------------------------------------------------
  1 | <!DOCTYPE html>
  2 | <html>
  3 | <head>
  4 |     <meta charset="UTF-8">
  5 |     
  6 |     <!-- Global site tag (gtag.js) - Google Analytics -->
  7 |     <script async src="https://www.googletagmanager.com/gtag/js?id=UA-115758488-2"></script>
  8 |     <script>
  9 |         window.dataLayer = window.dataLayer || [];
 10 |         function gtag(){dataLayer.push(arguments);}
 11 |         gtag('js', new Date());
 12 |         gtag('config', 'UA-115758488-2');
 13 |     </script>
 14 | 
 15 | </head>
 16 | 
 17 | <style>
 18 | 
 19 | #dragAndDrop {
 20 |     min-height: 200px;
 21 |     white-space: pre-line;
 22 |     border: 1px dashed black;
 23 |     border-radius: 10px;
 24 | }
 25 | 
 26 | #status_span {
 27 |     margin-left: 25px;
 28 | }
 29 | 
 30 | #info {
 31 |     width:100%
 32 |     display: -webkit-flex;
 33 |     display: flex;
 34 | }
 35 | 
 36 | #simplify_box {
 37 |     height: 20vh;
 38 |     overflow: auto;
 39 |     border: 1px solid black;
 40 | }
 41 | 
 42 | #note_box {
 43 |     margin-right: 30px;
 44 | }
 45 | 
 46 | #download_container {
 47 |     display: flex;
 48 |     justify-content: center;
 49 | }
 50 | 
 51 | #download {
 52 |     width:50%;
 53 |     height:40px;
 54 |     margin:0 auto;
 55 |     margin-top:10px;
 56 |     cursor: pointer;
 57 |     border: 0px;
 58 |     border-radius: 10px;
 59 |     background: rgb(240, 240, 240);
 60 | }
 61 | 
 62 | .download_ready_signal {
 63 |   animation: download_ready_signal 2s linear infinite;
 64 | }
 65 | 
 66 | @keyframes download_ready_signal {
 67 |   50% {
 68 |     opacity: 0.1;
 69 |     color: green;
 70 |     width: 60%;
 71 |   }
 72 | }
 73 | 
 74 | button {
 75 |     cursor: pointer;
 76 | }
 77 | 
 78 | </style>
 79 | 
 80 | <script type="text/javascript">
 81 |     `use strict`;
 82 | 
 83 |     var worker = new Worker('worker.js');
 84 | 
 85 |     worker.onmessage = function(e) {
 86 |         const log = e.data.log;
 87 |         if (log !== undefined) {
 88 |             var logger = document.getElementById('simplify_log');
 89 |             logger.innerHTML += '<li>' + log + '</li>';
 90 |             var box = document.getElementById('simplify_box');
 91 |             box.scrollTop = box.scrollHeight;
 92 |             return;
 93 |         }
 94 | 
 95 |         const file = e.data.blob;
 96 |         if (file !== undefined) {
 97 |             const s_name = SIMPLIFY_FILE.simplify_name;
 98 | 
 99 |             let url = window.URL.createObjectURL(file);
100 | 
101 |             var download = document.getElementById('download_a');
102 |             download.href = url;
103 |             download.download = s_name;
104 | 
105 |             var download_button = document.getElementById('download');
106 |             download_button.innerHTML = "Click to download "+"<br>" +s_name;
107 |             download_button.disabled = false;
108 |             download_button.className = "download_ready_signal";
109 |             setTimeout(function (){
110 |                 download_button.className = "";
111 |             }, 2000);
112 |             put_status("Ready to download");
113 |             return;
114 |         }
115 | 
116 |         console.error("Unknown Message from WebWorker", e.data);
117 |     }
118 | 
119 |     worker.onerror = function(e) {
120 |         console.log(e);
121 |         console.log(
122 |             'ERROR: Line ', e.lineno, ' in ', e.filename, ': ', e.message
123 |         );
124 |     }
125 | 
126 |     let SIMPLIFY_FILE = {
127 |         'blob': undefined,
128 |         get name(){
129 |             if (this.exists)
130 |                 return this.blob.name;
131 |         },
132 |         get simplify_name() {
133 |             if (this.exists)
134 |                 return "simplify_"+this.name;
135 |         },
136 |         get size() {
137 |             if (this.exists)
138 |                 return this.blob.size;
139 |         },
140 |         get exists() {
141 |             if (this.blob)
142 |                 return true;
143 |             else {
144 |                 return false;
145 |             }
146 |         }
147 |     }
148 | 
149 | 
150 |     function update_simplify_to() {
151 |         const slider_value = get_value_from_slider();
152 |         document.getElementById('percentage').innerHTML = slider_value;
153 |         if (SIMPLIFY_FILE.exists) {
154 |             document.getElementById('size').innerHTML = Math.ceil(slider_value * SIMPLIFY_FILE.size/(100*1024*1024));
155 |         }
156 |     }
157 | 
158 |     function get_value_from_slider() {
159 |         return document.getElementById('slider').valueAsNumber;
160 |     }
161 | 
162 |     function uploaded(file) {
163 |         if (file === undefined) { // via upload button
164 |             var uploadform = document.getElementById("fileuploadform");
165 |             file = uploadform.files[0];
166 | 
167 |             // this helps to force trigger even if user upload the same file
168 |             // https://stackoverflow.com/a/12102992/5260518
169 |             uploadform.value = null;
170 |         }
171 | 
172 |         if (file === undefined) // not via upload button defined otherwise
173 |             file = SIMPLIFY_FILE.blob;
174 | 
175 |         var uploadform = document.getElementById("fileuploadform");
176 |         console.log(uploadform);
177 |         uploadform.files[0] = file;
178 | 
179 |         SIMPLIFY_FILE.blob = file;
180 |         check_file(post_to_worker);
181 | 
182 |     }
183 | 
184 |     function check_file(success_cb) {
185 |         put_status("Checking file");
186 |         const filename = SIMPLIFY_FILE.name;
187 |         const extension = filename.toLowerCase().slice(filename.lastIndexOf(".")+1, filename.length);
188 |         if (extension!=="stl" && extension!=="obj") {
189 |             put_status("Please upload an stl or obj file not "+ extension);
190 |             return;
191 |         }
192 |         success_cb();
193 |     }
194 | 
195 |     function post_to_worker() {
196 |         update_simplify_to();
197 |         put_status("Simplifying by your browser...See log below");
198 |         worker.postMessage(
199 |             {"blob":SIMPLIFY_FILE.blob,
200 |              "percentage": get_value_from_slider()/100,
201 |              "simplify_name": SIMPLIFY_FILE.simplify_name
202 |             }
203 |         );
204 |     }
205 |     function dodrop(event) {
206 |         var dt = event.dataTransfer;
207 |         var file = dt.files[0];
208 |         var filename = file.name;
209 |         uploaded(file);
210 |     }
211 | 
212 |     function put_status(text)
213 |     {
214 |         document.getElementById("status").textContent = text;
215 |     }
216 | 
217 |     function init() {
218 |         update_simplify_to();
219 |     }
220 | 
221 |     window.onload = init;
222 | </script>
223 | 
224 | <body>
225 |     <h1>Mesh Simplification</h1>
226 |     <div id="dragAndDrop"
227 |          ondragenter="event.stopPropagation(); event.preventDefault();"
228 |          ondragover="event.stopPropagation(); event.preventDefault();"
229 |          ondrop="event.stopPropagation(); event.preventDefault();
230 |          dodrop(event);">
231 |         <ol>
232 |             <li>Drag and drop STL/OBJ files here or <input type="file" onChange="uploaded()" id="fileuploadform"/></li>
233 |             <li>Select simplify to <input type="range" min="1" max="99" value="50" id="slider" step="1" oninput="update_simplify_to()"><span id="percentage">50</span>% <span id="size">??</span>Mb</li>
234 |             <li>Click <button onclick="uploaded()">Simplify</button></li>
235 |         </ol>
236 |         <span id="status_span">Status: <span id="status">Waiting for upload</span>
237 |     </div>
238 |     <div id="download_container">
239 |         <a id="download_a" hidden></a>
240 |         <button id="download" disabled onclick="document.getElementById('download_a').click()">Nothing to download yet...</button>
241 |     </div>
242 |     <div id="info">
243 |         <div id="note_box">
244 |             <h3>Note</h3>
245 |             <ul>
246 |                 <li>WebAssembly Version of <a href="https://github.com/sp4cerat/Fast-Quadric-Mesh-Simplification">Fast Quadric Mesh Simplification by spe4cerat</a></li>
247 |             </ul>
248 |         </div>
249 |         <div>
250 |             <h3>Simplification Process Log</h3>
251 |                 <div id="simplify_box">
252 |                     <ul id="simplify_log">
253 |                     </ul>
254 |                 </div>
255 |             </div>
256 |         </div>
257 |     </div>
258 | 
259 |     <div id="log"></div>
260 | 
261 |     <a href="https://github.com/MyMiniFactory/Fast-Quadric-Mesh-Simplification"><img style="position: absolute; top: 0; right: 0; border: 0;" src="https://s3.amazonaws.com/github/ribbons/forkme_right_orange_ff7600.png" alt="Fork me on GitHub"></a>
262 | </body>
263 | 
264 | </html>
265 | 


--------------------------------------------------------------------------------
/vertex.h:
--------------------------------------------------------------------------------
1 | // https://github.com/mkeeter/fstl/blob/master/src/vertex.h
2 | 
3 | #ifndef VEC3_H
4 | #define VEC3_H
5 | 
6 | 
7 | #endif
8 | 


--------------------------------------------------------------------------------
/worker.js:
--------------------------------------------------------------------------------
 1 | self.addEventListener('message', function(e) {
 2 |     const file = e.data.blob;
 3 |     const percentage = e.data.percentage;
 4 |     const simplify_name = e.data.simplify_name;
 5 | 
 6 |     prepare_and_simplify(file, percentage, simplify_name);
 7 | }, false);
 8 | 
 9 | var Module = {
10 |     'print': function(text) {
11 |         self.postMessage({"log":text});
12 |     }
13 | };
14 | 
15 | self.importScripts("a.out.js");
16 | 
17 | let last_file_name = undefined;
18 | 
19 | function prepare_and_simplify(file, percentage, simplify_name) {
20 | 
21 |     var filename = file.name;
22 | 
23 |     // if simplify on the same file, don't even read the file
24 |     if (filename === last_file_name) {
25 |         console.log("skipping load and create data file");
26 |         simplify(filename, percentage, simplify_name);
27 |         return;
28 |     } else { // remove last file in memory
29 |         if (last_file_name !== undefined)
30 |             Module.FS_unlink(last_file_name);
31 |     }
32 | 
33 |     last_file_name = filename;
34 |     var fr = new FileReader();
35 |     fr.readAsArrayBuffer(file);
36 |     fr. onloadend = function (e) {
37 |         var data = new Uint8Array(fr.result);
38 |         Module.FS_createDataFile(".", filename, data, true, true);
39 |         simplify(filename, percentage, simplify_name);
40 |     }
41 | }
42 | 
43 | function simplify(filename, percentage, simplify_name) {
44 |     Module.ccall("simplify", // c function name
45 |         undefined, // return
46 |         ["string", "number", "string"], // param
47 |         [filename, percentage, simplify_name]
48 |     );
49 |     let out_bin = Module.FS_readFile(simplify_name);
50 |     // sla should work for binary stl
51 |     let file = new Blob([out_bin], {type: 'application/sla'});
52 |     self.postMessage({"blob":file});
53 | }
54 | 


--------------------------------------------------------------------------------