def test_CCodePrinter():
with warnings.catch_warnings():
warnings.filterwarnings("error", category=SymPyDeprecationWarning)
with raises(SymPyDeprecationWarning):
CCodePrinter()
with warnings.catch_warnings():
warnings.filterwarnings("ignore", category=SymPyDeprecationWarning)
assert CCodePrinter().language == 'C'
def test_ccode_Indexed():
from sympy.tensor import IndexedBase, Idx
from sympy import symbols
s, n, m, o = symbols('s n m o', integer=True)
i, j, k = Idx('i', n), Idx('j', m), Idx('k', o)
x = IndexedBase('x')[j]
A = IndexedBase('A')[i, j]
B = IndexedBase('B')[i, j, k]
with warnings.catch_warnings():
warnings.filterwarnings("ignore", category=SymPyDeprecationWarning)
p = CCodePrinter()
p._not_c = set()
assert p._print_Indexed(x) == 'x[j]'
assert p._print_Indexed(A) == 'A[%s]' % (m * i + j)
assert p._print_Indexed(B) == 'B[%s]' % (i * o * m + j * o + k)
assert p._not_c == set()
A = IndexedBase('A', shape=(5, 3))[i, j]
assert p._print_Indexed(A) == 'A[%s]' % (3 * i + j)
A = IndexedBase('A', shape=(5, 3), strides='F')[i, j]
assert ccode(A) == 'A[%s]' % (i + 5 * j)
A = IndexedBase('A', shape=(29, 29), strides=(1, s), offset=o)[i, j]
assert ccode(A) == 'A[o + s*j + i]'
Abase = IndexedBase('A', strides=(s, m, n), offset=o)
assert ccode(Abase[i, j, k]) == 'A[m*j + n*k + o + s*i]'
assert ccode(Abase[2, 3, k]) == 'A[3*m + n*k + o + 2*s]'
def convert_to(self, language, symbols={}, namespace={}):
'''
Converts expression into a string for the given ``language`` using the
given set of ``symbols``. Replaces each :class:`Symbol` appearing in the
expression with ``sym.read()``, and if the language is C++ or GPU then
uses ``sympy.CCodePrinter().doprint()`` to convert the syntax, e.g.
``x**y`` becomes ``pow(x,y)``.
'''
if language.name == 'python':
return substitute_symbols(self.expr, symbols)
elif language.name == 'c' or language.name == 'gpu':
return substitute_symbols(CCodePrinter().doprint(self.sympy_expr),
symbols)
def test_ccode_Indexed():
from sympy.tensor import IndexedBase, Idx
from sympy import symbols
i, j, k, n, m, o = symbols('i j k n m o', integer=True)
p = CCodePrinter()
p._not_c = set()
x = IndexedBase('x')[Idx(j, n)]
assert p._print_Indexed(x) == 'x[j]'
A = IndexedBase('A')[Idx(i, m), Idx(j, n)]
assert p._print_Indexed(A) == 'A[%s]' % str(j + n * i)
B = IndexedBase('B')[Idx(i, m), Idx(j, n), Idx(k, o)]
assert p._print_Indexed(B) == 'B[%s]' % str(k + i * n * o + j * o)
assert p._not_c == set()
def test_ccode_Indexed():
from sympy.tensor import IndexedBase, Idx
from sympy import symbols
n, m, o = symbols('n m o', integer=True)
i, j, k = Idx('i', n), Idx('j', m), Idx('k', o)
p = CCodePrinter()
p._not_c = set()
x = IndexedBase('x')[j]
assert p._print_Indexed(x) == 'x[j]'
A = IndexedBase('A')[i, j]
assert p._print_Indexed(A) == 'A[%s]' % (m * i + j)
B = IndexedBase('B')[i, j, k]
assert p._print_Indexed(B) == 'B[%s]' % (i * o * m + j * o + k)
assert p._not_c == set()
def __init__(self, **kwargs):
"""Initialize a C code printer.
The printer class, the name of the template, the line continuation
symbol, and the statement ending will be set automatically.
"""
super().__init__(
CCodePrinter(),
lambda base, indices: ''.join(
[base] + ['[{}]'.format(i.index) for i in indices]),
line_cont='\\',
stmt_end=';',
add_filters={
'form_loop_beg': _form_c_loop_beg,
'form_loop_end': _form_c_loop_end,
},
add_globals={'zero_literal': '0.0'},
**kwargs)
def test_ccode_Indexed():
from sympy.tensor import IndexedBase, Idx
from sympy import symbols
n, m, o = symbols('n m o', integer=True)
i, j, k = Idx('i', n), Idx('j', m), Idx('k', o)
x = IndexedBase('x')[j]
A = IndexedBase('A')[i, j]
B = IndexedBase('B')[i, j, k]
with warnings.catch_warnings():
warnings.filterwarnings("ignore", category=SymPyDeprecationWarning)
p = CCodePrinter()
p._not_c = set()
assert p._print_Indexed(x) == 'x[j]'
assert p._print_Indexed(A) == 'A[%s]' % (m * i + j)
assert p._print_Indexed(B) == 'B[%s]' % (i * o * m + j * o + k)
assert p._not_c == set()
def _indent_code(self, codelines):
p = CCodePrinter()
return p.indent_code(codelines)
def test_CCodePrinter():
with warns_deprecated_sympy():
CCodePrinter()
with warns_deprecated_sympy():
assert CCodePrinter().language == 'C'
def convert_to(self, language, symbols={}):
if language.name == 'python':
return self.substitute_symbols(self.expr, symbols)
elif language.name == 'c':
return self.substitute_symbols(
CCodePrinter().doprint(self.sympy_expr), symbols)
