Source code for sfepy.discrete.parse_equations

from pyparsing import Combine, Literal, Word, delimitedList, Group, Optional,\
     ZeroOrMore, OneOrMore, nums, alphas, alphanums,\
     StringStart, StringEnd, CaselessLiteral, Forward, oneOf

[docs]class TermParse(object): def __str__(self): ss = "%s\n" % self.__class__ for key, val in self.__dict__.iteritems(): ss += " %s:\n %s\n" % (key, self.__dict__[key]) return ss
[docs]def collect_term(term_descs, lc): signs = {'+': 1.0, '-': -1.0} def append(str, loc, toks): sign = signs[toks.sign] * signs[lc[0]] tp = TermParse() tp.integral = toks.term_desc.integral if not tp.integral: tp.integral = 'a' tp.region = toks.term_desc.region tp.flag = toks.term_desc.flag tp.sign = sign * eval(''.join(toks.mul)) tp.name = toks.term_desc.name tp.args = ', '.join(toks.args[0]) term_descs.append(tp) return append
[docs]def rhs(lc): def aux(str, loc, toks): if toks: lc[0] = '-' return aux
[docs]def create_bnf(term_descs): """term_descs .. list of TermParse objects (sign, term_name, term_arg_names), where sign can be real or complex multiplier""" lc = ['+'] # Linear combination context. equal = Literal("=").setParseAction(rhs(lc)) zero = Literal("0").suppress() point = Literal(".") e = CaselessLiteral("E") inumber = Word("+-" + nums, nums) fnumber = Combine(Word("+-" + nums, nums) + Optional(point + Optional(Word(nums))) + Optional(e + Word("+-" + nums, nums))) number = fnumber + Optional(Literal('j'), default='') add_op = oneOf('+ -') number_expr = Forward() number_expr << ZeroOrMore('(') + number \ + ZeroOrMore(add_op + number_expr) \ + ZeroOrMore(')') ident = Word(alphas, alphanums + "_") integral = Combine((Literal('i') + Word(alphanums)) | Literal('i') | Literal('a') | Word(nums))("integral") history = Optional('[' + inumber + ']', default='')("history") variable = Combine(Word(alphas, alphanums + '._') + history) derivative = Combine(Literal('d') + variable \ + Literal('/') + Literal('dt')) trace = Combine(Literal('tr') + '(' + variable + ')') generalized_var = derivative | trace | variable args = Group(delimitedList(generalized_var)) flag = Literal('a') term = Optional(Literal('+') | Literal('-'), default='+')("sign") \ + Optional(number_expr + Literal('*').suppress(), default=['1.0', ''])("mul") \ + Combine(ident("name") \ + Optional("." + (integral + "." + ident("region") + "." + flag("flag") | integral + "." + ident("region") | ident("region") )))("term_desc") + "(" \ + Optional(args, default=[''])("args") + ")" term.setParseAction(collect_term(term_descs, lc)) rhs1 = equal + OneOrMore(term) rhs2 = equal + zero equation = StringStart() + OneOrMore(term) \ + Optional(rhs1 | rhs2) + StringEnd() ## term.setDebug() return equation
if __name__ == "__main__": test_str = """d_term1.Y(fluid, u, w, Nu, dcf, mode) + 5.0 * d_term2.Omega(u, w, Nu, dcf, mode) - d_another_term.Elsewhere(w, p[-1], Nu, dcf, mode) = - dw_rhs.a.Y3(u, q, Nu, dcf, mode)""" term_descs = [] bnf = create_bnf(term_descs) out = bnf.parseString(test_str) print 'out:', out, '\n' for tp in term_descs: print tp