diff --git a/interp/baker/__init__.py b/interp/baker/__init__.py index 26392c3..0d8c5e5 100644 --- a/interp/baker/__init__.py +++ b/interp/baker/__init__.py @@ -103,13 +103,13 @@ def get_error(phi, R, S, order = 2): B = [] # baker eq 9 w = [] # baker eq 11 - p = pattern(len(phi), order) - log.info("pattern: %s" % p) + cur_pattern = pattern(len(phi), order) + log.info("pattern: %s" % cur_pattern) for (s,q) in zip(S.verts, S.q): cur_phi, cur_qlin = qlinear(s, R) l = [] - for i in p: + for i in cur_pattern: cur_sum = cur_phi[i[0]] for j in i[1:]: cur_sum *= cur_phi[j] @@ -136,7 +136,7 @@ def get_error(phi, R, S, order = 2): abc = np.dot(np.linalg.pinv(A), b) error_term = 0.0 - for (a, i) in zip(abc, p): + for (a, i) in zip(abc, cur_pattern): cur_sum = a for j in i: cur_sum *= phi[j]