@woosh/meep-engine 2.110.13 → 2.110.15

package/src/core/geom/3d/sphere/harmonics/sh3_vector_accumulate.spec.js

@@ -1,27 +0,0 @@
-import { sh3_vector_accumulate } from "./sh3_vector_accumulate.js";
-
-test("basics", () => {
-    const out = [
-        1, 2,
-        3,4, 5,6,7,8,
-        9,10, 11,12, 13,14, 15,16, 17,18
-    ];
-    const basis = [
-        2,
-        3, 5, 7,
-        11, 13, 17, 19, 23
-    ];
-
-    sh3_vector_accumulate(
-        out, 0,
-        basis, 0,
-        [-3, 5], 0,
-        2
-    );
-
-    expect(out).toEqual([
-        -5, 12,
-        -6, 19, -10, 31, -14, 43,
-        -24, 65, -28, 77, -38, 99, -42, 111, -52, 133
-    ]);
-});

package/src/core/geom/3d/sphere/sphere_intersects_point.spec.js

@@ -1,134 +0,0 @@
-import { sphere_intersects_point } from "./sphere_intersects_point.js";
-
-test("zero volume sphere", () => {
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 0,
-        0, 0, 0
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 0,
-        -1, 0, 0
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 0,
-        1, 0, 0
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 0,
-        0, -1, 0
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 0,
-        0, 1, 0
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 0,
-        0, 0, -1
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 0,
-        0, 0, 1
-    )).toBe(false)
-
-});
-
-test("out of bounds, centered on origin", () => {
-
-    // X
-    expect(sphere_intersects_point(
-        0, 0, 0, 1,
-        -2, 0, 0
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 1,
-        2, 0, 0
-    )).toBe(false)
-
-
-    // Y
-    expect(sphere_intersects_point(
-        0, 0, 0, 1,
-        0, -2, 0
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 1,
-        0, 2, 0
-    )).toBe(false)
-
-    // Z
-    expect(sphere_intersects_point(
-        0, 0, 0, 1,
-        0, 0, -2
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 1,
-        0, 0, 2
-    )).toBe(false)
-
-});
-
-test("out of bounds with offset origin", () => {
-
-    expect(sphere_intersects_point(
-        3, -5, 7, 1,
-        1, -5, 7
-    )).toBe(false)
-
-    expect(sphere_intersects_point(
-        3, -5, 7, 1,
-        5, -5, 7
-    )).toBe(false)
-
-});
-
-test("in bounds", () => {
-
-    expect(sphere_intersects_point(
-        0, 0, 0, 2,
-        0, 0, 0
-    )).toBe(true)
-
-    // X
-    expect(sphere_intersects_point(
-        0, 0, 0, 2,
-        -1, 0, 0
-    )).toBe(true)
-
-    expect(sphere_intersects_point(
-        5, 5, 5, 2,
-        4, 5, 5
-    )).toBe(true)
-
-    // Y
-    expect(sphere_intersects_point(
-        0, 0, 0, 2,
-        0, -1, 0
-    )).toBe(true)
-
-    expect(sphere_intersects_point(
-        5, 5, 5, 2,
-        5, 4, 5
-    )).toBe(true)
-
-    // Z
-    expect(sphere_intersects_point(
-        0, 0, 0, 2,
-        0, 0, -1
-    )).toBe(true)
-
-    expect(sphere_intersects_point(
-        5, 5, 5, 2,
-        5, 5, 4
-    )).toBe(true)
-
-});

package/src/core/geom/3d/sphere/sphere_intersects_ray.spec.js

@@ -1,60 +0,0 @@
-import { sphere_intersects_ray } from "./sphere_intersects_ray.js";
-
-test("ray from -X to +X intersects sphere", () => {
-    expect(sphere_intersects_ray(
-        0, 0, 0, 1,
-        -2, 0, 0,
-        1, 0, 0
-    )).toBe(true);
-});
-
-test("ray from +X to -X intersects sphere", () => {
-    expect(sphere_intersects_ray(
-        0, 0, 0, 1,
-        2, 0, 0,
-        -1, 0, 0
-    )).toBe(true);
-});
-
-test("ray from -Y to +Y intersects sphere", () => {
-    expect(sphere_intersects_ray(
-        0, 0, 0, 1,
-        0, -2, 0,
-        0, 1, 0
-    )).toBe(true);
-});
-
-test("ray from +Y to -Y intersects sphere", () => {
-    expect(sphere_intersects_ray(
-        0, 0, 0, 1,
-        0, 2, 0,
-        0, -1, 0
-    )).toBe(true);
-});
-
-test("ray from -Z to +Z intersects sphere", () => {
-    expect(sphere_intersects_ray(
-        0, 0, 0, 1,
-        0, 0, -2,
-        0, 0, 1
-    )).toBe(true);
-});
-
-test("ray from +Z to -Z intersects sphere", () => {
-    expect(sphere_intersects_ray(
-        0, 0, 0, 1,
-        0, 0, 2,
-        0, 0, -1
-    )).toBe(true);
-});
-
-
-// negative cases
-
-test("ray from -X to -X", () => {
-    expect(sphere_intersects_ray(
-        0, 0, 0, 1,
-        -2, 0, 0,
-        -1, 0, 0
-    )).toBe(false);
-});

package/src/core/geom/3d/sphere/sphere_radius_sqr_from_v3_array_transformed.spec.js

@@ -1,8 +0,0 @@
-import { MATRIX_4_IDENTITY } from "../mat4/MATRIX_4_IDENTITY.js";
-import { sphere_radius_sqr_from_v3_array_transformed } from "./sphere_radius_sqr_from_v3_array_transformed.js";
-
-test("base case", () => {
-    expect(sphere_radius_sqr_from_v3_array_transformed([0, 0, 0], 3, MATRIX_4_IDENTITY)).toBe(0);
-    expect(sphere_radius_sqr_from_v3_array_transformed([1, 0, 0], 3, MATRIX_4_IDENTITY)).toBe(1);
-    expect(sphere_radius_sqr_from_v3_array_transformed([-1, 0, 0], 3, MATRIX_4_IDENTITY)).toBe(1);
-});

package/src/core/geom/3d/tetrahedra/TetrahedralMesh.spec.js

@@ -1,233 +0,0 @@
-import { jest } from "@jest/globals";
-import { TetrahedralMesh } from "./TetrahedralMesh.js";
-
-test("constructor does not throw", () => {
-    expect(() => new TetrahedralMesh(1)).not.toThrow();
-});
-
-test("capacity set in constructor is respected", () => {
-    const mesh = new TetrahedralMesh(7);
-
-    expect(mesh.getCapacity()).toBe(7);
-});
-
-test("attempting to set capacity below what's used will throw", () => {
-    const mesh = new TetrahedralMesh(1);
-
-    mesh.allocate();
-    mesh.allocate();
-
-    expect(() => mesh.setCapacity(1)).toThrow();
-});
-
-test("set/get capacity consistence", () => {
-    const mesh = new TetrahedralMesh(2);
-
-    mesh.setCapacity(1);
-
-    expect(mesh.getCapacity()).toBe(1);
-
-    mesh.setCapacity(1); // idempotency check
-
-    expect(mesh.getCapacity()).toBe(1);
-
-    mesh.setCapacity(2);
-
-    expect(mesh.getCapacity()).toBe(2);
-});
-
-
-test("set/get neighbour consistency", () => {
-    const mesh = new TetrahedralMesh(7);
-
-    const tet = mesh.allocate();
-
-    mesh.setNeighbour(tet, 0, 7);
-    mesh.setNeighbour(tet, 1, 3);
-    mesh.setNeighbour(tet, 2, 11);
-    mesh.setNeighbour(tet, 3, 13);
-
-    //
-    expect(mesh.getNeighbour(tet, 0)).toBe(7);
-    expect(mesh.getNeighbour(tet, 1)).toBe(3);
-    expect(mesh.getNeighbour(tet, 2)).toBe(11);
-    expect(mesh.getNeighbour(tet, 3)).toBe(13);
-});
-
-test("set/get vertex index consistency", () => {
-    const mesh = new TetrahedralMesh(7);
-    const tet = mesh.allocate();
-
-    mesh.setVertexIndex(tet, 0, 3);
-    mesh.setVertexIndex(tet, 1, 7);
-    mesh.setVertexIndex(tet, 2, 11);
-    mesh.setVertexIndex(tet, 3, 13);
-
-    //
-    expect(mesh.getVertexIndex(tet, 0)).toBe(3);
-    expect(mesh.getVertexIndex(tet, 1)).toBe(7);
-    expect(mesh.getVertexIndex(tet, 2)).toBe(11);
-    expect(mesh.getVertexIndex(tet, 3)).toBe(13);
-});
-
-test("allocation within capacity doesn't resize capacity", () => {
-    const mesh = new TetrahedralMesh(1);
-
-    expect(mesh.getCapacity()).toBe(1);
-
-    mesh.allocate();
-
-    expect(mesh.getCapacity()).toBe(1); // allocation fits in the capacity, so no growth is expected
-});
-
-test("allocation beyond capacity does not destroy existing data", () => {
-    const mesh = new TetrahedralMesh(1);
-
-    const t0 = mesh.allocate();
-
-    expect(mesh.getCapacity()).toBe(1);
-
-    for (let i = 0; i < 4; i++) {
-        mesh.setNeighbour(t0, i, 1);
-        mesh.setVertexIndex(t0, i, 3);
-    }
-
-    mesh.allocate(); // allocation over capacity
-
-    expect(mesh.getCapacity()).toBeGreaterThan(1);
-
-    for (let i = 0; i < 4; i++) {
-        expect(mesh.getNeighbour(t0, i)).toBe(1);
-        expect(mesh.getVertexIndex(t0, i)).toBe(3);
-    }
-
-});
-
-test("allocation increases size", () => {
-    const mesh = new TetrahedralMesh(1);
-
-    expect(mesh.size()).toBe(0);
-
-    mesh.allocate();
-
-    expect(mesh.size()).toBe(1);
-});
-
-test("remove method reduces size", () => {
-    const mesh = new TetrahedralMesh(1);
-
-    const tet = mesh.allocate();
-
-    expect(mesh.size()).toBe(1);
-
-    mesh.delete(tet);
-
-    expect(mesh.size()).toBe(0);
-});
-
-
-test("forEach visits all allocated tets", () => {
-    const mesh = new TetrahedralMesh();
-
-    const t0 = mesh.allocate();
-    const t1 = mesh.allocate();
-
-    const mock = jest.fn();
-
-    mesh.forEach(mock);
-
-    expect(mock).toHaveBeenCalledTimes(2);
-    expect(mock).toHaveBeenCalledWith(t0, mesh);
-    expect(mock).toHaveBeenCalledWith(t1, mesh);
-});
-
-test("'exists' method", () => {
-    const mesh = new TetrahedralMesh();
-
-    expect(mesh.exists(0)).toBe(false);
-
-    const tet = mesh.allocate();
-
-    expect(mesh.exists(tet)).toBe(true);
-    expect(mesh.exists(tet + 1)).toBe(false);
-
-    mesh.delete(tet);
-
-    expect(mesh.exists(tet)).toBe(false);
-    expect(mesh.exists(tet + 1)).toBe(false);
-});
-
-
-test("after calling 'removeTetrasConnectedToPoints' matching points are no logner indexed by live tetrahedrons", () => {
-
-
-    const tetrahedra = new TetrahedralMesh();
-
-    /*
-     * This "magic" mesh is a result of tetrahedral mesh generator with following input:  0, 0, 0, 10, 0, 0, 10, 0, 10, 0, 0, 10, 0, 10, 0,
-     * The mesh was generated on 27/10/2022
-     * It's a very small example with 21 total tets that was failing this test
-     *
-     */
-    tetrahedra.deserialize_base64('AAAAARUAAQAAAAAAAAAFAAAABgAAAAsAAAAUAAAAHwAAADQAAAADAAAABwAAAAUAAAAIAAAA/////yoAAABAAAAALgAAAAAAAAAFAAAABgAAAAgAAAD/////HAAAACgAAAAAAAAAAgAAAAUAAAAHAAAABgAAAP////8TAAAAFQAAACwAAAACAAAABgAAAAcAAAAIAAAA/////zwAAAAwAAAADQAAAAIAAAABAAAABQAAAAYAAAABAAAADgAAACAAAAA4AAAABAAAAAEAAAAGAAAACAAAAB0AAAAyAAAASQAAACEAAAABAAAAAAAAAAYAAAAIAAAACQAAABgAAABIAAAAAgAAAAQAAAACAAAABgAAAAEAAAAWAAAAGwAAACUAAAAzAAAABAAAAAMAAAACAAAAAQAAADsAAAAiAAAAUwAAAEYAAAADAAAABQAAAAAAAAAIAAAACgAAAEwAAAAFAAAANQAAAAMAAAAHAAAAAgAAAAUAAAAPAAAAOgAAAAcAAABEAAAABAAAAAYAAAACAAAACAAAABIAAAA+AAAAGQAAACMAAAADAAAAAQAAAAAAAAAFAAAAAwAAACsAAAA5AAAAUAAAAAMAAAACAAAAAQAAAAUAAAAXAAAANgAAAC0AAAAkAAAABAAAAAIAAAAHAAAACAAAABEAAABCAAAAMQAAAEUAAAAEAAAABwAAAAMAAAAIAAAABgAAAE4AAAA9AAAARwAAAAQAAAADAAAABwAAAAIAAAAvAAAAPwAAACcAAABDAAAABAAAAAAAAAABAAAACAAAAB4AAAAaAAAATQAAAFEAAAAEAAAAAwAAAAAAAAAIAAAAKQAAAEoAAABBAAAAUgAAAAQAAAADAAAAAQAAAAAAAAA3AAAASwAAAE8AAAAmAAAA');
-
-    tetrahedra.removeTetrasConnectedToPoints(5, 8);
-
-    //The mesh should contain only 5 vertices
-
-    tetrahedra.forEach((tet, mesh) => {
-
-        for (let i = 0; i < 4; i++) {
-            expect(mesh.getVertexIndex(tet, i)).toBeLessThan(5);
-        }
-
-    });
-});
-
-test('relocate operation copies tet data and updates back-links correctly', () => {
-    const mesh = new TetrahedralMesh();
-
-    const a = mesh.allocate();
-    const b = mesh.allocate();
-    const dest = mesh.allocate();
-
-    mesh.setVertexIndex(a, 0, 3);
-    mesh.setVertexIndex(a, 1, 7);
-    mesh.setVertexIndex(a, 2, 11);
-    mesh.setVertexIndex(a, 3, 13);
-
-    mesh.setNeighbour(a, 0, b << 2);
-    mesh.setNeighbour(b, 0, a << 2);
-
-    mesh.relocate(a, dest);
-
-    expect(mesh.getVertexIndex(dest, 0)).toBe(3);
-    expect(mesh.getVertexIndex(dest, 1)).toBe(7);
-    expect(mesh.getVertexIndex(dest, 2)).toBe(11);
-    expect(mesh.getVertexIndex(dest, 3)).toBe(13);
-
-    expect(mesh.getNeighbour(dest, 0)).toBe(b << 2);
-
-    // check back-link
-    expect(mesh.getNeighbour(b, 0)).toBe(dest << 2);
-});
-
-test('Compaction is non-destructive when moving resident data from the very end of the allocated region', () => {
-    const mesh = new TetrahedralMesh();
-
-    const a = mesh.allocate();
-    const b = mesh.allocate();
-
-    // force move "a" to free set
-    mesh.delete(a);
-
-    expect(mesh.exists(b)).toBe(true);
-    expect(mesh.exists(a)).toBe(false);
-
-    mesh.compact();
-
-    expect(mesh.size()).toBe(1);
-
-    expect(mesh.exists(0)).toBe(true);
-    expect(mesh.exists(1)).toBe(false);
-});

package/src/core/geom/3d/tetrahedra/compute_bounding_simplex_3d.spec.js

@@ -1,66 +0,0 @@
-import { compute_bounding_simplex_3d } from "./compute_bounding_simplex_3d.js";
-import { v3_distance } from "../../vec3/v3_distance.js";
-import { orient3d_fast } from "../plane/orient3d_fast.js";
-
-test('one point', () => {
-    const r = [];
-    compute_bounding_simplex_3d(r, 0, [1, 2, 3], 1, 0);
-
-    for (let i = 0; i < 4; i++) {
-        const i3 = i * 3;
-
-        expect(r[i3]).toBeCloseTo(1);
-        expect(r[i3 + 1]).toBeCloseTo(2);
-        expect(r[i3 + 2]).toBeCloseTo(3);
-    }
-});
-
-test('planar containment, one point', () => {
-    const r = [];
-
-    compute_bounding_simplex_3d(r, 0, [0, 0, 0], 1, 1);
-
-    r.push(0, 0, 0);
-
-    expect(orient3d_fast(r, 2, 3, 1, 4)).toBeLessThan(0);
-    expect(orient3d_fast(r, 0, 3, 2, 4)).toBeLessThan(0);
-    expect(orient3d_fast(r, 0, 1, 3, 4)).toBeLessThan(0);
-    expect(orient3d_fast(r, 2, 1, 0, 4)).toBeLessThan(0);
-});
-
-test('positive padding results in larger simplex', () => {
-    const r_0 = [];
-    const r_n = [];
-    const r_p = [];
-
-    const volume = [
-        -100, 0, 0,
-        100, 0, 0,
-        0, -100, 0,
-        0, 100, 0,
-        0, 0, -100,
-        0, 0, 100,
-    ];
-
-    compute_bounding_simplex_3d(r_0, 0, volume, 6, 0);
-    compute_bounding_simplex_3d(r_n, 0, volume, 6, -1);
-    compute_bounding_simplex_3d(r_p, 0, volume, 6, 1);
-
-    function compute_distance_sum(ref, points, point_count) {
-        let r = 0;
-        for (let i = 0; i < point_count; i++) {
-            const i3 = i * 3;
-
-            r += v3_distance(points[i3], points[i3 + 1], points[i3 + 2], ref[0], ref[1], ref[2]);
-        }
-
-        return r;
-    }
-
-    const d_0 = compute_distance_sum([0, 0, 0], r_0, 4);
-    const d_n = compute_distance_sum([0, 0, 0], r_n, 4);
-    const d_p = compute_distance_sum([0, 0, 0], r_p, 4);
-
-    expect(d_0).toBeLessThan(d_p);
-    expect(d_0).toBeGreaterThan(d_n);
-});

package/src/core/geom/3d/tetrahedra/delaunay/compute_delaunay_tetrahedral_mesh.spec.js

@@ -1,36 +0,0 @@
-/**
- * @jest-environment jsdom
- */
-
-import { TetrahedralMesh } from "../TetrahedralMesh.js";
-import { validate_tetrahedral_mesh } from "../validate_tetrahedral_mesh.js";
-import { compute_delaunay_tetrahedral_mesh } from "./compute_delaunay_tetrahedral_mesh.js";
-
-test('single point', () => {
-    const mesh = new TetrahedralMesh();
-
-    compute_delaunay_tetrahedral_mesh(mesh, [1, 2, 3], 1);
-});
-
-test('exactly 4 points', () => {
-    const mesh = new TetrahedralMesh();
-    const points = [
-        -1, 0, 0,
-        1, 0, 0,
-        0.5, 1, 0,
-        0.5, 0, 1
-    ];
-    compute_delaunay_tetrahedral_mesh(mesh, points, 4);
-
-
-    // do full validation
-    const problems = [];
-
-    const is_valid = validate_tetrahedral_mesh(mesh, points, (problem) => problems.push(problem));
-
-    if (!is_valid) {
-        console.error(problems.join('\n'));
-    }
-
-    expect(is_valid).toBe(true);
-});

package/src/core/geom/3d/tetrahedra/delaunay/tetrahedral_mesh_build_from_grid.spec.js

@@ -1,61 +0,0 @@
-import { AABB3 } from "../../aabb/AABB3.js";
-import { compute_tetrahedron_volume } from "../compute_tetrahedron_volume.js";
-import { TetrahedralMesh } from "../TetrahedralMesh.js";
-import { validate_tetrahedral_mesh } from "../validate_tetrahedral_mesh.js";
-import { tetrahedral_mesh_build_from_grid } from "./tetrahedral_mesh_build_from_grid.js";
-
-test("Single cube", () => {
-
-    const mesh = new TetrahedralMesh();
-
-    /**
-     *
-     * @type {number[]}
-     */
-    const points = [];
-
-    tetrahedral_mesh_build_from_grid(
-        mesh,
-        points,
-        new AABB3(0, 0, 0, 1, 1, 1),
-        2, 2, 2
-    );
-
-    // number of points
-    expect(points.length).toBe(8 * 3);
-
-    expect(mesh.count).toBeGreaterThanOrEqual(5);
-    expect(mesh.count).toBeLessThanOrEqual(6);
-
-    // mesh should only reference created points
-    const tetrahedrons = mesh.getLive();
-
-    for (let i = 0; i < tetrahedrons.length; i++) {
-        const tet_index = tetrahedrons[i];
-
-        for (let vertex_index = 0; vertex_index < 4; vertex_index++) {
-
-            const vert = mesh.getVertexIndex(tet_index, vertex_index);
-
-            expect(vert).toBeGreaterThanOrEqual(0);
-            expect(vert).toBeLessThan(8);
-        }
-
-        // tetrahedrons should have positive volume
-        const volume = compute_tetrahedron_volume(mesh, points, tet_index);
-
-        expect(volume).toBeGreaterThan(0);
-    }
-
-    // do full validation
-    const problems = [];
-
-    const is_valid = validate_tetrahedral_mesh(mesh, points, (problem) => problems.push(problem));
-
-    if (!is_valid) {
-        console.error(problems.join('\n'));
-    }
-
-    expect(is_valid).toBe(true);
-
-});

package/src/core/geom/3d/tetrahedra/delaunay/tetrahedral_mesh_compute_sub_determinant.spec.js

@@ -1,30 +0,0 @@
-import { TetrahedralMesh } from "../TetrahedralMesh.js";
-import { tetrahedral_mesh_compute_sub_determinant } from "./tetrahedral_mesh_compute_sub_determinant.js";
-
-test("sanity test", () => {
-    const mesh = new TetrahedralMesh(1);
-
-    const tet = mesh.allocate();
-
-    mesh.setVertexIndex(tet, 0, 0);
-    mesh.setVertexIndex(tet, 1, 1);
-    mesh.setVertexIndex(tet, 2, 2);
-    mesh.setVertexIndex(tet, 3, 3);
-
-    const result = [];
-
-    const points = [
-        1, 2, 3,
-        4, 5, 6,
-        7, 8, 9,
-        10, 11, 12
-    ];
-
-    tetrahedral_mesh_compute_sub_determinant(result, 0, mesh, points, tet);
-
-    for (let i = 0; i < 4; i++) {
-        expect(typeof result[i]).toBe("number");
-        expect(result[i]).not.toBeNaN();
-    }
-
-});

package/src/core/geom/3d/tetrahedra/tetrahedron_compute_signed_volume.spec.js

@@ -1,24 +0,0 @@
-import { tetrahedron_compute_signed_volume } from "./tetrahedron_compute_signed_volume.js";
-
-test("0 volume", () => {
-    expect(tetrahedron_compute_signed_volume([
-        0, 0, 0
-    ], 0, 0, 0, 0)).toBe(0);
-});
-
-test("positive volume", () => {
-
-    expect(tetrahedron_compute_signed_volume([
-        1, 0, 1,
-        1, 0, 0,
-        0, 0, 0,
-        1, 1, 1
-    ], 1, 2, 3, 0)).toBeGreaterThan(0);
-
-    expect(tetrahedron_compute_signed_volume([
-        0, 0, 0,
-        2, 0, 0,
-        0, 2, 0,
-        0, 0, 2
-    ], 1, 2, 3, 0)).toBeGreaterThan(0);
-})