From cc3d16cf2c0696a7aadbc6440f5b3b25155275bc Mon Sep 17 00:00:00 2001 From: numzero Date: Thu, 7 Nov 2024 00:37:30 +0300 Subject: [PATCH] oops, these were still 2D... --- src/riemann.rs | 2 +- src/tube/metric.rs | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/riemann.rs b/src/riemann.rs index dfeea72..6725986 100644 --- a/src/riemann.rs +++ b/src/riemann.rs @@ -58,7 +58,7 @@ pub fn krist(space: &impl Metric, pos: Vec3) -> Tens3 { let d = space.part_derivs_at(pos); // ret[i][l][k] = sum((m) => .5f * g[m][i] * (d[k][l][m] + d[l][k][m] - d[m][k][l])) make_tens3(|i, l, k| { - 0.5 * (0..2) + 0.5 * (0..3) .map(|m| g.col(m)[i] * (d[l].col(k)[m] + d[k].col(m)[l] - d[m].col(k)[l])) .sum::() }) diff --git a/src/tube/metric.rs b/src/tube/metric.rs index 74efe53..ba65a81 100644 --- a/src/tube/metric.rs +++ b/src/tube/metric.rs @@ -107,7 +107,7 @@ mod test { let pos = vec3(x, y, z); let computed = testee.part_derivs_at(pos); let reference = approx.part_derivs_at(pos); - let eq = (0..2).all(|coord| computed[coord].abs_diff_eq(reference[coord], epsilon)); + let eq = (0..3).all(|coord| computed[coord].abs_diff_eq(reference[coord], epsilon)); assert!( eq, "Bad derivative computation at {pos}:\n explicit: {computed:?}\n numerical: {reference:?}\n"