[4c87de0] | 1 | r""" |
---|
| 2 | py2c |
---|
| 3 | ~~~~ |
---|
| 4 | |
---|
| 5 | Convert simple numeric python code into C code. |
---|
| 6 | |
---|
| 7 | This code is intended to translate direct algorithms for scientific code |
---|
| 8 | (mostly if statements and for loops operating on double precision values) |
---|
| 9 | into C code. Unlike projects like numba, cython, pypy and nuitka, the |
---|
| 10 | :func:`translate` function returns the corresponding C which can then be |
---|
| 11 | compiled with tinycc or sent to the GPU using CUDA or OpenCL. |
---|
| 12 | |
---|
| 13 | There is special handling certain constructs, such as *for i in range* and |
---|
| 14 | small integer powers. |
---|
| 15 | |
---|
| 16 | **TODO: make a nice list of supported constructs*** |
---|
| 17 | |
---|
| 18 | Imports are not supported, but they are at least ignored so that properly |
---|
| 19 | constructed code can be run via python or translated to C without change. |
---|
| 20 | |
---|
| 21 | Most other python constructs are **not** supported: |
---|
| 22 | * classes |
---|
| 23 | * builtin types (dict, set, list) |
---|
| 24 | * exceptions |
---|
| 25 | * with context |
---|
| 26 | * del |
---|
| 27 | * yield |
---|
| 28 | * async |
---|
| 29 | * list slicing |
---|
| 30 | * multiple return values |
---|
| 31 | * "is/is not", "in/not in" conditionals |
---|
| 32 | |
---|
| 33 | There is limited support for list and list comprehensions, so long as they |
---|
| 34 | can be represented by a fixed array whose size is known at compile time, and |
---|
| 35 | they are small enough to be stored on the stack. |
---|
| 36 | |
---|
| 37 | Variables definition in C |
---|
| 38 | ------------------------- |
---|
| 39 | Defining variables within the translate function is a bit of a guess work, |
---|
| 40 | using following rules: |
---|
| 41 | * By default, a variable is a 'double'. |
---|
| 42 | * Variable in a for loop is an int. |
---|
| 43 | * Variable that is references with brackets is an array of doubles. The |
---|
| 44 | variable within the brackets is integer. For example, in the |
---|
| 45 | reference 'var1[var2]', var1 is a double array, and var2 is an integer. |
---|
| 46 | * Assignment to an argument makes that argument an array, and the index |
---|
| 47 | in that assignment is 0. |
---|
| 48 | For example, the following python code:: |
---|
| 49 | def func(arg1, arg2): |
---|
| 50 | arg2 = 17. |
---|
| 51 | is translated to the following C code:: |
---|
| 52 | double func(double arg1) |
---|
| 53 | { |
---|
| 54 | arg2[0] = 17.0; |
---|
| 55 | } |
---|
| 56 | For example, the following python code is translated to the |
---|
| 57 | following C code:: |
---|
| 58 | |
---|
| 59 | def func(arg1, arg2): double func(double arg1) { |
---|
| 60 | arg2 = 17. arg2[0] = 17.0; |
---|
| 61 | } |
---|
| 62 | * All functions are defined as double, even if there is no |
---|
| 63 | return statement. |
---|
| 64 | |
---|
| 65 | Debugging |
---|
| 66 | --------- |
---|
| 67 | |
---|
| 68 | *print* is partially supported using a simple regular expression. This |
---|
| 69 | requires a stylized form. Be sure to use print as a function instead of |
---|
| 70 | the print statement. If you are including substition variables, use the |
---|
| 71 | % string substitution style. Include parentheses around the substitution |
---|
| 72 | tuple, even if there is only one item; do not include the final comma even |
---|
| 73 | if it is a single item (yes, it won't be a tuple, but it makes the regexp |
---|
| 74 | much simpler). Keep the item on a single line. Here are three forms that work:: |
---|
| 75 | |
---|
| 76 | print("x") => printf("x\n"); |
---|
| 77 | print("x %g"%(a)) => printf("x %g\n", a); |
---|
| 78 | print("x %g %g %g"%(a, b, c)) => printf("x %g %g %g\n", a, b, c); |
---|
| 79 | |
---|
| 80 | You can generate *main* using the *if __name__ == "__main__":* construct. |
---|
| 81 | This does a simple substitution with "def main():" before translation and |
---|
| 82 | a substitution with "int main(int argc, double *argv[])" after translation. |
---|
| 83 | The result is that the content of the *if* block becomes the content of *main*. |
---|
| 84 | Along with the print statement, you can run and test a translation standalone |
---|
| 85 | using:: |
---|
| 86 | |
---|
| 87 | python py2c.py source.py |
---|
| 88 | cc source.c |
---|
| 89 | ./a.out |
---|
| 90 | |
---|
| 91 | Known issues |
---|
| 92 | ------------ |
---|
| 93 | The following constructs may cause problems: |
---|
| 94 | |
---|
| 95 | * implicit arrays: possible namespace collision for variable "vec#" |
---|
| 96 | * swap fails: "x,y = y,x" will set x==y |
---|
| 97 | * top-level statements: code outside a function body causes errors |
---|
| 98 | * line number skew: each statement should be tagged with its own #line |
---|
| 99 | to avoid skew as comments are skipped and loop bodies are wrapped with |
---|
| 100 | braces, etc. |
---|
[765d025] | 101 | * doesn't support docstrings |
---|
[4c87de0] | 102 | |
---|
| 103 | References |
---|
| 104 | ---------- |
---|
| 105 | |
---|
| 106 | Based on a variant of codegen.py: |
---|
| 107 | |
---|
| 108 | https://github.com/andreif/codegen |
---|
[c01ed3e] | 109 | :copyright: Copyright 2008 by Armin Ronacher. |
---|
| 110 | :license: BSD. |
---|
| 111 | """ |
---|
[4c87de0] | 112 | |
---|
[6a37819] | 113 | # Update Notes |
---|
| 114 | # ============ |
---|
[d5014e4] | 115 | # 2017-11-22, OE: Each 'visit_*' method is to build a C statement string. It |
---|
| 116 | # shold insert 4 blanks per indentation level. The 'body' |
---|
| 117 | # method will combine all the strings, by adding the |
---|
| 118 | # 'current_statement' to the c_proc string list |
---|
| 119 | # 2017-11-22, OE: variables, argument definition implemented. Note: An |
---|
| 120 | # argument is considered an array if it is the target of an |
---|
| 121 | # assignment. In that case it is translated to <var>[0] |
---|
| 122 | # 2017-11-27, OE: 'pow' basicly working |
---|
| 123 | # 2017-12-07, OE: Multiple assignment: a1,a2,...,an=b1,b2,...bn implemented |
---|
| 124 | # 2017-12-07, OE: Power function, including special cases of |
---|
[6a37819] | 125 | # square(x)(pow(x,2)) and cube(x)(pow(x,3)), implemented in |
---|
| 126 | # translate_power, called from visit_BinOp |
---|
[d5014e4] | 127 | # 2017-12-07, OE: Translation of integer division, '\\' in python, implemented |
---|
[6a37819] | 128 | # in translate_integer_divide, called from visit_BinOp |
---|
[d5014e4] | 129 | # 2017-12-07, OE: C variable definition handled in 'define_c_vars' |
---|
[6a37819] | 130 | # : Python integer division, '//', translated to C in |
---|
| 131 | # 'translate_integer_divide' |
---|
[d5014e4] | 132 | # 2017-12-15, OE: Precedence maintained by writing opening and closing |
---|
[6a37819] | 133 | # parenthesesm '(',')', in procedure 'visit_BinOp'. |
---|
[d5014e4] | 134 | # 2017-12-18, OE: Added call to 'add_current_line()' at the beginning |
---|
[6a37819] | 135 | # of visit_Return |
---|
| 136 | # 2018-01-03, PK: Update interface for use in sasmodels |
---|
| 137 | # 2018-01-03, PK: support "expr if cond else expr" syntax |
---|
| 138 | # 2018-01-03, PK: x//y => (int)((x)/(y)) and x/y => ((double)(x)/(double)(y)) |
---|
| 139 | # 2018-01-03, PK: True/False => true/false |
---|
| 140 | # 2018-01-03, PK: f(x) was introducing an extra semicolon |
---|
| 141 | # 2018-01-03, PK: simplistic print function, for debugging |
---|
| 142 | # 2018-01-03, PK: while expr: ... => while (expr) { ... } |
---|
[d5014e4] | 143 | # 2018-01-04, OE: Fixed bug in 'visit_If': visiting node.orelse in case else exists. |
---|
[6a37819] | 144 | |
---|
| 145 | from __future__ import print_function |
---|
[4c87de0] | 146 | |
---|
[71779b2] | 147 | import sys |
---|
[6a37819] | 148 | import ast |
---|
[71779b2] | 149 | from ast import NodeVisitor |
---|
[6f91c91] | 150 | from inspect import currentframe, getframeinfo |
---|
| 151 | |
---|
[5dd7cfb] | 152 | try: # for debugging, astor lets us print out the node as python |
---|
| 153 | import astor |
---|
| 154 | except ImportError: |
---|
| 155 | pass |
---|
[71779b2] | 156 | |
---|
| 157 | BINOP_SYMBOLS = {} |
---|
| 158 | BINOP_SYMBOLS[ast.Add] = '+' |
---|
| 159 | BINOP_SYMBOLS[ast.Sub] = '-' |
---|
| 160 | BINOP_SYMBOLS[ast.Mult] = '*' |
---|
| 161 | BINOP_SYMBOLS[ast.Div] = '/' |
---|
| 162 | BINOP_SYMBOLS[ast.Mod] = '%' |
---|
| 163 | BINOP_SYMBOLS[ast.Pow] = '**' |
---|
| 164 | BINOP_SYMBOLS[ast.LShift] = '<<' |
---|
| 165 | BINOP_SYMBOLS[ast.RShift] = '>>' |
---|
| 166 | BINOP_SYMBOLS[ast.BitOr] = '|' |
---|
| 167 | BINOP_SYMBOLS[ast.BitXor] = '^' |
---|
| 168 | BINOP_SYMBOLS[ast.BitAnd] = '&' |
---|
| 169 | BINOP_SYMBOLS[ast.FloorDiv] = '//' |
---|
| 170 | |
---|
| 171 | BOOLOP_SYMBOLS = {} |
---|
| 172 | BOOLOP_SYMBOLS[ast.And] = '&&' |
---|
[7b1dcf9] | 173 | BOOLOP_SYMBOLS[ast.Or] = '||' |
---|
[71779b2] | 174 | |
---|
| 175 | CMPOP_SYMBOLS = {} |
---|
[7b1dcf9] | 176 | CMPOP_SYMBOLS[ast.Eq] = '==' |
---|
[71779b2] | 177 | CMPOP_SYMBOLS[ast.NotEq] = '!=' |
---|
| 178 | CMPOP_SYMBOLS[ast.Lt] = '<' |
---|
| 179 | CMPOP_SYMBOLS[ast.LtE] = '<=' |
---|
| 180 | CMPOP_SYMBOLS[ast.Gt] = '>' |
---|
| 181 | CMPOP_SYMBOLS[ast.GtE] = '>=' |
---|
| 182 | CMPOP_SYMBOLS[ast.Is] = 'is' |
---|
| 183 | CMPOP_SYMBOLS[ast.IsNot] = 'is not' |
---|
| 184 | CMPOP_SYMBOLS[ast.In] = 'in' |
---|
| 185 | CMPOP_SYMBOLS[ast.NotIn] = 'not in' |
---|
| 186 | |
---|
| 187 | UNARYOP_SYMBOLS = {} |
---|
| 188 | UNARYOP_SYMBOLS[ast.Invert] = '~' |
---|
| 189 | UNARYOP_SYMBOLS[ast.Not] = 'not' |
---|
| 190 | UNARYOP_SYMBOLS[ast.UAdd] = '+' |
---|
| 191 | UNARYOP_SYMBOLS[ast.USub] = '-' |
---|
| 192 | |
---|
| 193 | |
---|
[d5014e4] | 194 | # TODO: should not allow eval of arbitrary python |
---|
[3f9db6e] | 195 | def isevaluable(s): |
---|
| 196 | try: |
---|
| 197 | eval(s) |
---|
| 198 | return True |
---|
[c01ed3e] | 199 | except Exception: |
---|
[3f9db6e] | 200 | return False |
---|
[fa74acf] | 201 | |
---|
[5dd7cfb] | 202 | def render_expression(tree): |
---|
| 203 | generator = SourceGenerator() |
---|
| 204 | generator.visit(tree) |
---|
| 205 | c_code = "".join(generator.current_statement) |
---|
| 206 | return c_code |
---|
| 207 | |
---|
[71779b2] | 208 | class SourceGenerator(NodeVisitor): |
---|
| 209 | """This visitor is able to transform a well formed syntax tree into python |
---|
| 210 | sourcecode. For more details have a look at the docstring of the |
---|
| 211 | `node_to_source` function. |
---|
| 212 | """ |
---|
| 213 | |
---|
[5dd7cfb] | 214 | def __init__(self, indent_with=" ", constants=None, fname=None, lineno=0): |
---|
[71779b2] | 215 | self.indent_with = indent_with |
---|
| 216 | self.indentation = 0 |
---|
[c01ed3e] | 217 | |
---|
[7b1dcf9] | 218 | # for C |
---|
[c01ed3e] | 219 | self.c_proc = [] |
---|
[71779b2] | 220 | self.signature_line = 0 |
---|
| 221 | self.arguments = [] |
---|
[c01ed3e] | 222 | self.current_function = "" |
---|
| 223 | self.fname = fname |
---|
| 224 | self.lineno_offset = lineno |
---|
[71779b2] | 225 | self.warnings = [] |
---|
| 226 | self.current_statement = "" |
---|
[c01ed3e] | 227 | # TODO: use set rather than list for c_vars, ... |
---|
| 228 | self.c_vars = [] |
---|
| 229 | self.c_int_vars = [] |
---|
| 230 | self.c_pointers = [] |
---|
| 231 | self.c_dcl_pointers = [] |
---|
| 232 | self.c_functions = [] |
---|
| 233 | self.c_vectors = [] |
---|
| 234 | self.c_constants = constants if constants is not None else {} |
---|
[4c87de0] | 235 | self.in_expr = False |
---|
[c01ed3e] | 236 | self.in_subref = False |
---|
| 237 | self.in_subscript = False |
---|
| 238 | self.tuples = [] |
---|
[71779b2] | 239 | self.required_functions = [] |
---|
| 240 | self.visited_args = False |
---|
| 241 | |
---|
[c01ed3e] | 242 | def write_c(self, statement): |
---|
| 243 | # TODO: build up as a list rather than adding to string |
---|
| 244 | self.current_statement += statement |
---|
[fa74acf] | 245 | |
---|
[5dd7cfb] | 246 | def write_python(self, x): |
---|
| 247 | raise NotImplementedError("shouldn't be trying to write pythnon") |
---|
| 248 | |
---|
[c01ed3e] | 249 | def add_c_line(self, line): |
---|
| 250 | indentation = self.indent_with * self.indentation |
---|
| 251 | self.c_proc.append("".join((indentation, line, "\n"))) |
---|
[71779b2] | 252 | |
---|
[fa74acf] | 253 | def add_current_line(self): |
---|
[7b1dcf9] | 254 | if self.current_statement: |
---|
[fa74acf] | 255 | self.add_c_line(self.current_statement) |
---|
[71779b2] | 256 | self.current_statement = '' |
---|
| 257 | |
---|
[c01ed3e] | 258 | def add_unique_var(self, new_var): |
---|
| 259 | if new_var not in self.c_vars: |
---|
| 260 | self.c_vars.append(str(new_var)) |
---|
[71779b2] | 261 | |
---|
[c01ed3e] | 262 | def write_sincos(self, node): |
---|
[71779b2] | 263 | angle = str(node.args[0].id) |
---|
[fa74acf] | 264 | self.write_c(node.args[1].id + " = sin(" + angle + ");") |
---|
[71779b2] | 265 | self.add_current_line() |
---|
[fa74acf] | 266 | self.write_c(node.args[2].id + " = cos(" + angle + ");") |
---|
[71779b2] | 267 | self.add_current_line() |
---|
| 268 | for arg in node.args: |
---|
[c01ed3e] | 269 | self.add_unique_var(arg.id) |
---|
[71779b2] | 270 | |
---|
[5dd7cfb] | 271 | def track_lineno(self, node): |
---|
| 272 | #print("newline", node, [s for s in dir(node) if not s.startswith('_')]) |
---|
| 273 | if hasattr(node, 'lineno'): |
---|
| 274 | line = '#line %d "%s"\n' % (node.lineno+self.lineno_offset-1, self.fname) |
---|
| 275 | self.c_proc.append(line) |
---|
[71779b2] | 276 | |
---|
| 277 | def body(self, statements): |
---|
[7b1dcf9] | 278 | if self.current_statement: |
---|
[fa74acf] | 279 | self.add_current_line() |
---|
[71779b2] | 280 | self.new_line = True |
---|
| 281 | self.indentation += 1 |
---|
| 282 | for stmt in statements: |
---|
[7b1dcf9] | 283 | #if hasattr(stmt, 'targets') and hasattr(stmt.targets[0], 'id'): |
---|
| 284 | # target_name = stmt.targets[0].id # target name needed for debug only |
---|
[71779b2] | 285 | self.visit(stmt) |
---|
[fa74acf] | 286 | self.add_current_line() # just for breaking point. to be deleted. |
---|
[71779b2] | 287 | self.indentation -= 1 |
---|
| 288 | |
---|
| 289 | def body_or_else(self, node): |
---|
| 290 | self.body(node.body) |
---|
| 291 | if node.orelse: |
---|
[15be191] | 292 | self.unsupported(node, "for...else/while...else not supported") |
---|
| 293 | |
---|
[5dd7cfb] | 294 | self.track_lineno(node) |
---|
[71779b2] | 295 | self.write_c('else:') |
---|
| 296 | self.body(node.orelse) |
---|
| 297 | |
---|
| 298 | def signature(self, node): |
---|
| 299 | want_comma = [] |
---|
| 300 | def write_comma(): |
---|
| 301 | if want_comma: |
---|
| 302 | self.write_c(', ') |
---|
| 303 | else: |
---|
| 304 | want_comma.append(True) |
---|
[7b1dcf9] | 305 | |
---|
| 306 | # for C |
---|
[71779b2] | 307 | for arg in node.args: |
---|
[1ddb794] | 308 | # CRUFT: 2.7 uses arg.id, 3.x uses arg.arg |
---|
| 309 | try: |
---|
| 310 | arg_name = arg.arg |
---|
| 311 | except AttributeError: |
---|
| 312 | arg_name = arg.id |
---|
| 313 | self.arguments.append(arg_name) |
---|
[71779b2] | 314 | |
---|
[fa74acf] | 315 | padding = [None] *(len(node.args) - len(node.defaults)) |
---|
[71779b2] | 316 | for arg, default in zip(node.args, padding + node.defaults): |
---|
| 317 | if default is not None: |
---|
[1ddb794] | 318 | # CRUFT: 2.7 uses arg.id, 3.x uses arg.arg |
---|
| 319 | try: |
---|
| 320 | arg_name = arg.arg |
---|
| 321 | except AttributeError: |
---|
| 322 | arg_name = arg.id |
---|
[d5014e4] | 323 | w_str = ("C does not support default parameters: %s=%s" |
---|
[15be191] | 324 | % (arg_name, str(default.n))) |
---|
[fa74acf] | 325 | self.warnings.append(w_str) |
---|
[71779b2] | 326 | |
---|
| 327 | def decorators(self, node): |
---|
[5dd7cfb] | 328 | if node.decorator_list: |
---|
| 329 | self.unsupported(node.decorator_list[0]) |
---|
[71779b2] | 330 | for decorator in node.decorator_list: |
---|
[5dd7cfb] | 331 | self.trac_lineno(decorator) |
---|
[71779b2] | 332 | self.write_python('@') |
---|
| 333 | self.visit(decorator) |
---|
| 334 | |
---|
| 335 | # Statements |
---|
| 336 | |
---|
| 337 | def visit_Assert(self, node): |
---|
[5dd7cfb] | 338 | self.unsupported(node) |
---|
| 339 | |
---|
| 340 | self.track_lineno(node) |
---|
[71779b2] | 341 | self.write_c('assert ') |
---|
| 342 | self.visit(node.test) |
---|
| 343 | if node.msg is not None: |
---|
[fa74acf] | 344 | self.write_python(', ') |
---|
| 345 | self.visit(node.msg) |
---|
[71779b2] | 346 | |
---|
[c01ed3e] | 347 | def define_c_vars(self, target): |
---|
[7b1dcf9] | 348 | if hasattr(target, 'id'): |
---|
| 349 | # a variable is considered an array if it apears in the agrument list |
---|
| 350 | # and being assigned to. For example, the variable p in the following |
---|
| 351 | # sniplet is a pointer, while q is not |
---|
| 352 | # def somefunc(p, q): |
---|
| 353 | # p = q + 1 |
---|
| 354 | # return |
---|
| 355 | # |
---|
[c01ed3e] | 356 | if target.id not in self.c_vars: |
---|
[7b1dcf9] | 357 | if target.id in self.arguments: |
---|
[fa74acf] | 358 | idx = self.arguments.index(target.id) |
---|
[71779b2] | 359 | new_target = self.arguments[idx] + "[0]" |
---|
[c01ed3e] | 360 | if new_target not in self.c_pointers: |
---|
[71779b2] | 361 | target.id = new_target |
---|
[c01ed3e] | 362 | self.c_pointers.append(self.arguments[idx]) |
---|
[71779b2] | 363 | else: |
---|
[c01ed3e] | 364 | self.c_vars.append(target.id) |
---|
[71779b2] | 365 | |
---|
[fa74acf] | 366 | def add_semi_colon(self): |
---|
[4c87de0] | 367 | #semi_pos = self.current_statement.find(';') |
---|
| 368 | #if semi_pos >= 0: |
---|
| 369 | # self.current_statement = self.current_statement.replace(';', '') |
---|
[98a4f14] | 370 | self.write_c(';') |
---|
[937afef] | 371 | |
---|
[71779b2] | 372 | def visit_Assign(self, node): |
---|
[fa74acf] | 373 | self.add_current_line() |
---|
[5dd7cfb] | 374 | self.track_lineno(node) |
---|
[4c87de0] | 375 | self.in_expr = True |
---|
[71779b2] | 376 | for idx, target in enumerate(node.targets): # multi assign, as in 'a = b = c = 7' |
---|
| 377 | if idx: |
---|
| 378 | self.write_c(' = ') |
---|
[c01ed3e] | 379 | self.define_c_vars(target) |
---|
[71779b2] | 380 | self.visit(target) |
---|
[c01ed3e] | 381 | # Capture assigned tuple names, if any |
---|
| 382 | targets = self.tuples[:] |
---|
| 383 | del self.tuples[:] |
---|
[71779b2] | 384 | self.write_c(' = ') |
---|
| 385 | self.visited_args = False |
---|
| 386 | self.visit(node.value) |
---|
[fa74acf] | 387 | self.add_semi_colon() |
---|
| 388 | self.add_current_line() |
---|
[c01ed3e] | 389 | # Assign tuples to tuples, if any |
---|
| 390 | # TODO: doesn't handle swap: a,b = b,a |
---|
| 391 | for target, item in zip(targets, self.tuples): |
---|
| 392 | self.visit(target) |
---|
[71779b2] | 393 | self.write_c(' = ') |
---|
| 394 | self.visit(item) |
---|
[fa74acf] | 395 | self.add_semi_colon() |
---|
| 396 | self.add_current_line() |
---|
[5dd7cfb] | 397 | #if self.is_sequence and not self.visited_args: |
---|
| 398 | # for target in node.targets: |
---|
| 399 | # if hasattr(target, 'id'): |
---|
| 400 | # if target.id in self.c_vars and target.id not in self.c_dcl_pointers: |
---|
| 401 | # if target.id not in self.c_dcl_pointers: |
---|
| 402 | # self.c_dcl_pointers.append(target.id) |
---|
| 403 | # if target.id in self.c_vars: |
---|
| 404 | # self.c_vars.remove(target.id) |
---|
[71779b2] | 405 | self.current_statement = '' |
---|
[4c87de0] | 406 | self.in_expr = False |
---|
[71779b2] | 407 | |
---|
| 408 | def visit_AugAssign(self, node): |
---|
[c01ed3e] | 409 | if node.target.id not in self.c_vars: |
---|
[7b1dcf9] | 410 | if node.target.id not in self.arguments: |
---|
[c01ed3e] | 411 | self.c_vars.append(node.target.id) |
---|
[4c87de0] | 412 | self.in_expr = True |
---|
[71779b2] | 413 | self.visit(node.target) |
---|
| 414 | self.write_c(' ' + BINOP_SYMBOLS[type(node.op)] + '= ') |
---|
| 415 | self.visit(node.value) |
---|
[fa74acf] | 416 | self.add_semi_colon() |
---|
[4c87de0] | 417 | self.in_expr = False |
---|
[fa74acf] | 418 | self.add_current_line() |
---|
[71779b2] | 419 | |
---|
| 420 | def visit_ImportFrom(self, node): |
---|
[c01ed3e] | 421 | return # import ignored |
---|
[5dd7cfb] | 422 | self.track_lineno(node) |
---|
[fa74acf] | 423 | self.write_python('from %s%s import ' %('.' * node.level, node.module)) |
---|
[71779b2] | 424 | for idx, item in enumerate(node.names): |
---|
| 425 | if idx: |
---|
| 426 | self.write_python(', ') |
---|
| 427 | self.write_python(item) |
---|
| 428 | |
---|
| 429 | def visit_Import(self, node): |
---|
[c01ed3e] | 430 | return # import ignored |
---|
[5dd7cfb] | 431 | self.track_lineno(node) |
---|
[71779b2] | 432 | for item in node.names: |
---|
| 433 | self.write_python('import ') |
---|
| 434 | self.visit(item) |
---|
| 435 | |
---|
| 436 | def visit_Expr(self, node): |
---|
[4c87de0] | 437 | #self.in_expr = True |
---|
[5dd7cfb] | 438 | #self.track_lineno(node) |
---|
[71779b2] | 439 | self.generic_visit(node) |
---|
[4c87de0] | 440 | #self.in_expr = False |
---|
[71779b2] | 441 | |
---|
[c01ed3e] | 442 | def write_c_pointers(self, start_var): |
---|
| 443 | if self.c_dcl_pointers: |
---|
[7b1dcf9] | 444 | var_list = [] |
---|
[c01ed3e] | 445 | for c_ptr in self.c_dcl_pointers: |
---|
[7b1dcf9] | 446 | if c_ptr not in self.arguments: |
---|
| 447 | var_list.append("*" + c_ptr) |
---|
[c01ed3e] | 448 | if c_ptr in self.c_vars: |
---|
| 449 | self.c_vars.remove(c_ptr) |
---|
[7b1dcf9] | 450 | if var_list: |
---|
| 451 | c_dcl = " double " + ", ".join(var_list) + ";\n" |
---|
| 452 | self.c_proc.insert(start_var, c_dcl) |
---|
[71c5f4d] | 453 | start_var += 1 |
---|
| 454 | return start_var |
---|
| 455 | |
---|
[c01ed3e] | 456 | def insert_c_vars(self, start_var): |
---|
| 457 | have_decls = False |
---|
| 458 | start_var = self.write_c_pointers(start_var) |
---|
| 459 | if self.c_int_vars: |
---|
| 460 | for var in self.c_int_vars: |
---|
| 461 | if var in self.c_vars: |
---|
| 462 | self.c_vars.remove(var) |
---|
| 463 | decls = ", ".join(self.c_int_vars) |
---|
| 464 | self.c_proc.insert(start_var, " int " + decls + ";\n") |
---|
| 465 | have_decls = True |
---|
[71779b2] | 466 | start_var += 1 |
---|
[fa74acf] | 467 | |
---|
[c01ed3e] | 468 | if self.c_vars: |
---|
| 469 | decls = ", ".join(self.c_vars) |
---|
| 470 | self.c_proc.insert(start_var, " double " + decls + ";\n") |
---|
| 471 | have_decls = True |
---|
[937afef] | 472 | start_var += 1 |
---|
[7b1dcf9] | 473 | |
---|
[c01ed3e] | 474 | if self.c_vectors: |
---|
| 475 | for vec_number, vec_value in enumerate(self.c_vectors): |
---|
| 476 | name = "vec" + str(vec_number + 1) |
---|
| 477 | decl = " double " + name + "[] = {" + vec_value + "};" |
---|
| 478 | self.c_proc.insert(start_var, decl + "\n") |
---|
[71779b2] | 479 | start_var += 1 |
---|
[7b1dcf9] | 480 | |
---|
[c01ed3e] | 481 | del self.c_vars[:] |
---|
| 482 | del self.c_int_vars[:] |
---|
| 483 | del self.c_vectors[:] |
---|
| 484 | del self.c_pointers[:] |
---|
| 485 | del self.c_dcl_pointers[:] |
---|
| 486 | if have_decls: |
---|
[fa74acf] | 487 | self.c_proc.insert(start_var, "\n") |
---|
[71779b2] | 488 | |
---|
[c01ed3e] | 489 | def insert_signature(self): |
---|
[7b1dcf9] | 490 | arg_decls = [] |
---|
| 491 | for arg in self.arguments: |
---|
| 492 | decl = "double " + arg |
---|
[c01ed3e] | 493 | if arg in self.c_pointers: |
---|
[7b1dcf9] | 494 | decl += "[]" |
---|
| 495 | arg_decls.append(decl) |
---|
| 496 | args_str = ", ".join(arg_decls) |
---|
[c01ed3e] | 497 | method_sig = 'double ' + self.current_function + '(' + args_str + ")" |
---|
[7b1dcf9] | 498 | if self.signature_line >= 0: |
---|
[c01ed3e] | 499 | self.c_proc.insert(self.signature_line, method_sig) |
---|
[71779b2] | 500 | |
---|
| 501 | def visit_FunctionDef(self, node): |
---|
[c01ed3e] | 502 | if self.current_function: |
---|
| 503 | self.unsupported(node, "function within a function") |
---|
| 504 | self.current_function = node.name |
---|
| 505 | |
---|
[d5014e4] | 506 | # remember the location of the next warning that will be inserted |
---|
| 507 | # so that we can stuff the function name ahead of the warning list |
---|
| 508 | # if any warnings are generated by the function. |
---|
[4339764] | 509 | warning_index = len(self.warnings) |
---|
[d5014e4] | 510 | |
---|
[71779b2] | 511 | self.decorators(node) |
---|
[5dd7cfb] | 512 | self.track_lineno(node) |
---|
[71779b2] | 513 | self.arguments = [] |
---|
| 514 | self.visit(node.args) |
---|
[7b1dcf9] | 515 | # for C |
---|
[71779b2] | 516 | self.signature_line = len(self.c_proc) |
---|
| 517 | self.add_c_line("\n{") |
---|
| 518 | start_vars = len(self.c_proc) + 1 |
---|
| 519 | self.body(node.body) |
---|
| 520 | self.add_c_line("}\n") |
---|
[c01ed3e] | 521 | self.insert_signature() |
---|
| 522 | self.insert_c_vars(start_vars) |
---|
[15be191] | 523 | del self.c_pointers[:] |
---|
[c01ed3e] | 524 | self.current_function = "" |
---|
[d5014e4] | 525 | if warning_index != len(self.warnings): |
---|
[4339764] | 526 | self.warnings.insert(warning_index, "Warning in function '" + node.name + "':") |
---|
[71779b2] | 527 | |
---|
| 528 | def visit_ClassDef(self, node): |
---|
[5dd7cfb] | 529 | self.unsupported(node) |
---|
| 530 | |
---|
[71779b2] | 531 | have_args = [] |
---|
| 532 | def paren_or_comma(): |
---|
| 533 | if have_args: |
---|
| 534 | self.write_python(', ') |
---|
| 535 | else: |
---|
| 536 | have_args.append(True) |
---|
| 537 | self.write_python('(') |
---|
| 538 | |
---|
| 539 | self.decorators(node) |
---|
[5dd7cfb] | 540 | self.track_lineno(node) |
---|
[71779b2] | 541 | self.write_python('class %s' % node.name) |
---|
| 542 | for base in node.bases: |
---|
| 543 | paren_or_comma() |
---|
| 544 | self.visit(base) |
---|
[7b1dcf9] | 545 | # CRUFT: python 2.6 does not have "keywords" attribute |
---|
[71779b2] | 546 | if hasattr(node, 'keywords'): |
---|
| 547 | for keyword in node.keywords: |
---|
| 548 | paren_or_comma() |
---|
| 549 | self.write_python(keyword.arg + '=') |
---|
| 550 | self.visit(keyword.value) |
---|
| 551 | if node.starargs is not None: |
---|
| 552 | paren_or_comma() |
---|
| 553 | self.write_python('*') |
---|
| 554 | self.visit(node.starargs) |
---|
| 555 | if node.kwargs is not None: |
---|
| 556 | paren_or_comma() |
---|
| 557 | self.write_python('**') |
---|
| 558 | self.visit(node.kwargs) |
---|
| 559 | self.write_python(have_args and '):' or ':') |
---|
| 560 | self.body(node.body) |
---|
| 561 | |
---|
| 562 | def visit_If(self, node): |
---|
[4c87de0] | 563 | |
---|
[5dd7cfb] | 564 | self.track_lineno(node) |
---|
[71779b2] | 565 | self.write_c('if ') |
---|
[4c87de0] | 566 | self.in_expr = True |
---|
[71779b2] | 567 | self.visit(node.test) |
---|
[4c87de0] | 568 | self.in_expr = False |
---|
[71779b2] | 569 | self.write_c(' {') |
---|
| 570 | self.body(node.body) |
---|
| 571 | self.add_c_line('}') |
---|
| 572 | while True: |
---|
| 573 | else_ = node.orelse |
---|
| 574 | if len(else_) == 0: |
---|
| 575 | break |
---|
[7b1dcf9] | 576 | #elif hasattr(else_, 'orelse'): |
---|
[71779b2] | 577 | elif len(else_) == 1 and isinstance(else_[0], ast.If): |
---|
| 578 | node = else_[0] |
---|
[5dd7cfb] | 579 | self.track_lineno(node) |
---|
[71779b2] | 580 | self.write_c('else if ') |
---|
[4c87de0] | 581 | self.in_expr = True |
---|
[71779b2] | 582 | self.visit(node.test) |
---|
[4c87de0] | 583 | self.in_expr = False |
---|
[71779b2] | 584 | self.write_c(' {') |
---|
| 585 | self.body(node.body) |
---|
[fa74acf] | 586 | self.add_current_line() |
---|
[71779b2] | 587 | self.add_c_line('}') |
---|
[7b1dcf9] | 588 | #break |
---|
[71779b2] | 589 | else: |
---|
[5dd7cfb] | 590 | self.track_lineno(else_) |
---|
[71779b2] | 591 | self.write_c('else {') |
---|
[d5014e4] | 592 | self.body(else_) |
---|
[71779b2] | 593 | self.add_c_line('}') |
---|
| 594 | break |
---|
| 595 | |
---|
[c01ed3e] | 596 | def get_for_range(self, node): |
---|
[937afef] | 597 | stop = "" |
---|
| 598 | start = '0' |
---|
| 599 | step = '1' |
---|
| 600 | for_args = [] |
---|
| 601 | temp_statement = self.current_statement |
---|
| 602 | self.current_statement = '' |
---|
| 603 | for arg in node.iter.args: |
---|
| 604 | self.visit(arg) |
---|
| 605 | for_args.append(self.current_statement) |
---|
[71779b2] | 606 | self.current_statement = '' |
---|
[937afef] | 607 | self.current_statement = temp_statement |
---|
[7b1dcf9] | 608 | if len(for_args) == 1: |
---|
[937afef] | 609 | stop = for_args[0] |
---|
[7b1dcf9] | 610 | elif len(for_args) == 2: |
---|
[937afef] | 611 | start = for_args[0] |
---|
| 612 | stop = for_args[1] |
---|
[7b1dcf9] | 613 | elif len(for_args) == 3: |
---|
[937afef] | 614 | start = for_args[0] |
---|
| 615 | stop = for_args[1] |
---|
| 616 | start = for_args[2] |
---|
| 617 | else: |
---|
| 618 | raise("Ilegal for loop parameters") |
---|
[7b1dcf9] | 619 | return start, stop, step |
---|
[937afef] | 620 | |
---|
[d5014e4] | 621 | def add_c_int_var(self, name): |
---|
| 622 | if name not in self.c_int_vars: |
---|
| 623 | self.c_int_vars.append(name) |
---|
[4339764] | 624 | |
---|
[937afef] | 625 | def visit_For(self, node): |
---|
[7b1dcf9] | 626 | # node: for iterator is stored in node.target. |
---|
| 627 | # Iterator name is in node.target.id. |
---|
[937afef] | 628 | self.add_current_line() |
---|
[71779b2] | 629 | fForDone = False |
---|
[937afef] | 630 | self.current_statement = '' |
---|
[7b1dcf9] | 631 | if hasattr(node.iter, 'func'): |
---|
| 632 | if hasattr(node.iter.func, 'id'): |
---|
| 633 | if node.iter.func.id == 'range': |
---|
[937afef] | 634 | self.visit(node.target) |
---|
| 635 | iterator = self.current_statement |
---|
| 636 | self.current_statement = '' |
---|
[d5014e4] | 637 | self.add_c_int_var(iterator) |
---|
[c01ed3e] | 638 | start, stop, step = self.get_for_range(node) |
---|
[4c87de0] | 639 | self.write_c("for (" + iterator + "=" + str(start) + |
---|
[7b1dcf9] | 640 | " ; " + iterator + " < " + str(stop) + |
---|
| 641 | " ; " + iterator + " += " + str(step) + ") {") |
---|
[71779b2] | 642 | self.body_or_else(node) |
---|
[fa74acf] | 643 | self.write_c("}") |
---|
[71779b2] | 644 | fForDone = True |
---|
[7b1dcf9] | 645 | if not fForDone: |
---|
[c01ed3e] | 646 | # Generate the statement that is causing the error |
---|
[71779b2] | 647 | self.current_statement = '' |
---|
| 648 | self.write_c('for ') |
---|
| 649 | self.visit(node.target) |
---|
| 650 | self.write_c(' in ') |
---|
| 651 | self.visit(node.iter) |
---|
| 652 | self.write_c(':') |
---|
[c01ed3e] | 653 | # report the error |
---|
[15be191] | 654 | self.unsupported(node, "unsupported " + self.current_statement) |
---|
[71779b2] | 655 | |
---|
| 656 | def visit_While(self, node): |
---|
[5dd7cfb] | 657 | self.track_lineno(node) |
---|
[71779b2] | 658 | self.write_c('while ') |
---|
| 659 | self.visit(node.test) |
---|
[4c87de0] | 660 | self.write_c(' {') |
---|
[71779b2] | 661 | self.body_or_else(node) |
---|
[4c87de0] | 662 | self.write_c('}') |
---|
| 663 | self.add_current_line() |
---|
[71779b2] | 664 | |
---|
| 665 | def visit_With(self, node): |
---|
[c01ed3e] | 666 | self.unsupported(node) |
---|
[15be191] | 667 | |
---|
[5dd7cfb] | 668 | self.track_lineno(node) |
---|
[71779b2] | 669 | self.write_python('with ') |
---|
| 670 | self.visit(node.context_expr) |
---|
| 671 | if node.optional_vars is not None: |
---|
| 672 | self.write_python(' as ') |
---|
| 673 | self.visit(node.optional_vars) |
---|
| 674 | self.write_python(':') |
---|
| 675 | self.body(node.body) |
---|
| 676 | |
---|
| 677 | def visit_Pass(self, node): |
---|
[5dd7cfb] | 678 | #self.track_lineno(node) |
---|
[c01ed3e] | 679 | #self.write_python('pass') |
---|
[5dd7cfb] | 680 | pass |
---|
[71779b2] | 681 | |
---|
| 682 | def visit_Print(self, node): |
---|
[c01ed3e] | 683 | self.unsupported(node) |
---|
[15be191] | 684 | |
---|
[7b1dcf9] | 685 | # CRUFT: python 2.6 only |
---|
[5dd7cfb] | 686 | self.track_lineno(node) |
---|
[71779b2] | 687 | self.write_c('print ') |
---|
| 688 | want_comma = False |
---|
| 689 | if node.dest is not None: |
---|
| 690 | self.write_c(' >> ') |
---|
| 691 | self.visit(node.dest) |
---|
| 692 | want_comma = True |
---|
| 693 | for value in node.values: |
---|
| 694 | if want_comma: |
---|
| 695 | self.write_c(', ') |
---|
| 696 | self.visit(value) |
---|
| 697 | want_comma = True |
---|
| 698 | if not node.nl: |
---|
| 699 | self.write_c(',') |
---|
| 700 | |
---|
| 701 | def visit_Delete(self, node): |
---|
[c01ed3e] | 702 | self.unsupported(node) |
---|
[15be191] | 703 | |
---|
[5dd7cfb] | 704 | self.track_lineno(node) |
---|
[71779b2] | 705 | self.write_python('del ') |
---|
| 706 | for idx, target in enumerate(node): |
---|
| 707 | if idx: |
---|
| 708 | self.write_python(', ') |
---|
| 709 | self.visit(target) |
---|
| 710 | |
---|
| 711 | def visit_TryExcept(self, node): |
---|
[c01ed3e] | 712 | self.unsupported(node) |
---|
[15be191] | 713 | |
---|
[5dd7cfb] | 714 | self.track_linno(node) |
---|
[71779b2] | 715 | self.write_python('try:') |
---|
| 716 | self.body(node.body) |
---|
| 717 | for handler in node.handlers: |
---|
| 718 | self.visit(handler) |
---|
| 719 | |
---|
| 720 | def visit_TryFinally(self, node): |
---|
[c01ed3e] | 721 | self.unsupported(node) |
---|
[15be191] | 722 | |
---|
[5dd7cfb] | 723 | self.track_lineno(node) |
---|
[71779b2] | 724 | self.write_python('try:') |
---|
| 725 | self.body(node.body) |
---|
[5dd7cfb] | 726 | self.track_lineno(node) |
---|
[71779b2] | 727 | self.write_python('finally:') |
---|
| 728 | self.body(node.finalbody) |
---|
| 729 | |
---|
| 730 | def visit_Global(self, node): |
---|
[c01ed3e] | 731 | self.unsupported(node) |
---|
[15be191] | 732 | |
---|
[5dd7cfb] | 733 | self.track_lineno(node) |
---|
[71779b2] | 734 | self.write_python('global ' + ', '.join(node.names)) |
---|
| 735 | |
---|
| 736 | def visit_Nonlocal(self, node): |
---|
[5dd7cfb] | 737 | self.track_lineno(node) |
---|
[71779b2] | 738 | self.write_python('nonlocal ' + ', '.join(node.names)) |
---|
| 739 | |
---|
| 740 | def visit_Return(self, node): |
---|
[4c72117] | 741 | self.add_current_line() |
---|
[5dd7cfb] | 742 | self.track_lineno(node) |
---|
[4c87de0] | 743 | self.in_expr = True |
---|
[71779b2] | 744 | if node.value is None: |
---|
| 745 | self.write_c('return') |
---|
| 746 | else: |
---|
[d5014e4] | 747 | self.write_c('return ') |
---|
[71779b2] | 748 | self.visit(node.value) |
---|
[98a4f14] | 749 | self.add_semi_colon() |
---|
[4c87de0] | 750 | self.in_expr = False |
---|
[fa74acf] | 751 | self.add_c_line(self.current_statement) |
---|
[71779b2] | 752 | self.current_statement = '' |
---|
| 753 | |
---|
| 754 | def visit_Break(self, node): |
---|
[5dd7cfb] | 755 | self.track_lineno(node) |
---|
[71779b2] | 756 | self.write_c('break') |
---|
| 757 | |
---|
| 758 | def visit_Continue(self, node): |
---|
[5dd7cfb] | 759 | self.track_lineno(node) |
---|
[71779b2] | 760 | self.write_c('continue') |
---|
| 761 | |
---|
| 762 | def visit_Raise(self, node): |
---|
[c01ed3e] | 763 | self.unsupported(node) |
---|
[15be191] | 764 | |
---|
[7b1dcf9] | 765 | # CRUFT: Python 2.6 / 3.0 compatibility |
---|
[5dd7cfb] | 766 | self.track_lineno(node) |
---|
[71779b2] | 767 | self.write_python('raise') |
---|
| 768 | if hasattr(node, 'exc') and node.exc is not None: |
---|
| 769 | self.write_python(' ') |
---|
| 770 | self.visit(node.exc) |
---|
| 771 | if node.cause is not None: |
---|
| 772 | self.write_python(' from ') |
---|
| 773 | self.visit(node.cause) |
---|
| 774 | elif hasattr(node, 'type') and node.type is not None: |
---|
| 775 | self.visit(node.type) |
---|
| 776 | if node.inst is not None: |
---|
| 777 | self.write_python(', ') |
---|
| 778 | self.visit(node.inst) |
---|
| 779 | if node.tback is not None: |
---|
| 780 | self.write_python(', ') |
---|
| 781 | self.visit(node.tback) |
---|
| 782 | |
---|
| 783 | # Expressions |
---|
| 784 | |
---|
| 785 | def visit_Attribute(self, node): |
---|
[c01ed3e] | 786 | self.unsupported(node, "attribute reference a.b not supported") |
---|
[15be191] | 787 | |
---|
[71779b2] | 788 | self.visit(node.value) |
---|
| 789 | self.write_python('.' + node.attr) |
---|
| 790 | |
---|
| 791 | def visit_Call(self, node): |
---|
| 792 | want_comma = [] |
---|
| 793 | def write_comma(): |
---|
| 794 | if want_comma: |
---|
| 795 | self.write_c(', ') |
---|
| 796 | else: |
---|
| 797 | want_comma.append(True) |
---|
[7b1dcf9] | 798 | if hasattr(node.func, 'id'): |
---|
[c01ed3e] | 799 | if node.func.id not in self.c_functions: |
---|
| 800 | self.c_functions.append(node.func.id) |
---|
[7b1dcf9] | 801 | if node.func.id == 'abs': |
---|
[fa74acf] | 802 | self.write_c("fabs ") |
---|
[7b1dcf9] | 803 | elif node.func.id == 'int': |
---|
[71779b2] | 804 | self.write_c('(int) ') |
---|
[7b1dcf9] | 805 | elif node.func.id == "SINCOS": |
---|
[c01ed3e] | 806 | self.write_sincos(node) |
---|
[71779b2] | 807 | return |
---|
| 808 | else: |
---|
| 809 | self.visit(node.func) |
---|
| 810 | else: |
---|
| 811 | self.visit(node.func) |
---|
| 812 | self.write_c('(') |
---|
| 813 | for arg in node.args: |
---|
| 814 | write_comma() |
---|
[fa74acf] | 815 | self.visited_args = True |
---|
[71779b2] | 816 | self.visit(arg) |
---|
| 817 | for keyword in node.keywords: |
---|
| 818 | write_comma() |
---|
| 819 | self.write_c(keyword.arg + '=') |
---|
| 820 | self.visit(keyword.value) |
---|
[fa74acf] | 821 | if hasattr(node, 'starargs'): |
---|
[71779b2] | 822 | if node.starargs is not None: |
---|
| 823 | write_comma() |
---|
| 824 | self.write_c('*') |
---|
| 825 | self.visit(node.starargs) |
---|
[fa74acf] | 826 | if hasattr(node, 'kwargs'): |
---|
[71779b2] | 827 | if node.kwargs is not None: |
---|
| 828 | write_comma() |
---|
| 829 | self.write_c('**') |
---|
| 830 | self.visit(node.kwargs) |
---|
[4c87de0] | 831 | self.write_c(')') |
---|
| 832 | if not self.in_expr: |
---|
| 833 | self.add_semi_colon() |
---|
[71779b2] | 834 | |
---|
[4c87de0] | 835 | TRANSLATE_CONSTANTS = { |
---|
[6a37819] | 836 | # python 2 uses normal name references through vist_Name |
---|
[4c87de0] | 837 | 'True': 'true', |
---|
| 838 | 'False': 'false', |
---|
| 839 | 'None': 'NULL', # "None" will probably fail for other reasons |
---|
[6a37819] | 840 | # python 3 uses NameConstant |
---|
| 841 | True: 'true', |
---|
| 842 | False: 'false', |
---|
| 843 | None: 'NULL', # "None" will probably fail for other reasons |
---|
[4c87de0] | 844 | } |
---|
[6a37819] | 845 | |
---|
[71779b2] | 846 | def visit_Name(self, node): |
---|
[4c87de0] | 847 | translation = self.TRANSLATE_CONSTANTS.get(node.id, None) |
---|
| 848 | if translation: |
---|
| 849 | self.write_c(translation) |
---|
| 850 | return |
---|
[71779b2] | 851 | self.write_c(node.id) |
---|
[c01ed3e] | 852 | if node.id in self.c_pointers and not self.in_subref: |
---|
[71779b2] | 853 | self.write_c("[0]") |
---|
| 854 | name = "" |
---|
| 855 | sub = node.id.find("[") |
---|
[7b1dcf9] | 856 | if sub > 0: |
---|
[71779b2] | 857 | name = node.id[0:sub].strip() |
---|
| 858 | else: |
---|
| 859 | name = node.id |
---|
[c01ed3e] | 860 | # add variable to c_vars if it ins't there yet, not an argument and not a number |
---|
| 861 | if (name not in self.c_functions and name not in self.c_vars and |
---|
| 862 | name not in self.c_int_vars and name not in self.arguments and |
---|
| 863 | name not in self.c_constants and not name.isdigit()): |
---|
| 864 | if self.in_subscript: |
---|
[d5014e4] | 865 | self.add_c_int_var(node.id) |
---|
[937afef] | 866 | else: |
---|
[c01ed3e] | 867 | self.c_vars.append(node.id) |
---|
[71779b2] | 868 | |
---|
[6a37819] | 869 | def visit_NameConstant(self, node): |
---|
| 870 | translation = self.TRANSLATE_CONSTANTS.get(node.value, None) |
---|
| 871 | if translation is not None: |
---|
| 872 | self.write_c(translation) |
---|
| 873 | else: |
---|
| 874 | self.unsupported(node, "don't know how to translate %r"%node.value) |
---|
| 875 | |
---|
[71779b2] | 876 | def visit_Str(self, node): |
---|
[4c87de0] | 877 | s = node.s |
---|
[5dd7cfb] | 878 | s = s.replace('\\', '\\\\').replace('"', '\\"').replace('\n', '\\n') |
---|
[4c87de0] | 879 | self.write_c('"') |
---|
| 880 | self.write_c(s) |
---|
| 881 | self.write_c('"') |
---|
[71779b2] | 882 | |
---|
| 883 | def visit_Bytes(self, node): |
---|
[4c87de0] | 884 | s = node.s |
---|
[5dd7cfb] | 885 | s = s.replace('\\', '\\\\').replace('"', '\\"').replace('\n', '\\n') |
---|
[4c87de0] | 886 | self.write_c('"') |
---|
| 887 | self.write_c(s) |
---|
| 888 | self.write_c('"') |
---|
[71779b2] | 889 | |
---|
| 890 | def visit_Num(self, node): |
---|
| 891 | self.write_c(repr(node.n)) |
---|
| 892 | |
---|
| 893 | def visit_Tuple(self, node): |
---|
| 894 | for idx, item in enumerate(node.elts): |
---|
| 895 | if idx: |
---|
[c01ed3e] | 896 | self.tuples.append(item) |
---|
[71779b2] | 897 | else: |
---|
| 898 | self.visit(item) |
---|
| 899 | |
---|
[5dd7cfb] | 900 | def visit_List(self, node): |
---|
| 901 | #self.unsupported(node) |
---|
| 902 | #print("visiting", node) |
---|
| 903 | #print(astor.to_source(node)) |
---|
| 904 | #print(node.elts) |
---|
| 905 | exprs = [render_expression(item) for item in node.elts] |
---|
| 906 | if exprs: |
---|
| 907 | self.c_vectors.append(', '.join(exprs)) |
---|
| 908 | vec_name = "vec" + str(len(self.c_vectors)) |
---|
| 909 | self.write_c(vec_name) |
---|
| 910 | |
---|
| 911 | def visit_Set(self, node): |
---|
| 912 | self.unsupported(node) |
---|
[71779b2] | 913 | |
---|
| 914 | def visit_Dict(self, node): |
---|
[c01ed3e] | 915 | self.unsupported(node) |
---|
[15be191] | 916 | |
---|
[71779b2] | 917 | self.write_python('{') |
---|
| 918 | for idx, (key, value) in enumerate(zip(node.keys, node.values)): |
---|
| 919 | if idx: |
---|
| 920 | self.write_python(', ') |
---|
| 921 | self.visit(key) |
---|
| 922 | self.write_python(': ') |
---|
| 923 | self.visit(value) |
---|
| 924 | self.write_python('}') |
---|
| 925 | |
---|
[fa74acf] | 926 | def get_special_power(self, string): |
---|
[3f9db6e] | 927 | function_name = '' |
---|
| 928 | is_negative_exp = False |
---|
[7b1dcf9] | 929 | if isevaluable(str(self.current_statement)): |
---|
[3f9db6e] | 930 | exponent = eval(string) |
---|
| 931 | is_negative_exp = exponent < 0 |
---|
| 932 | abs_exponent = abs(exponent) |
---|
[7b1dcf9] | 933 | if abs_exponent == 2: |
---|
[3f9db6e] | 934 | function_name = "square" |
---|
[7b1dcf9] | 935 | elif abs_exponent == 3: |
---|
[3f9db6e] | 936 | function_name = "cube" |
---|
[7b1dcf9] | 937 | elif abs_exponent == 0.5: |
---|
[3f9db6e] | 938 | function_name = "sqrt" |
---|
[7b1dcf9] | 939 | elif abs_exponent == 1.0/3.0: |
---|
[3f9db6e] | 940 | function_name = "cbrt" |
---|
[7b1dcf9] | 941 | if function_name == '': |
---|
[3f9db6e] | 942 | function_name = "pow" |
---|
| 943 | return function_name, is_negative_exp |
---|
| 944 | |
---|
[fa74acf] | 945 | def translate_power(self, node): |
---|
[7b1dcf9] | 946 | # get exponent by visiting the right hand argument. |
---|
[71779b2] | 947 | function_name = "pow" |
---|
| 948 | temp_statement = self.current_statement |
---|
[7b1dcf9] | 949 | # 'visit' functions write the results to the 'current_statement' class memnber |
---|
| 950 | # Here, a temporary variable, 'temp_statement', is used, that enables the |
---|
| 951 | # use of the 'visit' function |
---|
[71779b2] | 952 | self.current_statement = '' |
---|
| 953 | self.visit(node.right) |
---|
[fa74acf] | 954 | exponent = self.current_statement.replace(' ', '') |
---|
| 955 | function_name, is_negative_exp = self.get_special_power(self.current_statement) |
---|
[71779b2] | 956 | self.current_statement = temp_statement |
---|
[7b1dcf9] | 957 | if is_negative_exp: |
---|
[fa74acf] | 958 | self.write_c("1.0 /(") |
---|
| 959 | self.write_c(function_name + "(") |
---|
[71779b2] | 960 | self.visit(node.left) |
---|
[7b1dcf9] | 961 | if function_name == "pow": |
---|
[71779b2] | 962 | self.write_c(", ") |
---|
| 963 | self.visit(node.right) |
---|
| 964 | self.write_c(")") |
---|
[7b1dcf9] | 965 | if is_negative_exp: |
---|
[3f9db6e] | 966 | self.write_c(")") |
---|
[4c87de0] | 967 | #self.write_c(" ") |
---|
[71779b2] | 968 | |
---|
[fa74acf] | 969 | def translate_integer_divide(self, node): |
---|
[4c87de0] | 970 | self.write_c("(int)((") |
---|
[71779b2] | 971 | self.visit(node.left) |
---|
[4c87de0] | 972 | self.write_c(")/(") |
---|
[71779b2] | 973 | self.visit(node.right) |
---|
[4c87de0] | 974 | self.write_c("))") |
---|
| 975 | |
---|
| 976 | def translate_float_divide(self, node): |
---|
| 977 | self.write_c("((double)(") |
---|
| 978 | self.visit(node.left) |
---|
| 979 | self.write_c(")/(double)(") |
---|
| 980 | self.visit(node.right) |
---|
| 981 | self.write_c("))") |
---|
[71779b2] | 982 | |
---|
| 983 | def visit_BinOp(self, node): |
---|
[fb5c8c7] | 984 | self.write_c("(") |
---|
[7b1dcf9] | 985 | if '%s' % BINOP_SYMBOLS[type(node.op)] == BINOP_SYMBOLS[ast.Pow]: |
---|
[fa74acf] | 986 | self.translate_power(node) |
---|
[7b1dcf9] | 987 | elif '%s' % BINOP_SYMBOLS[type(node.op)] == BINOP_SYMBOLS[ast.FloorDiv]: |
---|
[fa74acf] | 988 | self.translate_integer_divide(node) |
---|
[4c87de0] | 989 | elif '%s' % BINOP_SYMBOLS[type(node.op)] == BINOP_SYMBOLS[ast.Div]: |
---|
| 990 | self.translate_float_divide(node) |
---|
[71779b2] | 991 | else: |
---|
| 992 | self.visit(node.left) |
---|
| 993 | self.write_c(' %s ' % BINOP_SYMBOLS[type(node.op)]) |
---|
| 994 | self.visit(node.right) |
---|
[fb5c8c7] | 995 | self.write_c(")") |
---|
| 996 | |
---|
[7b1dcf9] | 997 | # for C |
---|
[71779b2] | 998 | def visit_BoolOp(self, node): |
---|
| 999 | self.write_c('(') |
---|
| 1000 | for idx, value in enumerate(node.values): |
---|
| 1001 | if idx: |
---|
| 1002 | self.write_c(' %s ' % BOOLOP_SYMBOLS[type(node.op)]) |
---|
| 1003 | self.visit(value) |
---|
| 1004 | self.write_c(')') |
---|
| 1005 | |
---|
| 1006 | def visit_Compare(self, node): |
---|
| 1007 | self.write_c('(') |
---|
| 1008 | self.visit(node.left) |
---|
| 1009 | for op, right in zip(node.ops, node.comparators): |
---|
| 1010 | self.write_c(' %s ' % CMPOP_SYMBOLS[type(op)]) |
---|
| 1011 | self.visit(right) |
---|
| 1012 | self.write_c(')') |
---|
| 1013 | |
---|
| 1014 | def visit_UnaryOp(self, node): |
---|
| 1015 | self.write_c('(') |
---|
| 1016 | op = UNARYOP_SYMBOLS[type(node.op)] |
---|
| 1017 | self.write_c(op) |
---|
| 1018 | if op == 'not': |
---|
| 1019 | self.write_c(' ') |
---|
| 1020 | self.visit(node.operand) |
---|
| 1021 | self.write_c(')') |
---|
| 1022 | |
---|
| 1023 | def visit_Subscript(self, node): |
---|
[c01ed3e] | 1024 | if node.value.id not in self.c_constants: |
---|
| 1025 | if node.value.id not in self.c_pointers: |
---|
| 1026 | self.c_pointers.append(node.value.id) |
---|
| 1027 | self.in_subref = True |
---|
[71779b2] | 1028 | self.visit(node.value) |
---|
[c01ed3e] | 1029 | self.in_subref = False |
---|
[71779b2] | 1030 | self.write_c('[') |
---|
[c01ed3e] | 1031 | self.in_subscript = True |
---|
[71779b2] | 1032 | self.visit(node.slice) |
---|
[c01ed3e] | 1033 | self.in_subscript = False |
---|
[71779b2] | 1034 | self.write_c(']') |
---|
| 1035 | |
---|
| 1036 | def visit_Slice(self, node): |
---|
| 1037 | if node.lower is not None: |
---|
| 1038 | self.visit(node.lower) |
---|
| 1039 | self.write_python(':') |
---|
| 1040 | if node.upper is not None: |
---|
| 1041 | self.visit(node.upper) |
---|
| 1042 | if node.step is not None: |
---|
| 1043 | self.write_python(':') |
---|
[fa74acf] | 1044 | if not(isinstance(node.step, Name) and node.step.id == 'None'): |
---|
[71779b2] | 1045 | self.visit(node.step) |
---|
| 1046 | |
---|
| 1047 | def visit_ExtSlice(self, node): |
---|
| 1048 | for idx, item in node.dims: |
---|
| 1049 | if idx: |
---|
| 1050 | self.write_python(', ') |
---|
| 1051 | self.visit(item) |
---|
| 1052 | |
---|
| 1053 | def visit_Yield(self, node): |
---|
[c01ed3e] | 1054 | self.unsupported(node) |
---|
[15be191] | 1055 | |
---|
[71779b2] | 1056 | self.write_python('yield ') |
---|
| 1057 | self.visit(node.value) |
---|
| 1058 | |
---|
| 1059 | def visit_Lambda(self, node): |
---|
[c01ed3e] | 1060 | self.unsupported(node) |
---|
[15be191] | 1061 | |
---|
[71779b2] | 1062 | self.write_python('lambda ') |
---|
| 1063 | self.visit(node.args) |
---|
| 1064 | self.write_python(': ') |
---|
| 1065 | self.visit(node.body) |
---|
| 1066 | |
---|
| 1067 | def visit_Ellipsis(self, node): |
---|
[c01ed3e] | 1068 | self.unsupported(node) |
---|
[15be191] | 1069 | |
---|
[71779b2] | 1070 | self.write_python('Ellipsis') |
---|
| 1071 | |
---|
| 1072 | def generator_visit(left, right): |
---|
| 1073 | def visit(self, node): |
---|
| 1074 | self.write_python(left) |
---|
[937afef] | 1075 | self.write_c(left) |
---|
[71779b2] | 1076 | self.visit(node.elt) |
---|
| 1077 | for comprehension in node.generators: |
---|
| 1078 | self.visit(comprehension) |
---|
[937afef] | 1079 | self.write_c(right) |
---|
[7b1dcf9] | 1080 | #self.write_python(right) |
---|
[71779b2] | 1081 | return visit |
---|
| 1082 | |
---|
| 1083 | visit_ListComp = generator_visit('[', ']') |
---|
| 1084 | visit_GeneratorExp = generator_visit('(', ')') |
---|
| 1085 | visit_SetComp = generator_visit('{', '}') |
---|
| 1086 | del generator_visit |
---|
| 1087 | |
---|
| 1088 | def visit_DictComp(self, node): |
---|
[c01ed3e] | 1089 | self.unsupported(node) |
---|
[15be191] | 1090 | |
---|
[71779b2] | 1091 | self.write_python('{') |
---|
| 1092 | self.visit(node.key) |
---|
| 1093 | self.write_python(': ') |
---|
| 1094 | self.visit(node.value) |
---|
| 1095 | for comprehension in node.generators: |
---|
| 1096 | self.visit(comprehension) |
---|
| 1097 | self.write_python('}') |
---|
| 1098 | |
---|
| 1099 | def visit_IfExp(self, node): |
---|
[4c87de0] | 1100 | self.write_c('((') |
---|
[71779b2] | 1101 | self.visit(node.test) |
---|
[4c87de0] | 1102 | self.write_c(')?(') |
---|
| 1103 | self.visit(node.body) |
---|
| 1104 | self.write_c('):(') |
---|
[71779b2] | 1105 | self.visit(node.orelse) |
---|
[4c87de0] | 1106 | self.write_c('))') |
---|
[71779b2] | 1107 | |
---|
| 1108 | def visit_Starred(self, node): |
---|
| 1109 | self.write_c('*') |
---|
| 1110 | self.visit(node.value) |
---|
| 1111 | |
---|
| 1112 | def visit_Repr(self, node): |
---|
[7b1dcf9] | 1113 | # CRUFT: python 2.6 only |
---|
[71779b2] | 1114 | self.write_c('`') |
---|
| 1115 | self.visit(node.value) |
---|
| 1116 | self.write_python('`') |
---|
| 1117 | |
---|
| 1118 | # Helper Nodes |
---|
| 1119 | |
---|
| 1120 | def visit_alias(self, node): |
---|
[c01ed3e] | 1121 | self.unsupported(node) |
---|
[15be191] | 1122 | |
---|
[71779b2] | 1123 | self.write_python(node.name) |
---|
| 1124 | if node.asname is not None: |
---|
| 1125 | self.write_python(' as ' + node.asname) |
---|
| 1126 | |
---|
| 1127 | def visit_comprehension(self, node): |
---|
[5dd7cfb] | 1128 | self.unsupported(node) |
---|
| 1129 | |
---|
[71779b2] | 1130 | self.write_c(' for ') |
---|
| 1131 | self.visit(node.target) |
---|
[5dd7cfb] | 1132 | self.write_c(' in ') |
---|
[7b1dcf9] | 1133 | #self.write_python(' in ') |
---|
[71779b2] | 1134 | self.visit(node.iter) |
---|
| 1135 | if node.ifs: |
---|
| 1136 | for if_ in node.ifs: |
---|
| 1137 | self.write_python(' if ') |
---|
| 1138 | self.visit(if_) |
---|
| 1139 | |
---|
| 1140 | def visit_arguments(self, node): |
---|
| 1141 | self.signature(node) |
---|
| 1142 | |
---|
[c01ed3e] | 1143 | def unsupported(self, node, message=None): |
---|
| 1144 | if hasattr(node, "value"): |
---|
| 1145 | lineno = node.value.lineno |
---|
| 1146 | elif hasattr(node, "iter"): |
---|
| 1147 | lineno = node.iter.lineno |
---|
| 1148 | else: |
---|
| 1149 | #print(dir(node)) |
---|
| 1150 | lineno = 0 |
---|
| 1151 | |
---|
| 1152 | lineno += self.lineno_offset |
---|
| 1153 | if self.fname: |
---|
| 1154 | location = "%s(%d)" % (self.fname, lineno) |
---|
| 1155 | else: |
---|
| 1156 | location = "%d" % (self.fname, lineno) |
---|
| 1157 | if self.current_function: |
---|
| 1158 | location += ", function %s" % self.current_function |
---|
| 1159 | if message is None: |
---|
| 1160 | message = node.__class__.__name__ + " syntax not supported" |
---|
| 1161 | raise SyntaxError("[%s] %s" % (location, message)) |
---|
| 1162 | |
---|
[71779b2] | 1163 | def print_function(f=None): |
---|
| 1164 | """ |
---|
| 1165 | Print out the code for the function |
---|
| 1166 | """ |
---|
| 1167 | # Include some comments to see if they get printed |
---|
| 1168 | import ast |
---|
| 1169 | import inspect |
---|
| 1170 | if f is not None: |
---|
| 1171 | tree = ast.parse(inspect.getsource(f)) |
---|
[fa74acf] | 1172 | tree_source = to_source(tree) |
---|
[71779b2] | 1173 | print(tree_source) |
---|
| 1174 | |
---|
[c01ed3e] | 1175 | def define_constant(name, value, block_size=1): |
---|
| 1176 | # type: (str, any, int) -> str |
---|
| 1177 | """ |
---|
| 1178 | Convert a python constant into a C constant of the same name. |
---|
| 1179 | |
---|
| 1180 | Returns the C declaration of the constant as a string, possibly containing |
---|
| 1181 | line feeds. The string will not be indented. |
---|
| 1182 | |
---|
| 1183 | Supports int, double and sequences of double. |
---|
| 1184 | """ |
---|
| 1185 | const = "constant " # OpenCL needs globals to be constant |
---|
| 1186 | if isinstance(value, int): |
---|
[15be191] | 1187 | parts = [const + "int ", name, " = ", "%d"%value, ";\n"] |
---|
[c01ed3e] | 1188 | elif isinstance(value, float): |
---|
[15be191] | 1189 | parts = [const + "double ", name, " = ", "%.15g"%value, ";\n"] |
---|
[c01ed3e] | 1190 | else: |
---|
| 1191 | try: |
---|
| 1192 | len(value) |
---|
| 1193 | except TypeError: |
---|
| 1194 | raise TypeError("constant %s must be int, float or [float, ...]"%name) |
---|
| 1195 | # extend constant arrays to a multiple of 4; not sure if this |
---|
| 1196 | # is necessary, but some OpenCL targets broke if the number |
---|
| 1197 | # of parameters in the parameter table was not a multiple of 4, |
---|
| 1198 | # so do it for all constant arrays to be safe. |
---|
| 1199 | if len(value)%block_size != 0: |
---|
| 1200 | value = list(value) + [0.]*(block_size - len(value)%block_size) |
---|
| 1201 | elements = ["%.15g"%v for v in value] |
---|
| 1202 | parts = [const + "double ", name, "[]", " = ", |
---|
[15be191] | 1203 | "{\n ", ", ".join(elements), "\n};\n"] |
---|
[c01ed3e] | 1204 | |
---|
| 1205 | return "".join(parts) |
---|
| 1206 | |
---|
| 1207 | |
---|
| 1208 | # Modified from the following: |
---|
| 1209 | # |
---|
| 1210 | # http://code.activestate.com/recipes/578272-topological-sort/ |
---|
| 1211 | # Copyright (C) 2012 Sam Denton |
---|
| 1212 | # License: MIT |
---|
| 1213 | def ordered_dag(dag): |
---|
| 1214 | # type: (Dict[T, Set[T]]) -> Iterator[T] |
---|
| 1215 | """ |
---|
| 1216 | Given a dag defined by a dictionary of {k1: [k2, ...]} yield keys |
---|
| 1217 | in order such that every key occurs after the keys it depends upon. |
---|
| 1218 | |
---|
| 1219 | This is an iterator not a sequence. To reverse it use:: |
---|
| 1220 | |
---|
| 1221 | reversed(tuple(ordered_dag(dag))) |
---|
| 1222 | |
---|
| 1223 | Raise an error if there are any cycles. |
---|
| 1224 | |
---|
| 1225 | Keys are arbitrary hashable values. |
---|
| 1226 | """ |
---|
| 1227 | # Local import to make the function stand-alone, and easier to borrow |
---|
| 1228 | from functools import reduce |
---|
| 1229 | |
---|
| 1230 | dag = dag.copy() |
---|
| 1231 | |
---|
| 1232 | # make leaves depend on the empty set |
---|
| 1233 | leaves = reduce(set.union, dag.values()) - set(dag.keys()) |
---|
| 1234 | dag.update({node: set() for node in leaves}) |
---|
| 1235 | while True: |
---|
| 1236 | leaves = set(node for node, links in dag.items() if not links) |
---|
| 1237 | if not leaves: |
---|
| 1238 | break |
---|
| 1239 | for node in leaves: |
---|
| 1240 | yield node |
---|
| 1241 | dag = {node: (links-leaves) |
---|
| 1242 | for node, links in dag.items() if node not in leaves} |
---|
| 1243 | if dag: |
---|
| 1244 | raise ValueError("Cyclic dependes exists amongst these items:\n%s" |
---|
| 1245 | % ", ".join(str(node) for node in dag.keys())) |
---|
| 1246 | |
---|
[4c87de0] | 1247 | import re |
---|
| 1248 | PRINT_ARGS = re.compile(r'print[(]"(?P<template>[^"]*)" *% *[(](?P<args>[^\n]*)[)] *[)] *\n') |
---|
| 1249 | SUBST_ARGS = r'printf("\g<template>\\n", \g<args>)\n' |
---|
| 1250 | PRINT_STR = re.compile(r'print[(]"(?P<template>[^"]*)" *[)] *\n') |
---|
| 1251 | SUBST_STR = r'printf("\g<template>\n")' |
---|
[c01ed3e] | 1252 | def translate(functions, constants=None): |
---|
[4c87de0] | 1253 | # type: (Sequence[(str, str, int)], Dict[str, any]) -> List[str] |
---|
[c01ed3e] | 1254 | """ |
---|
[4c87de0] | 1255 | Convert a list of functions to a list of C code strings. |
---|
| 1256 | |
---|
[d5014e4] | 1257 | Returns list of corresponding code snippets (with trailing lines in |
---|
| 1258 | each block) and a list of warnings generated by the translator. |
---|
| 1259 | |
---|
[4c87de0] | 1260 | A function is given by the tuple (source, filename, line number). |
---|
| 1261 | |
---|
| 1262 | Global constants are given in a dictionary of {name: value}. The |
---|
| 1263 | constants are used for name space resolution and type inferencing. |
---|
| 1264 | Constants are not translated by this code. Instead, call |
---|
| 1265 | :func:`define_constant` with name and value, and maybe block_size |
---|
| 1266 | if arrays need to be padded to the next block boundary. |
---|
| 1267 | |
---|
| 1268 | Function prototypes are not generated. Use :func:`ordered_dag` |
---|
| 1269 | to list the functions in reverse order of dependency before calling |
---|
| 1270 | translate. [Maybe a future revision will return the function prototypes |
---|
| 1271 | so that a suitable "*.h" file can be generated. |
---|
[c01ed3e] | 1272 | """ |
---|
[1ddb794] | 1273 | snippets = [] |
---|
[4339764] | 1274 | warnings = [] |
---|
[c01ed3e] | 1275 | for source, fname, lineno in functions: |
---|
[15be191] | 1276 | line_directive = '#line %d "%s"\n'%(lineno, fname.replace('\\', '\\\\')) |
---|
[1ddb794] | 1277 | snippets.append(line_directive) |
---|
[4c87de0] | 1278 | # Replace simple print function calls with printf statements |
---|
| 1279 | source = PRINT_ARGS.sub(SUBST_ARGS, source) |
---|
| 1280 | source = PRINT_STR.sub(SUBST_STR, source) |
---|
[98a4f14] | 1281 | tree = ast.parse(source) |
---|
[d5014e4] | 1282 | generator = SourceGenerator(constants=constants, fname=fname, lineno=lineno) |
---|
| 1283 | generator.visit(tree) |
---|
| 1284 | c_code = "".join(generator.c_proc) |
---|
[1ddb794] | 1285 | snippets.append(c_code) |
---|
[d5014e4] | 1286 | warnings.extend(generator.warnings) |
---|
| 1287 | return snippets, warnings |
---|
[71779b2] | 1288 | |
---|
[4339764] | 1289 | |
---|
[6f91c91] | 1290 | C_HEADER_LINENO = getframeinfo(currentframe()).lineno + 2 |
---|
[0bd0877] | 1291 | C_HEADER = """ |
---|
[6f91c91] | 1292 | #line %d "%s" |
---|
[0bd0877] | 1293 | #include <stdio.h> |
---|
| 1294 | #include <stdbool.h> |
---|
| 1295 | #include <math.h> |
---|
| 1296 | #define constant const |
---|
| 1297 | double square(double x) { return x*x; } |
---|
| 1298 | double cube(double x) { return x*x*x; } |
---|
| 1299 | double polyval(constant double *coef, double x, int N) |
---|
| 1300 | { |
---|
| 1301 | int i = 0; |
---|
| 1302 | double ans = coef[0]; |
---|
| 1303 | |
---|
| 1304 | while (i < N) { |
---|
| 1305 | ans = ans * x + coef[i++]; |
---|
| 1306 | } |
---|
| 1307 | |
---|
| 1308 | return ans; |
---|
| 1309 | } |
---|
| 1310 | """ |
---|
[d5014e4] | 1311 | |
---|
| 1312 | USAGE = """\ |
---|
| 1313 | Usage: python py2c.py <infile> [<outfile>] |
---|
| 1314 | |
---|
| 1315 | if outfile is omitted, output file is '<infile>.c' |
---|
| 1316 | """ |
---|
[71779b2] | 1317 | |
---|
[7b1dcf9] | 1318 | def main(): |
---|
[71779b2] | 1319 | import os |
---|
[d7f33e5] | 1320 | #print("Parsing...using Python" + sys.version) |
---|
[7b1dcf9] | 1321 | if len(sys.argv) == 1: |
---|
[d5014e4] | 1322 | print(USAGE) |
---|
[7b1dcf9] | 1323 | return |
---|
| 1324 | |
---|
| 1325 | fname_in = sys.argv[1] |
---|
| 1326 | if len(sys.argv) == 2: |
---|
| 1327 | fname_base = os.path.splitext(fname_in)[0] |
---|
| 1328 | fname_out = str(fname_base) + '.c' |
---|
| 1329 | else: |
---|
| 1330 | fname_out = sys.argv[2] |
---|
| 1331 | |
---|
| 1332 | with open(fname_in, "r") as python_file: |
---|
| 1333 | code = python_file.read() |
---|
[c01ed3e] | 1334 | name = "gauss" |
---|
| 1335 | code = (code |
---|
| 1336 | .replace(name+'.n', 'GAUSS_N') |
---|
| 1337 | .replace(name+'.z', 'GAUSS_Z') |
---|
[4c87de0] | 1338 | .replace(name+'.w', 'GAUSS_W') |
---|
| 1339 | .replace('if __name__ == "__main__"', "def main()") |
---|
[0bd0877] | 1340 | ) |
---|
[7b1dcf9] | 1341 | |
---|
[d5014e4] | 1342 | translation, warnings = translate([(code, fname_in, 1)]) |
---|
| 1343 | c_code = "".join(translation) |
---|
[4c87de0] | 1344 | c_code = c_code.replace("double main()", "int main(int argc, char *argv[])") |
---|
[7b1dcf9] | 1345 | |
---|
| 1346 | with open(fname_out, "w") as file_out: |
---|
[6f91c91] | 1347 | file_out.write(C_HEADER%(C_HEADER_LINENO, __file__)) |
---|
[4c87de0] | 1348 | file_out.write(c_code) |
---|
[d5014e4] | 1349 | |
---|
[4339764] | 1350 | if warnings: |
---|
[d5014e4] | 1351 | print("\n".join(warnings)) |
---|
[d7f33e5] | 1352 | #print("...Done") |
---|
[7b1dcf9] | 1353 | |
---|
| 1354 | if __name__ == "__main__": |
---|
| 1355 | main() |
---|