Commit f5103a6c authored by Radek Fučík's avatar Radek Fučík
Browse files

added txt outputs from supp materials

Changes to be committed:
	new file:   supp_d1q3_ade_txt/output_d1q3_ade_clbm1_symbolic_pde_00.txt
	new file:   supp_d1q3_ade_txt/output_d1q3_ade_mrt1_symbolic_pde_00.txt
	new file:   supp_d1q3_ade_txt/output_d1q3_ade_srt_symbolic_pde_00.txt
	new file:   supp_d1q3_nse_txt/output_d1q3_nse_clbm1_symbolic_pde_00.txt
	new file:   supp_d1q3_nse_txt/output_d1q3_nse_clbm1_symbolic_pde_01.txt
	new file:   supp_d1q3_nse_txt/output_d1q3_nse_mrt1_symbolic_pde_00.txt
	new file:   supp_d1q3_nse_txt/output_d1q3_nse_mrt1_symbolic_pde_01.txt
	new file:   supp_d1q3_nse_txt/output_d1q3_nse_srt_symbolic_pde_00.txt
	new file:   supp_d1q3_nse_txt/output_d1q3_nse_srt_symbolic_pde_01.txt
	new file:   supp_d2q5_ade_txt/output_d2q5_ade_clbm1_symbolic_pde_00.txt
	new file:   supp_d2q5_ade_txt/output_d2q5_ade_clbm2_symbolic_pde_00.txt
	new file:   supp_d2q5_ade_txt/output_d2q5_ade_mrt1_symbolic_pde_00.txt
	new file:   supp_d2q5_ade_txt/output_d2q5_ade_mrt2_symbolic_pde_00.txt
	new file:   supp_d2q5_ade_txt/output_d2q5_ade_srt_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_clbm1_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_clbm1_symbolic_pde_01.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_clbm1_symbolic_pde_02.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_clbm2_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_clbm2_symbolic_pde_01.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_clbm2_symbolic_pde_02.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_culbm1_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_culbm1_symbolic_pde_01.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_culbm1_symbolic_pde_02.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_culbm2_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_culbm2_symbolic_pde_01.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_culbm2_symbolic_pde_02.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_mrt1_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_mrt1_symbolic_pde_01.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_mrt1_symbolic_pde_02.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_mrt2_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_mrt2_symbolic_pde_01.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_mrt2_symbolic_pde_02.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_srt_symbolic_pde_00.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_srt_symbolic_pde_01.txt
	new file:   supp_d2q9_nse_txt/output_d2q9_nse_srt_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm1_symbolic_pde_00.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm1_symbolic_pde_01.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm1_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm1_symbolic_pde_03.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm2_symbolic_pde_00.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm2_symbolic_pde_01.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm2_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_clbm2_symbolic_pde_03.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm1_symbolic_pde_00.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm1_symbolic_pde_01.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm1_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm1_symbolic_pde_03.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm2_symbolic_pde_00.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm2_symbolic_pde_01.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm2_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_culbm2_symbolic_pde_03.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt1_symbolic_pde_00.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt1_symbolic_pde_01.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt1_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt1_symbolic_pde_03.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt2_symbolic_pde_00.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt2_symbolic_pde_01.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt2_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_mrt2_symbolic_pde_03.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_srt_symbolic_pde_00.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_srt_symbolic_pde_01.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_srt_symbolic_pde_02.txt
	new file:   supp_d3q27_nse_txt/output_d3q27_nse_srt_symbolic_pde_03.txt
	new file:   supp_d3q7_ade_txt/output_d3q7_ade_clbm1_symbolic_pde_00.txt
	new file:   supp_d3q7_ade_txt/output_d3q7_ade_clbm2_symbolic_pde_00.txt
	new file:   supp_d3q7_ade_txt/output_d3q7_ade_mrt1_symbolic_pde_00.txt
	new file:   supp_d3q7_ade_txt/output_d3q7_ade_mrt2_symbolic_pde_00.txt
	new file:   supp_d3q7_ade_txt/output_d3q7_ade_srt_symbolic_pde_00.txt
parent da3f1cdf
Loading
Loading
Loading
Loading
+43 −0
Original line number Original line Diff line number Diff line
Complete expression of PDE that can be directly imported to Maxima:
D[0](rho)(t,x,y,z)-1/2*D[1](u)(t,x,y,z)^2*(2*dl^2*dt^(-1)*rho(t,x,y,z)*omega_2^(-1)-dl^2*dt^(-1)*rho(t,x,y,z))-1/2*D[0,0,0,1](u)(t,x,y,z)*(dl*dt^2*rho(t,x,y,z)*omega_2^(-1)+2*dl*dt^2*rho(t,x,y,z)*omega_2^(-3)-3*dl*dt^2*rho(t,x,y,z)*omega_2^(-2))+dl*D[1](u)(t,x,y,z)*dt^(-1)*rho(t,x,y,z)+1/6*D[0,1,1](u)(t,x,y,z)*(12*dl^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-2)-12*dl^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)+dl^2*rho(t,x,y,z)*u(t,x,y,z))-1/2*D[1,1](u)(t,x,y,z)*(2*dl^2*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)-dl^2*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z))+3/2*D[0,0,1,1](u)(t,x,y,z)*(3*dl^2*dt*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-2)-2*dl^2*dt*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-3)-dl^2*dt*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1))+1/12*D[1,1,1,1](u)(t,x,y,z)*(2*dl^4*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3+42*dl^4*omega_3^(-2)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3-84*dl^4*omega_3^(-2)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3*omega_2^(-1)+2*dl^4*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3*omega_2^(-1)-12*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-2)+24*dl^4*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3*omega_2^(-1)-30*dl^4*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-2)+6*dl^4*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-2)-18*dl^4*omega_3^(-2)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)-dl^4*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)+24*dl^4*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-3)-12*dl^4*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-2)-24*dl^4*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)-24*dl^4*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3+2*dl^4*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)+36*dl^4*omega_3^(-2)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)-60*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)-2*dl^4*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)+12*dl^4*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)-12*dl^4*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3*omega_2^(-3)+72*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)*omega_2^(-1)+30*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)+dl^4*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)+60*dl^4*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3*omega_2^(-2)-24*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*rho(t,x,y,z)*u(t,x,y,z)-12*dl^4*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^3*omega_2^(-2))+1/12*D[0,0,1](u)(t,x,y,z)*(12*dl*dt*rho(t,x,y,z)*omega_2^(-2)-12*dl*dt*rho(t,x,y,z)*omega_2^(-1)+dl*dt*rho(t,x,y,z))-1/2*D[0,1](u)(t,x,y,z)*(2*dl*rho(t,x,y,z)*omega_2^(-1)-dl*rho(t,x,y,z))-1/24*D[1,1,1,1](rho)(t,x,y,z)*(30*dl^4*omega_3^(-1)*dt^(-1)*u(t,x,y,z)^4-24*dl^4*dt^(-1)*c_s^4*omega_2^(-3)+24*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*omega_2^(-2)-72*dl^4*omega_3^(-1)*dt^(-1)*u(t,x,y,z)^4*omega_2^(-1)+36*dl^4*omega_3^(-2)*dt^(-1)*u(t,x,y,z)^2-12*dl^4*dt^(-1)*u(t,x,y,z)^2*omega_2^(-1)-72*dl^4*omega_3^(-2)*dt^(-1)*u(t,x,y,z)^2*omega_2^(-1)+48*dl^4*dt^(-1)*c_s^4*omega_2^(-2)-dl^4*dt^(-1)*c_s^2+24*dl^4*omega_3^(-1)*dt^(-1)*c_s^4*omega_2^(-1)+72*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*u(t,x,y,z)^2-3*dl^4*dt^(-1)*u(t,x,y,z)^4+36*dl^4*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-2)-72*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-2)+8*dl^4*dt^(-1)*c_s^2*omega_2^(-1)+6*dl^4*omega_3^(-1)*dt^(-1)*c_s^2+72*dl^4*omega_3^(-2)*dt^(-1)*u(t,x,y,z)^4*omega_2^(-1)+12*dl^4*dt^(-1)*u(t,x,y,z)^4*omega_2^(-1)-12*dl^4*dt^(-1)*c_s^2*omega_2^(-2)-144*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-1)+3*dl^4*dt^(-1)*u(t,x,y,z)^2-6*dl^4*omega_3^(-1)*dt^(-1)*c_s^4+216*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-1)+12*dl^4*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-1)-24*dl^4*omega_3^(-1)*dt^(-1)*c_s^4*omega_2^(-2)+72*dl^4*omega_3^(-1)*dt^(-1)*u(t,x,y,z)^2*omega_2^(-1)-36*dl^4*omega_3^(-2)*dt^(-1)*u(t,x,y,z)^4-30*dl^4*omega_3^(-1)*dt^(-1)*u(t,x,y,z)^2-24*dl^4*dt^(-1)*c_s^4*omega_2^(-1)-6*dl^4*dt^(-1)*c_s^2*u(t,x,y,z)^2-108*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*u(t,x,y,z)^2-24*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*omega_2^(-1)+3*dl^4*dt^(-1)*c_s^4)-1/2*D[1](rho)(t,x,y,z)*(2*dl^2*dt^(-1)*u(t,x,y,z)*omega_2^(-1)-dl^2*dt^(-1)*u(t,x,y,z))*D[1](u)(t,x,y,z)+1/6*D[1,1,1](rho)(t,x,y,z)*(dl^3*dt^(-1)*u(t,x,y,z)+3*dl^3*omega_3^(-1)*dt^(-1)*u(t,x,y,z)^3-3*dl^3*dt^(-1)*u(t,x,y,z)*omega_2^(-1)-6*dl^3*omega_3^(-1)*dt^(-1)*u(t,x,y,z)^3*omega_2^(-1)-3*dl^3*omega_3^(-1)*dt^(-1)*u(t,x,y,z)-dl^3*dt^(-1)*u(t,x,y,z)^3-3*dl^3*dt^(-1)*c_s^2*u(t,x,y,z)+9*dl^3*dt^(-1)*c_s^2*u(t,x,y,z)*omega_2^(-1)+6*dl^3*omega_3^(-1)*dt^(-1)*u(t,x,y,z)*omega_2^(-1)+3*dl^3*dt^(-1)*u(t,x,y,z)^3*omega_2^(-1)-18*dl^3*omega_3^(-1)*dt^(-1)*c_s^2*u(t,x,y,z)*omega_2^(-1)+9*dl^3*omega_3^(-1)*dt^(-1)*c_s^2*u(t,x,y,z))+D[1](rho)(t,x,y,z)*dl*dt^(-1)*u(t,x,y,z)+1/2*D[1](rho)(t,x,y,z)*(dl-2*dl*omega_2^(-1))*D[0](u)(t,x,y,z)-1/12*D[1,1,1](u)(t,x,y,z)*(12*dl^3*dt^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-2)-2*dl^3*dt^(-1)*rho(t,x,y,z)+6*dl^3*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)+6*dl^3*dt^(-1)*rho(t,x,y,z)*omega_2^(-1)+3*dl^3*dt^(-1)*c_s^2*rho(t,x,y,z)-12*dl^3*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-2)-18*dl^3*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2-6*dl^3*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-1)+12*dl^3*omega_3^(-1)*dt^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-1)-18*dl^3*dt^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-1)+5*dl^3*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2-12*dl^3*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*omega_2^(-1)-6*dl^3*omega_3^(-1)*dt^(-1)*c_s^2*rho(t,x,y,z)+36*dl^3*omega_3^(-1)*dt^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-1))+1/2*D[1,1](rho)(t,x,y,z)*(dl^2*dt^(-1)*c_s^2-2*dl^2*dt^(-1)*c_s^2*omega_2^(-1))-1/12*D[0,1,1,1](u)(t,x,y,z)*(3*dl^3*rho(t,x,y,z)*u(t,x,y,z)^2+24*dl^3*omega_3^(-1)*rho(t,x,y,z)*omega_2^(-2)-15*dl^3*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-1)-24*dl^3*omega_3^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-2)+12*dl^3*omega_3^(-2)*rho(t,x,y,z)*omega_2^(-1)+6*dl^3*omega_3^(-2)*c_s^2*rho(t,x,y,z)-36*dl^3*omega_3^(-2)*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-1)-72*dl^3*omega_3^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-2)-25*dl^3*c_s^2*rho(t,x,y,z)*omega_2^(-1)-12*dl^3*rho(t,x,y,z)*omega_2^(-2)-dl^3*rho(t,x,y,z)-27*dl^3*omega_3^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2+9*dl^3*omega_3^(-1)*rho(t,x,y,z)-12*dl^3*omega_3^(-2)*c_s^2*rho(t,x,y,z)*omega_2^(-1)-24*dl^3*c_s^2*rho(t,x,y,z)*omega_2^(-3)-6*dl^3*omega_3^(-2)*rho(t,x,y,z)+18*dl^3*omega_3^(-2)*rho(t,x,y,z)*u(t,x,y,z)^2+108*dl^3*omega_3^(-1)*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-1)+2*dl^3*c_s^2*rho(t,x,y,z)+11*dl^3*rho(t,x,y,z)*omega_2^(-1)+48*dl^3*c_s^2*rho(t,x,y,z)*omega_2^(-2)+36*dl^3*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-3)-18*dl^3*rho(t,x,y,z)*u(t,x,y,z)^2*omega_2^(-2)+36*dl^3*omega_3^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-1)-9*dl^3*omega_3^(-1)*c_s^2*rho(t,x,y,z)-36*dl^3*omega_3^(-1)*rho(t,x,y,z)*omega_2^(-1))=0


List of all coefficients of each partial derivatives:
D[1](rho)
	dl*dt^(-1)*u
D[1](u)
	dl*dt^(-1)*rho
D[0](rho)
	1
D[1,1](rho)
	1/2*dl^2*dt^(-1)*c_s^2-dl^2*dt^(-1)*c_s^2*omega_2^(-1)
D[1,1](u)
	-dl^2*dt^(-1)*rho*u*omega_2^(-1)+1/2*dl^2*dt^(-1)*rho*u
D[0,1](u)
	-dl*rho*omega_2^(-1)+1/2*dl*rho
D[1,1,1](rho)
	1/6*dl^3*dt^(-1)*u+1/2*dl^3*omega_3^(-1)*dt^(-1)*u^3-1/2*dl^3*dt^(-1)*u*omega_2^(-1)-dl^3*omega_3^(-1)*dt^(-1)*u^3*omega_2^(-1)-1/2*dl^3*omega_3^(-1)*dt^(-1)*u-1/6*dl^3*dt^(-1)*u^3-1/2*dl^3*dt^(-1)*c_s^2*u+3/2*dl^3*dt^(-1)*c_s^2*u*omega_2^(-1)+dl^3*omega_3^(-1)*dt^(-1)*u*omega_2^(-1)+1/2*dl^3*dt^(-1)*u^3*omega_2^(-1)-3*dl^3*omega_3^(-1)*dt^(-1)*c_s^2*u*omega_2^(-1)+3/2*dl^3*omega_3^(-1)*dt^(-1)*c_s^2*u
D[1,1,1](u)
	-dl^3*dt^(-1)*c_s^2*rho*omega_2^(-2)+1/6*dl^3*dt^(-1)*rho-1/2*dl^3*omega_3^(-1)*dt^(-1)*rho-1/2*dl^3*dt^(-1)*rho*omega_2^(-1)-1/4*dl^3*dt^(-1)*c_s^2*rho+dl^3*dt^(-1)*rho*u^2*omega_2^(-2)+3/2*dl^3*omega_3^(-1)*dt^(-1)*rho*u^2+1/2*dl^3*dt^(-1)*rho*u^2*omega_2^(-1)-dl^3*omega_3^(-1)*dt^(-1)*c_s^2*rho*omega_2^(-1)+3/2*dl^3*dt^(-1)*c_s^2*rho*omega_2^(-1)-5/12*dl^3*dt^(-1)*rho*u^2+dl^3*omega_3^(-1)*dt^(-1)*rho*omega_2^(-1)+1/2*dl^3*omega_3^(-1)*dt^(-1)*c_s^2*rho-3*dl^3*omega_3^(-1)*dt^(-1)*rho*u^2*omega_2^(-1)
D[0,1,1](u)
	2*dl^2*rho*u*omega_2^(-2)-2*dl^2*rho*u*omega_2^(-1)+1/6*dl^2*rho*u
D[0,0,1](u)
	dl*dt*rho*omega_2^(-2)-dl*dt*rho*omega_2^(-1)+1/12*dl*dt*rho
D[1,1,1,1](rho)
	-5/4*dl^4*omega_3^(-1)*dt^(-1)*u^4+dl^4*dt^(-1)*c_s^4*omega_2^(-3)-dl^4*omega_3^(-1)*dt^(-1)*c_s^2*omega_2^(-2)+3*dl^4*omega_3^(-1)*dt^(-1)*u^4*omega_2^(-1)-3/2*dl^4*omega_3^(-2)*dt^(-1)*u^2+1/2*dl^4*dt^(-1)*u^2*omega_2^(-1)+3*dl^4*omega_3^(-2)*dt^(-1)*u^2*omega_2^(-1)-2*dl^4*dt^(-1)*c_s^4*omega_2^(-2)+1/24*dl^4*dt^(-1)*c_s^2-dl^4*omega_3^(-1)*dt^(-1)*c_s^4*omega_2^(-1)-3*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*u^2+1/8*dl^4*dt^(-1)*u^4-3/2*dl^4*dt^(-1)*c_s^2*u^2*omega_2^(-2)+3*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*u^2*omega_2^(-2)-1/3*dl^4*dt^(-1)*c_s^2*omega_2^(-1)-1/4*dl^4*omega_3^(-1)*dt^(-1)*c_s^2-3*dl^4*omega_3^(-2)*dt^(-1)*u^4*omega_2^(-1)-1/2*dl^4*dt^(-1)*u^4*omega_2^(-1)+1/2*dl^4*dt^(-1)*c_s^2*omega_2^(-2)+6*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*u^2*omega_2^(-1)-1/8*dl^4*dt^(-1)*u^2+1/4*dl^4*omega_3^(-1)*dt^(-1)*c_s^4-9*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*u^2*omega_2^(-1)-1/2*dl^4*dt^(-1)*c_s^2*u^2*omega_2^(-1)+dl^4*omega_3^(-1)*dt^(-1)*c_s^4*omega_2^(-2)-3*dl^4*omega_3^(-1)*dt^(-1)*u^2*omega_2^(-1)+3/2*dl^4*omega_3^(-2)*dt^(-1)*u^4+5/4*dl^4*omega_3^(-1)*dt^(-1)*u^2+dl^4*dt^(-1)*c_s^4*omega_2^(-1)+1/4*dl^4*dt^(-1)*c_s^2*u^2+9/2*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*u^2+dl^4*omega_3^(-1)*dt^(-1)*c_s^2*omega_2^(-1)-1/8*dl^4*dt^(-1)*c_s^4
D[1,1,1,1](u)
	1/6*dl^4*dt^(-1)*rho*u^3+7/2*dl^4*omega_3^(-2)*dt^(-1)*rho*u^3-7*dl^4*omega_3^(-2)*dt^(-1)*rho*u^3*omega_2^(-1)+1/6*dl^4*dt^(-1)*rho*u^3*omega_2^(-1)-dl^4*omega_3^(-1)*dt^(-1)*c_s^2*rho*u*omega_2^(-2)+2*dl^4*omega_3^(-1)*dt^(-1)*rho*u^3*omega_2^(-1)-5/2*dl^4*dt^(-1)*c_s^2*rho*u*omega_2^(-2)+1/2*dl^4*dt^(-1)*rho*u*omega_2^(-2)-3/2*dl^4*omega_3^(-2)*dt^(-1)*rho*u-1/12*dl^4*dt^(-1)*rho*u+2*dl^4*dt^(-1)*c_s^2*rho*u*omega_2^(-3)-dl^4*omega_3^(-1)*dt^(-1)*rho*u*omega_2^(-2)-2*dl^4*omega_3^(-1)*dt^(-1)*rho*u*omega_2^(-1)-2*dl^4*omega_3^(-1)*dt^(-1)*rho*u^3+1/6*dl^4*dt^(-1)*rho*u*omega_2^(-1)+3*dl^4*omega_3^(-2)*dt^(-1)*rho*u*omega_2^(-1)-5*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*rho*u*omega_2^(-1)-1/6*dl^4*dt^(-1)*c_s^2*rho*u*omega_2^(-1)+dl^4*omega_3^(-1)*dt^(-1)*rho*u-dl^4*dt^(-1)*rho*u^3*omega_2^(-3)+6*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*rho*u*omega_2^(-1)+5/2*dl^4*omega_3^(-2)*dt^(-1)*c_s^2*rho*u+1/12*dl^4*dt^(-1)*c_s^2*rho*u+5*dl^4*omega_3^(-1)*dt^(-1)*rho*u^3*omega_2^(-2)-2*dl^4*omega_3^(-1)*dt^(-1)*c_s^2*rho*u-dl^4*dt^(-1)*rho*u^3*omega_2^(-2)
D[0,1,1,1](u)
	-1/4*dl^3*rho*u^2-2*dl^3*omega_3^(-1)*rho*omega_2^(-2)+5/4*dl^3*rho*u^2*omega_2^(-1)+2*dl^3*omega_3^(-1)*c_s^2*rho*omega_2^(-2)-dl^3*omega_3^(-2)*rho*omega_2^(-1)-1/2*dl^3*omega_3^(-2)*c_s^2*rho+3*dl^3*omega_3^(-2)*rho*u^2*omega_2^(-1)+6*dl^3*omega_3^(-1)*rho*u^2*omega_2^(-2)+25/12*dl^3*c_s^2*rho*omega_2^(-1)+dl^3*rho*omega_2^(-2)+1/12*dl^3*rho+9/4*dl^3*omega_3^(-1)*rho*u^2-3/4*dl^3*omega_3^(-1)*rho+dl^3*omega_3^(-2)*c_s^2*rho*omega_2^(-1)+2*dl^3*c_s^2*rho*omega_2^(-3)+1/2*dl^3*omega_3^(-2)*rho-3/2*dl^3*omega_3^(-2)*rho*u^2-9*dl^3*omega_3^(-1)*rho*u^2*omega_2^(-1)-1/6*dl^3*c_s^2*rho-11/12*dl^3*rho*omega_2^(-1)-4*dl^3*c_s^2*rho*omega_2^(-2)-3*dl^3*rho*u^2*omega_2^(-3)+3/2*dl^3*rho*u^2*omega_2^(-2)-3*dl^3*omega_3^(-1)*c_s^2*rho*omega_2^(-1)+3/4*dl^3*omega_3^(-1)*c_s^2*rho+3*dl^3*omega_3^(-1)*rho*omega_2^(-1)
D[0,0,1,1](u)
	9/2*dl^2*dt*rho*u*omega_2^(-2)-3*dl^2*dt*rho*u*omega_2^(-3)-3/2*dl^2*dt*rho*u*omega_2^(-1)
D[0,0,0,1](u)
	-1/2*dl*dt^2*rho*omega_2^(-1)-dl*dt^2*rho*omega_2^(-3)+3/2*dl*dt^2*rho*omega_2^(-2)
D[1](rho)  *  D[1](u)
	-dl^2*dt^(-1)*u*omega_2^(-1)+1/2*dl^2*dt^(-1)*u
D[1](u)  *  D[1](u)
	-dl^2*dt^(-1)*rho*omega_2^(-1)+1/2*dl^2*dt^(-1)*rho
D[1](rho)  *  D[0](u)
	1/2*dl-dl*omega_2^(-1)

+43 −0
Original line number Original line Diff line number Diff line
Complete expression of PDE that can be directly imported to Maxima:
1/2*(dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*dl^2-2*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^2)*D[1,1](u)(t,x,y,z)-1/2*D[0](u)(t,x,y,z)*D[1](rho)(t,x,y,z)*(2*omega_2^(-1)*dl-dl)-1/2*D[1,1](rho)(t,x,y,z)*(2*dt^(-1)*c_s^2*omega_2^(-1)*dl^2-dt^(-1)*c_s^2*dl^2)+1/12*D[0,0,1](u)(t,x,y,z)*(dt*rho(t,x,y,z)*dl+12*dt*rho(t,x,y,z)*omega_2^(-2)*dl-12*dt*rho(t,x,y,z)*omega_2^(-1)*dl)+1/12*D[1,1,1](u)(t,x,y,z)*(6*dt^(-1)*u(t,x,y,z)^2*rho(t,x,y,z)*dl^3*omega_3^(-1)+12*dt^(-1)*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)+18*dt^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-1)*dl^3+6*dt^(-1)*c_s^2*rho(t,x,y,z)*dl^3*omega_3^(-1)+18*dt^(-1)*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-1)*dl^3-12*dt^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)+2*dt^(-1)*rho(t,x,y,z)*dl^3-5*dt^(-1)*u(t,x,y,z)^2*rho(t,x,y,z)*dl^3-6*dt^(-1)*rho(t,x,y,z)*omega_2^(-1)*dl^3-12*dt^(-1)*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-2)*dl^3-3*dt^(-1)*c_s^2*rho(t,x,y,z)*dl^3-12*dt^(-1)*c_s^2*rho(t,x,y,z)*omega_2^(-2)*dl^3-12*dt^(-1)*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)-6*dt^(-1)*rho(t,x,y,z)*dl^3*omega_3^(-1))+1/2*(dt^(-1)*u(t,x,y,z)*dl^2-2*dt^(-1)*u(t,x,y,z)*omega_2^(-1)*dl^2)*D[1](u)(t,x,y,z)*D[1](rho)(t,x,y,z)-1/12*D[1,1,1,1](u)(t,x,y,z)*(12*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*omega_2^(-1)*dl^4*omega_3^(-2)-dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*dl^4-8*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^4+12*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^4*omega_3^(-2)+6*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*dl^4*omega_3^(-1)+24*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^4*omega_3^(-1)-24*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*omega_2^(-2)*dl^4+12*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-2)*dl^4*omega_3^(-1)+12*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*omega_2^(-2)*dl^4*omega_3^(-1)+12*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*omega_2^(-3)*dl^4-6*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*dl^4*omega_3^(-1)-6*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*dl^4*omega_3^(-2)+20*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^4-42*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-2)*dl^4+6*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*dl^4*omega_3^(-2)+6*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*dl^4*omega_3^(-1)+dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*dl^4+24*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-3)*dl^4-2*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*dl^4-12*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-2)*dl^4*omega_3^(-1)-12*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^4*omega_3^(-2)-6*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*dl^4*omega_3^(-2)-24*dt^(-1)*c_s^2*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^4*omega_3^(-1)+6*dt^(-1)*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-2)*dl^4+16*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*omega_2^(-1)*dl^4-24*dt^(-1)*u(t,x,y,z)^3*rho(t,x,y,z)*omega_2^(-1)*dl^4*omega_3^(-1))-1/24*D[1,1,1,1](rho)(t,x,y,z)*(8*dt^(-1)*c_s^2*omega_2^(-1)*dl^4-12*dt^(-1)*c_s^2*u(t,x,y,z)^2*dl^4*omega_3^(-2)-24*dt^(-1)*c_s^2*omega_2^(-1)*dl^4*omega_3^(-1)+18*dt^(-1)*u(t,x,y,z)^4*dl^4*omega_3^(-1)+48*dt^(-1)*u(t,x,y,z)^4*omega_2^(-2)*dl^4*omega_3^(-1)+3*dt^(-1)*u(t,x,y,z)^2*dl^4-6*dt^(-1)*c_s^4*dl^4*omega_3^(-1)-48*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-1)*dl^4*omega_3^(-1)+24*dt^(-1)*u(t,x,y,z)^4*omega_2^(-1)*dl^4*omega_3^(-2)-24*dt^(-1)*u(t,x,y,z)^2*omega_2^(-1)*dl^4+96*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-3)*dl^4-dt^(-1)*c_s^2*dl^4-24*dt^(-1)*c_s^4*omega_2^(-1)*dl^4-18*dt^(-1)*u(t,x,y,z)^2*dl^4*omega_3^(-1)-24*dt^(-1)*c_s^4*omega_2^(-2)*dl^4*omega_3^(-1)+72*dt^(-1)*u(t,x,y,z)^2*omega_2^(-1)*dl^4*omega_3^(-1)+6*dt^(-1)*c_s^2*dl^4*omega_3^(-1)+24*dt^(-1)*u(t,x,y,z)^4*omega_2^(-1)*dl^4-156*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-2)*dl^4-3*dt^(-1)*u(t,x,y,z)^4*dl^4+24*dt^(-1)*c_s^4*omega_2^(-1)*dl^4*omega_3^(-1)-24*dt^(-1)*u(t,x,y,z)^4*omega_2^(-2)*dl^4+72*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-1)*dl^4+3*dt^(-1)*c_s^4*dl^4-48*dt^(-1)*u(t,x,y,z)^2*omega_2^(-2)*dl^4*omega_3^(-1)-24*dt^(-1)*u(t,x,y,z)^2*omega_2^(-1)*dl^4*omega_3^(-2)+48*dt^(-1)*c_s^4*omega_2^(-2)*dl^4+12*dt^(-1)*u(t,x,y,z)^2*dl^4*omega_3^(-2)+24*dt^(-1)*c_s^2*omega_2^(-2)*dl^4*omega_3^(-1)+24*dt^(-1)*u(t,x,y,z)^2*omega_2^(-2)*dl^4+24*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-1)*dl^4*omega_3^(-2)-72*dt^(-1)*u(t,x,y,z)^4*omega_2^(-1)*dl^4*omega_3^(-1)-12*dt^(-1)*c_s^2*omega_2^(-2)*dl^4+24*dt^(-1)*c_s^2*u(t,x,y,z)^2*omega_2^(-2)*dl^4*omega_3^(-1)+12*dt^(-1)*c_s^2*u(t,x,y,z)^2*dl^4*omega_3^(-1)-24*dt^(-1)*c_s^4*omega_2^(-3)*dl^4-12*dt^(-1)*u(t,x,y,z)^4*dl^4*omega_3^(-2)-6*dt^(-1)*c_s^2*u(t,x,y,z)^2*dl^4)+1/6*(u(t,x,y,z)*rho(t,x,y,z)*dl^2+12*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^2*omega_3^(-1)-6*u(t,x,y,z)*rho(t,x,y,z)*dl^2*omega_3^(-1)-6*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^2)*D[0,1,1](u)(t,x,y,z)+dt^(-1)*rho(t,x,y,z)*D[1](u)(t,x,y,z)*dl-1/2*D[0,0,0,1](u)(t,x,y,z)*(dt^2*rho(t,x,y,z)*omega_2^(-1)*dl-3*dt^2*rho(t,x,y,z)*omega_2^(-2)*dl+2*dt^2*rho(t,x,y,z)*omega_2^(-3)*dl)+1/2*(dt^(-1)*rho(t,x,y,z)*dl^2-2*dt^(-1)*rho(t,x,y,z)*omega_2^(-1)*dl^2)*D[1](u)(t,x,y,z)^2+dt^(-1)*u(t,x,y,z)*D[1](rho)(t,x,y,z)*dl+1/2*(rho(t,x,y,z)*dl-2*rho(t,x,y,z)*omega_2^(-1)*dl)*D[0,1](u)(t,x,y,z)+D[0](rho)(t,x,y,z)+1/2*(2*dt*u(t,x,y,z)*rho(t,x,y,z)*dl^2*omega_3^(-2)-dt*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-2)*dl^2-4*dt*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^2*omega_3^(-2)+2*dt*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-3)*dl^2-4*dt*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-2)*dl^2*omega_3^(-1)+8*dt*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^2*omega_3^(-1)-2*dt*u(t,x,y,z)*rho(t,x,y,z)*dl^2*omega_3^(-1)-dt*u(t,x,y,z)*rho(t,x,y,z)*omega_2^(-1)*dl^2)*D[0,0,1,1](u)(t,x,y,z)+1/6*(3*dt^(-1)*u(t,x,y,z)^3*dl^3*omega_3^(-1)-3*dt^(-1)*u(t,x,y,z)*omega_2^(-1)*dl^3+6*dt^(-1)*u(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)+dt^(-1)*u(t,x,y,z)*dl^3-12*dt^(-1)*c_s^2*u(t,x,y,z)*omega_2^(-2)*dl^3-6*dt^(-1)*c_s^2*u(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)+3*dt^(-1)*u(t,x,y,z)^3*omega_2^(-1)*dl^3-3*dt^(-1)*u(t,x,y,z)*dl^3*omega_3^(-1)-dt^(-1)*u(t,x,y,z)^3*dl^3+15*dt^(-1)*c_s^2*u(t,x,y,z)*omega_2^(-1)*dl^3+3*dt^(-1)*c_s^2*u(t,x,y,z)*dl^3*omega_3^(-1)-6*dt^(-1)*u(t,x,y,z)^3*omega_2^(-1)*dl^3*omega_3^(-1)-3*dt^(-1)*c_s^2*u(t,x,y,z)*dl^3)*D[1,1,1](rho)(t,x,y,z)+1/12*D[0,1,1,1](u)(t,x,y,z)*(15*u(t,x,y,z)^2*rho(t,x,y,z)*dl^3*omega_3^(-1)+36*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)+48*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-2)*dl^3*omega_3^(-1)+6*rho(t,x,y,z)*dl^3*omega_3^(-2)+12*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-2)+25*c_s^2*rho(t,x,y,z)*omega_2^(-1)*dl^3+9*c_s^2*rho(t,x,y,z)*dl^3*omega_3^(-1)-36*c_s^2*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)+12*rho(t,x,y,z)*omega_2^(-2)*dl^3+27*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-1)*dl^3+rho(t,x,y,z)*dl^3-3*u(t,x,y,z)^2*rho(t,x,y,z)*dl^3+24*c_s^2*rho(t,x,y,z)*omega_2^(-3)*dl^3-42*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-2)*dl^3-11*rho(t,x,y,z)*omega_2^(-1)*dl^3+12*c_s^2*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-2)-6*c_s^2*rho(t,x,y,z)*dl^3*omega_3^(-2)+24*c_s^2*rho(t,x,y,z)*omega_2^(-2)*dl^3*omega_3^(-1)-2*c_s^2*rho(t,x,y,z)*dl^3-48*c_s^2*rho(t,x,y,z)*omega_2^(-2)*dl^3+12*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-3)*dl^3-60*u(t,x,y,z)^2*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-1)-24*rho(t,x,y,z)*omega_2^(-2)*dl^3*omega_3^(-1)-12*rho(t,x,y,z)*omega_2^(-1)*dl^3*omega_3^(-2)-9*rho(t,x,y,z)*dl^3*omega_3^(-1)-6*u(t,x,y,z)^2*rho(t,x,y,z)*dl^3*omega_3^(-2))=0


List of all coefficients of each partial derivatives:
D[1](rho)
	dt^(-1)*u*dl
D[1](u)
	dt^(-1)*rho*dl
D[0](rho)
	1
D[1,1](rho)
	-dt^(-1)*c_s^2*omega_2^(-1)*dl^2+1/2*dt^(-1)*c_s^2*dl^2
D[1,1](u)
	1/2*dt^(-1)*u*rho*dl^2-dt^(-1)*u*rho*omega_2^(-1)*dl^2
D[0,1](u)
	1/2*rho*dl-rho*omega_2^(-1)*dl
D[1,1,1](rho)
	1/2*dt^(-1)*u^3*dl^3*omega_3^(-1)-1/2*dt^(-1)*u*omega_2^(-1)*dl^3+dt^(-1)*u*omega_2^(-1)*dl^3*omega_3^(-1)+1/6*dt^(-1)*u*dl^3-2*dt^(-1)*c_s^2*u*omega_2^(-2)*dl^3-dt^(-1)*c_s^2*u*omega_2^(-1)*dl^3*omega_3^(-1)+1/2*dt^(-1)*u^3*omega_2^(-1)*dl^3-1/2*dt^(-1)*u*dl^3*omega_3^(-1)-1/6*dt^(-1)*u^3*dl^3+5/2*dt^(-1)*c_s^2*u*omega_2^(-1)*dl^3+1/2*dt^(-1)*c_s^2*u*dl^3*omega_3^(-1)-dt^(-1)*u^3*omega_2^(-1)*dl^3*omega_3^(-1)-1/2*dt^(-1)*c_s^2*u*dl^3
D[1,1,1](u)
	1/2*dt^(-1)*u^2*rho*dl^3*omega_3^(-1)+dt^(-1)*rho*omega_2^(-1)*dl^3*omega_3^(-1)+3/2*dt^(-1)*c_s^2*rho*omega_2^(-1)*dl^3+1/2*dt^(-1)*c_s^2*rho*dl^3*omega_3^(-1)+3/2*dt^(-1)*u^2*rho*omega_2^(-1)*dl^3-dt^(-1)*c_s^2*rho*omega_2^(-1)*dl^3*omega_3^(-1)+1/6*dt^(-1)*rho*dl^3-5/12*dt^(-1)*u^2*rho*dl^3-1/2*dt^(-1)*rho*omega_2^(-1)*dl^3-dt^(-1)*u^2*rho*omega_2^(-2)*dl^3-1/4*dt^(-1)*c_s^2*rho*dl^3-dt^(-1)*c_s^2*rho*omega_2^(-2)*dl^3-dt^(-1)*u^2*rho*omega_2^(-1)*dl^3*omega_3^(-1)-1/2*dt^(-1)*rho*dl^3*omega_3^(-1)
D[0,1,1](u)
	1/6*u*rho*dl^2+2*u*rho*omega_2^(-1)*dl^2*omega_3^(-1)-u*rho*dl^2*omega_3^(-1)-u*rho*omega_2^(-1)*dl^2
D[0,0,1](u)
	1/12*dt*rho*dl+dt*rho*omega_2^(-2)*dl-dt*rho*omega_2^(-1)*dl
D[1,1,1,1](rho)
	-1/3*dt^(-1)*c_s^2*omega_2^(-1)*dl^4+1/2*dt^(-1)*c_s^2*u^2*dl^4*omega_3^(-2)+dt^(-1)*c_s^2*omega_2^(-1)*dl^4*omega_3^(-1)-3/4*dt^(-1)*u^4*dl^4*omega_3^(-1)-2*dt^(-1)*u^4*omega_2^(-2)*dl^4*omega_3^(-1)-1/8*dt^(-1)*u^2*dl^4+1/4*dt^(-1)*c_s^4*dl^4*omega_3^(-1)+2*dt^(-1)*c_s^2*u^2*omega_2^(-1)*dl^4*omega_3^(-1)-dt^(-1)*u^4*omega_2^(-1)*dl^4*omega_3^(-2)+dt^(-1)*u^2*omega_2^(-1)*dl^4-4*dt^(-1)*c_s^2*u^2*omega_2^(-3)*dl^4+1/24*dt^(-1)*c_s^2*dl^4+dt^(-1)*c_s^4*omega_2^(-1)*dl^4+3/4*dt^(-1)*u^2*dl^4*omega_3^(-1)+dt^(-1)*c_s^4*omega_2^(-2)*dl^4*omega_3^(-1)-3*dt^(-1)*u^2*omega_2^(-1)*dl^4*omega_3^(-1)-1/4*dt^(-1)*c_s^2*dl^4*omega_3^(-1)-dt^(-1)*u^4*omega_2^(-1)*dl^4+13/2*dt^(-1)*c_s^2*u^2*omega_2^(-2)*dl^4+1/8*dt^(-1)*u^4*dl^4-dt^(-1)*c_s^4*omega_2^(-1)*dl^4*omega_3^(-1)+dt^(-1)*u^4*omega_2^(-2)*dl^4-3*dt^(-1)*c_s^2*u^2*omega_2^(-1)*dl^4-1/8*dt^(-1)*c_s^4*dl^4+2*dt^(-1)*u^2*omega_2^(-2)*dl^4*omega_3^(-1)+dt^(-1)*u^2*omega_2^(-1)*dl^4*omega_3^(-2)-2*dt^(-1)*c_s^4*omega_2^(-2)*dl^4-1/2*dt^(-1)*u^2*dl^4*omega_3^(-2)-dt^(-1)*c_s^2*omega_2^(-2)*dl^4*omega_3^(-1)-dt^(-1)*u^2*omega_2^(-2)*dl^4-dt^(-1)*c_s^2*u^2*omega_2^(-1)*dl^4*omega_3^(-2)+3*dt^(-1)*u^4*omega_2^(-1)*dl^4*omega_3^(-1)+1/2*dt^(-1)*c_s^2*omega_2^(-2)*dl^4-dt^(-1)*c_s^2*u^2*omega_2^(-2)*dl^4*omega_3^(-1)-1/2*dt^(-1)*c_s^2*u^2*dl^4*omega_3^(-1)+dt^(-1)*c_s^4*omega_2^(-3)*dl^4+1/2*dt^(-1)*u^4*dl^4*omega_3^(-2)+1/4*dt^(-1)*c_s^2*u^2*dl^4
D[1,1,1,1](u)
	-dt^(-1)*u^3*rho*omega_2^(-1)*dl^4*omega_3^(-2)+1/12*dt^(-1)*c_s^2*u*rho*dl^4+2/3*dt^(-1)*u*rho*omega_2^(-1)*dl^4-dt^(-1)*c_s^2*u*rho*omega_2^(-1)*dl^4*omega_3^(-2)-1/2*dt^(-1)*c_s^2*u*rho*dl^4*omega_3^(-1)-2*dt^(-1)*u*rho*omega_2^(-1)*dl^4*omega_3^(-1)+2*dt^(-1)*u^3*rho*omega_2^(-2)*dl^4-dt^(-1)*c_s^2*u*rho*omega_2^(-2)*dl^4*omega_3^(-1)-dt^(-1)*u^3*rho*omega_2^(-2)*dl^4*omega_3^(-1)-dt^(-1)*u^3*rho*omega_2^(-3)*dl^4+1/2*dt^(-1)*u*rho*dl^4*omega_3^(-1)+1/2*dt^(-1)*u^3*rho*dl^4*omega_3^(-2)-5/3*dt^(-1)*c_s^2*u*rho*omega_2^(-1)*dl^4+7/2*dt^(-1)*c_s^2*u*rho*omega_2^(-2)*dl^4-1/2*dt^(-1)*u*rho*dl^4*omega_3^(-2)-1/2*dt^(-1)*u^3*rho*dl^4*omega_3^(-1)-1/12*dt^(-1)*u*rho*dl^4-2*dt^(-1)*c_s^2*u*rho*omega_2^(-3)*dl^4+1/6*dt^(-1)*u^3*rho*dl^4+dt^(-1)*u*rho*omega_2^(-2)*dl^4*omega_3^(-1)+dt^(-1)*u*rho*omega_2^(-1)*dl^4*omega_3^(-2)+1/2*dt^(-1)*c_s^2*u*rho*dl^4*omega_3^(-2)+2*dt^(-1)*c_s^2*u*rho*omega_2^(-1)*dl^4*omega_3^(-1)-1/2*dt^(-1)*u*rho*omega_2^(-2)*dl^4-4/3*dt^(-1)*u^3*rho*omega_2^(-1)*dl^4+2*dt^(-1)*u^3*rho*omega_2^(-1)*dl^4*omega_3^(-1)
D[0,1,1,1](u)
	5/4*u^2*rho*dl^3*omega_3^(-1)+3*rho*omega_2^(-1)*dl^3*omega_3^(-1)+4*u^2*rho*omega_2^(-2)*dl^3*omega_3^(-1)+1/2*rho*dl^3*omega_3^(-2)+u^2*rho*omega_2^(-1)*dl^3*omega_3^(-2)+25/12*c_s^2*rho*omega_2^(-1)*dl^3+3/4*c_s^2*rho*dl^3*omega_3^(-1)-3*c_s^2*rho*omega_2^(-1)*dl^3*omega_3^(-1)+rho*omega_2^(-2)*dl^3+9/4*u^2*rho*omega_2^(-1)*dl^3+1/12*rho*dl^3-1/4*u^2*rho*dl^3+2*c_s^2*rho*omega_2^(-3)*dl^3-7/2*u^2*rho*omega_2^(-2)*dl^3-11/12*rho*omega_2^(-1)*dl^3+c_s^2*rho*omega_2^(-1)*dl^3*omega_3^(-2)-1/2*c_s^2*rho*dl^3*omega_3^(-2)+2*c_s^2*rho*omega_2^(-2)*dl^3*omega_3^(-1)-1/6*c_s^2*rho*dl^3-4*c_s^2*rho*omega_2^(-2)*dl^3+u^2*rho*omega_2^(-3)*dl^3-5*u^2*rho*omega_2^(-1)*dl^3*omega_3^(-1)-2*rho*omega_2^(-2)*dl^3*omega_3^(-1)-rho*omega_2^(-1)*dl^3*omega_3^(-2)-3/4*rho*dl^3*omega_3^(-1)-1/2*u^2*rho*dl^3*omega_3^(-2)
D[0,0,1,1](u)
	dt*u*rho*dl^2*omega_3^(-2)-1/2*dt*u*rho*omega_2^(-2)*dl^2-2*dt*u*rho*omega_2^(-1)*dl^2*omega_3^(-2)+dt*u*rho*omega_2^(-3)*dl^2-2*dt*u*rho*omega_2^(-2)*dl^2*omega_3^(-1)+4*dt*u*rho*omega_2^(-1)*dl^2*omega_3^(-1)-dt*u*rho*dl^2*omega_3^(-1)-1/2*dt*u*rho*omega_2^(-1)*dl^2
D[0,0,0,1](u)
	-1/2*dt^2*rho*omega_2^(-1)*dl+3/2*dt^2*rho*omega_2^(-2)*dl-dt^2*rho*omega_2^(-3)*dl
D[1](rho)  *  D[1](u)
	1/2*dt^(-1)*u*dl^2-dt^(-1)*u*omega_2^(-1)*dl^2
D[1](u)  *  D[1](u)
	1/2*dt^(-1)*rho*dl^2-dt^(-1)*rho*omega_2^(-1)*dl^2
D[1](rho)  *  D[0](u)
	-omega_2^(-1)*dl+1/2*dl

+43 −0

File added.

Preview size limit exceeded, changes collapsed.

+21 −0

File added.

Preview size limit exceeded, changes collapsed.

+31 −0

File added.

Preview size limit exceeded, changes collapsed.

Loading