blob: acdbb077754309b02168ac84a7cef600b1a34da2 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
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
|
