chiark / gitweb /
curveopt: in some comments, treat rad as a dimension (nfc)
authorIan Jackson <ijackson@chiark.greenend.org.uk>
Mon, 16 Apr 2018 20:23:52 +0000 (21:23 +0100)
committerIan Jackson <ijackson@chiark.greenend.org.uk>
Mon, 16 Apr 2018 20:23:52 +0000 (21:23 +0100)
Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
symbolic.py

index 8680378ea06d906cbf58b4d230f09451e5ec281d..21ae25691c990b345000acb573fdce7a93786ec2 100644 (file)
@@ -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]