/// @description Insert description here
#region params
	domain           = noone;
	particlePosBuff  = noone;

	width            = 0;
	height           = 0;
	particleSize     = 0;
	density          = 0;
	viscosity        = 0;
	friction         = 0;
	maxParticles     = 0;
	
	numParticles     = 0;
	velocityDamping  = 0.9;

	dt               = 0.1;
	iteration        = 8;
	numPressureIters = 2;
	numParticleIters = 2;
	
	g                = 1;
	flipRatio        = 0.8;
	overRelaxation   = 1.5;
	
	obstracles       = [];
	wallCollide      = true;
	wallElasticity   = 1;
	
	skip_incompressible = false;
	
	particleRadius   = 0;
#endregion

function init(width, height, particleSize, density, maxParticles) { #region domain init
	particlePos  = array_create(maxParticles * 2);
	particleHist = array_create(maxParticles * 2);
	particleLife = array_create(maxParticles);
	obstracles   = [];
	numParticles = 0;
	
	if(domain            != noone        && 
	   self.width        == width        && 
	   self.height       == height       && 
	   self.particleSize == particleSize && 
	   self.density      == density      && 
	   self.maxParticles == maxParticles) {
		
		FLIP_resetDomain(domain);
		return;
	}
	
	self.width        = width       ;
	self.height       = height      ;
	self.particleSize = particleSize;
	self.density      = density     ;
	self.maxParticles = maxParticles;
	
	particlePosBuff   = buffer_create(maxParticles * 2 * 8, buffer_grow, 8);
	particleLifeBuff  = buffer_create(maxParticles * 8, buffer_grow, 8);
	
	aPosBuff  = buffer_get_address(particlePosBuff);
	aLifeBuff = buffer_get_address(particleLifeBuff);
	
	domain            = FLIP_initDomain(width, height, particleSize, density, maxParticles);
	particleRadius    = FLIP_getParticleRadius(domain);
	
	cellX = floor(width  / particleSize) + 1;
	cellY = floor(height / particleSize) + 1;
} #endregion

function update() { #region
	FLIP_setQuality(		 domain, iteration, numPressureIters, numParticleIters);
	FLIP_setGravity(		 domain, g);
	FLIP_setViscosity(		 domain, viscosity);
	FLIP_setFriction(		 domain, power(1 - friction, 0.025));
	FLIP_setFlipRatio(		 domain, flipRatio);
	FLIP_setVelocityDamping( domain, velocityDamping);
	FLIP_setOverRelaxation(	 domain, overRelaxation);
	FLIP_setWallCollisions(  domain, wallCollide, wallElasticity);
} #endregion

function step() { #region
	FLIP_resetDensity(domain);
	
	if(skip_incompressible) { 
		FLIP_setTimeStep(domain, dt);
		repeat(iteration) {
			FLIP_simulate_integrateParticles(domain);
			FLIP_simulate_pushParticlesApart(domain);
			FLIP_simulate_handleParticleCollisions(domain);
			//FLIP_simulate_transferVelocities(domain, 1);
			//FLIP_simulate_updateParticleDensity(domain);
			//FLIP_simulate_solveIncompressibility(domain);
			//FLIP_simulate_transferVelocities(domain, 0);
		}
	} else 
		FLIP_simulate(domain, dt);
	
	FLIP_setParticleBuffer(domain, aPosBuff, aLifeBuff);
	
	buffer_seek(particlePosBuff,  buffer_seek_start, 0);
	buffer_seek(particleLifeBuff, buffer_seek_start, 0);
	
	for(var i = 0; i < maxParticles * 2; i++) {
		particleHist[i] = particlePos[i];
		particlePos[i]  = buffer_read(particlePosBuff, buffer_f64);
	}
	
	for(var i = 0; i < maxParticles; i++)
		particleLife[i] = buffer_read(particleLifeBuff, buffer_f64);
} #endregion