From: Ian Jackson Date: Mon, 16 Apr 2018 20:23:52 +0000 (+0100) Subject: curveopt: in some comments, treat rad as a dimension (nfc) X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=ee4500220b0c3b5376b86c6cb77f6d206559ae0f;p=moebius3.git curveopt: in some comments, treat rad as a dimension (nfc) Signed-off-by: Ian Jackson --- diff --git a/symbolic.py b/symbolic.py index 8680378..21ae256 100644 --- a/symbolic.py +++ b/symbolic.py @@ -132,12 +132,12 @@ def calculate(): al, bl = ScalarArray('al bl', 'NP-1', a.magnitude() ).s() # [mm] au, bu = CoordArray ('au bu', 'NP-1', a / al ).s() # [1] - tan_theta = (au ^ bu) / (au & bu) # [1] bending + tan_theta = (au ^ bu) / (au & bu) # [rad] bending global mu, nu - mu, nu = CoordArray ('mu nu', 'NP-2', tan_theta ).s() # [1] + mu, nu = CoordArray ('mu nu', 'NP-2', tan_theta ).s() # [rad] - CostComponent('NP-3', sqnorm(mu - nu)) # [1] + CostComponent('NP-3', sqnorm(mu - nu)) # [rad^2] dl2 = pow(al - bl, 2) # [mm^2] CostComponent('NP-2', dl2 / (al*bl)) # [1]