From 3795ff9d16ebf4d30ee7e33c687b6c26a67c025e Mon Sep 17 00:00:00 2001
From: Bryan Newbold <bnewbold@robocracy.org>
Date: Thu, 15 Sep 2022 17:00:28 -0700
Subject: 2022-09-15: manifolds chapter, started ch3

---
 chapters/2-Manifolds.scm | 190 +++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 190 insertions(+)
 create mode 100644 chapters/2-Manifolds.scm

(limited to 'chapters')

diff --git a/chapters/2-Manifolds.scm b/chapters/2-Manifolds.scm
new file mode 100644
index 0000000..e929510
--- /dev/null
+++ b/chapters/2-Manifolds.scm
@@ -0,0 +1,190 @@
+
+(define R2 (make-manifold R^n 2))
+
+(define U (patch 'origin R2))
+
+; these are coordinate systems
+(define R2-rect (coordinate-system 'rectangular U))
+(define R2-polar (coordinate-system 'polar/cylindrical U))
+
+; these are charts, and their inverses
+(define R2-rect-chi (chart R2-rect))
+(define R2-rect-chi-inverse (point R2-rect))
+(define R2-polar-chi (chart R2-polar))
+(define R2-polar-chi-inverse (point R2-polar))
+
+((compose R2-polar-chi R2-rect-chi-inverse) (up 'x0 'y0))
+; (up 
+;   (sqrt (+ (expt x0 2) (expt y0 2)))  ; radius
+;   (atan y0 x0))                       ; angle (theta)
+
+((compose R2-rect-chi R2-polar-chi-inverse) (up 'r0 'theta0))
+; (up
+;   (* r0 (cos theta0))     ; x
+;   (* r0 (sin theta0)))    ; y
+
+
+((D (compose R2-rect-chi R2-polar-chi-inverse)) (up 'r0 'theta0))
+; (down
+;   (up
+;       (cos theta0)
+;       (sin theta0))
+;   (up
+;       (* -1 r0 (sin theta0))
+;       (* r0 (cos theta0))))
+
+(define R2->R (-> (UP Real Real) Real))
+(define f (compose (literal-function 'f-rect R2->R) R2-rect-chi))
+
+(define R2-rect-point (R2-rect-chi-inverse (up 'x0 'y0)))
+
+(define corresponding-polar-point
+  (R2-polar-chi-inverse
+    (up (sqrt (+ (square 'x0) (square 'y0)))
+        (atan 'y0 'x0))))
+
+(f R2-rect-point)
+; (f-rect (up x0 y0))
+
+(f corresponding-polar-point)
+; (f-rect (up x0 y0))
+
+; confirms that the CAS simplifies to the same point
+
+(define-coordinates (up x y) R2-rect)
+(define-coordinates (up r theta) R2-polar)
+
+(x (R2-rect-chi-inverse (up 'x0 'y0)))
+; x0
+
+; expect r0 * cos(theta0)
+(x (R2-polar-chi-inverse (up 'r0 'theta0)))
+; (* r0 (cos theta0))
+
+; expect r0 * sin(theta0)
+(y (R2-polar-chi-inverse (up 'r0 'theta0)))
+; (* r0 (sin theta0))
+
+
+;(r (R2-polar-chi-inverse (up 'r0 'theta0)))
+; r0
+
+; expect sqrt(x0^2 + y0^2)
+(r (R2-rect-chi-inverse (up 'x0 'y0)))
+; (sqrt (+ (expt x0 2) (expt y0 2)))
+
+; expect atan(y0, x0)
+(theta (R2-rect-chi-inverse (up 'x0 'y0)))
+; (atan y0 x0)
+
+; h: x * r^2 + y^3
+(define h (+ (* x (square r)) (cube y)))
+
+(h R2-rect-point)
+; (+ (expt x0 3) (* x0 (expt y0 2)) (expt y0 3))
+; aka: x0^3 + x0 y0^2 + y0^3 
+;      x0 (x0^2 + y0^2) + y0^3
+
+(h (R2-polar-chi-inverse (up 'r0 'theta0)))
+; (+ (* (expt r0 3) (expt (sin theta0) 3)) (* (expt r0 3) (cos theta0)))
+
+
+;;; Exersize 2.1a
+
+((- r (* 2 'a (+ 1 (cos theta))))
+ ((point R2-rect) (up 'x 'y)))
+; (/ 
+;   (+ (* -2 a x)
+;      (* -2 a (sqrt (+ (expt x 2) (expt y 2))))
+;      (expt x 2) (expt y 2))
+;   (sqrt (+ (expt x 2) (expt y 2))))
+
+
+; "Lemniscate of Bernoulli"
+; (x^2 + y^2)^2 = 2 a^2 (x^2 - y^2)
+
+((-
+   (square (+ (square x) (square y)))
+   (* 2
+      (square 'a)
+      (- (square x) (square y))))
+ ((point R2-polar) (up 'r0 'theta0)))
+; (+ (* -4 (expt a 2) (expt r0 2) (expt (cos theta0) 2))
+;    (* 2 (expt a 2) (expt r0 2))
+;    (expt r0 4))
+;
+; r^2 / a^2 + 2 = 4 cos^2 (theta)
+; => r^2 = 2 a^2 cos(2 theta)
+
+; this matches Wikipedia, with substitution of a^2 = 2 c^2 (in the rectangular version)
+
+
+;;; Exersize 2.1b
+
+(define R3-rect-chi (chart R3-rect))
+(define R3-rect-chi-inverse (point R3-rect))
+(define R3-cyl-chi (chart R3-cyl))
+(define R3-cyl-chi-inverse (point R3-cyl))
+
+(define-coordinates (up x y z) R3-rect)
+(define-coordinates (up r theta z-cyl) R3-cyl)
+
+; R->R3
+(define (helix-rect t)
+  (up (* 'R (cos t))
+      (* 'R (sin t))
+      (* 's t)))
+
+; R->R3
+(define (helix-cyl t)
+  (up 'R
+      t
+      (* 's t)))
+
+((compose R3-rect-chi R3-cyl-chi-inverse helix-cyl) 't)
+; (up
+;   (* R (cos t))
+;   (* R (sin t))
+;   (* s t))
+
+
+((compose R3-cyl-chi R3-rect-chi-inverse helix-rect) 't)
+; (up R
+;     t
+;     (* s t))
+
+((- helix-rect (compose R3-rect-chi R3-cyl-chi-inverse helix-cyl)) 't)
+; (up 0 0 0)
+
+;;; Exersize 2.2
+
+; given polar coordinates of point on plane, get to spherical coordinates of point on sphere with:
+((compose
+   (chart S2-spherical)
+   (point S2-Riemann)
+   (chart R2-rect)
+   (point R2-polar))
+ (up 'rho 'theta))
+; (up (acos (/ (+ -1 (expt rho 2))
+;              (+ 1 (expt rho 2))))
+;     theta)
+
+; given spherical coordinates of point on sphere, what are the polar
+; coordinates of corresponding point on the plane?
+((compose
+   (chart R2-polar)
+   (point R2-rect)
+   (chart S2-Riemann)
+   (point S2-spherical))
+ (up 'phi 'lambda))
+; (up (/ (sin phi)
+;        (+ -1 (cos phi)))
+;     (atan (* -1 (sin lambda))
+;           (* -1 (cos lambda))))
+
+; this is just the inverse of the above, right?
+
+; huh, I would expect the second term to just be lambda.
+; atan (-sin(lambda), -cos(lambda)) -> tan^-1(tan(lambda)) -> lambda
+; seems like an identity, ok
+
-- 
cgit v1.2.3