diff options
Diffstat (limited to 'sys/mwcs/mwmkidmd.x')
| -rw-r--r-- | sys/mwcs/mwmkidmd.x | 18 |
1 files changed, 18 insertions, 0 deletions
diff --git a/sys/mwcs/mwmkidmd.x b/sys/mwcs/mwmkidmd.x new file mode 100644 index 00000000..acdbb077 --- /dev/null +++ b/sys/mwcs/mwmkidmd.x @@ -0,0 +1,18 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# MW_MKIDMD -- Make the identity matrix. + +procedure mw_mkidmd (ltm, ndim) + +double 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 |