From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Wed, 8 Jul 2015 20:46:52 -0400
Subject: Initial commit

---
 sys/mwcs/mwmkidmr.x | 18 ++++++++++++++++++
 1 file changed, 18 insertions(+)
 create mode 100644 sys/mwcs/mwmkidmr.x

(limited to 'sys/mwcs/mwmkidmr.x')

diff --git a/sys/mwcs/mwmkidmr.x b/sys/mwcs/mwmkidmr.x
new file mode 100644
index 00000000..f4771217
--- /dev/null
+++ b/sys/mwcs/mwmkidmr.x
@@ -0,0 +1,18 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+# MW_MKIDMR -- Make the identity matrix.
+
+procedure mw_mkidmr (ltm, ndim)
+
+real	ltm[ndim,ndim]		#O set to the identity matrix
+int	ndim			#I dimension of (square) matrix
+
+int	i, j
+
+begin
+	do j = 1, ndim {
+	    do i = 1, ndim
+		ltm[i,j] = 0.0
+	    ltm[j,j] = 1.0
+	}
+end
-- 
cgit