Python converters from SymPy math expression to C code optimized for speed (SymPy CSE, GCC -Ofast)

gistfile1.py

from sympy import Symbol cse, ccode
from sympy.utilities import numbered_symbols

# Inspired by http://stackoverflow.com/questions/22665990/optimize-code-generated-by-sympy
def sympyToC(symname, symfunc):
    tmpsyms = numbered_symbols("tmp")
    symbols, simple = cse(symfunc, symbols=tmpsyms)
    symbolslist = map(lambda x:str(x), list(symfunc.atoms(Symbol)) )
    symbolslist.sort()
    varstring=",".join( " double "+x for x in symbolslist )[1:]
    c_code = "\n// Needed GCC 4.7+\n\n"
    c_code += "#pragma GCC push_options\n"
    c_code += "#pragma GCC optimize (\"Ofast\")\n\n"
    c_code += "#include \"" + symname + ".h\"\n"
    c_code += "#include <math.h>\n\n"
    c_code += "double "+str(symname)+"("+varstring+") {\n"
    for s in symbols:
        c_code +=  "  double " +ccode(s[0]) + " = " + ccode(s[1]) + ";\n"
    c_code +=  "  double r;\n"
    c_code +=  "  r = " + ccode(simple[0])+";\n"
    c_code +=  "  return r;\n"
    c_code += "}\n\n"
    c_code += "#pragma GCC pop_options\n"
    return c_code
    
def sympyToH(symname, symfunc):
    h_code = "\n"
    h_code += "#ifndef PROJECT__" + symname.upper() + "__H\n"
    h_code += "#define PROJECT__" + symname.upper() + "__H\n"
    h_code += "#ifdef __cplusplus\n"
    h_code += "extern \"C\" {\n"
    h_code += "#endif\n"
    symbolslist = map(lambda x:str(x), list(symfunc.atoms(Symbol)) )
    symbolslist.sort()
    varstring=",".join( " double "+x for x in symbolslist )[1:]
    h_code += "double "+str(symname)+"("+varstring+");\n"
    h_code += "#ifdef __cplusplus\n"
    h_code += "}\n"
    h_code += "#endif\n"
    h_code += "#endif\n"
    return h_code