diff options
Diffstat (limited to 'sys/mwcs/mwmkidmr.x')
| -rw-r--r-- | sys/mwcs/mwmkidmr.x | 18 |
1 files changed, 18 insertions, 0 deletions
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 |