ipwr py and for skeleton

.gitignore

@@ -35,4 +35,5 @@ lang-sentences*
 /venv310/
 /cmake-build-debug-docker/
 /build-docs/
-.vscode
+/a.out.dSYM/
+a.out

.vscode/launch.json

@@ -0,0 +1,21 @@
+{
+    // Use IntelliSense to learn about possible attributes.
+    // Hover to view descriptions of existing attributes.
+    // For more information, visit: https://go.microsoft.com/fwlink/?linkid=830387
+    "version": "0.2.0",
+    "configurations": [
+      {
+        "name": "for",
+        "type": "cppdbg",
+        "request": "launch",
+        "preLaunchTask": "ipwr",
+        "program": "${workspaceFolder}/a.out",
+        "args": [],
+        "stopAtEntry": false,
+        "cwd": "${workspaceFolder}",
+        "environment": [],
+        "externalConsole": false,
+        "MIMode": "gdb"
+      },
+    ]
+}

.vscode/tasks.json

@@ -0,0 +1,10 @@
+{
+    "version": "2.0.0",
+    "tasks": [
+        {
+            "type": "shell",
+            "label": "ipwr",
+            "command": "gfortran pdp10/algebra/ipwr.for -g -O0 -fbacktrace -fcheck=all -ffpe-trap=zero,overflow,underflow -Wall -Wextra -Warray-temporaries -Wconversion -ffree-line-length-0",
+        }
+    ]
+}

pdp10/docs/algebra/ipwr.for

@@ -0,0 +1,11 @@
+
+
+      integer i
+      real pi
+      pi = 3.1415927
+      i = 0
+      i = i + 1
+10    format (f10.7)
+      write (*, 10) pi
+      end
+

pdp10/docs/algebra/ipwr.py

@@ -0,0 +1,21 @@
+""" imaginary powers of 10, starting with 10**(i/1024), and squaring successively ten times. this matches with 
+feynman's table 22-3. the rounded version gives the same x value in the last iteration as feynman's, -0.66928. note 
+the initial y = .00225 should be .0022486, just as in the text. the same rounding is used here as there."""
+from math import sqrt
+
+print('rounded version')
+y = .00225
+x = round(sqrt(1. - y**2), 7)
+for _ in range(11):
+    print('%10.5f %10.5f' % (x, y))
+    x2 = round(x**2 - y**2, 7)
+    y2 = round(2 * x * y, 7)
+    x, y = x2, y2    
+print('\nnon rounded version')
+y = .00225
+x = sqrt(1. - y**2)
+for _ in range(11):
+    print('%10.5f %10.5f' % (x, y))
+    x2 = x**2 - y**2
+    y2 = 2 * x * y
+    x, y = x2, y2

algebra/simplex.py → pdp10/docs/algebra/simplex.py

@@ -1,28 +1,25 @@
-"""simplest possible simplex implementation"""
-import numpy as np
-
-
-def main():
-    A = np.array([[-2, 1, 1, 0, 0], [-1, 2, 0, 1, 0], [1, 0, 0, 0, 1]])
-    b = np.array([2, 7, 3])
-    c = np.array([-1, -2, 0, 0, 0])
-    x = simplex(A, b, c)
-    return x
-
-
-def simplex(A, b, c):
-    """revised simplex method. """
-    m, n, o = A.shape[0], A.shape[1], A.shape[1] - A.shape[0]
-    ns, v1, v2 = [i for i in range(n)], np.array(c[o:]), np.array(c[:o])
-    Binv = np.linalg.inv(A[:, ns[o:]])
-    while not np.min(v2 - v1 @ Binv @ A[:, ns[:o]]) > 0:
-        n1 = np.argmin(v2 - v1 @ Binv @ A[:, ns[:o]])
-        t1, t2 = Binv @ b, Binv @ A[:, ns[n1]]
-        n2 = np.argmin([t1[i] / t2[i] if t2[i] > 0 else np.inf for i in range(m)])
-        ns[n1], ns[n2 + o], v1[n2], v2[n1] = ns[n2 + o], ns[n1], v2[n1], v1[n2]
-        Binv = np.linalg.inv(A[:, ns[o:]])
-    return Binv @ b
-
-
-if __name__ == "__main__":
-    main()
+"""simplest possible simplex implementation"""
+import numpy as np
+
+def main():
+    A = np.array([[-2, 1, 1, 0, 0], [-1, 2, 0, 1, 0], [1, 0, 0, 0, 1]])
+    b = np.array([2, 7, 3])
+    c = np.array([-1, -2, 0, 0, 0])
+    x = simplex(A, b, c)
+    return x
+
+def simplex(A, b, c):
+    """revised simplex method. """
+    m, n, o = A.shape[0], A.shape[1], A.shape[1] - A.shape[0]
+    ns, v1, v2 = [i for i in range(n)], np.array(c[o:]), np.array(c[:o])
+    Binv = np.linalg.inv(A[:, ns[o:]])
+    while not np.min(v2 - v1 @ Binv @ A[:, ns[:o]]) > 0:
+        n1 = np.argmin(v2 - v1 @ Binv @ A[:, ns[:o]])
+        t1, t2 = Binv @ b, Binv @ A[:, ns[n1]]
+        n2 = np.argmin([t1[i] / t2[i] if t2[i] > 0 else np.inf for i in range(m)])
+        ns[n1], ns[n2 + o], v1[n2], v2[n1] = ns[n2 + o], ns[n1], v2[n1], v1[n2]
+        Binv = np.linalg.inv(A[:, ns[o:]])
+    return Binv @ b
+
+if __name__ == "__main__":
+    main()

pdp10/docs/sec5-minimalist-walkthrough.md

@@ -1,47 +1,47 @@
-assume you've started simh pdp10 and it's listening on 2010 for connections.
-
-# on raspi side, start kermit
-
-    ./kermit/wermit
-    C-Kermit>set host localhost 2010 /raw-socket
-    C-Kermit>c
-
-# on tops10 side, start kermit
-
-assume you've created a subfolder [,,decwar] in your home folder [,].
-
-    .r setsrc
-    *cp [,,decwar]
-    *^C
-    .r kermit
-    Kermit-10>set file byte-size 36-bit
-    Kermit-10>server
-
-ctrl-\c back over to raspi kermit
-
-# in raspi kermit
-
-    C-Kermit>set transfer mode manual
-    C-Kermit>set file type binary
-    C-Kermit>send utexas/*.*
-    C-Kermit>send compuserve/*.*
-    C-Kermit>send scripts/COM1.CMD
-    C-Kermit>send scripts/COM2.CMD
-    C-Kermit>send scripts/CAN1.CMD
-
-can leave this connection in this state for further usage and get on tops10 via another connection, but if you like to use just this connection
-
-    C-Kermit>c
-    Kermit-10>^c^c^c
-    .
-
-# on tops10 side
-
-    .compile @com1
-    .compile @com2
-    .r link
-    *@can1
-    .get decwar
-    .ssave
-
-if it's not already there, create p,pn [1,27] using react, then copy DECWAR.EXE, DECWAR.GRP, DECWAR.HLP, DECWAR.NWS to its home folder to 'deploy/install' it. [1,27] is where decwar was deployed for compuserve. for utexas it was logical device name 'gam:', assigned to [5,30].
+assume you've started simh pdp10 and it's listening on 2010 for connections.
+
+# on raspi side, start kermit
+
+    ./kermit/wermit
+    C-Kermit>set host localhost 2010 /raw-socket
+    C-Kermit>c
+
+# on tops10 side, start kermit
+
+assume you've created a subfolder [,,decwar] in your home folder [,].
+
+    .r setsrc
+    *cp [,,decwar]
+    *^C
+    .r kermit
+    Kermit-10>set file byte-size 36-bit
+    Kermit-10>server
+
+ctrl-\c back over to raspi kermit
+
+# in raspi kermit
+
+    C-Kermit>set transfer mode manual
+    C-Kermit>set file type binary
+    C-Kermit>send utexas/*.*
+    C-Kermit>send compuserve/*.*
+    C-Kermit>send scripts/COM1.CMD
+    C-Kermit>send scripts/COM2.CMD
+    C-Kermit>send scripts/CAN1.CMD
+
+can leave this connection in this state for further usage and get on tops10 via another connection, but if you like to use just this connection
+
+    C-Kermit>c
+    Kermit-10>^c^c^c
+    .
+
+# on tops10 side
+
+    .compile @com1
+    .compile @com2
+    .r link
+    *@can1
+    .get decwar
+    .ssave
+
+if it's not already there, create p,pn [1,27] using react, then copy DECWAR.EXE, DECWAR.GRP, DECWAR.HLP, DECWAR.NWS to its home folder to 'deploy/install' it. [1,27] is where decwar was deployed for compuserve. for utexas it was logical device name 'gam:', assigned to [5,30].