File: transpose.asy

package info (click to toggle)
asymptote 1.43-1
  • links: PTS
  • area: main
  • in suites: lenny
  • size: 6,156 kB
  • ctags: 6,328
  • sloc: cpp: 26,241; python: 4,327; ansic: 1,895; perl: 1,405; lisp: 1,329; yacc: 547; sh: 431; lex: 406; makefile: 359
file content (55 lines) | stat: -rw-r--r-- 1,096 bytes parent folder | download | duplicates (13)
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
53
54
55
import TestLib;
import math;

StartTest("transpose");
int n=3;
real[][] a=new real[n][n]; 
real[][] b=new real[n][n]; 
for(int i=0; i < n; ++i) {
  for(int j=0; j < n; ++j) {
    a[i][j]=b[j][i]=rand();
  }
}

bool operator == (real[][] a, real[][] b)
{
  int n=a.length;
  for(int i=0; i < n; ++i)
    if(!all(a[i] == b[i])) return false;
  return true;
}

bool operator == (real[][][] a, real[][][] b)
{
  int n=a.length;
  for(int i=0; i < n; ++i) {
    real[][] ai=a[i];
    real[][] bi=b[i];
    int m=ai.length;
    for(int j=0; j < m; ++j) {
      if(!all(ai[j] == bi[j])) return false;
    }
  }
  return true;
}

assert(a == transpose(b));

int n=3;
real[][][] a=new real[n][n][n]; 
real[][][] b=new real[n][n][n]; 
real[][][] c=new real[n][n][n]; 
real[][][] d=new real[n][n][n]; 
for(int i=0; i < n; ++i) {
  for(int j=0; j < n; ++j) {
    for(int k=0; k < n; ++k) {
      a[i][j][k]=b[j][i][k]=c[i][k][j]=d[k][j][i]=rand();
    }
  }
}

assert(a == transpose(b,new int[] {1,0,2}));
assert(a == transpose(c,new int[] {0,2,1}));
assert(a == transpose(d,new int[] {2,1,0}));

EndTest();