mirror of
https://github.com/firestar5683/StarPilot.git
synced 2026-10-28 00:53:45 +08:00
770 lines
26 KiB
Python
770 lines
26 KiB
Python
#!/usr/bin/env python
|
|
import unittest
|
|
import numpy as np
|
|
from test.helpers import assert_jit_cache_len, not_support_multi_device, needs_second_gpu
|
|
from tinygrad import Tensor, TinyJit, Device, dtypes, Context, UOp
|
|
from tinygrad.helpers import DEV, GlobalCounters
|
|
from tinygrad.engine.jit import JitError
|
|
from tinygrad.renderer.isa.x86 import X86Renderer
|
|
|
|
def _simple_test(add, extract=lambda x: x, N=10):
|
|
for _ in range(5):
|
|
a = Tensor.randn(N, N)
|
|
b = Tensor.randn(N, N)
|
|
c = add(a, b)
|
|
np.testing.assert_allclose(extract(c).numpy(), a.numpy()+b.numpy(), atol=1e-4, rtol=1e-5)
|
|
assert_jit_cache_len(add, 1)
|
|
|
|
class TestJit(unittest.TestCase):
|
|
def test_jit_input_view(self):
|
|
@TinyJit
|
|
def f(x): return (x[2:5].contiguous() + 1).realize()
|
|
for i in range(5):
|
|
x = (Tensor.arange(10).float() + i * 10).clone().realize()
|
|
np.testing.assert_allclose(f(x).numpy(), x.numpy()[2:5] + 1)
|
|
|
|
@unittest.skipIf(isinstance(Device[Device.DEFAULT].renderer, X86Renderer), "estimates are wrong for x86")
|
|
def test_global_counters_jit(self):
|
|
@TinyJit
|
|
def f(a, b):
|
|
c = (a + b).realize()
|
|
d = (c * 2).realize()
|
|
return (d - a).realize()
|
|
a, b = Tensor.randn(64, 64).realize(), Tensor.randn(64, 64).realize()
|
|
for _ in range(4):
|
|
GlobalCounters.reset()
|
|
f(a, b)
|
|
Device[a.device].synchronize()
|
|
self.assertGreater(GlobalCounters.global_mem, 0)
|
|
self.assertGreater(GlobalCounters.global_ops, 0)
|
|
|
|
def test_jit_assign(self, dtype=dtypes.float32):
|
|
@TinyJit
|
|
def add(a):
|
|
a += 1
|
|
a.realize()
|
|
a = Tensor.zeros(1, dtype=dtype).contiguous().realize()
|
|
for _ in range(5): add(a)
|
|
self.assertEqual(a.item(), 5)
|
|
|
|
def test_jit_assign_int8(self): self.test_jit_assign(dtypes.int8)
|
|
|
|
def test_jit_copyin(self):
|
|
@TinyJit
|
|
def f(a):
|
|
return a + Tensor([1,2,3])
|
|
for _ in range(5):
|
|
b = Tensor.randn(3)
|
|
c = f(b)
|
|
np.testing.assert_allclose(c.numpy(), b.numpy()+[1,2,3], atol=1e-4, rtol=1e-5)
|
|
|
|
def test_jitted_clone(self):
|
|
def f(a): return a.clone().realize()
|
|
jf = TinyJit(f)
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 10, device=Device.DEFAULT).realize()
|
|
ja = jf(a)
|
|
np.testing.assert_allclose(a.numpy(), ja.numpy(), atol=1e-4, rtol=1e-5)
|
|
|
|
@needs_second_gpu
|
|
@unittest.skipIf(not_support_multi_device(), "no multi")
|
|
def test_jitted_transfers(self):
|
|
d0, d1 = f"{Device.DEFAULT}:0", f"{Device.DEFAULT}:1"
|
|
|
|
def f(a, b):
|
|
x = a.to(d1)
|
|
y = b.to(d1)
|
|
return x.realize(), y.realize()
|
|
|
|
jf = TinyJit(f)
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 10, device=d0).realize()
|
|
b = Tensor.randn(10, 10, device=d0).realize()
|
|
xc, yc = jf(a, b)
|
|
np.testing.assert_allclose(a.numpy(), xc.numpy(), atol=1e-4, rtol=1e-5)
|
|
np.testing.assert_allclose(b.numpy(), yc.numpy(), atol=1e-4, rtol=1e-5)
|
|
|
|
def test_jit_several_devs(self):
|
|
d0, d1 = f"{Device.DEFAULT}:0", "CPU"
|
|
|
|
def f(a, b):
|
|
x = a.to(d0).realize()
|
|
y = b.to(d0).realize()
|
|
return x+y.realize(), x*y.realize()
|
|
|
|
jf = TinyJit(f)
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 10, device=d1).realize()
|
|
b = Tensor.randn(10, 10, device=d1).realize()
|
|
zc, wc = jf(a, b)
|
|
np.testing.assert_allclose((a.numpy()+b.numpy()), zc.numpy(), atol=1e-4, rtol=1e-5)
|
|
np.testing.assert_allclose((a.numpy()*b.numpy()), wc.numpy(), atol=1e-4, rtol=1e-5)
|
|
|
|
@needs_second_gpu
|
|
@unittest.skipIf(not_support_multi_device(), "no multi")
|
|
def test_jitted_view(self):
|
|
d0, d1 = f"{Device.DEFAULT}:0", f"{Device.DEFAULT}:1"
|
|
|
|
def f(a):
|
|
x1 = a.sum(axis=(1,))
|
|
x = (x1 + 5).bitcast(dtypes.int32)
|
|
y = x.to(d1)
|
|
return y.realize()
|
|
|
|
jf = TinyJit(f)
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 1000, device=d0).realize()
|
|
xc = jf(a)
|
|
np.testing.assert_allclose(a.numpy().sum(axis=(1,)) + 5, xc.numpy().view(np.float32), atol=1e-4, rtol=5e-5)
|
|
|
|
def test_jitbeam_triggers_beam(self):
|
|
from unittest.mock import patch
|
|
from tinygrad.helpers import getenv as _getenv
|
|
@TinyJit
|
|
def add(a, b): return (a+b).realize()
|
|
a, b = Tensor.ones(10, 10).contiguous().realize(), Tensor.ones(10, 10).contiguous().realize()
|
|
with patch("tinygrad.codegen.opt.search.beam_search", wraps=lambda k,*a,**kw: k) as mock_beam:
|
|
add(a, b)
|
|
assert mock_beam.call_count == 0
|
|
with patch("tinygrad.engine.jit.getenv", side_effect=lambda k, d=0: 1 if k == "JITBEAM" else _getenv(k, d)): add(a, b)
|
|
assert mock_beam.call_count == 1
|
|
|
|
def test_simple_jit_reset(self):
|
|
@TinyJit
|
|
def add(a, b): return (a+b).realize()
|
|
_simple_test(add)
|
|
add.reset()
|
|
_simple_test(add, N=20)
|
|
|
|
def test_simple_jit_norealize(self):
|
|
@TinyJit
|
|
def add(a, b): return (a+b)
|
|
_simple_test(add)
|
|
|
|
def test_simple_jit_norealize_list(self):
|
|
@TinyJit
|
|
def add(a, b): return [a+b]
|
|
_simple_test(add, extract=lambda x: x[0])
|
|
|
|
def test_simple_jit_norealize_dict(self):
|
|
@TinyJit
|
|
def add(a, b): return {"billy": a+b}
|
|
_simple_test(add, extract=lambda x: x["billy"])
|
|
|
|
def test_jit_multiple_outputs(self):
|
|
@TinyJit
|
|
def f(a, b): return (a+b).realize(), (a-b).realize(), (a*b).realize()
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 10)
|
|
b = Tensor.randn(10, 10)
|
|
c, d, e = f(a, b)
|
|
np.testing.assert_allclose(c.numpy(), a.numpy()+b.numpy(), atol=1e-4, rtol=1e-5)
|
|
np.testing.assert_allclose(d.numpy(), a.numpy()-b.numpy(), atol=1e-4, rtol=1e-5)
|
|
np.testing.assert_allclose(e.numpy(), a.numpy()*b.numpy(), atol=1e-4, rtol=1e-5)
|
|
assert_jit_cache_len(f, 3)
|
|
|
|
def test_nothing_jitted(self):
|
|
@TinyJit
|
|
def add(a, b): return None
|
|
with self.assertRaises(JitError):
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 10)
|
|
b = Tensor.randn(10, 10)
|
|
add(a, b)
|
|
|
|
def test_jit_zero_does_not_jit(self):
|
|
@TinyJit
|
|
def add(a, b): return (a+b).realize()
|
|
with Context(JIT=0):
|
|
for i in range(5):
|
|
a = Tensor([i])
|
|
b = Tensor([i])
|
|
c = add(a, b)
|
|
np.testing.assert_allclose(c.numpy(), 2*i)
|
|
assert_jit_cache_len(add, 0)
|
|
|
|
def test_jit_not_capturing(self):
|
|
@TinyJit
|
|
def add(a, b):
|
|
Tensor.zeros(4, 4).contiguous().realize() # no-op kernel is captured
|
|
return (a+b).realize()
|
|
for i in range(5):
|
|
a = Tensor([i])
|
|
b = Tensor([i])
|
|
c = add(a, b)
|
|
np.testing.assert_allclose(c.numpy(), 2*i)
|
|
assert_jit_cache_len(add, 2)
|
|
|
|
@TinyJit
|
|
def add2(a, b):
|
|
with Context(CAPTURING=0): # not captured
|
|
Tensor.zeros(4, 4).contiguous().realize()
|
|
return (a+b).realize()
|
|
for i in range(5):
|
|
a = Tensor([i])
|
|
b = Tensor([i])
|
|
c = add2(a, b)
|
|
np.testing.assert_allclose(c.numpy(), 2*i)
|
|
assert_jit_cache_len(add2, 1)
|
|
|
|
def test_jit_shape_mismatch(self):
|
|
@TinyJit
|
|
def add(a, b): return (a+b).realize()
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 10)
|
|
b = Tensor.randn(10, 10)
|
|
add(a, b)
|
|
bad = Tensor.randn(20, 20)
|
|
with self.assertRaises(JitError):
|
|
add(a, bad)
|
|
|
|
def test_jit_shape_views_mismatch(self):
|
|
@TinyJit
|
|
def add(a): return (a+1).realize()
|
|
with self.assertRaises(JitError):
|
|
for i in range(1,5):
|
|
# a has an offset that the kernel doesn't know about
|
|
a = Tensor.randn(10, 10).realize()[:, i:i+2]
|
|
add(a)
|
|
|
|
def test_jit_duplicate_fail(self):
|
|
# the jit doesn't support duplicate arguments
|
|
@TinyJit
|
|
def add(a, b): return (a+b).realize()
|
|
a = Tensor.randn(10, 10)
|
|
with self.assertRaises(JitError):
|
|
add(a, a)
|
|
|
|
def test_kwargs_jit(self):
|
|
@TinyJit
|
|
def add_kwargs(first, second): return (first+second).realize()
|
|
for _ in range(5):
|
|
a = Tensor.randn(10, 10)
|
|
b = Tensor.randn(10, 10)
|
|
c = add_kwargs(first=a, second=b)
|
|
np.testing.assert_allclose(c.numpy(), a.numpy()+b.numpy(), atol=1e-4, rtol=1e-5)
|
|
assert_jit_cache_len(add_kwargs, 1)
|
|
|
|
def test_array_jit(self):
|
|
@TinyJit
|
|
def add_array(a, arr): return (a+arr[0]).realize()
|
|
for _ in range(5):
|
|
a, b = Tensor.randn(10, 10).realize(), Tensor.randn(10, 10).realize()
|
|
np.testing.assert_allclose(add_array(a, [b]).numpy(), a.numpy()+b.numpy(), atol=1e-4, rtol=1e-5)
|
|
assert_jit_cache_len(add_array, 1)
|
|
|
|
def test_method_jit(self):
|
|
class Fun:
|
|
def __init__(self):
|
|
self.a = Tensor.randn(10, 10)
|
|
@TinyJit
|
|
def __call__(self, b:Tensor) -> Tensor:
|
|
return (self.a+b).realize()
|
|
fun = Fun()
|
|
for _ in range(5):
|
|
b = Tensor.randn(10, 10)
|
|
c = fun(b)
|
|
np.testing.assert_allclose(c.numpy(), fun.a.numpy()+b.numpy(), atol=1e-4, rtol=1e-5)
|
|
assert_jit_cache_len(fun.__call__.func.__self__, 1)
|
|
|
|
def test_jit_size1_input(self):
|
|
@TinyJit
|
|
def f(a, b): return (a+b).realize()
|
|
a = Tensor([1, 2, 3])
|
|
for i in range(5):
|
|
np.testing.assert_allclose(f(a, Tensor([i])).numpy(), (a+i).numpy(), atol=1e-4, rtol=1e-5)
|
|
assert_jit_cache_len(f, 1)
|
|
|
|
def test_jit_output_non_tensor_fail(self):
|
|
@TinyJit
|
|
def f(a, b, i): return (a+b).realize(), i
|
|
with self.assertRaises(JitError):
|
|
for i in range(3):
|
|
f(Tensor.randn(10, 10), Tensor.randn(10, 10), i)
|
|
|
|
def test_jit_random_regen(self):
|
|
def f(a, b):
|
|
rn = Tensor.randn(*a.shape)
|
|
return ((a+b)*rn).realize()
|
|
a = Tensor.randn(10, 10).realize() # realize these before resetting the random seed
|
|
b = Tensor.randn(10, 10).realize()
|
|
|
|
Tensor.manual_seed(1234)
|
|
jf = TinyJit(f)
|
|
res = set()
|
|
for _ in range(5):
|
|
o1 = jf(a, b)
|
|
res.add(o1.numpy()[0][0])
|
|
assert len(res) == 5, "All values should be different, rand works in jit."
|
|
|
|
Tensor.manual_seed(1234)
|
|
jf2 = TinyJit(f)
|
|
res2 = set()
|
|
for _ in range(5):
|
|
o1 = jf2(a, b)
|
|
res2.add(o1.numpy()[0][0])
|
|
assert len(res2) == 5, "All values should be different, rand works in jit."
|
|
assert res == res2, "Jit rand is not reproducible with the same seed"
|
|
|
|
Tensor.manual_seed(3421)
|
|
jf3 = TinyJit(f)
|
|
res3 = set()
|
|
for _ in range(5):
|
|
o1 = jf3(a, b)
|
|
res3.add(o1.numpy()[0][0])
|
|
assert len(res3) == 5, "All values should be different, rand works in jit."
|
|
assert res3 != res2, "Jit rand is diff with diff seeds"
|
|
|
|
def test_jit_v_nojit_random_regen(self):
|
|
def f(a, b):
|
|
rn = Tensor.randn(*a.shape)
|
|
rn = rn * a
|
|
rn2 = Tensor.randn(*a.shape)
|
|
rn2 = rn2 * b
|
|
rn = rn + rn2
|
|
rn2 = rn2 + Tensor.randn(*a.shape)
|
|
return ((a+b)*rn).realize(), ((a+b)*rn2).realize()
|
|
Tensor.manual_seed(0)
|
|
a = Tensor.randn(10, 10).realize() # realize these before resetting the random seed
|
|
b = Tensor.randn(10, 10).realize()
|
|
|
|
Tensor.manual_seed(1234)
|
|
without_jit = set()
|
|
for _ in range(5):
|
|
o1, o2 = f(a, b)
|
|
without_jit.add(o1.numpy()[0][0])
|
|
without_jit.add(o2.numpy()[0][0])
|
|
assert len(without_jit) == 10, "All values should be different."
|
|
|
|
Tensor.manual_seed(1234)
|
|
jf = TinyJit(f)
|
|
with_jit = set()
|
|
for _ in range(5):
|
|
o1, o2 = jf(a, b)
|
|
with_jit.add(o1.numpy()[0][0])
|
|
with_jit.add(o2.numpy()[0][0])
|
|
assert len(with_jit) == 10, "All values should be different."
|
|
assert with_jit == without_jit, "jit and non-jit should produce the same random values with the same seed"
|
|
|
|
def test_jit_multiple_random_regen(self):
|
|
def f(a, b):
|
|
rn = Tensor.randn(*a.shape)
|
|
rn = rn * a
|
|
rn2 = Tensor.randn(*a.shape)
|
|
rn2 = rn2 * b
|
|
rn = rn + rn2
|
|
rn2 = rn2 + Tensor.randn(*a.shape)
|
|
return ((a+b)*rn).realize(), ((a+b)*rn2).realize()
|
|
a = Tensor.randn(10, 10).realize() # realize these before resetting the random seed
|
|
b = Tensor.randn(10, 10).realize()
|
|
|
|
Tensor.manual_seed(1234)
|
|
jf = TinyJit(f)
|
|
res = set()
|
|
for _ in range(5):
|
|
o1, o2 = jf(a, b)
|
|
res.add(o1.numpy()[0][0])
|
|
res.add(o2.numpy()[0][0])
|
|
assert len(res) == 10, "All values should be different, rand works in jit."
|
|
|
|
Tensor.manual_seed(1234)
|
|
jf2 = TinyJit(f)
|
|
res2 = set()
|
|
for _ in range(5):
|
|
o1, o2 = jf2(a, b)
|
|
res2.add(o1.numpy()[0][0])
|
|
res2.add(o2.numpy()[0][0])
|
|
assert len(res2) == 10, "All values should be different, rand works in jit."
|
|
assert res == res2, "Jit rand is not reproducible with the same seed"
|
|
|
|
Tensor.manual_seed(3421)
|
|
jf3 = TinyJit(f)
|
|
res3 = set()
|
|
for _ in range(5):
|
|
o1, o2 = jf3(a, b)
|
|
res3.add(o1.numpy()[0][0])
|
|
res3.add(o2.numpy()[0][0])
|
|
assert len(res3) == 10, "All values should be different, rand works in jit."
|
|
assert res3 != res2, "Jit rand is diff with diff seeds"
|
|
|
|
def test_jit_random_after_unrealized_random(self):
|
|
@TinyJit
|
|
def f(): return Tensor.rand()
|
|
Tensor.manual_seed(1234)
|
|
Tensor.rand()
|
|
res = [f().numpy() for _ in range(3)]
|
|
assert res[1] != res[2]
|
|
|
|
def test_jit_realization_and_sampling(self):
|
|
w = Tensor.eye(5)
|
|
|
|
@TinyJit
|
|
def foo (x): return w.dot(x).realize()
|
|
|
|
arg = [
|
|
Tensor([1,2,3,4,5]),
|
|
Tensor([1,3,3,4,6]),
|
|
Tensor([1,2,5,4,7]),
|
|
Tensor([0,2,3,1,0]),
|
|
]
|
|
|
|
Y = [foo(e).numpy() for e in arg]
|
|
|
|
foo(Tensor([7,7,7,7,7]))
|
|
want = [[1., 2., 3., 4., 5.],
|
|
[1., 3., 3., 4., 6.],
|
|
[1., 2., 5., 4., 7.],
|
|
[0., 2., 3., 1., 0.]]
|
|
np.testing.assert_allclose(want, Y)
|
|
|
|
def test_jit_output_clone(self):
|
|
@TinyJit
|
|
def f(x:Tensor) -> Tensor: return (x + 1).realize()
|
|
|
|
f(Tensor([0.0]))
|
|
f(Tensor([0.0]))
|
|
|
|
a = f(Tensor([1.0])).clone().realize()
|
|
b = f(Tensor([2.0]))
|
|
assert abs((a - b).item()) > 0.5
|
|
|
|
def test_jit_init_empty(self):
|
|
@TinyJit
|
|
def f(x:Tensor) -> Tensor: return (x + 1).realize()
|
|
|
|
f(Tensor.empty(1))
|
|
f(Tensor.empty(1))
|
|
# scalar const input is not allowed
|
|
with self.assertRaises(JitError):
|
|
f(Tensor(2.0)).item()
|
|
# self.assertEqual(f(Tensor([2.0])).item(), 1.0) # TODO: wrong output, should be 3.0. currently depends on empty value
|
|
|
|
def test_jit_const_input(self):
|
|
@TinyJit
|
|
def f(x:Tensor) -> Tensor: return (x + 1).realize()
|
|
with self.assertRaises(JitError):
|
|
f(Tensor(UOp.const(2.0).cast(dtypes.float))).item()
|
|
|
|
def test_jit_deviceless_compute_input(self):
|
|
@TinyJit
|
|
def f(x:Tensor) -> Tensor: return (x + 1).realize()
|
|
with self.assertRaises(JitError):
|
|
f(Tensor(UOp.const(2.0).cast(dtypes.float) + UOp.const(1.0).cast(dtypes.float))).item()
|
|
|
|
def test_jit_init_empty_alt(self):
|
|
@TinyJit
|
|
def f(a:Tensor, b:Tensor) -> Tensor: return b.assign(a+1)
|
|
for i in range(4):
|
|
a = Tensor([i])
|
|
b = Tensor.empty_like(a)
|
|
c = f(a, b)
|
|
self.assertEqual(c.item(), i+1)
|
|
|
|
def test_jit_lazy_grad_after_replay(self):
|
|
# the lazy .grad created during capture is read outside the JIT, the memory planner must not suballocate its buffers (issue #16571)
|
|
from tinygrad import nn
|
|
def step(conv, x, y):
|
|
out = conv(x.permute(0, 3, 1, 2).contiguous()).relu().flatten(1)
|
|
loss = (out * y).sum(axis=1) # per-example loss
|
|
loss.sum().backward()
|
|
conv.weight.grad = None
|
|
(loss * 0.5).sum().backward()
|
|
return loss.mean().realize()
|
|
|
|
Tensor.manual_seed(42)
|
|
conv = nn.Conv2d(3, 4, kernel_size=3, padding=1)
|
|
x, y = Tensor.randn(4, 8, 8, 3).realize(), Tensor.randn(4, 4*8*8).realize()
|
|
|
|
step(conv, x, y)
|
|
ref = conv.weight.grad.numpy()
|
|
jit_step = TinyJit(step)
|
|
for _ in range(4):
|
|
jit_step(conv, x, y)
|
|
np.testing.assert_allclose(conv.weight.grad.numpy(), ref, atol=1e-4, rtol=1e-5)
|
|
|
|
@unittest.skip("Pending multioutput implementation #3607")
|
|
class TestMultioutputJit(unittest.TestCase):
|
|
def _test(self, f):
|
|
for _ in range(5):
|
|
a, b = Tensor.randn(10, 10), Tensor.randn(10, 10)
|
|
out0, out1, out2 = f(a, b)
|
|
np.testing.assert_allclose(out0.numpy(), a.numpy()+b.numpy(), atol=1e-4, rtol=1e-5)
|
|
np.testing.assert_allclose(out1.numpy(), a.numpy()-b.numpy(), atol=1e-4, rtol=1e-5)
|
|
np.testing.assert_allclose(out2.numpy(), a.numpy()*b.numpy(), atol=1e-4, rtol=1e-5)
|
|
|
|
def test_jit_multioutput_realize(self):
|
|
@TinyJit
|
|
def fxn(a, b): return (a+b).realize(), (a-b).realize(), (a*b).realize()
|
|
self._test(fxn)
|
|
assert_jit_cache_len(fxn, 3)
|
|
|
|
def test_jit_multioutput_norealize(self):
|
|
@TinyJit
|
|
def fxn(a, b): return a+b, a-b, a*b
|
|
self._test(fxn)
|
|
assert_jit_cache_len(fxn, 1)
|
|
|
|
def test_jit_multioutput_mix(self):
|
|
@TinyJit
|
|
def fxn(a, b): return a+b, a-b, (a*b).realize()
|
|
self._test(fxn)
|
|
assert_jit_cache_len(fxn, 2)
|
|
|
|
class TestCopyInsideJit(unittest.TestCase):
|
|
def test_copy_inside_jit(self):
|
|
@TinyJit
|
|
def add(x,y) -> Tensor: return x.to(Device.DEFAULT)+y
|
|
for _ in range(5):
|
|
# create a Tensor on CPU
|
|
a = Tensor.rand(16,16,device="CPU").realize()
|
|
b = Tensor.rand(16,16).realize()
|
|
out = add(a,b)
|
|
np.testing.assert_allclose(out.flatten().tolist(), [x+y for x,y in zip(a.flatten().tolist(), b.flatten().tolist())])
|
|
|
|
class TestJitPrune(unittest.TestCase):
|
|
def test_prune_w_copy_correct(self):
|
|
weights = Tensor.rand(16).realize()
|
|
def w2(x) -> Tensor: return (weights*2).contiguous() + x.to(Device.DEFAULT)
|
|
w2_noprune = TinyJit(w2)
|
|
w2_prune = TinyJit(w2, prune=True)
|
|
|
|
for _ in range(3):
|
|
a = Tensor.rand(16, device="CPU").realize()
|
|
out = w2_noprune(a)
|
|
np.testing.assert_allclose(out.tolist(), [x*2+y for x,y in zip(weights.tolist(), a.tolist())])
|
|
|
|
for _ in range(3):
|
|
a = Tensor.rand(16, device="CPU").realize()
|
|
out = w2_prune(a)
|
|
np.testing.assert_allclose(out.tolist(), [x*2+y for x,y in zip(weights.tolist(), a.tolist())])
|
|
|
|
def test_prune_w_independent_copy_correct(self):
|
|
weights = Tensor.rand(16, device="CPU").realize()
|
|
def w2(x) -> Tensor: return (weights*2).contiguous().to(Device.DEFAULT) + x
|
|
w2_noprune = TinyJit(w2)
|
|
w2_prune = TinyJit(w2, prune=True)
|
|
|
|
for _ in range(3):
|
|
a = Tensor.rand(16).realize()
|
|
out = w2_noprune(a)
|
|
np.testing.assert_allclose(out.tolist(), [x*2+y for x,y in zip(weights.tolist(), a.tolist())])
|
|
|
|
for _ in range(3):
|
|
a = Tensor.rand(16).realize()
|
|
out = w2_prune(a)
|
|
np.testing.assert_allclose(out.tolist(), [x*2+y for x,y in zip(weights.tolist(), a.tolist())])
|
|
|
|
assert_jit_cache_len(w2_prune, 1)
|
|
|
|
def test_simple_prune(self):
|
|
weights = Tensor.rand(16).realize()
|
|
def w2(x) -> Tensor: return (weights*2).contiguous() + x
|
|
w2_noprune = TinyJit(w2)
|
|
w2_prune = TinyJit(w2, prune=True)
|
|
|
|
for _ in range(3):
|
|
a = Tensor.rand(16).realize()
|
|
out = w2_noprune(a)
|
|
np.testing.assert_allclose(out.tolist(), [x*2+y for x,y in zip(weights.tolist(), a.tolist())])
|
|
assert_jit_cache_len(w2_noprune, 2)
|
|
|
|
for _ in range(3):
|
|
a = Tensor.rand(16).realize()
|
|
out = w2_prune(a)
|
|
np.testing.assert_allclose(out.tolist(), [x*2+y for x,y in zip(weights.tolist(), a.tolist())])
|
|
assert_jit_cache_len(w2_prune, 1)
|
|
|
|
class TestJitFree(unittest.TestCase):
|
|
def test_free_intermediates(self):
|
|
ext_tensor = Tensor([1,24,23,45,1])
|
|
@TinyJit
|
|
def fxn(x:Tensor):
|
|
t1 = (x * 2).contiguous().realize()
|
|
t2 = (t1 + ext_tensor).contiguous().realize()
|
|
out = (t2.sum()).contiguous().realize()
|
|
return out
|
|
for i in range(5):
|
|
out = fxn(inp:=Tensor([i,1,2,3,4]))
|
|
self.assertEqual(out.item(), 114+2*i)
|
|
pre_free = GlobalCounters.mem_used
|
|
fxn.captured.free_intermediates()
|
|
savings_after_free = pre_free - GlobalCounters.mem_used
|
|
|
|
expected_savings = (len(inp) * inp.dtype.itemsize * 2) + dtypes.float32.itemsize # (t1 and t2) + out
|
|
|
|
self.assertGreaterEqual(savings_after_free, expected_savings)
|
|
out = fxn(Tensor([11,1,2,3,4]))
|
|
self.assertEqual(out.item(), 136)
|
|
|
|
# Try one more time...
|
|
pre_free = GlobalCounters.mem_used
|
|
fxn.captured.free_intermediates()
|
|
fxn.captured.free_intermediates() # 2nd time to validate
|
|
savings_after_free = pre_free - GlobalCounters.mem_used
|
|
|
|
self.assertGreaterEqual(savings_after_free, expected_savings)
|
|
out = fxn(Tensor([11,1,2,3,4]))
|
|
self.assertEqual(out.item(), 136)
|
|
|
|
def test_updated_not_freed(self):
|
|
x = Tensor([1]).realize()
|
|
@TinyJit
|
|
def fxn(y):
|
|
nonlocal x
|
|
x += y
|
|
return x
|
|
for _ in range(5): fxn(Tensor([1]))
|
|
self.assertEqual(x.item(), 6)
|
|
pre_free = GlobalCounters.mem_used
|
|
fxn.captured.free_intermediates()
|
|
savings_after_free = pre_free - GlobalCounters.mem_used
|
|
self.assertEqual(savings_after_free, 0)
|
|
fxn(Tensor([2]))
|
|
self.assertEqual(x.item(), 8)
|
|
|
|
class TestJitGraphSplit(unittest.TestCase):
|
|
def compute(self, device, inp):
|
|
assert inp.device == device, f"Input device {inp.device} does not match expected {device}"
|
|
return (inp + 1.0).contiguous().realize()
|
|
|
|
def copy(self, device, to_device, inp):
|
|
assert inp.device == device, f"Input device {inp.device} does not match expected {device}"
|
|
return inp.to(to_device).realize()
|
|
|
|
def expect(self, f, *args):
|
|
def _numpies(tpl): return tpl.numpy() if tpl.__class__ is Tensor else tuple([t.numpy() for t in tpl])
|
|
|
|
expected = _numpies(f(*args))
|
|
for i in range(4):
|
|
res = _numpies(f(*args))
|
|
np.testing.assert_allclose(res, expected, atol=1e-4, rtol=1e-5)
|
|
|
|
def test_jit_split_simple(self):
|
|
@TinyJit
|
|
def f(inp):
|
|
op0 = self.compute(Device.DEFAULT, inp)
|
|
op1 = self.compute(Device.DEFAULT, op0)
|
|
op2 = self.compute(Device.DEFAULT, op1)
|
|
return op2
|
|
|
|
inp = Tensor.randn(10, 10, device=Device.DEFAULT).realize()
|
|
self.expect(f, inp)
|
|
|
|
def test_jit_cpu_simple(self):
|
|
if Device.DEFAULT == "CPU": raise unittest.SkipTest("CPU is not a valid default device for this test")
|
|
|
|
@TinyJit
|
|
def f(inp, inp_cpu):
|
|
op0 = self.compute(Device.DEFAULT, inp)
|
|
op1 = self.compute(Device.DEFAULT, op0)
|
|
op2 = self.compute("CPU", inp_cpu)
|
|
op3 = self.compute(Device.DEFAULT, op1)
|
|
return op2, op3
|
|
|
|
inp = Tensor.randn(10, 10, device=Device.DEFAULT).realize()
|
|
inp_cpu = Tensor.randn(10, 10, device="CPU").realize()
|
|
self.expect(f, inp, inp_cpu)
|
|
|
|
def test_jit_cpu_several(self):
|
|
if Device.DEFAULT == "CPU": raise unittest.SkipTest("CPU is not a valid default device for this test")
|
|
|
|
@TinyJit
|
|
def f(inp, inp_cpu):
|
|
op0 = self.compute(Device.DEFAULT, inp)
|
|
op1 = self.compute(Device.DEFAULT, op0)
|
|
op2 = self.compute("CPU", inp_cpu)
|
|
op3 = self.compute("CPU", op2)
|
|
op4 = self.compute(Device.DEFAULT, op1)
|
|
return op3, op4
|
|
|
|
inp = Tensor.randn(10, 10, device=Device.DEFAULT).realize()
|
|
inp_cpu = Tensor.randn(10, 10, device="CPU").realize()
|
|
self.expect(f, inp, inp_cpu)
|
|
|
|
def test_jit_multidev(self):
|
|
if Device.DEFAULT == "CPU": raise unittest.SkipTest("CPU is not a valid default device for this test")
|
|
|
|
try: Device[f"{Device.DEFAULT}:1"]
|
|
except Exception: raise unittest.SkipTest("no multidevice")
|
|
|
|
@TinyJit
|
|
def f(inp, inp_d1):
|
|
op0 = self.compute(Device.DEFAULT, inp)
|
|
op1 = self.compute(Device.DEFAULT, op0)
|
|
op2 = self.compute(f"{Device.DEFAULT}:1", inp_d1)
|
|
op3 = self.compute(f"{Device.DEFAULT}:1", op2)
|
|
op4 = self.compute(Device.DEFAULT, op1)
|
|
return op3, op4
|
|
|
|
inp = Tensor.randn(10, 10, device=Device.DEFAULT).realize()
|
|
inp_d1 = Tensor.randn(10, 10, device=f"{Device.DEFAULT}:1").realize()
|
|
self.expect(f, inp, inp_d1)
|
|
|
|
def test_jit_multidev_xfer(self):
|
|
if Device.DEFAULT in {"CPU"}: raise unittest.SkipTest("CPU is not a valid default device for this test (zero-copies)")
|
|
if Device.DEFAULT == "METAL": raise unittest.SkipTest("Metal is flaky, with multidevice (same as metal llama 4gpu?)")
|
|
|
|
try: Device[f"{Device.DEFAULT}:1"]
|
|
except Exception: raise unittest.SkipTest("no multidevice")
|
|
|
|
@TinyJit
|
|
def f(inp, inp_d1):
|
|
op0 = self.compute(Device.DEFAULT, inp)
|
|
op1 = self.compute(Device.DEFAULT, op0)
|
|
op2 = self.compute(f"{Device.DEFAULT}:1", inp_d1)
|
|
op3 = self.copy(f"{Device.DEFAULT}:1", Device.DEFAULT, op2)
|
|
op4 = self.compute(f"{Device.DEFAULT}:1", op2)
|
|
op5 = self.compute(Device.DEFAULT, op3)
|
|
return op1, op4, op5
|
|
|
|
inp = Tensor.randn(10, 10, device=Device.DEFAULT).realize()
|
|
inp_d1 = Tensor.randn(10, 10, device=f"{Device.DEFAULT}:1").realize()
|
|
self.expect(f, inp, inp_d1)
|
|
|
|
@unittest.skip("this fails if you don't have SDMA or are using AMD_DISABLE_SDMA=1")
|
|
@unittest.skipIf(DEV.interface.startswith("MOCK"), "MockGPU does not support parallel copies")
|
|
def test_jit_multidev_copy(self):
|
|
if Device.DEFAULT in {"CPU"}: raise unittest.SkipTest("CPU/LLVM is not a valid default device for this test (zero-copies)")
|
|
|
|
@TinyJit
|
|
def f(inp):
|
|
op0 = self.compute(Device.DEFAULT, inp)
|
|
op1 = self.compute(Device.DEFAULT, op0)
|
|
op2 = self.copy(Device.DEFAULT, "CPU", op1)
|
|
op3 = self.compute("CPU", op2)
|
|
return op3
|
|
|
|
inp = Tensor.randn(10, 10, device=Device.DEFAULT).realize()
|
|
self.expect(f, inp)
|
|
|
|
class TestJitInsideJit(unittest.TestCase):
|
|
def test_jit_jit_error(self):
|
|
@TinyJit
|
|
def f(t): return t + 1
|
|
|
|
@TinyJit
|
|
def g(t): return f(t) * 3
|
|
|
|
# NOTE: first does not raise
|
|
self.assertEqual(g(Tensor([1])).realize().item(), 6)
|
|
with self.assertRaisesRegex(RuntimeError, "having TinyJit inside another TinyJit is not supported"):
|
|
g(Tensor([2])).realize()
|
|
|
|
class TestJitRandom(unittest.TestCase):
|
|
def test_jit_rangeify(self):
|
|
tst = {0:[], 1:[]}
|
|
for r in [0,1]:
|
|
Tensor.manual_seed(1337)
|
|
with Context(JIT=r):
|
|
_ = Tensor.randint(4, high=3)
|
|
# this second one makes the behavior different
|
|
_ = Tensor.randint(4, high=3)
|
|
@TinyJit
|
|
def f(): return Tensor.randint(20, high=5)
|
|
for _ in range(5): tst[r].append(f().tolist())
|
|
for i, (t0, t1) in enumerate(zip(tst[0], tst[1])):
|
|
self.assertListEqual(t0, t1, msg=f"mismatch at list {i}")
|
|
|
|
if __name__ == '__main__':
|
|
unittest.main()
|