def test_readout(): engine = ParticleMeshEngine(pm) code = CodeSegment(engine) s = pm.BoxSize / pm.Nmesh * 0.001 + 0.99 * engine.q / pm.Nmesh # sample all positions. field = pm.generate_whitenoise(seed=1234, mode='real') code.get_x(s='s', x='x') code.decompose(x='x', layout='layout') code.readout(x='x', mesh='density', layout='layout', value='value') check_grad(code, 'value', 'density', init={ 'density': field, 's': s }, eps=1e-4, rtol=1e-8) check_grad(code, 'value', 's', init={ 'density': field, 's': s }, eps=1e-4, rtol=1e-8)
def test_paint(): engine = ParticleMeshEngine(pm) code = CodeSegment(engine) s = pm.BoxSize / pm.Nmesh * 0.001 + 0.99 * engine.q / pm.Nmesh # sample all positions. code.get_x(s='s', x='x') code.decompose(x='x', layout='layout') code.paint(x='x', mesh='density', layout='layout') check_grad(code, 'density', 's', init={'s': s}, eps=1e-4, rtol=1e-8)
def test_paintdirect(): engine = ParticleMeshEngine(pm) code = CodeSegment(engine) s = pm.BoxSize / pm.Nmesh * 0.001 + 0.99 * engine.q / pm.Nmesh # sample all positions. m = numpy.ones(len(engine.q)) * 3 code.get_x(s='s', x='x') code.decompose(x='x', layout='layout') code.paintdirect(x='x', mesh='density', layout='layout', mass='m') check_grad(code, 'density', 's', init={'s': s, 'm' : m}, eps=1e-4, rtol=1e-8, atol=1e-11) check_grad(code, 'density', 'm', init={'s': s, 'm' : m}, eps=1e-4, rtol=1e-8, atol=1e-11)
def test_readout(): engine = ParticleMeshEngine(pm) code = CodeSegment(engine) s = pm.BoxSize / pm.Nmesh * 0.001 + 0.99 * engine.q / pm.Nmesh # sample all positions. field = pm.generate_whitenoise(seed=1234, type='real') code.get_x(s='s', x='x') code.decompose(x='x', layout='layout') code.readout(x='x', mesh='density', layout='layout', value='value') check_grad(code, 'value', 'density', init={'density' : field, 's': s}, eps=1e-4, rtol=1e-8) check_grad(code, 'value', 's', init={'density' : field, 's': s}, eps=1e-4, rtol=1e-8)