seed!(42)
L = 3
D = 10
N = 20
nSamples = 10
b = 5
x = rand(N, D)
y = rand(N)
@testset "MCIntegrationVI" begin
    C = 10.0
    i = AGP.tuple_inference(MCIntegrationVI(clipping = C), 1, D, N, N, [], [])
    l = GaussianLikelihood()
    # @test AGP.grad_quad(l, 200.0, 0.0, 1.0, i) == (C, -C)
end
