source: doc/theses/mike_brooks_MMath/programs/hello-md.cfa @ 558999f

ADTast-experimentalpthread-emulation
Last change on this file since 558999f was bbf6a180, checked in by Peter A. Buhr <pabuhr@…>, 3 years ago

change to uw-ethesis and restructure source

  • Property mode set to 100644
File size: 1.6 KB
Line 
1#include "array.hfa"
2
3
4trait ix( C &, E &, ztype(N) ) {
5    E & ?[?]( C &, ptrdiff_t );
6    void __taglen( tag(C), tag(N) );
7};
8
9forall( ztype(Zn), ztype(S), Timmed &, Tbase & )
10void __taglen( tag(arpk(Zn, S, Timmed, Tbase)), tag(Zn) ) {}
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
40forall( ztype( N ) )
41void print1d_cstyle( array(float, N) & c );
42
43forall( C &, ztype( N ) | ix( C, float, N ) )
44void print1d( C & c );
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60forall( ztype( N ) )
61void print1d_cstyle( array(float, N) & c ) {
62    for( i; z(N) ) {
63        printf("%.1f  ", c[i]);
64    }
65    printf("\n");
66}
67
68
69
70
71
72
73
74
75
76
77
78
79
80forall( C &, ztype( N ) | ix( C, float, N ) )
81void print1d( C & c ) {
82    for( i; z(N) ) {
83        printf("%.1f  ", c[i]);
84    }
85    printf("\n");
86}
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101void fill( array(float, Z(5), Z(7)) & a ) {
102    for ( i; (ptrdiff_t) 5 ) {
103        for ( j; 7 ) {
104            a[[i,j]] = 1.0 * i + 0.1 * j;
105            printf("%.1f  ", a[[i,j]]);
106        }
107        printf("\n");
108    }
109    printf("\n");
110}
111
112int main() {
113
114
115
116
117
118
119
120array( float, Z(5), Z(7) ) a;
121fill(a);
122/*
1230.0  0.1  0.2  0.3  0.4  0.5  0.6 
1241.0  1.1  1.2  1.3  1.4  1.5  1.6 
1252.0  2.1  2.2  2.3  2.4  2.5  2.6 
1263.0  3.1  3.2  3.3  3.4  3.5  3.6 
1274.0  4.1  4.2  4.3  4.4  4.5  4.6
128*/
129   
130
131
132
133
134
135
136
137
138
139
140print1d_cstyle( a[ 2 ] );  // 2.0  2.1  2.2  2.3  2.4  2.5  2.6
141
142
143
144
145print1d( a[ 2 ] );  // 2.0  2.1  2.2  2.3  2.4  2.5  2.6
146
147
148
149
150print1d( a[[ 2, all ]] );  // 2.0  2.1  2.2  2.3  2.4  2.5  2.6
151print1d( a[[ all, 3 ]] );  // 0.3  1.3  2.3  3.3  4.3
152
153
154
155print1d_cstyle( a[[ 2, all ]] );
156
157
158
159
160
161
162
163#ifdef SHOWERR1
164
165print1d_cstyle( a[[ all, 2 ]] );  // bad
166
167#endif
168
169}
170
171
172
Note: See TracBrowser for help on using the repository browser.