/******************************************************************************
 * Spine Runtimes License Agreement
 * Last updated January 1, 2020. Replaces all prior versions.
 *
 * Copyright (c) 2013-2020, Esoteric Software LLC
 *
 * Integration of the Spine Runtimes into software or otherwise creating
 * derivative works of the Spine Runtimes is permitted under the terms and
 * conditions of Section 2 of the Spine Editor License Agreement:
 * http://esotericsoftware.com/spine-editor-license
 *
 * Otherwise, it is permitted to integrate the Spine Runtimes into software
 * or otherwise create derivative works of the Spine Runtimes (collectively,
 * "Products"), provided that each user of the Products must obtain their own
 * Spine Editor license and redistribution of the Products in any form must
 * include this license and copyright notice.
 *
 * THE SPINE RUNTIMES ARE PROVIDED BY ESOTERIC SOFTWARE LLC "AS IS" AND ANY
 * EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
 * DISCLAIMED. IN NO EVENT SHALL ESOTERIC SOFTWARE LLC BE LIABLE FOR ANY
 * DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES,
 * BUSINESS INTERRUPTION, OR LOSS OF USE, DATA, OR PROFITS) HOWEVER CAUSED AND
 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
 * THE SPINE RUNTIMES, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 *****************************************************************************/

module spine.webgl {
	export const M00 = 0;
	export const M01 = 4;
	export const M02 = 8;
	export const M03 = 12;
	export const M10 = 1;
	export const M11 = 5;
	export const M12 = 9;
	export const M13 = 13;
	export const M20 = 2;
	export const M21 = 6;
	export const M22 = 10;
	export const M23 = 14;
	export const M30 = 3;
	export const M31 = 7;
	export const M32 = 11;
	export const M33 = 15;

	export class Matrix4 {
		temp: Float32Array = new Float32Array(16);
		values: Float32Array = new Float32Array(16);

		private static xAxis: Vector3 = null;
		private static yAxis: Vector3 = null;
		private static zAxis: Vector3 = null;
		private static tmpMatrix = new Matrix4();

		constructor () {
			let v = this.values;
			v[M00] = 1;
			v[M11] = 1;
			v[M22] = 1;
			v[M33] = 1;
		}

		set (values: ArrayLike<number>): Matrix4 {
			this.values.set(values);
			return this;
		}

		transpose (): Matrix4 {
			let t = this.temp;
			let v = this.values;
			t[M00] = v[M00];
			t[M01] = v[M10];
			t[M02] = v[M20];
			t[M03] = v[M30];
			t[M10] = v[M01];
			t[M11] = v[M11];
			t[M12] = v[M21];
			t[M13] = v[M31];
			t[M20] = v[M02];
			t[M21] = v[M12];
			t[M22] = v[M22];
			t[M23] = v[M32];
			t[M30] = v[M03];
			t[M31] = v[M13];
			t[M32] = v[M23];
			t[M33] = v[M33];
			return this.set(t);
		}

		identity (): Matrix4 {
			let v = this.values;
			v[M00] = 1;
			v[M01] = 0;
			v[M02] = 0;
			v[M03] = 0;
			v[M10] = 0;
			v[M11] = 1;
			v[M12] = 0;
			v[M13] = 0;
			v[M20] = 0;
			v[M21] = 0;
			v[M22] = 1;
			v[M23] = 0;
			v[M30] = 0;
			v[M31] = 0;
			v[M32] = 0;
			v[M33] = 1;
			return this;
		}

		invert (): Matrix4 {
			let v = this.values;
			let t = this.temp;
			let l_det = v[M30] * v[M21] * v[M12] * v[M03] - v[M20] * v[M31] * v[M12] * v[M03] - v[M30] * v[M11] * v[M22] * v[M03]
				+ v[M10] * v[M31] * v[M22] * v[M03] + v[M20] * v[M11] * v[M32] * v[M03] - v[M10] * v[M21] * v[M32] * v[M03]
				- v[M30] * v[M21] * v[M02] * v[M13] + v[M20] * v[M31] * v[M02] * v[M13] + v[M30] * v[M01] * v[M22] * v[M13]
				- v[M00] * v[M31] * v[M22] * v[M13] - v[M20] * v[M01] * v[M32] * v[M13] + v[M00] * v[M21] * v[M32] * v[M13]
				+ v[M30] * v[M11] * v[M02] * v[M23] - v[M10] * v[M31] * v[M02] * v[M23] - v[M30] * v[M01] * v[M12] * v[M23]
				+ v[M00] * v[M31] * v[M12] * v[M23] + v[M10] * v[M01] * v[M32] * v[M23] - v[M00] * v[M11] * v[M32] * v[M23]
				- v[M20] * v[M11] * v[M02] * v[M33] + v[M10] * v[M21] * v[M02] * v[M33] + v[M20] * v[M01] * v[M12] * v[M33]
				- v[M00] * v[M21] * v[M12] * v[M33] - v[M10] * v[M01] * v[M22] * v[M33] + v[M00] * v[M11] * v[M22] * v[M33];
			if (l_det == 0) throw new Error("non-invertible matrix");
			let inv_det = 1.0 / l_det;
			t[M00] = v[M12] * v[M23] * v[M31] - v[M13] * v[M22] * v[M31] + v[M13] * v[M21] * v[M32]
				- v[M11] * v[M23] * v[M32] - v[M12] * v[M21] * v[M33] + v[M11] * v[M22] * v[M33];
			t[M01] = v[M03] * v[M22] * v[M31] - v[M02] * v[M23] * v[M31] - v[M03] * v[M21] * v[M32]
				+ v[M01] * v[M23] * v[M32] + v[M02] * v[M21] * v[M33] - v[M01] * v[M22] * v[M33];
			t[M02] = v[M02] * v[M13] * v[M31] - v[M03] * v[M12] * v[M31] + v[M03] * v[M11] * v[M32]
				- v[M01] * v[M13] * v[M32] - v[M02] * v[M11] * v[M33] + v[M01] * v[M12] * v[M33];
			t[M03] = v[M03] * v[M12] * v[M21] - v[M02] * v[M13] * v[M21] - v[M03] * v[M11] * v[M22]
				+ v[M01] * v[M13] * v[M22] + v[M02] * v[M11] * v[M23] - v[M01] * v[M12] * v[M23];
			t[M10] = v[M13] * v[M22] * v[M30] - v[M12] * v[M23] * v[M30] - v[M13] * v[M20] * v[M32]
				+ v[M10] * v[M23] * v[M32] + v[M12] * v[M20] * v[M33] - v[M10] * v[M22] * v[M33];
			t[M11] = v[M02] * v[M23] * v[M30] - v[M03] * v[M22] * v[M30] + v[M03] * v[M20] * v[M32]
				- v[M00] * v[M23] * v[M32] - v[M02] * v[M20] * v[M33] + v[M00] * v[M22] * v[M33];
			t[M12] = v[M03] * v[M12] * v[M30] - v[M02] * v[M13] * v[M30] - v[M03] * v[M10] * v[M32]
				+ v[M00] * v[M13] * v[M32] + v[M02] * v[M10] * v[M33] - v[M00] * v[M12] * v[M33];
			t[M13] = v[M02] * v[M13] * v[M20] - v[M03] * v[M12] * v[M20] + v[M03] * v[M10] * v[M22]
				- v[M00] * v[M13] * v[M22] - v[M02] * v[M10] * v[M23] + v[M00] * v[M12] * v[M23];
			t[M20] = v[M11] * v[M23] * v[M30] - v[M13] * v[M21] * v[M30] + v[M13] * v[M20] * v[M31]
				- v[M10] * v[M23] * v[M31] - v[M11] * v[M20] * v[M33] + v[M10] * v[M21] * v[M33];
			t[M21] = v[M03] * v[M21] * v[M30] - v[M01] * v[M23] * v[M30] - v[M03] * v[M20] * v[M31]
				+ v[M00] * v[M23] * v[M31] + v[M01] * v[M20] * v[M33] - v[M00] * v[M21] * v[M33];
			t[M22] = v[M01] * v[M13] * v[M30] - v[M03] * v[M11] * v[M30] + v[M03] * v[M10] * v[M31]
				- v[M00] * v[M13] * v[M31] - v[M01] * v[M10] * v[M33] + v[M00] * v[M11] * v[M33];
			t[M23] = v[M03] * v[M11] * v[M20] - v[M01] * v[M13] * v[M20] - v[M03] * v[M10] * v[M21]
				+ v[M00] * v[M13] * v[M21] + v[M01] * v[M10] * v[M23] - v[M00] * v[M11] * v[M23];
			t[M30] = v[M12] * v[M21] * v[M30] - v[M11] * v[M22] * v[M30] - v[M12] * v[M20] * v[M31]
				+ v[M10] * v[M22] * v[M31] + v[M11] * v[M20] * v[M32] - v[M10] * v[M21] * v[M32];
			t[M31] = v[M01] * v[M22] * v[M30] - v[M02] * v[M21] * v[M30] + v[M02] * v[M20] * v[M31]
				- v[M00] * v[M22] * v[M31] - v[M01] * v[M20] * v[M32] + v[M00] * v[M21] * v[M32];
			t[M32] = v[M02] * v[M11] * v[M30] - v[M01] * v[M12] * v[M30] - v[M02] * v[M10] * v[M31]
				+ v[M00] * v[M12] * v[M31] + v[M01] * v[M10] * v[M32] - v[M00] * v[M11] * v[M32];
			t[M33] = v[M01] * v[M12] * v[M20] - v[M02] * v[M11] * v[M20] + v[M02] * v[M10] * v[M21]
				- v[M00] * v[M12] * v[M21] - v[M01] * v[M10] * v[M22] + v[M00] * v[M11] * v[M22];
			v[M00] = t[M00] * inv_det;
			v[M01] = t[M01] * inv_det;
			v[M02] = t[M02] * inv_det;
			v[M03] = t[M03] * inv_det;
			v[M10] = t[M10] * inv_det;
			v[M11] = t[M11] * inv_det;
			v[M12] = t[M12] * inv_det;
			v[M13] = t[M13] * inv_det;
			v[M20] = t[M20] * inv_det;
			v[M21] = t[M21] * inv_det;
			v[M22] = t[M22] * inv_det;
			v[M23] = t[M23] * inv_det;
			v[M30] = t[M30] * inv_det;
			v[M31] = t[M31] * inv_det;
			v[M32] = t[M32] * inv_det;
			v[M33] = t[M33] * inv_det;
			return this;
		}

		determinant (): number {
			let v = this.values;
			return v[M30] * v[M21] * v[M12] * v[M03] - v[M20] * v[M31] * v[M12] * v[M03] - v[M30] * v[M11] * v[M22] * v[M03]
				+ v[M10] * v[M31] * v[M22] * v[M03] + v[M20] * v[M11] * v[M32] * v[M03] - v[M10] * v[M21] * v[M32] * v[M03]
				- v[M30] * v[M21] * v[M02] * v[M13] + v[M20] * v[M31] * v[M02] * v[M13] + v[M30] * v[M01] * v[M22] * v[M13]
				- v[M00] * v[M31] * v[M22] * v[M13] - v[M20] * v[M01] * v[M32] * v[M13] + v[M00] * v[M21] * v[M32] * v[M13]
				+ v[M30] * v[M11] * v[M02] * v[M23] - v[M10] * v[M31] * v[M02] * v[M23] - v[M30] * v[M01] * v[M12] * v[M23]
				+ v[M00] * v[M31] * v[M12] * v[M23] + v[M10] * v[M01] * v[M32] * v[M23] - v[M00] * v[M11] * v[M32] * v[M23]
				- v[M20] * v[M11] * v[M02] * v[M33] + v[M10] * v[M21] * v[M02] * v[M33] + v[M20] * v[M01] * v[M12] * v[M33]
				- v[M00] * v[M21] * v[M12] * v[M33] - v[M10] * v[M01] * v[M22] * v[M33] + v[M00] * v[M11] * v[M22] * v[M33];
		}

		translate (x: number, y: number, z: number): Matrix4 {
			let v = this.values;
			v[M03] += x;
			v[M13] += y;
			v[M23] += z;
			return this;
		}

		copy (): Matrix4 {
			return new Matrix4().set(this.values);
		}

		projection (near: number, far: number, fovy: number, aspectRatio: number): Matrix4 {
			this.identity();
			let l_fd = (1.0 / Math.tan((fovy * (Math.PI / 180)) / 2.0));
			let l_a1 = (far + near) / (near - far);
			let l_a2 = (2 * far * near) / (near - far);
			let v = this.values;
			v[M00] = l_fd / aspectRatio;
			v[M10] = 0;
			v[M20] = 0;
			v[M30] = 0;
			v[M01] = 0;
			v[M11] = l_fd;
			v[M21] = 0;
			v[M31] = 0;
			v[M02] = 0;
			v[M12] = 0;
			v[M22] = l_a1;
			v[M32] = -1;
			v[M03] = 0;
			v[M13] = 0;
			v[M23] = l_a2;
			v[M33] = 0;
			return this;
		}

		ortho2d (x: number, y: number, width: number, height: number): Matrix4 {
			return this.ortho(x, x + width, y, y + height, 0, 1);
		}

		ortho (left: number, right: number, bottom: number, top: number, near: number, far: number): Matrix4 {
			this.identity();
			let x_orth = 2 / (right - left);
			let y_orth = 2 / (top - bottom);
			let z_orth = -2 / (far - near);

			let tx = -(right + left) / (right - left);
			let ty = -(top + bottom) / (top - bottom);
			let tz = -(far + near) / (far - near);

			let v = this.values;
			v[M00] = x_orth;
			v[M10] = 0;
			v[M20] = 0;
			v[M30] = 0;
			v[M01] = 0;
			v[M11] = y_orth;
			v[M21] = 0;
			v[M31] = 0;
			v[M02] = 0;
			v[M12] = 0;
			v[M22] = z_orth;
			v[M32] = 0;
			v[M03] = tx;
			v[M13] = ty;
			v[M23] = tz;
			v[M33] = 1;
			return this;
		}

		multiply (matrix: Matrix4): Matrix4 {
			let t = this.temp;
			let v = this.values;
			let m = matrix.values;
			t[M00] = v[M00] * m[M00] + v[M01] * m[M10] + v[M02] * m[M20] + v[M03] * m[M30];
			t[M01] = v[M00] * m[M01] + v[M01] * m[M11] + v[M02] * m[M21] + v[M03] * m[M31];
			t[M02] = v[M00] * m[M02] + v[M01] * m[M12] + v[M02] * m[M22] + v[M03] * m[M32];
			t[M03] = v[M00] * m[M03] + v[M01] * m[M13] + v[M02] * m[M23] + v[M03] * m[M33];
			t[M10] = v[M10] * m[M00] + v[M11] * m[M10] + v[M12] * m[M20] + v[M13] * m[M30];
			t[M11] = v[M10] * m[M01] + v[M11] * m[M11] + v[M12] * m[M21] + v[M13] * m[M31];
			t[M12] = v[M10] * m[M02] + v[M11] * m[M12] + v[M12] * m[M22] + v[M13] * m[M32];
			t[M13] = v[M10] * m[M03] + v[M11] * m[M13] + v[M12] * m[M23] + v[M13] * m[M33];
			t[M20] = v[M20] * m[M00] + v[M21] * m[M10] + v[M22] * m[M20] + v[M23] * m[M30];
			t[M21] = v[M20] * m[M01] + v[M21] * m[M11] + v[M22] * m[M21] + v[M23] * m[M31];
			t[M22] = v[M20] * m[M02] + v[M21] * m[M12] + v[M22] * m[M22] + v[M23] * m[M32];
			t[M23] = v[M20] * m[M03] + v[M21] * m[M13] + v[M22] * m[M23] + v[M23] * m[M33];
			t[M30] = v[M30] * m[M00] + v[M31] * m[M10] + v[M32] * m[M20] + v[M33] * m[M30];
			t[M31] = v[M30] * m[M01] + v[M31] * m[M11] + v[M32] * m[M21] + v[M33] * m[M31];
			t[M32] = v[M30] * m[M02] + v[M31] * m[M12] + v[M32] * m[M22] + v[M33] * m[M32];
			t[M33] = v[M30] * m[M03] + v[M31] * m[M13] + v[M32] * m[M23] + v[M33] * m[M33];
			return this.set(this.temp);
		}

		multiplyLeft (matrix: Matrix4): Matrix4 {
			let t = this.temp;
			let v = this.values;
			let m = matrix.values;
			t[M00] = m[M00] * v[M00] + m[M01] * v[M10] + m[M02] * v[M20] + m[M03] * v[M30];
			t[M01] = m[M00] * v[M01] + m[M01] * v[M11] + m[M02] * v[M21] + m[M03] * v[M31];
			t[M02] = m[M00] * v[M02] + m[M01] * v[M12] + m[M02] * v[M22] + m[M03] * v[M32];
			t[M03] = m[M00] * v[M03] + m[M01] * v[M13] + m[M02] * v[M23] + m[M03] * v[M33];
			t[M10] = m[M10] * v[M00] + m[M11] * v[M10] + m[M12] * v[M20] + m[M13] * v[M30];
			t[M11] = m[M10] * v[M01] + m[M11] * v[M11] + m[M12] * v[M21] + m[M13] * v[M31];
			t[M12] = m[M10] * v[M02] + m[M11] * v[M12] + m[M12] * v[M22] + m[M13] * v[M32];
			t[M13] = m[M10] * v[M03] + m[M11] * v[M13] + m[M12] * v[M23] + m[M13] * v[M33];
			t[M20] = m[M20] * v[M00] + m[M21] * v[M10] + m[M22] * v[M20] + m[M23] * v[M30];
			t[M21] = m[M20] * v[M01] + m[M21] * v[M11] + m[M22] * v[M21] + m[M23] * v[M31];
			t[M22] = m[M20] * v[M02] + m[M21] * v[M12] + m[M22] * v[M22] + m[M23] * v[M32];
			t[M23] = m[M20] * v[M03] + m[M21] * v[M13] + m[M22] * v[M23] + m[M23] * v[M33];
			t[M30] = m[M30] * v[M00] + m[M31] * v[M10] + m[M32] * v[M20] + m[M33] * v[M30];
			t[M31] = m[M30] * v[M01] + m[M31] * v[M11] + m[M32] * v[M21] + m[M33] * v[M31];
			t[M32] = m[M30] * v[M02] + m[M31] * v[M12] + m[M32] * v[M22] + m[M33] * v[M32];
			t[M33] = m[M30] * v[M03] + m[M31] * v[M13] + m[M32] * v[M23] + m[M33] * v[M33];
			return this.set(this.temp);
		}

		lookAt (position: Vector3, direction: Vector3, up: Vector3) {
			Matrix4.initTemps();
			let xAxis = Matrix4.xAxis, yAxis = Matrix4.yAxis, zAxis = Matrix4.zAxis;
			zAxis.setFrom(direction).normalize();
			xAxis.setFrom(direction).normalize();
			xAxis.cross(up).normalize();
			yAxis.setFrom(xAxis).cross(zAxis).normalize();
			this.identity();
			let val = this.values;
			val[M00] = xAxis.x;
			val[M01] = xAxis.y;
			val[M02] = xAxis.z;
			val[M10] = yAxis.x;
			val[M11] = yAxis.y;
			val[M12] = yAxis.z;
			val[M20] = -zAxis.x;
			val[M21] = -zAxis.y;
			val[M22] = -zAxis.z;

			Matrix4.tmpMatrix.identity();
			Matrix4.tmpMatrix.values[M03] = -position.x;
			Matrix4.tmpMatrix.values[M13] = -position.y;
			Matrix4.tmpMatrix.values[M23] = -position.z;
			this.multiply(Matrix4.tmpMatrix)

			return this;
		}

		static initTemps() {
			if (Matrix4.xAxis === null) Matrix4.xAxis = new Vector3();
			if (Matrix4.yAxis === null) Matrix4.yAxis = new Vector3();
			if (Matrix4.zAxis === null) Matrix4.zAxis = new Vector3();
		}
	}
}