File: solve.asy

package info (click to toggle)
texlive-bin 2009-8
  • links: PTS
  • area: main
  • in suites: squeeze
  • size: 231,976 kB
  • ctags: 165,663
  • sloc: ansic: 1,113,212; cpp: 454,597; sh: 49,189; perl: 43,889; makefile: 15,564; python: 11,825; xml: 7,446; lisp: 6,494; asm: 4,361; lex: 3,571; java: 3,569; pascal: 2,487; yacc: 2,451; exp: 2,031; ruby: 2,012; ada: 1,681; objc: 1,362; cs: 879; sed: 519; csh: 47; tcl: 45; awk: 35
file content (35 lines) | stat: -rw-r--r-- 966 bytes parent folder | download | duplicates (14) 
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
 
import TestLib;

StartTest("solve");
real[][] a=new real[][] {{1,1,1},{1,2,2},{0,0,1}};
real[][] b=new real[][] {{3,9},{5,11},{1,3}};
real[][] x=new real[][] {{1,7},{1,-1,},{1,3}};
real[][] c=solve(a,b);
for(int i=0; i < c.length; ++i)
  for(int j=0; j < c[i].length; ++j)
    assert(close(c[i][j],x[i][j]));

real[][] a={{1,-2,3,0},{4,-5,6,2},{-7,-8,10,5},{1,50,1,-2}};
real[] b={7,19,33,3};
real[] x=solve(a,b);
real[] c=a*x;
for(int i=0; i < c.length; ++i)
  assert(close(c[i],b[i]));
EndTest();

StartTest("inverse");
real[][] a=new real[][] {{1,1,1},{1,2,2},{0,0,1}};
real[][] ainverse=new real[][] {{2,-1,0},{-1,1,-1},{0,0,1}};
real[][] d=inverse(a);
real[][] l=d*a;
real[][] r=a*d;
real[][] I=identity(a.length);
for(int i=0; i < d.length; ++i) {
  for(int j=0; j < d[i].length; ++j) {
    assert(close(d[i][j],ainverse[i][j]));
    assert(I[i][j] == (i == j ? 1 : 0));
    assert(close(l[i][j],I[i][j]));
    assert(close(r[i][j],I[i][j]));
  }
}
EndTest();