From 3c41733a4cb46d74d3391e6dec33496cc2a7d0d7 Mon Sep 17 00:00:00 2001
From: Stephan Hilb <stephan@ecshi.net>
Date: Thu, 17 Feb 2022 08:54:02 +0100
Subject: [PATCH] add norm calculations
---
scripts/symbolic_calculations.jl | 18 ++++++++++++++++++
1 file changed, 18 insertions(+)
diff --git a/scripts/symbolic_calculations.jl b/scripts/symbolic_calculations.jl
index 2b98fc4..a4cc08f 100644
--- a/scripts/symbolic_calculations.jl
+++ b/scripts/symbolic_calculations.jl
@@ -17,3 +17,21 @@ function l1triangle()
return result
end
+
+function norm_grad()
+ function ref_int(f)
+ inty_f = Algebra.int(f, :y, 0, :(1 - x))
+ return Algebra.int(inty_f, :x, 0, 1)
+ end
+
+ # \|f\|
+ f = :(a * x + b * y + c * (1 - x - y))
+ n2_f = ref_int(:(($f)^2))
+
+ # \|grad f\|
+ dx_f = Algebra.df(f, :x)
+ dy_f = Algebra.df(f, :y)
+ n2_gradf = ref_int(:(($dx_f)^2 + ($dy_f)^2))
+
+ return :($n2_gradf / $n2_f)
+end
--
GitLab