From 01f44324c12108fd04ff33577a32eb3d770dc4a5 Mon Sep 17 00:00:00 2001 From: Aaron Fenyes Date: Wed, 17 Jul 2024 22:45:17 -0700 Subject: [PATCH] Tetrahedron radius ratio: find radius ratio --- engine-proto/gram-test/tetrahedron-radius-ratio.jl | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/engine-proto/gram-test/tetrahedron-radius-ratio.jl b/engine-proto/gram-test/tetrahedron-radius-ratio.jl index 9e79c05..4218cb7 100644 --- a/engine-proto/gram-test/tetrahedron-radius-ratio.jl +++ b/engine-proto/gram-test/tetrahedron-radius-ratio.jl @@ -1,5 +1,6 @@ include("Engine.jl") +using LinearAlgebra using SparseArrays using Random @@ -87,4 +88,9 @@ else println("\nFailed to reach target accuracy") end println("Steps: ", size(history.scaled_loss, 1)) -println("Loss: ", history.scaled_loss[end], "\n") \ No newline at end of file +println("Loss: ", history.scaled_loss[end]) +if success + infty = BigFloat[0, 0, 0, 0, 1] + radius_ratio = dot(infty, Engine.Q * L[:,5]) / dot(infty, Engine.Q * L[:,6]) + println("\nCircumradius / inradius: ", radius_ratio) +end \ No newline at end of file