summaryrefslogtreecommitdiffstats
path: root/tests/autodiff/vector-arithmetic-jvp.slang
blob: a1abd12b49e9c11eceb4666a75fd2d1f0e876bb6 (plain)
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
//TEST(compute):COMPARE_COMPUTE_EX:-slang -compute -shaderobj -output-using-type
//TEST(compute, vulkan):COMPARE_COMPUTE_EX:-vk -compute -shaderobj -output-using-type
//TEST(compute):COMPARE_COMPUTE_EX:-cuda -compute -shaderobj -output-using-type

//TEST_INPUT:ubuffer(data=[0 0 0 0], stride=4):out,name=outputBuffer
RWStructuredBuffer<float> outputBuffer;

typedef DifferentialPair<float2> dpfloat2;
typedef DifferentialPair<float3> dpfloat3;
typedef DifferentialPair<float4> dpfloat4;

[ForwardDifferentiable]
float3 f(float3 x)
{
    return x;
}

[ForwardDifferentiable]
float3 g(float3 x, float3 y)
{
    float3 a = x + y;
    float3 b = x - y;
    return a * b + 2 * x * y;
}

[ForwardDifferentiable]
float2 h(float2 x, float2 y)
{
    float2 a = x + y;
    float2 b = x - y;
    return a * b + 2 * x * y;
}

[ForwardDifferentiable]
float4 j(float4 x, float4 y)
{
    float4 a = x + y;
    float4 b = x - y;
    return a * b + 2 * x * y;
}

[numthreads(1, 1, 1)]
void computeMain(uint3 dispatchThreadID : SV_DispatchThreadID)
{
    {
        float3 a = float3(2.0, 2.0, 2.0);
        float3 b = float3(1.5, 1.5, 1.5);
        float3 da = float3(1.0, 1.0, 1.0);

        float2 a2 = float2(2.0, 1.0);
        float2 b2 = float2(1.5, -2.0);

        float4 a4 = float4(2.0, 1.0, 0.0, 2.0);
        float4 b4 = float4(1.5, -2.0, 1.0, 1.5);

        outputBuffer[0] = __fwd_diff(f)(dpfloat3(a, da)).d.z;    // Expect: 1

        outputBuffer[1] = __fwd_diff(g)(
            dpfloat3(a, da),
            dpfloat3(b, float3(2.0, 1.0, 0.0))).d.y;        // Expect: 8

        outputBuffer[2] = __fwd_diff(h)(
            dpfloat2(a2, float2(1.0, 0.0)),
            dpfloat2(b2, float2(1.0, 1.0))).d.x;            // Expect: 8

        outputBuffer[3] = __fwd_diff(j)(
            dpfloat4(a4, float4(1.0)),
            dpfloat4(b4, float4(2.0))).d.w;                 // Expect: 9
    }
}