@chemistry/math 2.9.0 → 3.0.1

package/README.md

@@ -1,27 +1,29 @@
 # @chemistry/math
+
 [![npm version](https://badge.fury.io/js/%40chemistry%2Fmath.svg)](https://badge.fury.io/js/%40chemistry%2Fmath)
-[![Build Status](https://travis-ci.com/chemistry/chemical-libraries.svg?branch=master)](https://travis-ci.org/chemistry/chemical-libraries)
-[![codecov](https://codecov.io/gh/chemistry/chemical-libraries/branch/master/graph/badge.svg)](https://codecov.io/gh/chemistry/chemical-libraries)
 [![License: MIT](https://img.shields.io/badge/License-MIT-gren.svg)](https://opensource.org/licenses/MIT)
 
 Simple linear algebra Math library to support [chemistry js project](http://vreshch.com/chemistry-js.html)
 
 ## Install
+
 ```bash
 npm install @chemistry/math
 ```
 
 ## Include following clases:
-  * Vec3
-  * Matrix3x3
-  * Matrix3x4
-  * Transform3d
-  * Quaternion
+
+- Vec3
+- Matrix3x3
+- Matrix3x4
+- Transform3d
+- Quaternion
 
 ## Getting started:
+
 ```javascript
 import { V, ChemElementData } from '@chemistry/math';
-const matrix = new Matrix3x3([1, 0, 0,  0, 2, 0,  0, 0, 3]);
+const matrix = new Matrix3x3([1, 0, 0, 0, 2, 0, 0, 0, 3]);
 const vector = new Vec3(2, 3, 1);
 
 const res = matrix.project(vector);
@@ -30,4 +32,5 @@ console.log(res);
 ```
 
 ## Commands:
-  * Build project: `npm run build`
+
+- Build project: `npm run build`

package/dist/index.d.ts

@@ -1,5 +1,6 @@
-export * from "./matrix3x3";
-export * from "./matrix3x4";
-export * from "./transform3d";
-export * from "./vec3";
-export * from "./vec2";
+export * from './matrix3x3';
+export * from './matrix3x4';
+export * from './transform3d';
+export * from './vec3';
+export * from './vec2';
+//# sourceMappingURL=index.d.ts.map

package/dist/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,cAAc,aAAa,CAAC;AAC5B,cAAc,aAAa,CAAC;AAC5B,cAAc,eAAe,CAAC;AAC9B,cAAc,QAAQ,CAAC;AACvB,cAAc,QAAQ,CAAC"}

package/dist/index.js

@@ -0,0 +1,6 @@
+export * from './matrix3x3';
+export * from './matrix3x4';
+export * from './transform3d';
+export * from './vec3';
+export * from './vec2';
+//# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,cAAc,aAAa,CAAC;AAC5B,cAAc,aAAa,CAAC;AAC5B,cAAc,eAAe,CAAC;AAC9B,cAAc,QAAQ,CAAC;AACvB,cAAc,QAAQ,CAAC"}

package/dist/matrix3x3.d.ts

@@ -1,5 +1,5 @@
-import { ICloneable, IEquatable } from "@chemistry/common";
-import { Vec3 } from "./vec3";
+import { ICloneable, IEquatable } from '@chemistry/common';
+import { Vec3 } from './vec3';
 export declare class Matrix3x3 implements ICloneable<Matrix3x3>, IEquatable<Matrix3x3> {
     static equals(m1: Matrix3x3, m2: Matrix3x3): boolean;
     static inverse(matrix: Matrix3x3): Matrix3x3;
@@ -36,3 +36,4 @@ export declare class Matrix3x3 implements ICloneable<Matrix3x3>, IEquatable<Matr
      */
     toString(): string;
 }
+//# sourceMappingURL=matrix3x3.d.ts.map

package/dist/matrix3x3.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"matrix3x3.d.ts","sourceRoot":"","sources":["../src/matrix3x3.ts"],"names":[],"mappings":"AAAA,OAAO,EAAW,UAAU,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACpE,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAE9B,qBAAa,SAAU,YAAW,UAAU,CAAC,SAAS,CAAC,EAAE,UAAU,CAAC,SAAS,CAAC;WAC9D,MAAM,CAAC,EAAE,EAAE,SAAS,EAAE,EAAE,EAAE,SAAS,GAAG,OAAO;WAc7C,OAAO,CAAC,MAAM,EAAE,SAAS,GAAG,SAAS;IAsBnD,SAAS,CAAC,QAAQ,EAAE,MAAM,EAAE,CAAC;gBAEjB,QAAQ,CAAC,EAAE,MAAM,EAAE;IAI/B;;OAEG;IACI,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;IAO5C;;OAEG;IACI,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAOjC,KAAK,CAAC,IAAI,EAAE,MAAM,GAAG,SAAS;IAc9B,OAAO,CAAC,IAAI,EAAE,IAAI,GAAG,IAAI;IAQhC;;OAEG;IACI,QAAQ,CAAC,IAAI,EAAE,SAAS,GAAG,SAAS;IAgC3C;;OAEG;IACI,WAAW,IAAI,MAAM;IAQ5B;;OAEG;IACI,KAAK,IAAI,SAAS;IAIzB;;OAEG;IACI,MAAM,CAAC,MAAM,EAAE,SAAS,GAAG,OAAO;IAIzC;;OAEG;IACI,QAAQ,IAAI,MAAM;CAiB1B"}

package/dist/matrix3x3.js

@@ -0,0 +1,143 @@
+import { EPSILON } from '@chemistry/common';
+import { Vec3 } from './vec3';
+export class Matrix3x3 {
+    static equals(m1, m2) {
+        return (Math.abs(m1.elements[0] - m2.elements[0]) < EPSILON &&
+            Math.abs(m1.elements[1] - m2.elements[1]) < EPSILON &&
+            Math.abs(m1.elements[2] - m2.elements[2]) < EPSILON &&
+            Math.abs(m1.elements[3] - m2.elements[3]) < EPSILON &&
+            Math.abs(m1.elements[4] - m2.elements[4]) < EPSILON &&
+            Math.abs(m1.elements[5] - m2.elements[5]) < EPSILON &&
+            Math.abs(m1.elements[6] - m2.elements[6]) < EPSILON &&
+            Math.abs(m1.elements[7] - m2.elements[7]) < EPSILON &&
+            Math.abs(m1.elements[8] - m2.elements[8]) < EPSILON);
+    }
+    static inverse(matrix) {
+        const d = matrix.determinant();
+        if (d === 0) {
+            throw new Error('Matrix not invertable');
+        }
+        return new Matrix3x3([
+            (matrix.get(1, 1) * matrix.get(2, 2) - matrix.get(1, 2) * matrix.get(2, 1)) / d,
+            (-matrix.get(0, 1) * matrix.get(2, 2) + matrix.get(0, 2) * matrix.get(2, 1)) / d,
+            (matrix.get(0, 1) * matrix.get(1, 2) - matrix.get(0, 2) * matrix.get(1, 1)) / d,
+            (-matrix.get(1, 0) * matrix.get(2, 2) + matrix.get(1, 2) * matrix.get(2, 0)) / d,
+            (matrix.get(0, 0) * matrix.get(2, 2) - matrix.get(0, 2) * matrix.get(2, 0)) / d,
+            (-matrix.get(0, 0) * matrix.get(1, 2) + matrix.get(0, 2) * matrix.get(1, 0)) / d,
+            (matrix.get(1, 0) * matrix.get(2, 1) - matrix.get(1, 1) * matrix.get(2, 0)) / d,
+            (-matrix.get(0, 0) * matrix.get(2, 1) + matrix.get(0, 1) * matrix.get(2, 0)) / d,
+            (matrix.get(0, 0) * matrix.get(1, 1) - matrix.get(0, 1) * matrix.get(1, 0)) / d,
+        ]);
+    }
+    constructor(elements) {
+        this.elements = elements?.length === 9 ? elements.slice(0) : [1, 0, 0, 0, 1, 0, 0, 0, 1];
+    }
+    /**
+     * Set value to specified element of the matrix
+     */
+    set(i, j, val) {
+        if (i > 3 || i < 0 || j > 3 || j < 0) {
+            throw new Error('Wrong coefficients');
+        }
+        this.elements[i * 3 + j] = val;
+    }
+    /**
+     * return specified element of the matrix
+     */
+    get(i, j) {
+        if (i > 3 || i < 0 || j > 3 || j < 0) {
+            throw new Error('Wrong coefficients');
+        }
+        return this.elements[i * 3 + j];
+    }
+    scale(mult) {
+        return new Matrix3x3([
+            this.elements[0] * mult,
+            this.elements[1] * mult,
+            this.elements[2] * mult,
+            this.elements[3] * mult,
+            this.elements[4] * mult,
+            this.elements[5] * mult,
+            this.elements[6] * mult,
+            this.elements[7] * mult,
+            this.elements[8] * mult,
+        ]);
+    }
+    project(mult) {
+        return new Vec3(this.elements[0] * mult.x + this.elements[1] * mult.y + this.elements[2] * mult.z, this.elements[3] * mult.x + this.elements[4] * mult.y + this.elements[5] * mult.z, this.elements[6] * mult.x + this.elements[7] * mult.y + this.elements[8] * mult.z);
+    }
+    /**
+     * Returns matrix multiplication result to scalar, vector or another matrix
+     */
+    multiply(mult) {
+        return new Matrix3x3([
+            this.elements[0] * mult.elements[0] +
+                this.elements[3] * mult.elements[1] +
+                this.elements[6] * mult.elements[2],
+            this.elements[1] * mult.elements[0] +
+                this.elements[4] * mult.elements[1] +
+                this.elements[7] * mult.elements[2],
+            this.elements[2] * mult.elements[0] +
+                this.elements[5] * mult.elements[1] +
+                this.elements[8] * mult.elements[2],
+            this.elements[0] * mult.elements[3] +
+                this.elements[3] * mult.elements[4] +
+                this.elements[6] * mult.elements[5],
+            this.elements[1] * mult.elements[3] +
+                this.elements[4] * mult.elements[4] +
+                this.elements[7] * mult.elements[5],
+            this.elements[2] * mult.elements[3] +
+                this.elements[5] * mult.elements[4] +
+                this.elements[8] * mult.elements[5],
+            this.elements[0] * mult.elements[6] +
+                this.elements[3] * mult.elements[7] +
+                this.elements[6] * mult.elements[8],
+            this.elements[1] * mult.elements[6] +
+                this.elements[4] * mult.elements[7] +
+                this.elements[7] * mult.elements[8],
+            this.elements[2] * mult.elements[6] +
+                this.elements[5] * mult.elements[7] +
+                this.elements[8] * mult.elements[8],
+        ]);
+    }
+    /**
+     * Returns matrix determinant
+     */
+    determinant() {
+        let result = this.get(0, 0) * (this.get(1, 1) * this.get(2, 2) - this.get(2, 1) * this.get(1, 2));
+        result -= this.get(1, 0) * (this.get(0, 1) * this.get(2, 2) - this.get(2, 1) * this.get(0, 2));
+        result += this.get(2, 0) * (this.get(0, 1) * this.get(1, 2) - this.get(1, 1) * this.get(0, 2));
+        return result;
+    }
+    /**
+     * Clone the matrix
+     */
+    clone() {
+        return new Matrix3x3(this.elements.slice(0));
+    }
+    /**
+     * Compare 2 matrixes
+     */
+    equals(matrix) {
+        return Matrix3x3.equals(this, matrix);
+    }
+    /**
+     * Returns string representation of the matrix
+     */
+    toString() {
+        return ('(' +
+            [
+                this.elements[0].toFixed(3),
+                this.elements[1].toFixed(3),
+                this.elements[2].toFixed(3),
+                this.elements[3].toFixed(3),
+                this.elements[4].toFixed(3),
+                this.elements[5].toFixed(3),
+                this.elements[6].toFixed(3),
+                this.elements[7].toFixed(3),
+                this.elements[8].toFixed(3),
+            ].join(',') +
+            ')');
+    }
+}
+//# sourceMappingURL=matrix3x3.js.map

package/dist/matrix3x3.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"matrix3x3.js","sourceRoot":"","sources":["../src/matrix3x3.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAA0B,MAAM,mBAAmB,CAAC;AACpE,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAE9B,MAAM,OAAO,SAAS;IACb,MAAM,CAAC,MAAM,CAAC,EAAa,EAAE,EAAa;QAC/C,OAAO,CACL,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;YACnD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CACpD,CAAC;IACJ,CAAC;IAEM,MAAM,CAAC,OAAO,CAAC,MAAiB;QACrC,MAAM,CAAC,GAAG,MAAM,CAAC,WAAW,EAAE,CAAC;QAC/B,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACZ,MAAM,IAAI,KAAK,CAAC,uBAAuB,CAAC,CAAC;QAC3C,CAAC;QAED,OAAO,IAAI,SAAS,CAAC;YACnB,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAC/E,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAChF,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAC/E,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAChF,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAC/E,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAChF,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAC/E,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;YAChF,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;SAChF,CAAC,CAAC;IACL,CAAC;IAOD,YAAY,QAAmB;QAC7B,IAAI,CAAC,QAAQ,GAAG,QAAQ,EAAE,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3F,CAAC;IAED;;OAEG;IACI,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,GAAW;QAC1C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACrC,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;QACxC,CAAC;QACD,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;IACjC,CAAC;IAED;;OAEG;IACI,GAAG,CAAC,CAAS,EAAE,CAAS;QAC7B,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACrC,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;QACxC,CAAC;QACD,OAAO,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAClC,CAAC;IAEM,KAAK,CAAC,IAAY;QACvB,OAAO,IAAI,SAAS,CAAC;YACnB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;YACvB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI;SACxB,CAAC,CAAC;IACL,CAAC;IAEM,OAAO,CAAC,IAAU;QACvB,OAAO,IAAI,IAAI,CACb,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,EACjF,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,EACjF,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAClF,CAAC;IACJ,CAAC;IAED;;OAEG;IACI,QAAQ,CAAC,IAAe;QAC7B,OAAO,IAAI,SAAS,CAAC;YACnB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACjC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;SACtC,CAAC,CAAC;IACL,CAAC;IAED;;OAEG;IACI,WAAW;QAChB,IAAI,MAAM,GACR,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACvF,MAAM,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC/F,MAAM,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC/F,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;OAEG;IACI,KAAK;QACV,OAAO,IAAI,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/C,CAAC;IAED;;OAEG;IACI,MAAM,CAAC,MAAiB;QAC7B,OAAO,SAAS,CAAC,MAAM,CAAC,IAAI,EAAE,MAAM,CAAC,CAAC;IACxC,CAAC;IAED;;OAEG;IACI,QAAQ;QACb,OAAO,CACL,GAAG;YACH;gBACE,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;gBAC3B,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;aAC5B,CAAC,IAAI,CAAC,GAAG,CAAC;YACX,GAAG,CACJ,CAAC;IACJ,CAAC;CACF"}

package/dist/matrix3x4.d.ts

@@ -1,5 +1,5 @@
-import { ICloneable, IEquatable } from "@chemistry/common";
-import { Vec3 } from "./vec3";
+import { ICloneable, IEquatable } from '@chemistry/common';
+import { Vec3 } from './vec3';
 /**
  * Simple Class to work with special type of matrix.
  * | x1 x2  x3  x4 |
@@ -43,3 +43,4 @@ export declare class Matrix3x4 implements ICloneable<Matrix3x4>, IEquatable<Matr
      */
     toSymetryCode(): string;
 }
+//# sourceMappingURL=matrix3x4.d.ts.map

package/dist/matrix3x4.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"matrix3x4.d.ts","sourceRoot":"","sources":["../src/matrix3x4.ts"],"names":[],"mappings":"AAAA,OAAO,EAAW,UAAU,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACpE,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAE9B;;;;;;GAMG;AACH,qBAAa,SAAU,YAAW,UAAU,CAAC,SAAS,CAAC,EAAE,UAAU,CAAC,SAAS,CAAC;WAC9D,GAAG,CAAC,OAAO,EAAE,SAAS,EAAE,OAAO,EAAE,SAAS,GAAG,SAAS;WAmBtD,GAAG,CAAC,OAAO,EAAE,SAAS,EAAE,OAAO,EAAE,SAAS,GAAG,SAAS;WAmBtD,OAAO,CAAC,MAAM,EAAE,SAAS,GAAG,SAAS;WAgJrC,GAAG,CAAC,MAAM,EAAE,SAAS,GAAG,OAAO;WAiB/B,MAAM,CAAC,OAAO,EAAE,SAAS,EAAE,OAAO,EAAE,SAAS,GAAG,OAAO;IAiBrE,SAAS,CAAC,QAAQ,EAAE,MAAM,EAAE,CAAC;gBAEjB,QAAQ,CAAC,EAAE,MAAM,EAAE;IAK/B;;OAEG;IACI,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;IAO5C;;;;;OAKG;IACI,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;IAOxB,GAAG,CAAC,MAAM,EAAE,SAAS,GAAG,SAAS;IAIjC,GAAG,CAAC,MAAM,EAAE,SAAS,GAAG,SAAS;IAIjC,SAAS,CAAC,CAAC,EAAE,IAAI,GAAG,SAAS;IAQ7B,KAAK,CAAC,IAAI,EAAE,MAAM,GAAG,SAAS;IAiB9B,OAAO,CAAC,IAAI,EAAE,IAAI,GAAG,IAAI;IAiBzB,QAAQ,CAAC,IAAI,EAAE,SAAS;IA0DxB,OAAO,IAAI,SAAS;IAIpB,GAAG,IAAI,OAAO;IAId,MAAM,CAAC,MAAM,EAAE,SAAS,GAAG,OAAO;IAIlC,KAAK,IAAI,SAAS;IAIlB,QAAQ,IAAI,MAAM;IAIzB;;;OAGG;IACI,aAAa,IAAI,MAAM;CAmD/B"}

package/dist/matrix3x4.js

@@ -0,0 +1,393 @@
+import { EPSILON } from '@chemistry/common';
+import { Vec3 } from './vec3';
+/**
+ * Simple Class to work with special type of matrix.
+ * | x1 x2  x3  x4 |
+ * | x5 x6  x7  x8 |
+ * | x8 x10 x11 x12|
+ * | 0  0   0   1  |
+ */
+export class Matrix3x4 {
+    static add(matrix1, matrix2) {
+        const m1 = matrix1.elements;
+        const m2 = matrix2.elements;
+        return new Matrix3x4([
+            m1[0] + m2[0],
+            m1[1] + m2[1],
+            m1[2] + m2[2],
+            m1[3] + m2[3],
+            m1[4] + m2[4],
+            m1[5] + m2[5],
+            m1[6] + m2[6],
+            m1[7] + m2[7],
+            m1[8] + m2[8],
+            m1[9] + m2[9],
+            m1[10] + m2[10],
+            m1[11] + m2[11],
+        ]);
+    }
+    static sub(matrix1, matrix2) {
+        const m1 = matrix1.elements;
+        const m2 = matrix2.elements;
+        return new Matrix3x4([
+            m1[0] - m2[0],
+            m1[1] - m2[1],
+            m1[2] - m2[2],
+            m1[3] - m2[3],
+            m1[4] - m2[4],
+            m1[5] - m2[5],
+            m1[6] - m2[6],
+            m1[7] - m2[7],
+            m1[8] - m2[8],
+            m1[9] - m2[9],
+            m1[10] - m2[10],
+            m1[11] - m2[11],
+        ]);
+    }
+    static inverse(matrix) {
+        const r = [];
+        let i;
+        const m = matrix.elements.slice(0, 12);
+        m[12] = 0;
+        m[13] = 0;
+        m[14] = 0;
+        m[15] = 1;
+        r[0] =
+            m[5] * m[10] * m[15] -
+                m[5] * m[14] * m[11] -
+                m[6] * m[9] * m[15] +
+                m[6] * m[13] * m[11] +
+                m[7] * m[9] * m[14] -
+                m[7] * m[13] * m[10];
+        r[1] =
+            -m[1] * m[10] * m[15] +
+                m[1] * m[14] * m[11] +
+                m[2] * m[9] * m[15] -
+                m[2] * m[13] * m[11] -
+                m[3] * m[9] * m[14] +
+                m[3] * m[13] * m[10];
+        r[2] =
+            m[1] * m[6] * m[15] -
+                m[1] * m[14] * m[7] -
+                m[2] * m[5] * m[15] +
+                m[2] * m[13] * m[7] +
+                m[3] * m[5] * m[14] -
+                m[3] * m[13] * m[6];
+        r[3] =
+            -m[1] * m[6] * m[11] +
+                m[1] * m[10] * m[7] +
+                m[2] * m[5] * m[11] -
+                m[2] * m[9] * m[7] -
+                m[3] * m[5] * m[10] +
+                m[3] * m[9] * m[6];
+        r[4] =
+            -m[4] * m[10] * m[15] +
+                m[4] * m[14] * m[11] +
+                m[6] * m[8] * m[15] -
+                m[6] * m[12] * m[11] -
+                m[7] * m[8] * m[14] +
+                m[7] * m[12] * m[10];
+        r[5] =
+            m[0] * m[10] * m[15] -
+                m[0] * m[14] * m[11] -
+                m[2] * m[8] * m[15] +
+                m[2] * m[12] * m[11] +
+                m[3] * m[8] * m[14] -
+                m[3] * m[12] * m[10];
+        r[6] =
+            -m[0] * m[6] * m[15] +
+                m[0] * m[14] * m[7] +
+                m[2] * m[4] * m[15] -
+                m[2] * m[12] * m[7] -
+                m[3] * m[4] * m[14] +
+                m[3] * m[12] * m[6];
+        r[7] =
+            m[0] * m[6] * m[11] -
+                m[0] * m[10] * m[7] -
+                m[2] * m[4] * m[11] +
+                m[2] * m[8] * m[7] +
+                m[3] * m[4] * m[10] -
+                m[3] * m[8] * m[6];
+        r[8] =
+            m[4] * m[9] * m[15] -
+                m[4] * m[13] * m[11] -
+                m[5] * m[8] * m[15] +
+                m[5] * m[12] * m[11] +
+                m[7] * m[8] * m[13] -
+                m[7] * m[12] * m[9];
+        r[9] =
+            -m[0] * m[9] * m[15] +
+                m[0] * m[13] * m[11] +
+                m[1] * m[8] * m[15] -
+                m[1] * m[12] * m[11] -
+                m[3] * m[8] * m[13] +
+                m[3] * m[12] * m[9];
+        r[10] =
+            m[0] * m[5] * m[15] -
+                m[0] * m[13] * m[7] -
+                m[1] * m[4] * m[15] +
+                m[1] * m[12] * m[7] +
+                m[3] * m[4] * m[13] -
+                m[3] * m[12] * m[5];
+        r[11] =
+            -m[0] * m[5] * m[11] +
+                m[0] * m[9] * m[7] +
+                m[1] * m[4] * m[11] -
+                m[1] * m[8] * m[7] -
+                m[3] * m[4] * m[9] +
+                m[3] * m[8] * m[5];
+        r[12] =
+            -m[4] * m[9] * m[14] +
+                m[4] * m[13] * m[10] +
+                m[5] * m[8] * m[14] -
+                m[5] * m[12] * m[10] -
+                m[6] * m[8] * m[13] +
+                m[6] * m[12] * m[9];
+        r[13] =
+            m[0] * m[9] * m[14] -
+                m[0] * m[13] * m[10] -
+                m[1] * m[8] * m[14] +
+                m[1] * m[12] * m[10] +
+                m[2] * m[8] * m[13] -
+                m[2] * m[12] * m[9];
+        r[14] =
+            -m[0] * m[5] * m[14] +
+                m[0] * m[13] * m[6] +
+                m[1] * m[4] * m[14] -
+                m[1] * m[12] * m[6] -
+                m[2] * m[4] * m[13] +
+                m[2] * m[12] * m[5];
+        r[15] =
+            m[0] * m[5] * m[10] -
+                m[0] * m[9] * m[6] -
+                m[1] * m[4] * m[10] +
+                m[1] * m[8] * m[6] +
+                m[2] * m[4] * m[9] -
+                m[2] * m[8] * m[5];
+        const det = m[0] * r[0] + m[1] * r[4] + m[2] * r[8] + m[3] * r[12];
+        for (i = 0; i < 16; i++) {
+            r[i] /= det;
+        }
+        return new Matrix3x4(r.slice(0, 12));
+    }
+    static isE(matrix) {
+        return (Math.abs(matrix.elements[0] - 1) < EPSILON &&
+            Math.abs(matrix.elements[1] - 0) < EPSILON &&
+            Math.abs(matrix.elements[2] - 0) < EPSILON &&
+            Math.abs(matrix.elements[3] - 0) < EPSILON &&
+            Math.abs(matrix.elements[4] - 0) < EPSILON &&
+            Math.abs(matrix.elements[5] - 1) < EPSILON &&
+            Math.abs(matrix.elements[6] - 0) < EPSILON &&
+            Math.abs(matrix.elements[7] - 0) < EPSILON &&
+            Math.abs(matrix.elements[8] - 0) < EPSILON &&
+            Math.abs(matrix.elements[9] - 0) < EPSILON &&
+            Math.abs(matrix.elements[10] - 1) < EPSILON &&
+            Math.abs(matrix.elements[11] - 0) < EPSILON);
+    }
+    static equals(matrix1, matrix2) {
+        return (Math.abs(matrix1.elements[0] - matrix2.elements[0]) < EPSILON &&
+            Math.abs(matrix1.elements[1] - matrix2.elements[1]) < EPSILON &&
+            Math.abs(matrix1.elements[2] - matrix2.elements[2]) < EPSILON &&
+            Math.abs(matrix1.elements[3] - matrix2.elements[3]) < EPSILON &&
+            Math.abs(matrix1.elements[4] - matrix2.elements[4]) < EPSILON &&
+            Math.abs(matrix1.elements[5] - matrix2.elements[5]) < EPSILON &&
+            Math.abs(matrix1.elements[6] - matrix2.elements[6]) < EPSILON &&
+            Math.abs(matrix1.elements[7] - matrix2.elements[7]) < EPSILON &&
+            Math.abs(matrix1.elements[8] - matrix2.elements[8]) < EPSILON &&
+            Math.abs(matrix1.elements[9] - matrix2.elements[9]) < EPSILON &&
+            Math.abs(matrix1.elements[10] - matrix2.elements[10]) < EPSILON &&
+            Math.abs(matrix1.elements[11] - matrix2.elements[11]) < EPSILON);
+    }
+    constructor(elements) {
+        this.elements =
+            elements?.length === 12 ? elements.slice(0) : [1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0];
+    }
+    /**
+     * Set value to specified element of the matrix
+     */
+    set(i, j, val) {
+        if (i > 3 || i < 0 || j > 4 || j < 0) {
+            throw new Error('Wrong coefficients');
+        }
+        this.elements[i * 4 + j] = val;
+    }
+    /**
+     * return specified element of the matrix
+     * @param  {Number} i row
+     * @param  {Number} j column
+     * @return {Number}
+     */
+    get(i, j) {
+        if (i > 3 || i < 0 || j > 4 || j < 0) {
+            throw new Error('Incorrect  coeficients');
+        }
+        return this.elements[i * 4 + j];
+    }
+    add(matrix) {
+        return Matrix3x4.add(this, matrix);
+    }
+    sub(matrix) {
+        return Matrix3x4.sub(this, matrix);
+    }
+    translate(v) {
+        const elements = this.elements.slice(0);
+        elements[3] = elements[3] + v.x;
+        elements[7] = elements[7] + v.y;
+        elements[11] = elements[11] + v.z;
+        return new Matrix3x4(elements);
+    }
+    scale(mult) {
+        return new Matrix3x4([
+            this.elements[0] * mult,
+            this.elements[1] * mult,
+            this.elements[2] * mult,
+            this.elements[3] * mult,
+            this.elements[4] * mult,
+            this.elements[5] * mult,
+            this.elements[6] * mult,
+            this.elements[7] * mult,
+            this.elements[8] * mult,
+            this.elements[9] * mult,
+            this.elements[10] * mult,
+            this.elements[11] * mult,
+        ]);
+    }
+    project(mult) {
+        return new Vec3(this.elements[0] * mult.x +
+            this.elements[1] * mult.y +
+            this.elements[2] * mult.z +
+            this.elements[3], this.elements[4] * mult.x +
+            this.elements[5] * mult.y +
+            this.elements[6] * mult.z +
+            this.elements[7], this.elements[8] * mult.x +
+            this.elements[9] * mult.y +
+            this.elements[10] * mult.z +
+            this.elements[11]);
+    }
+    multiply(mult) {
+        const elements = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0];
+        elements[0] =
+            mult.elements[0] * this.elements[0] +
+                mult.elements[4] * this.elements[1] +
+                mult.elements[8] * this.elements[2];
+        elements[1] =
+            mult.elements[1] * this.elements[0] +
+                mult.elements[5] * this.elements[1] +
+                mult.elements[9] * this.elements[2];
+        elements[2] =
+            mult.elements[2] * this.elements[0] +
+                mult.elements[6] * this.elements[1] +
+                mult.elements[10] * this.elements[2];
+        elements[3] =
+            mult.elements[3] * this.elements[0] +
+                mult.elements[7] * this.elements[1] +
+                mult.elements[11] * this.elements[2] +
+                this.elements[3];
+        elements[4] =
+            mult.elements[0] * this.elements[4] +
+                mult.elements[4] * this.elements[5] +
+                mult.elements[8] * this.elements[6];
+        elements[5] =
+            mult.elements[1] * this.elements[4] +
+                mult.elements[5] * this.elements[5] +
+                mult.elements[9] * this.elements[6];
+        elements[6] =
+            mult.elements[2] * this.elements[4] +
+                mult.elements[6] * this.elements[5] +
+                mult.elements[10] * this.elements[6];
+        elements[7] =
+            mult.elements[3] * this.elements[4] +
+                mult.elements[7] * this.elements[5] +
+                mult.elements[11] * this.elements[6] +
+                this.elements[7];
+        elements[8] =
+            mult.elements[0] * this.elements[8] +
+                mult.elements[4] * this.elements[9] +
+                mult.elements[8] * this.elements[10];
+        elements[9] =
+            mult.elements[1] * this.elements[8] +
+                mult.elements[5] * this.elements[9] +
+                mult.elements[9] * this.elements[10];
+        elements[10] =
+            mult.elements[2] * this.elements[8] +
+                mult.elements[6] * this.elements[9] +
+                mult.elements[10] * this.elements[10];
+        elements[11] =
+            mult.elements[3] * this.elements[8] +
+                mult.elements[7] * this.elements[9] +
+                mult.elements[11] * this.elements[10] +
+                this.elements[11];
+        return new Matrix3x4(elements);
+    }
+    inverse() {
+        return Matrix3x4.inverse(this);
+    }
+    isE() {
+        return Matrix3x4.isE(this);
+    }
+    equals(matrix) {
+        return Matrix3x4.equals(this, matrix);
+    }
+    clone() {
+        return new Matrix3x4(this.elements);
+    }
+    toString() {
+        return '(' + this.elements.slice(0, 12).map((element) => element.toFixed(3)) + ')';
+    }
+    /**
+     * Return transformation matrix as symetry code (common in Crystallography)
+     * e.g. 1-x,1-y,1-z
+     */
+    toSymetryCode() {
+        function formatSymNumber(num) {
+            return (Math.round(num * 1000) / 1000).toString();
+        }
+        function closeToOne(num) {
+            return Math.abs(num - 1) < 0.001;
+        }
+        function codePart1(x, y, z, l) {
+            const arr = [];
+            if (x !== 0) {
+                if (closeToOne(Math.abs(x))) {
+                    arr.push(x > 0 ? 'x' : '-x');
+                }
+                else {
+                    arr.push(formatSymNumber(x) + 'x');
+                }
+            }
+            if (y !== 0) {
+                if (closeToOne(Math.abs(y))) {
+                    arr.push(y > 0 ? 'y' : '-y');
+                }
+                else {
+                    arr.push(formatSymNumber(y) + 'y');
+                }
+            }
+            if (z !== 0) {
+                if (closeToOne(Math.abs(z))) {
+                    arr.push(z > 0 ? 'z' : '-z');
+                }
+                else {
+                    arr.push(formatSymNumber(z) + 'z');
+                }
+            }
+            if (l !== 0) {
+                arr.push(formatSymNumber(l));
+            }
+            return arr.reduce((curr, next) => {
+                if (next?.startsWith('-')) {
+                    return curr + next;
+                }
+                if (curr === '') {
+                    return next;
+                }
+                return curr + '+' + next;
+            }, '');
+        }
+        const sx = codePart1(this.elements[0], this.elements[1], this.elements[2], this.elements[3]);
+        const sy = codePart1(this.elements[4], this.elements[5], this.elements[6], this.elements[7]);
+        const sz = codePart1(this.elements[8], this.elements[9], this.elements[10], this.elements[11]);
+        return sx + ',' + sy + ',' + sz;
+    }
+}
+//# sourceMappingURL=matrix3x4.js.map