diff --git a/code/code_gen/post.py b/code/code_gen/post.py index eacce7f..a6125dc 100644 --- a/code/code_gen/post.py +++ b/code/code_gen/post.py @@ -431,7 +431,13 @@ class Field (FieldBase): else: self.latex = self.latex_given - self.exporter = self.attr.get("export") + self.exporter = None + try: + if self.attr["export"]: + self.exporter = FieldExporter(self.name, self.attr, self) + except KeyError: + pass + ''' for a in self.dep: diff --git a/code/termspecs/termspec.vnsd_fluctuation b/code/termspecs/termspec.vnsd_fluctuation index 0255baa..25cace5 100644 --- a/code/termspecs/termspec.vnsd_fluctuation +++ b/code/termspecs/termspec.vnsd_fluctuation @@ -2,43 +2,42 @@ [u, v, w, y] -c = 1.0 - y -wrate = rxn_rate(c) -fsd = sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c))) +c(latex="c") = 1.0 - y +wrate (latex="\dot{\omega}") = rxn_rate(c) +fsd (latex="\Sigma_f'")= sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c))) -sd = ((d2dx(c) + d2dy(c) + d2dz(c)) * $rod + wrate) / fsd +lapc(latex="\Delta c") = d2dx(c) + d2dy(c) + d2dz(c) +sd(latex="S_d") = ($rod * lapc + wrate) / fsd -nx = - ddx(c) / fsd -ny = - ddy(c) / fsd -nz = - ddz(c) / fsd +nx(latex="n_x") = - ddx(c) / fsd +ny(latex="n_y") = - ddy(c) / fsd +nz(latex="n_z") = - ddz(c) / fsd unx = u * nx vny = v * ny wnz = w * nz -vn = unx + vny + wnz +vn(latex="v_n") = unx + vny + wnz -divn = ddx(nx) + ddy(ny) + ddz(nz) +divn(latex="\nabla \cdot \vec{n}") = ddx(nx) + ddy(ny) + ddz(nz) divn2 = sqr(divn) # grad(n) : grad(n) -gradn_gradn = sqr(ddx(nx)) + sqr(ddy(ny)) + sqr(ddz(nz)) + 2*(ddx(ny)*ddy(nx)) + 2*(ddx(nz)*ddz(nx)) + 2*(ddy(nz)*ddz(ny)) +gradn_gradn(latex="\nabla \vec{n} : \nabla \vec{n}") = sqr(ddx(nx)) + sqr(ddy(ny)) + sqr(ddz(nz)) + 2*(ddx(ny)*ddy(nx)) + 2*(ddx(nz)*ddz(nx)) + 2*(ddy(nz)*ddz(ny)) absk = abs(divn) -lapc = d2dx(c) + d2dy(c) + d2dz(c) - vnsd = vn+sd divn_nx = divn*nx -gradn_vnsd = ddx(vnsd) * nx + ddy(vnsd) * ny + ddz(vnsd) * nz +gradn_vnsd (latex="\nabla_{\mathrm{n}}(v_n+S_d)") = ddx(vnsd) * nx + ddy(vnsd) * ny + ddz(vnsd) * nz -gradt_vnsd1 = ddx(vnsd) - nx * gradn_vnsd +gradt_vnsd1 (latex="\nabla_{\top,x}(v_n+S_d)")= ddx(vnsd) - nx * gradn_vnsd divn2_dndn = divn2 - gradn_gradn -sign_divn = $dsign(1.0, divn) +sign_divn(latex="\mathrm{sign}(\nabla \cdot \vec{n})") = $dsign(1.0, divn) sign_divn2_dndn = sign_divn * divn2_dndn