-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathBDNK_full.txt
More file actions
5 lines (5 loc) · 9.66 KB
/
BDNK_full.txt
File metadata and controls
5 lines (5 loc) · 9.66 KB
1
2
3
4
5
Derivative(p(t, x, y, z), t) = (-(Gamma - 1)*(p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1) - (-Gamma*W(t, x, y, z) + W(t, x, y, z))*(rho(t, x, y, z) - p(t, x, y, z)/(1 - Gamma))*W(t, x, y, z)/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1))*Derivative(v1(t, x, y, z), x) + (-(Gamma - 1)*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1) - (-Gamma*W(t, x, y, z) + W(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1))*Derivative(rho(t, x, y, z), x) + (-(Gamma - 1)*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1) + (-Gamma*W(t, x, y, z) + W(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/((1 - Gamma)*(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1)))*Derivative(p(t, x, y, z), x)
Derivative(rho(t, x, y, z), t) = ((p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1) - (rho(t, x, y, z) - p(t, x, y, z)/(1 - Gamma))*(Gamma*W(t, x, y, z)**2 - Gamma - W(t, x, y, z)**2 + 1)*W(t, x, y, z)/(Gamma*W(t, x, y, z)**3 - Gamma*W(t, x, y, z) - 2*W(t, x, y, z)**3 + W(t, x, y, z)))*Derivative(v1(t, x, y, z), x) + (W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1) + (Gamma*W(t, x, y, z)**2 - Gamma - W(t, x, y, z)**2 + 1)*W(t, x, y, z)*v1(t, x, y, z)/((1 - Gamma)*(Gamma*W(t, x, y, z)**3 - Gamma*W(t, x, y, z) - 2*W(t, x, y, z)**3 + W(t, x, y, z))))*Derivative(p(t, x, y, z), x) + (-(Gamma*W(t, x, y, z)**2 - Gamma - W(t, x, y, z)**2 + 1)*W(t, x, y, z)*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**3 - Gamma*W(t, x, y, z) - 2*W(t, x, y, z)**3 + W(t, x, y, z)) + W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2 - Gamma - 2*W(t, x, y, z)**2 + 1))*Derivative(rho(t, x, y, z), x)
Derivative(v1(t, x, y, z), t) = (-(-Gamma*v1(t, x, y, z) + 2*v1(t, x, y, z))*(p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - (Gamma*v1(t, x, y, z) - v1(t, x, y, z))*(rho(t, x, y, z) - p(t, x, y, z)/(1 - Gamma))*W(t, x, y, z)/(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z)) - 2*(p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)))*Derivative(v1(t, x, y, z), x) + (-(-Gamma*v1(t, x, y, z) + 2*v1(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - (W(t, x, y, z)**2*v1(t, x, y, z)**2 + 1)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)) + (Gamma*v1(t, x, y, z) - v1(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/((1 - Gamma)*(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z))))*Derivative(p(t, x, y, z), x) + (-(-Gamma*v1(t, x, y, z) + 2*v1(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - (Gamma*v1(t, x, y, z) - v1(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z)) - W(t, x, y, z)**2*v1(t, x, y, z)**2/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)))*Derivative(rho(t, x, y, z), x)
Derivative(v2(t, x, y, z), t) = (-(-Gamma*v2(t, x, y, z) + 2*v2(t, x, y, z))*(p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - (Gamma*v2(t, x, y, z) - v2(t, x, y, z))*(rho(t, x, y, z) - p(t, x, y, z)/(1 - Gamma))*W(t, x, y, z)/(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z)) - (p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2*v2(t, x, y, z)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)))*Derivative(v1(t, x, y, z), x) + (-(-Gamma*v2(t, x, y, z) + 2*v2(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - (Gamma*v2(t, x, y, z) - v2(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z)) - W(t, x, y, z)**2*v1(t, x, y, z)*v2(t, x, y, z)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)))*Derivative(rho(t, x, y, z), x) + (-(-Gamma*v2(t, x, y, z) + 2*v2(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - W(t, x, y, z)**2*v1(t, x, y, z)*v2(t, x, y, z)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)) + (Gamma*v2(t, x, y, z) - v2(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/((1 - Gamma)*(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z))))*Derivative(p(t, x, y, z), x) - (p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)*Derivative(v2(t, x, y, z), x)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z))
Derivative(v3(t, x, y, z), t) = (-(-Gamma*v3(t, x, y, z) + 2*v3(t, x, y, z))*(p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - (Gamma*v3(t, x, y, z) - v3(t, x, y, z))*(rho(t, x, y, z) - p(t, x, y, z)/(1 - Gamma))*W(t, x, y, z)/(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z)) - (p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2*v3(t, x, y, z)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)))*Derivative(v1(t, x, y, z), x) + (-(-Gamma*v3(t, x, y, z) + 2*v3(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - (Gamma*v3(t, x, y, z) - v3(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z)) - W(t, x, y, z)**2*v1(t, x, y, z)*v3(t, x, y, z)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)))*Derivative(rho(t, x, y, z), x) + (-(-Gamma*v3(t, x, y, z) + 2*v3(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)/(Gamma*W(t, x, y, z)**2*p(t, x, y, z) + Gamma*W(t, x, y, z)**2*rho(t, x, y, z) - Gamma*p(t, x, y, z) - Gamma*rho(t, x, y, z) - 2*W(t, x, y, z)**2*p(t, x, y, z) - 2*W(t, x, y, z)**2*rho(t, x, y, z) + p(t, x, y, z) + rho(t, x, y, z)) - W(t, x, y, z)**2*v1(t, x, y, z)*v3(t, x, y, z)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z)) + (Gamma*v3(t, x, y, z) - v3(t, x, y, z))*W(t, x, y, z)*v1(t, x, y, z)/((1 - Gamma)*(Gamma*W(t, x, y, z)**3*p(t, x, y, z) + Gamma*W(t, x, y, z)**3*rho(t, x, y, z) - Gamma*W(t, x, y, z)*p(t, x, y, z) - Gamma*W(t, x, y, z)*rho(t, x, y, z) - 2*W(t, x, y, z)**3*p(t, x, y, z) - 2*W(t, x, y, z)**3*rho(t, x, y, z) + W(t, x, y, z)*p(t, x, y, z) + W(t, x, y, z)*rho(t, x, y, z))))*Derivative(p(t, x, y, z), x) - (p(t, x, y, z) + rho(t, x, y, z))*W(t, x, y, z)**2*v1(t, x, y, z)*Derivative(v3(t, x, y, z), x)/(W(t, x, y, z)**2*p(t, x, y, z) + W(t, x, y, z)**2*rho(t, x, y, z))