|
1 |
| -allapprox(x, y, tol = 1e-8) = all(isapprox.(x, y, atol = tol)) |
2 | 1 | FDM = central_fdm(5, 1)
|
3 | 2 |
|
4 | 3 | function gradient(::Val{:Zygote}, f::Function, args)
|
|
20 | 19 | function gradient(::Val{:FiniteDiff}, f::Function, args)
|
21 | 20 | first(FiniteDifferences.grad(FDM, f, args))
|
22 | 21 | end
|
23 |
| - |
24 |
| - |
25 |
| - |
26 |
| -function transform_AD(::Val{:Zygote}, t::Transform, A) |
27 |
| - ps = KernelFunctions.params(t) |
28 |
| - @test allapprox( |
29 |
| - first(Zygote.gradient(p -> transform_with_duplicate(p, t, A), ps)), |
30 |
| - first(FiniteDifferences.grad( |
31 |
| - FDM, |
32 |
| - p -> transform_with_duplicate(p, t, A), |
33 |
| - ps, |
34 |
| - )), |
35 |
| - ) |
36 |
| - @test allapprox( |
37 |
| - first(Zygote.gradient(X -> sum(transform(t, X, 2)), A)), |
38 |
| - first(FiniteDifferences.grad(FDM, X -> sum(transform(t, X, 2)), A)), |
39 |
| - ) |
40 |
| -end |
41 |
| - |
42 |
| -function transform_AD(::Val{:ForwardDiff}, t::Transform, A) |
43 |
| - ps = KernelFunctions.params(t) |
44 |
| - if t isa ScaleTransform |
45 |
| - @test allapprox( |
46 |
| - first(ForwardDiff.gradient( |
47 |
| - p -> transform_with_duplicate(first(p), t, A), |
48 |
| - [ps], |
49 |
| - )), |
50 |
| - first(FiniteDifferences.grad( |
51 |
| - FDM, |
52 |
| - p -> transform_with_duplicate(p, t, A), |
53 |
| - ps, |
54 |
| - )), |
55 |
| - ) |
56 |
| - else |
57 |
| - @test allapprox( |
58 |
| - ForwardDiff.gradient(p -> transform_with_duplicate(p, t, A), ps), |
59 |
| - first(FiniteDifferences.grad( |
60 |
| - FDM, |
61 |
| - p -> transform_with_duplicate(p, t, A), |
62 |
| - ps, |
63 |
| - )), |
64 |
| - ) |
65 |
| - end |
66 |
| - @test allapprox( |
67 |
| - ForwardDiff.gradient(X -> sum(transform(t, X, 2)), A), |
68 |
| - first(FiniteDifferences.grad(FDM, X -> sum(transform(t, X, 2)), A)), |
69 |
| - ) |
70 |
| -end |
0 commit comments