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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
// RUN: %clang -std=c11 -O0 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme -S | %lli -
// RUN: %clang -std=c11 -O1 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme -S | %lli -
// RUN: %clang -std=c11 -O2 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme -S | %lli -
// RUN: %clang -std=c11 -O3 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme -S | %lli -
// RUN: %clang -std=c11 -O0 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme --enzyme-inline=1 -S | %lli -
// RUN: %clang -std=c11 -O1 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme --enzyme-inline=1 -S | %lli -
// RUN: %clang -std=c11 -O2 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme --enzyme-inline=1 -S | %lli -
// RUN: %clang -std=c11 -O3 %s -S -emit-llvm -o - | %opt - %OPloadEnzyme %enzyme --enzyme-inline=1 -S | %lli -
#include "../test_utils.h"
double __enzyme_autodiff(void*, ...);
double alldiv(double* a, int N, double start) {
for (int i=0; i<N; i++) {
start /= a[i];
}
return start;
}
double alldiv2(double* a, int N) {
double start = 3.14;
for (int i=0; i<N; i++) {
start /= a[i];
}
return start;
}
int main(int argc, char** argv) {
int N = 5;
double a[N];
for(int i=0; i<N; i++) {
a[i] = i+1;
}
double d_a[N];
for(int i=0; i<N; i++)
d_a[i] = 0.0f;
double res = alldiv(a, N, 3.14);
//omp(*a, N);
printf("ran res=%e\n", res);
double dstart = __enzyme_autodiff((void*)alldiv, a, d_a, N, 3.14);
for(int i=0; i<N; i++) {
printf("a[%d]=%f d_a[%d]=%e cor=%e\n", i, a[i], i, d_a[i], - res / a[i]);
}
fflush(0);
for(int i=0; i<N; i++) {
printf("i=%d\n", i);
APPROX_EQ(d_a[i], - res / a[i], 1e-10);
}
APPROX_EQ(dstart, alldiv(a, N, 1.0), 1e-10);
for(int i=0; i<N; i++) {
a[i] = i+1;
}
for(int i=0; i<N; i++)
d_a[i] = 0.0f;
__enzyme_autodiff((void*)alldiv2, a, d_a, N);
for(int i=0; i<N; i++) {
printf("a[%d]=%f d_a[%d]=%e cor=%e\n", i, a[i], i, d_a[i], - res / a[i]);
}
fflush(0);
for(int i=0; i<N; i++) {
printf("i=%d\n", i);
APPROX_EQ(d_a[i], - res / a[i], 1e-10);
}
return 0;
}