Algebrite/sources/adj.coffee at 51d40342210d7b4dc574bc1a27fcf28613eb624c · davidedc/Algebrite · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
### adj =====================================================================

Tags
----
scripting, JS, internal, treenode, general concept

Parameters
----------
m

General description
-------------------
Returns the adjunct of matrix m. The inverse of m is equal to adj(m) divided by det(m).

###


Eval_adj = ->
	push(cadr(p1))
	Eval()
	adj()

adj = ->
	i = 0
	j = 0
	n = 0

	save()

	p1 = pop()

	if (istensor(p1) && p1.tensor.ndim == 2 && p1.tensor.dim[0] == p1.tensor.dim[1])
		doNothing = 1
	else
		stop("adj: square matrix expected")

	n = p1.tensor.dim[0]

	p2 = alloc_tensor(n * n)

	p2.tensor.ndim = 2
	p2.tensor.dim[0] = n
	p2.tensor.dim[1] = n

	for i in [0...n]
		for j in [0...n]
			cofactor(p1, n, i, j)
			p2.tensor.elem[n * j + i] = pop(); # transpose

	push(p2)

	restore()