Skip to content
Snippets Groups Projects
Commit 32d636c4 authored by Stephan Seitz's avatar Stephan Seitz
Browse files

Extend test in test_interpolation.py

parent fe94402a
Branches
No related merge requests found
...@@ -7,9 +7,11 @@ ...@@ -7,9 +7,11 @@
""" """
""" """
import itertools
from os.path import dirname, join from os.path import dirname, join
import numpy as np import numpy as np
import pytest
import sympy import sympy
import pycuda.autoinit # NOQA import pycuda.autoinit # NOQA
...@@ -215,19 +217,20 @@ def test_rotate_interpolation_size_change(): ...@@ -215,19 +217,20 @@ def test_rotate_interpolation_size_change():
pyconrad.imshow(out, "small out " + address_mode) pyconrad.imshow(out, "small out " + address_mode)
def test_field_interpolated(): @pytest.mark.parametrize('address_mode, target',
itertools.product(['border', 'wrap', 'clamp', 'mirror'], ['cpu', 'gpu']))
def test_field_interpolated(address_mode, target):
x_f, y_f = pystencils.fields('x,y: float64 [2d]') x_f, y_f = pystencils.fields('x,y: float64 [2d]')
for address_mode in ['border', 'wrap', 'clamp', 'mirror']: assignments = pystencils.AssignmentCollection({
assignments = pystencils.AssignmentCollection({ y_f.center(): x_f.interpolated_access([0.5 * x_ + 2.7, 0.25 * y_ + 7.2], address_mode=address_mode)
y_f.center(): x_f.interpolated_access([0.5 * x_ + 2.7, 0.25 * y_ + 7.2], address_mode=address_mode) })
}) print(assignments)
print(assignments) ast = pystencils.create_kernel(assignments)
ast = pystencils.create_kernel(assignments) print(ast)
print(ast) print(pystencils.show_code(ast))
print(pystencils.show_code(ast)) kernel = ast.compile()
kernel = ast.compile()
out = np.zeros_like(lenna) out = np.zeros_like(lenna)
kernel(x=lenna, y=out) kernel(x=lenna, y=out)
pyconrad.imshow(out, "out " + address_mode) pyconrad.imshow(out, "out " + address_mode)
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment