From 7e533976e1f48ce6832c653bf8bebac44f51932a Mon Sep 17 00:00:00 2001 From: dedner <a.s.dedner@warwick.ac.uk> Date: Fri, 27 Jul 2018 11:42:13 +0100 Subject: [PATCH] add a euler module --- pydemo/euler/euler.py | 9 +++++++++ 1 file changed, 9 insertions(+) create mode 100644 pydemo/euler/euler.py diff --git a/pydemo/euler/euler.py b/pydemo/euler/euler.py new file mode 100644 index 00000000..333932d5 --- /dev/null +++ b/pydemo/euler/euler.py @@ -0,0 +1,9 @@ +from ufl import * +exec(open("euler.ufl").read()) + +def RiemanProblem(UL,UR): + x = SpatialCoordinate(space.cell()) + return conditional(x[0]<0,as_vector(UL),as_vector(UR)) +def Sod(): + return RiemanProblem( Model.toCons([1,0,0,1]), + Model.toCons([0.125,0,0,0.1]) ) -- GitLab