noteuclid/refraction
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Remove useless code

Browse Source
This commit is contained in:
numzero 2024-04-28 23:15:19 +03:00
parent cbed57bd3a
commit 641b13149f
1 changed files with 0 additions and 28 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
src/bin/flat.rs
View File

@ -130,10 +130,6 @@ mod riemann {
} }
} }
fn trace(space: &impl Metric, base: Vec2, dir: Vec2, distance: f32, dt: f32) -> Vec<Vec2> {
trace_iter(space, base, dir, dt).take((distance / dt) as usize).collect()
}
pub fn trace_iter<M: Metric>(space: &M, base: Vec2, dir: Vec2, dt: f32) -> TraceIter<M> { pub fn trace_iter<M: Metric>(space: &M, base: Vec2, dir: Vec2, dt: f32) -> TraceIter<M> {
TraceIter { TraceIter {
space, space,
@ -143,30 +139,6 @@ mod riemann {
} }
} }
#[cfg(test)]
mod test {
use glm::*;
use crate::riemann::{trace, trace_iter};
use crate::Coil;
#[test]
fn t_iter() {
let space = Coil {
coil_scale: 2.0,
coil_r: 300.0,
coil_w: 50.0,
coil_m: 10.0,
};
let base = vec2(-500.0, 0.0);
let dir = vec2(1.0, 0.3);
let dt = 1.0;
let steps = 1000;
let a = trace(&space, base, dir, dt * (steps as f32), dt);
let b: Vec<Vec2> = trace_iter(&space, base, dir, dt).take(steps).collect();
assert_eq!(a, b);
}
}
fn krist(space: &impl Metric, pos: Vec2) -> Tens2 { fn krist(space: &impl Metric, pos: Vec2) -> Tens2 {
// Γ^i_k_l = .5 * g^i^m * (g_m_k,l + g_m_l,k - g_k_l,m) // Γ^i_k_l = .5 * g^i^m * (g_m_k,l + g_m_l,k - g_k_l,m)
let g = inverse(&space.metric(pos)); // с верхними индексами let g = inverse(&space.metric(pos)); // с верхними индексами