@woosh/meep-engine 2.131.3 → 2.131.6
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/package.json +1 -1
- package/src/core/collection/queue/Deque.d.ts.map +1 -1
- package/src/core/collection/queue/Deque.js +17 -12
- package/src/core/math/physics/mie/MIE_PARTICLES_STANDARD.d.ts +1 -0
- package/src/core/math/physics/mie/MIE_PARTICLES_STANDARD.d.ts.map +1 -1
- package/src/core/math/physics/mie/MIE_PARTICLES_STANDARD.js +1 -0
- package/src/core/math/physics/mie/MIE_PARTICLES_STANDARD_PRECOMPUTED.d.ts +39 -0
- package/src/core/math/physics/mie/MIE_PARTICLES_STANDARD_PRECOMPUTED.js +235 -214
- package/src/core/math/physics/mie/compute_lorenz_mie_optical_properties.d.ts.map +1 -1
- package/src/core/math/physics/mie/compute_lorenz_mie_optical_properties.js +3 -342
- package/src/core/math/physics/mie/compute_mie_particle_properties_rgb.js +1 -1
- package/src/core/math/physics/mie/compute_mie_phase.d.ts +16 -0
- package/src/core/math/physics/mie/compute_mie_phase.d.ts.map +1 -0
- package/src/core/math/physics/mie/compute_mie_phase.js +84 -0
- package/src/core/math/physics/mie/lorenz_mie_coefs.d.ts +28 -0
- package/src/core/math/physics/mie/lorenz_mie_coefs.d.ts.map +1 -0
- package/src/core/math/physics/mie/lorenz_mie_coefs.js +232 -0
- package/src/core/math/physics/mie/mie_ab_to_optical_properties.d.ts +14 -0
- package/src/core/math/physics/mie/mie_ab_to_optical_properties.d.ts.map +1 -0
- package/src/core/math/physics/mie/mie_ab_to_optical_properties.js +109 -0
- package/src/core/model/node-graph/node/NodeInstance.d.ts +2 -0
- package/src/core/model/node-graph/node/NodeInstance.d.ts.map +1 -1
- package/src/core/model/node-graph/node/NodeInstance.js +11 -5
- package/src/engine/animation/curve/AnimationCurve.d.ts +18 -3
- package/src/engine/animation/curve/AnimationCurve.d.ts.map +1 -1
- package/src/engine/animation/curve/AnimationCurve.js +69 -36
|
@@ -0,0 +1,232 @@
|
|
|
1
|
+
import { assert } from "../../../assert.js";
|
|
2
|
+
import { v2_length } from "../../../geom/vec2/v2_length.js";
|
|
3
|
+
import { complex_add } from "../../complex/complex_add.js";
|
|
4
|
+
import { complex_div } from "../../complex/complex_div.js";
|
|
5
|
+
import { complex_mul } from "../../complex/complex_mul.js";
|
|
6
|
+
import { complex_sub } from "../../complex/complex_sub.js";
|
|
7
|
+
|
|
8
|
+
export function vec2(x, y) {
|
|
9
|
+
return new Float64Array([x, y]);
|
|
10
|
+
}
|
|
11
|
+
|
|
12
|
+
function vec4(x, y, z, w) {
|
|
13
|
+
return new Float64Array([x, y, z, w]);
|
|
14
|
+
}
|
|
15
|
+
|
|
16
|
+
export function cadd(a, b) {
|
|
17
|
+
const r = vec2(0, 0);
|
|
18
|
+
|
|
19
|
+
complex_add(r, a, b);
|
|
20
|
+
|
|
21
|
+
return r;
|
|
22
|
+
}
|
|
23
|
+
|
|
24
|
+
function csub(a, b) {
|
|
25
|
+
const r = vec2(0, 0);
|
|
26
|
+
|
|
27
|
+
|
|
28
|
+
complex_sub(r, a, b);
|
|
29
|
+
|
|
30
|
+
return r;
|
|
31
|
+
}
|
|
32
|
+
|
|
33
|
+
export function cmult(a, b) {
|
|
34
|
+
const r = vec2(0, 0);
|
|
35
|
+
|
|
36
|
+
|
|
37
|
+
complex_mul(r, a, b);
|
|
38
|
+
|
|
39
|
+
return r;
|
|
40
|
+
}
|
|
41
|
+
|
|
42
|
+
export function cdiv(a, b) {
|
|
43
|
+
const r = vec2(0, 0);
|
|
44
|
+
|
|
45
|
+
|
|
46
|
+
complex_div(r, a, b);
|
|
47
|
+
|
|
48
|
+
return r;
|
|
49
|
+
}
|
|
50
|
+
|
|
51
|
+
// Number of terms to include in infinite sums (truncation)
|
|
52
|
+
function terms_to_sum(z) {
|
|
53
|
+
const size = v2_length(z[0], z[1]);
|
|
54
|
+
return Math.ceil(size + 4.3 * Math.cbrt(size) + 1.0);
|
|
55
|
+
}
|
|
56
|
+
|
|
57
|
+
// Downward recurrence for A
|
|
58
|
+
function A_all_n(z, M) {
|
|
59
|
+
const A = new Array(M + 2);
|
|
60
|
+
A[M + 1] = vec2(0.0, 0.0);
|
|
61
|
+
for (let n = M; n >= 0; --n) {
|
|
62
|
+
const tmp = cdiv(vec2(n + 1.0, 0.0), z);
|
|
63
|
+
A[n] = csub(tmp, cdiv(vec2(1.0, 0.0), cadd(tmp, A[n + 1])));
|
|
64
|
+
}
|
|
65
|
+
return A;
|
|
66
|
+
}
|
|
67
|
+
|
|
68
|
+
// Upward recurrences for B, psi_zeta, and R
|
|
69
|
+
function psi_zeta(n, z, A, old_B, old_psi_zeta) {
|
|
70
|
+
if (n > 0) {
|
|
71
|
+
const n_z = cdiv(vec2(n, 0.0), z);
|
|
72
|
+
const tmp = cmult(csub(n_z, A[n - 1]), csub(n_z, old_B));
|
|
73
|
+
old_psi_zeta = cmult(old_psi_zeta, tmp);
|
|
74
|
+
}
|
|
75
|
+
return old_psi_zeta;
|
|
76
|
+
}
|
|
77
|
+
|
|
78
|
+
/**
|
|
79
|
+
* Calculates the B_n(z) coefficient for term 'n' via upward recurrence.
|
|
80
|
+
* This is a helper function for LorenzMie_ab.
|
|
81
|
+
*
|
|
82
|
+
* @param {number} n - The current recurrence term index (n >= 0).
|
|
83
|
+
* @param {vec2} z - The complex size parameter (x or y).
|
|
84
|
+
* @param {Array<vec2>} A - Pre-computed log-derivatives A_n(z).
|
|
85
|
+
* @param {vec2} old_B - The previous B coefficient (B_{n-1}).
|
|
86
|
+
* @param {vec2} old_psi_zeta - The previous psi_zeta product term.
|
|
87
|
+
* @returns {vec2} The new B_n coefficient (complex number).
|
|
88
|
+
*/
|
|
89
|
+
function B(n, z, A, old_B, old_psi_zeta) {
|
|
90
|
+
|
|
91
|
+
if (n > 0) {
|
|
92
|
+
old_B.psi_zeta = psi_zeta(n, z, A, old_B, old_psi_zeta);
|
|
93
|
+
old_B = cadd(A[n], cdiv(vec2(0.0, 1.0), old_B.psi_zeta));
|
|
94
|
+
}
|
|
95
|
+
|
|
96
|
+
return old_B;
|
|
97
|
+
}
|
|
98
|
+
|
|
99
|
+
/**
|
|
100
|
+
* Calculates the R_n(z) coefficient (ratio psi_n / zeta_n) for term 'n'.
|
|
101
|
+
* This is a helper function for LorenzMie_ab.
|
|
102
|
+
*
|
|
103
|
+
* @param {number} n - The current recurrence term index (n >= 0).
|
|
104
|
+
* @param {vec2} z - The complex size parameter (x).
|
|
105
|
+
* @param {Array<vec2>} A - Pre-computed log-derivatives A_n(z).
|
|
106
|
+
* @param {vec2} B_n - The current B_n(z) coefficient.
|
|
107
|
+
* @param {vec2} old_R - The previous R coefficient (R_{n-1}).
|
|
108
|
+
* @returns {vec2} The new R_n coefficient (complex number).
|
|
109
|
+
*/
|
|
110
|
+
function R(n, z, A, B_n, old_R) {
|
|
111
|
+
if (n > 0) {
|
|
112
|
+
const n_z = cdiv(vec2(n, 0.0), z);
|
|
113
|
+
const tmp = cdiv(cadd(B_n, n_z), cadd(A[n], n_z));
|
|
114
|
+
old_R = cmult(old_R, tmp);
|
|
115
|
+
}
|
|
116
|
+
return old_R;
|
|
117
|
+
}
|
|
118
|
+
|
|
119
|
+
/**
|
|
120
|
+
* Calculates the Lorenz-Mie scattering coefficients a_n and b_n for a single term 'n'.
|
|
121
|
+
* This implements the formulas from Figure 3 of Frisvad et al. 2007.
|
|
122
|
+
*
|
|
123
|
+
* @param {number} n - The current recurrence term index (n >= 0).
|
|
124
|
+
* @param {number} size - Unitless size parameter x_vac = 2*pi*r / lambda_vac.
|
|
125
|
+
* @param {vec2} n_p - Complex refractive index of the particle.
|
|
126
|
+
* @param {vec2} n_med - Complex refractive index of the medium.
|
|
127
|
+
* @param {Array<vec2>} A_p - Pre-computed log-derivatives A_n(y) for the particle.
|
|
128
|
+
* @param {Array<vec2>} A_med - Pre-computed log-derivatives A_n(x) for the medium.
|
|
129
|
+
* @param {vec2} old_B - The previous B_{n-1}(x) coefficient.
|
|
130
|
+
* @param {vec2} old_R - The previous R_{n-1}(x) coefficient.
|
|
131
|
+
* @param {vec2} old_psi_zeta - The previous psi_zeta product.
|
|
132
|
+
* @returns {vec4} A vec4 [Re(a_n), Im(a_n), Re(b_n), Im(b_n)] with state properties
|
|
133
|
+
* (.old_B, .old_R, .old_psi_zeta) attached for the next iteration.
|
|
134
|
+
*/
|
|
135
|
+
function LorenzMie_ab(n, size, n_p, n_med, A_p, A_med, old_B, old_R, old_psi_zeta) {
|
|
136
|
+
|
|
137
|
+
const x = vec2(size * n_med[0], size * n_med[1]);
|
|
138
|
+
const B_n = B(n, x, A_med, old_B, old_psi_zeta);
|
|
139
|
+
const R_n = R(n, x, A_med, B_n, old_R);
|
|
140
|
+
const n_med_A_p = cmult(n_med, A_p[n]);
|
|
141
|
+
const n_p_A_med = cmult(n_p, A_med[n]);
|
|
142
|
+
const n_p_A_p = cmult(n_p, A_p[n]);
|
|
143
|
+
const n_med_A_med = cmult(n_med, A_med[n]);
|
|
144
|
+
const n_p_B_n = cmult(n_p, B_n);
|
|
145
|
+
const n_med_B_n = cmult(n_med, B_n);
|
|
146
|
+
const a = cmult(R_n, cdiv(csub(n_med_A_p, n_p_A_med), csub(n_med_A_p, n_p_B_n)));
|
|
147
|
+
const b = cmult(R_n, cdiv(csub(n_p_A_p, n_med_A_med), csub(n_p_A_p, n_med_B_n)));
|
|
148
|
+
|
|
149
|
+
|
|
150
|
+
const ab = vec4(a[0], a[1], b[0], b[1]);
|
|
151
|
+
|
|
152
|
+
ab.old_psi_zeta = n > 0 ? old_B.psi_zeta : old_psi_zeta;
|
|
153
|
+
ab.old_B = B_n;
|
|
154
|
+
ab.old_R = R_n;
|
|
155
|
+
|
|
156
|
+
return ab;
|
|
157
|
+
}
|
|
158
|
+
|
|
159
|
+
/**
|
|
160
|
+
* The following source code implements Lorenz-Mie theory using the formulas presented in the SIGGRAPH 2007 paper
|
|
161
|
+
*
|
|
162
|
+
* "Computing the Scattering Properties of Participating Media Using Lorenz-Mie Theory"
|
|
163
|
+
*
|
|
164
|
+
* By authors: Jeppe Revall Frisvad
|
|
165
|
+
* Niels Joergen Christensen
|
|
166
|
+
* Henrik Wann Jensen
|
|
167
|
+
*
|
|
168
|
+
* Code written by Jeppe Revall Frisvad, 2017.
|
|
169
|
+
* Copyright (c) Jeppe Revall Frisvad 2017
|
|
170
|
+
*
|
|
171
|
+
* Permission is granted to anyone to use this code as
|
|
172
|
+
* software for any purpose, including commercial applications.
|
|
173
|
+
* However, the software is provided 'as-is' without any warranty.
|
|
174
|
+
*
|
|
175
|
+
* @param {number} radius_m - Particle radius (in meters).
|
|
176
|
+
* @param {number} wavelength_m - Wavelength of light in vacuum (in meters).
|
|
177
|
+
* @param {vec2} n_p - Complex refractive index of the particle
|
|
178
|
+
* @param {vec2} n_med - Complex refractive index of the medium (e.g., [1.0, 0.0] for air).
|
|
179
|
+
* @returns {Float64Array} Mie coefficients a_n and b_n in form [a_n0.r, a_n0.i, b_n0.r, b_n0.i, ... a_nN.r, a_nN.i, b_nN.r, b_nN.i]
|
|
180
|
+
*/
|
|
181
|
+
export function lorenz_mie_coefs(
|
|
182
|
+
wavelength_m,
|
|
183
|
+
radius_m,
|
|
184
|
+
n_p,
|
|
185
|
+
n_med
|
|
186
|
+
) {
|
|
187
|
+
|
|
188
|
+
assert.isNumber(wavelength_m, 'wavelength_m');
|
|
189
|
+
assert.isFinite(wavelength_m, 'wavelength_m');
|
|
190
|
+
assert.greaterThan(wavelength_m, 0, 'wavelength_m');
|
|
191
|
+
|
|
192
|
+
assert.isNumber(radius_m, 'radius_m');
|
|
193
|
+
assert.isFinite(radius_m, 'radius_m');
|
|
194
|
+
assert.greaterThan(radius_m, 0, 'radius_m');
|
|
195
|
+
|
|
196
|
+
const psize = 2.0 * Math.PI * radius_m / wavelength_m;
|
|
197
|
+
|
|
198
|
+
const x = vec2(psize * n_med[0], psize * n_med[1]);
|
|
199
|
+
const y = vec2(psize * n_p[0], psize * n_p[1]);
|
|
200
|
+
const M = terms_to_sum(x);
|
|
201
|
+
const A_med = A_all_n(x, M);
|
|
202
|
+
const A_p = A_all_n(y, M);
|
|
203
|
+
|
|
204
|
+
const e_term = Math.exp(2.0 * x[1]);
|
|
205
|
+
|
|
206
|
+
let prev_psi_zeta = vec2(0.5 * (1.0 - Math.cos(2.0 * x[0]) / e_term), -0.5 * Math.sin(2.0 * x[0]) / e_term);
|
|
207
|
+
let prev_B = vec2(0.0, 1.0);
|
|
208
|
+
let prev_R = vec2(0.5 * (1.0 - Math.cos(2.0 * x[0]) * e_term), 0.5 * Math.sin(2.0 * x[0]) * e_term);
|
|
209
|
+
|
|
210
|
+
const ab = new Float64Array(M * 4);
|
|
211
|
+
|
|
212
|
+
for (let n = 0; n < M; ++n) {
|
|
213
|
+
|
|
214
|
+
const ab_n = LorenzMie_ab(
|
|
215
|
+
n, psize, n_p, n_med, A_p, A_med,
|
|
216
|
+
prev_B, prev_R, prev_psi_zeta
|
|
217
|
+
);
|
|
218
|
+
|
|
219
|
+
prev_psi_zeta = ab_n.old_psi_zeta;
|
|
220
|
+
prev_B = ab_n.old_B;
|
|
221
|
+
prev_R = ab_n.old_R;
|
|
222
|
+
|
|
223
|
+
// plug the data in
|
|
224
|
+
ab[4 * n] = ab_n[0];
|
|
225
|
+
ab[4 * n + 1] = ab_n[1];
|
|
226
|
+
ab[4 * n + 2] = ab_n[2];
|
|
227
|
+
ab[4 * n + 3] = ab_n[3];
|
|
228
|
+
|
|
229
|
+
}
|
|
230
|
+
|
|
231
|
+
return ab;
|
|
232
|
+
}
|
|
@@ -0,0 +1,14 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* Calculates the final optical properties by summing the Mie coefficients.
|
|
3
|
+
* This function implements Equations (2) and (3) from the 2007 paper
|
|
4
|
+
* "Computing the Scattering Properties of Participating Media..."
|
|
5
|
+
*
|
|
6
|
+
* @param {Float32Array} ab - The flat array of a_n, b_n coefficients from LorenzMie_coefs.
|
|
7
|
+
* @param {number} radius - Particle radius (in meters).
|
|
8
|
+
* @param {number} wavelength - Wavelength of light in vacuum (in meters).
|
|
9
|
+
* @param {vec2} n_med - Complex refractive index of the medium (e.g., [1.0, 0.0] for air).
|
|
10
|
+
* @returns {object} An object with Q_e, Q_s, albedo, C_ext, and C_sca.
|
|
11
|
+
*/
|
|
12
|
+
export function mie_ab_to_optical_properties(ab: Float32Array, wavelength: number, radius: number, n_med: typeof vec2): object;
|
|
13
|
+
import { vec2 } from "./lorenz_mie_coefs.js";
|
|
14
|
+
//# sourceMappingURL=mie_ab_to_optical_properties.d.ts.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"mie_ab_to_optical_properties.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/mie_ab_to_optical_properties.js"],"names":[],"mappings":"AAGA;;;;;;;;;;GAUG;AACH,iDANW,YAAY,cAEZ,MAAM,UADN,MAAM,uBAGJ,MAAM,CAgGlB;qBA3GuC,uBAAuB"}
|
|
@@ -0,0 +1,109 @@
|
|
|
1
|
+
import { v2_dot } from "../../../geom/vec2/v2_dot.js";
|
|
2
|
+
import { cadd, cdiv, cmult, vec2 } from "./lorenz_mie_coefs.js";
|
|
3
|
+
|
|
4
|
+
/**
|
|
5
|
+
* Calculates the final optical properties by summing the Mie coefficients.
|
|
6
|
+
* This function implements Equations (2) and (3) from the 2007 paper
|
|
7
|
+
* "Computing the Scattering Properties of Participating Media..."
|
|
8
|
+
*
|
|
9
|
+
* @param {Float32Array} ab - The flat array of a_n, b_n coefficients from LorenzMie_coefs.
|
|
10
|
+
* @param {number} radius - Particle radius (in meters).
|
|
11
|
+
* @param {number} wavelength - Wavelength of light in vacuum (in meters).
|
|
12
|
+
* @param {vec2} n_med - Complex refractive index of the medium (e.g., [1.0, 0.0] for air).
|
|
13
|
+
* @returns {object} An object with Q_e, Q_s, albedo, C_ext, and C_sca.
|
|
14
|
+
*/
|
|
15
|
+
export function mie_ab_to_optical_properties(
|
|
16
|
+
ab,
|
|
17
|
+
wavelength,
|
|
18
|
+
radius,
|
|
19
|
+
n_med
|
|
20
|
+
) {
|
|
21
|
+
const M = ab.length / 4;
|
|
22
|
+
|
|
23
|
+
let sum_t = vec2(0, 0); // Complex sum for Ct
|
|
24
|
+
let sum_s = 0; // Real sum for Cs
|
|
25
|
+
let sum_g = 0; // Real sum for Asymmetry (g)
|
|
26
|
+
|
|
27
|
+
for (let n = 1; n < M; ++n) {
|
|
28
|
+
// The 'n' in the paper's formulas
|
|
29
|
+
const multiplier = 2 * n + 1;
|
|
30
|
+
|
|
31
|
+
// Get an and bn from your flat array
|
|
32
|
+
const an = vec2(ab[4 * n], ab[4 * n + 1]);
|
|
33
|
+
const bn = vec2(ab[4 * n + 2], ab[4 * n + 3]);
|
|
34
|
+
|
|
35
|
+
// 1. For Extinction (Ct)
|
|
36
|
+
const an_plus_bn = cadd(an, bn);
|
|
37
|
+
// You need complex math functions cadd, cdiv
|
|
38
|
+
const n_med_sq = cmult(n_med, n_med);
|
|
39
|
+
const term_t = cdiv(an_plus_bn, n_med_sq);
|
|
40
|
+
sum_t = cadd(sum_t, cmult(vec2(multiplier, 0), term_t));
|
|
41
|
+
|
|
42
|
+
// 2. For Scattering (Cs)
|
|
43
|
+
const an_mag_sq = an[0] * an[0] + an[1] * an[1];
|
|
44
|
+
const bn_mag_sq = bn[0] * bn[0] + bn[1] * bn[1];
|
|
45
|
+
sum_s += multiplier * (an_mag_sq + bn_mag_sq);
|
|
46
|
+
|
|
47
|
+
// 3. Asymmetry Accumulation (New)
|
|
48
|
+
// Term A: Interference of an and bn
|
|
49
|
+
// Formula: (2n+1) / (n(n+1)) * Re(an * bn*)
|
|
50
|
+
const factor_A = (2 * n + 1) / (n * (n + 1));
|
|
51
|
+
sum_g += factor_A * v2_dot(an[0], an[1], bn[0], bn[1]);
|
|
52
|
+
|
|
53
|
+
// Term B: Interference of n and n+1
|
|
54
|
+
// Formula: (n(n+2)) / (n+1) * Re(an * an+1* + bn * bn+1*)
|
|
55
|
+
if (n < M - 1) {
|
|
56
|
+
const n_next = n + 1;
|
|
57
|
+
const an1_r = ab[4 * n_next];
|
|
58
|
+
const an1_i = ab[4 * n_next + 1];
|
|
59
|
+
const bn1_r = ab[4 * n_next + 2];
|
|
60
|
+
const bn1_i = ab[4 * n_next + 3];
|
|
61
|
+
|
|
62
|
+
const factor_B = (n * (n + 2)) / (n + 1);
|
|
63
|
+
|
|
64
|
+
const term_an_an1 = v2_dot(an[0], an[1], an1_r, an1_i);
|
|
65
|
+
const term_bn_bn1 = v2_dot(bn[0], bn[1], bn1_r, bn1_i);
|
|
66
|
+
|
|
67
|
+
sum_g += factor_B * (term_an_an1 + term_bn_bn1);
|
|
68
|
+
}
|
|
69
|
+
}
|
|
70
|
+
|
|
71
|
+
// --- Final Calculations ---
|
|
72
|
+
|
|
73
|
+
// Asymmetry Parameter (g)
|
|
74
|
+
// g = (4/x^2 * sum_g) / (2/x^2 * sum_s) => 2 * sum_g / sum_s
|
|
75
|
+
let g = 0;
|
|
76
|
+
if (sum_s > 1e-12) {
|
|
77
|
+
g = (2 * sum_g) / sum_s;
|
|
78
|
+
}
|
|
79
|
+
|
|
80
|
+
// Extinction (Ct from Eq. 22)
|
|
81
|
+
const C_t = (Math.pow(wavelength, 2) / (2 * Math.PI)) * sum_t[0]; //
|
|
82
|
+
|
|
83
|
+
// Scattering (Cs from Eq. 23)
|
|
84
|
+
const alpha = (4 * Math.PI * radius * n_med[1]) / wavelength;
|
|
85
|
+
let gamma;
|
|
86
|
+
if (Math.abs(alpha) < 1e-10) {
|
|
87
|
+
// avoid division by 0
|
|
88
|
+
gamma = 1;
|
|
89
|
+
} else {
|
|
90
|
+
gamma = (2 * (1 + (alpha - 1) * Math.exp(alpha))) / Math.pow(alpha, 2);
|
|
91
|
+
}
|
|
92
|
+
const n_med_mag_sq = n_med[0] * n_med[0] + n_med[1] * n_med[1];
|
|
93
|
+
const C_s_factor = (Math.pow(wavelength, 2) * Math.exp(-alpha)) / (2 * Math.PI * gamma * n_med_mag_sq);
|
|
94
|
+
const C_s = C_s_factor * sum_s;
|
|
95
|
+
|
|
96
|
+
// Bulk Coefficients (assuming number density N)
|
|
97
|
+
const sigma_a_med = (4 * Math.PI * n_med[1]) / wavelength;
|
|
98
|
+
|
|
99
|
+
// To get bulk coefficients sigma_t and sigma_s, use the following formulae, where N is numeric density (particles per cubic meter):
|
|
100
|
+
// sigma_t = sigma_a_med + C_t * N;
|
|
101
|
+
// sigma_s = C_s * N;
|
|
102
|
+
|
|
103
|
+
return {
|
|
104
|
+
sigma_a_med,
|
|
105
|
+
C_ext: C_t,
|
|
106
|
+
C_sca: C_s,
|
|
107
|
+
g, // Anisotropy parameter
|
|
108
|
+
}
|
|
109
|
+
}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"NodeInstance.d.ts","sourceRoot":"","sources":["../../../../../../src/core/model/node-graph/node/NodeInstance.js"],"names":[],"mappings":"AAcA;;;;GAIG;AACH;IAEI;;;;OAIG;IACH,aAFU,MAAM,CAEE;IAElB;;;OAGG;IACH,6BAAmB;IAEnB;;;OAGG;IACH,WAFU,yBAAyB,EAAE,CAEtB;IAEf;;;OAGG;IACH,gBAAgB;IAEhB;;;OAGG;IACH,qBAAmB;IAEnB;;;OAGG;IACH,sBAFU,gBAAgB,CAED;IAEzB;;;;;;OAMG;IACH,wBAFU,MAAM,CAEE;IAElB;;OAEG;IACH;QACI;;;WAGG;mCADO,OAAO,MAAM,WAAO;QAG9B;;;WAGG;iCADO,OAAO,MAAM,MAAI;QAG3B;;;WAGG;mCADO,OAAO,MAAM,MAAI;QAI3B;;;;WAIG;qCADO,yCAAyC,YAAY,CAAC;MAGlE;IAEF;;;;;;OAMG;IACH,8BALW,MAAM,aACN,aAAa,UACb,YAAY,GACV,MAAM,CAsBlB;IAED;;;OAGG;IACH,gDAEC;IAED;;;OAGG;IACH,+CAEC;IAGD;;;OAGG;IACH,mCAEC;IAED;;;;OAIG;IACH,kCAEC;IAED;;;;OAIG;IACH,mBAHW,MAAM,oBAKhB;IAED;;;;OAIG;IACH,mBAHW,MAAM,OAKhB;IAED;;;;OAIG;IACH,yBAHW,MAAM,KAOhB;IAED;;;;OAIG;IACH,2BAHW,MAAM,oBAmChB;IAED;;;;OAIG;IACH,iBAHW,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,oBAHW,MAAM,GACJ,OAAO,CAkBnB;IAED;;;OAGG;IACH,uCAQC;IAED,wBAIC;IAED;;;OAGG;IACH,mDAoCC;IAED;;;;OAIG;IACH,qBAHW,MAAM,GACJ,yBAAyB,GAAC,SAAS,CAoB/C;IAED;;;;OAIG;IACH,6BAHW,MAAM,GACJ,SAAS,GAAC,yBAAyB,
|
|
1
|
+
{"version":3,"file":"NodeInstance.d.ts","sourceRoot":"","sources":["../../../../../../src/core/model/node-graph/node/NodeInstance.js"],"names":[],"mappings":"AAcA;;;;GAIG;AACH;IAEI;;;;OAIG;IACH,aAFU,MAAM,CAEE;IAElB;;;OAGG;IACH,6BAAmB;IAEnB;;;OAGG;IACH,WAFU,yBAAyB,EAAE,CAEtB;IAEf;;;OAGG;IACH,gBAAgB;IAEhB;;;OAGG;IACH,qBAAmB;IAEnB;;;OAGG;IACH,sBAFU,gBAAgB,CAED;IAEzB;;;;;;OAMG;IACH,wBAFU,MAAM,CAEE;IAElB;;OAEG;IACH;QACI;;;WAGG;mCADO,OAAO,MAAM,WAAO;QAG9B;;;WAGG;iCADO,OAAO,MAAM,MAAI;QAG3B;;;WAGG;mCADO,OAAO,MAAM,MAAI;QAI3B;;;;WAIG;qCADO,yCAAyC,YAAY,CAAC;MAGlE;IAEF;;;;;;OAMG;IACH,8BALW,MAAM,aACN,aAAa,UACb,YAAY,GACV,MAAM,CAsBlB;IAED;;;OAGG;IACH,gDAEC;IAED;;;OAGG;IACH,+CAEC;IAGD;;;OAGG;IACH,mCAEC;IAED;;;;OAIG;IACH,kCAEC;IAED;;;;OAIG;IACH,mBAHW,MAAM,oBAKhB;IAED;;;;OAIG;IACH,mBAHW,MAAM,OAKhB;IAED;;;;OAIG;IACH,yBAHW,MAAM,KAOhB;IAED;;;;OAIG;IACH,2BAHW,MAAM,oBAmChB;IAED;;;;OAIG;IACH,iBAHW,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,oBAHW,MAAM,GACJ,OAAO,CAkBnB;IAED;;;OAGG;IACH,uCAQC;IAED,wBAIC;IAED;;;OAGG;IACH,mDAoCC;IAED;;;;OAIG;IACH,qBAHW,MAAM,GACJ,yBAAyB,GAAC,SAAS,CAoB/C;IAED;;;;OAIG;IACH,6BAHW,MAAM,GACJ,SAAS,GAAC,yBAAyB,CAmB/C;IAED;;;OAGG;IACH,QAFY,MAAM,CAIjB;IAED;;;;OAIG;IACH,cAHW,YAAY,GACV,OAAO,CAQnB;IAED,mBAYC;IAIL;;;OAGG;IACH,yBAFU,OAAO,CAEoB;CAPpC;0CA9YyC,gCAAgC;iBAHzD,kCAAkC;mBAChC,kCAAkC;8BAGvB,oBAAoB"}
|
|
@@ -351,13 +351,19 @@ export class NodeInstance {
|
|
|
351
351
|
getFirstEndpointByName(name) {
|
|
352
352
|
assert.isString(name, 'name');
|
|
353
353
|
|
|
354
|
-
const
|
|
354
|
+
const endpoints = this.endpoints;
|
|
355
|
+
const endpoint_count = endpoints.length;
|
|
356
|
+
|
|
357
|
+
for (let i = 0; i < endpoint_count; i++) {
|
|
358
|
+
const ref = endpoints[i];
|
|
355
359
|
|
|
356
|
-
|
|
357
|
-
|
|
360
|
+
if (ref.port.name === name) {
|
|
361
|
+
return ref;
|
|
362
|
+
}
|
|
358
363
|
}
|
|
359
364
|
|
|
360
|
-
|
|
365
|
+
// not found
|
|
366
|
+
return undefined;
|
|
361
367
|
|
|
362
368
|
}
|
|
363
369
|
|
|
@@ -387,7 +393,7 @@ export class NodeInstance {
|
|
|
387
393
|
|
|
388
394
|
// add description
|
|
389
395
|
const d = this.description;
|
|
390
|
-
if(d !== undefined && d !== null){
|
|
396
|
+
if (d !== undefined && d !== null) {
|
|
391
397
|
result += `, description = ${d.id}(${d.name})`
|
|
392
398
|
}
|
|
393
399
|
|
|
@@ -12,10 +12,17 @@
|
|
|
12
12
|
*
|
|
13
13
|
* jump_curve.evaluate(0.1); // what is the height at time 0.1?
|
|
14
14
|
*
|
|
15
|
+
* @example
|
|
16
|
+
* const curve = AnimationCurve.easeInOut();
|
|
17
|
+
*
|
|
18
|
+
* sprite.transparency = curve.evaluate(time); // smoothly animate transparency of the sprite
|
|
19
|
+
*
|
|
20
|
+
* @implements Iterable<Keyframe>
|
|
21
|
+
*
|
|
15
22
|
* @author Alex Goldring
|
|
16
23
|
* @copyright Company Named Limited (c) 2025
|
|
17
24
|
*/
|
|
18
|
-
export class AnimationCurve {
|
|
25
|
+
export class AnimationCurve implements Iterable<Keyframe> {
|
|
19
26
|
/**
|
|
20
27
|
* Utility constructor
|
|
21
28
|
* @param {Keyframe[]} keys
|
|
@@ -25,6 +32,7 @@ export class AnimationCurve {
|
|
|
25
32
|
/**
|
|
26
33
|
* S-shaped curve that starts slowly, ramps up and flattens out again.
|
|
27
34
|
* Useful for pleasing transitions where exit and entry should not be abrupt.
|
|
35
|
+
*
|
|
28
36
|
* @param {number} [timeStart]
|
|
29
37
|
* @param {number} [valueStart]
|
|
30
38
|
* @param {number} [timeEnd]
|
|
@@ -71,7 +79,7 @@ export class AnimationCurve {
|
|
|
71
79
|
/**
|
|
72
80
|
*
|
|
73
81
|
* @param {Keyframe} key
|
|
74
|
-
* @returns {boolean}
|
|
82
|
+
* @returns {boolean} true if the key was removed, false if the key was not found
|
|
75
83
|
*/
|
|
76
84
|
remove(key: Keyframe): boolean;
|
|
77
85
|
/**
|
|
@@ -108,10 +116,12 @@ export class AnimationCurve {
|
|
|
108
116
|
get duration(): number;
|
|
109
117
|
/**
|
|
110
118
|
* Returns index of a key that is just before or at the time specified.
|
|
119
|
+
* Useful for insertion and evaluation logic.
|
|
120
|
+
* Note: if time is past the end of last key - index of the last key will be returned instead
|
|
111
121
|
* @param {number} time
|
|
112
122
|
* @returns {number} index of the key
|
|
113
123
|
*/
|
|
114
|
-
|
|
124
|
+
getKeyIndexLow(time: number): number;
|
|
115
125
|
/**
|
|
116
126
|
* Evaluate interpolated value across the curve at a given time.
|
|
117
127
|
* @param {number} time time in seconds
|
|
@@ -164,6 +174,11 @@ export class AnimationCurve {
|
|
|
164
174
|
* @type {boolean}
|
|
165
175
|
*/
|
|
166
176
|
readonly isAnimationCurve: boolean;
|
|
177
|
+
/**
|
|
178
|
+
* @deprecated use `getKeyIndexLow` instead
|
|
179
|
+
*/
|
|
180
|
+
getKeyIndexByTime: (time: number) => number;
|
|
181
|
+
[Symbol.iterator](): Generator<Keyframe, void, unknown>;
|
|
167
182
|
}
|
|
168
183
|
import { Keyframe } from "./Keyframe.js";
|
|
169
184
|
//# sourceMappingURL=AnimationCurve.d.ts.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"AnimationCurve.d.ts","sourceRoot":"","sources":["../../../../../src/engine/animation/curve/AnimationCurve.js"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"AnimationCurve.d.ts","sourceRoot":"","sources":["../../../../../src/engine/animation/curve/AnimationCurve.js"],"names":[],"mappings":"AAUA;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,gDALwB,QAAQ;IAqZ5B;;;;OAIG;IACH,kBAHW,QAAQ,EAAE,GACR,cAAc,CAQ1B;IAED;;;;;;;;;OASG;IACH,6BANW,MAAM,eACN,MAAM,YACN,MAAM,aACN,MAAM,GACL,cAAc,CAUzB;IAED;;;;;;OAMG;IACH,4BALW,MAAM,YACN,MAAM,UACN,MAAM,GACL,cAAc,CAOzB;IAED;;;;;;;OAOG;IACH,0BANW,MAAM,eACN,MAAM,YACN,MAAM,aACN,MAAM,GACL,cAAc,CAezB;IAndD;;;;;OAKG;IACH,eAFU,QAAQ,EAAE,CAEV;IAEV;;;;;OAKG;IACH,SAHW,QAAQ,GACN,MAAM,CAuClB;IAED;;;OAGG;IACH,cAFW,QAAQ,EAAE,QAUpB;IAED;;;;OAIG;IACH,YAHW,QAAQ,GACN,OAAO,CAYnB;IAED;;OAEG;IACH,cAEC;IAGD;;;OAGG;IACH,WAFY,OAAO,CAIlB;IAED;;;OAGG;IACH,qBAEC;IAED;;;;OAIG;IACH,yBAUC;IAED;;;;OAIG;IACH,uBAUC;IAED;;;;OAIG;IACH,uBAcC;IAED;;;;;;OAMG;IACH,qBAHW,MAAM,GACJ,MAAM,CAmDlB;IAED;;;;OAIG;IACH,eAHW,MAAM,GACL,MAAM,CAoCjB;IAED;;;;OAIG;IACH,qBAFW,MAAM,QAmChB;IAED;;;;OAIG;IACH,sBAHW,MAAM,UACN,MAAM,QAWhB;IAED,0BAKC;IAGD;;;OAGG;IACH,YAFW,cAAc,QAIxB;IAED;;;OAGG;IACH,SAFY,cAAc,CAQzB;IAED;;;;OAIG;IACH,cAHW,cAAc,GACb,OAAO,CAIlB;IAED;;;OAGG;IACH,QAFY,MAAM,CAIjB;IAED;;MAIC;IAED;;aAkBC;IAoFL;;;;OAIG;IACH,2BAFU,OAAO,CAEwB;IAGzC;;OAEG;IACH,0BA3Te,MAAM,KACJ,MAAM,CA0TmB;IA7FtC,wDAQC;CAuEJ;yBA/ewB,eAAe"}
|