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

Added raw txt output that is included in the resulting supplementary pdfs

	modified:   d1q3_ade/supp_d1q3_ade.tex
	modified:   d1q3_nse/supp_d1q3_nse.tex
	modified:   d2q5_ade/supp_d2q5_ade.tex
	modified:   d2q9_nse/supp_d2q9_nse.tex
	modified:   d3q27_nse/supp_d3q27_nse.tex
	modified:   d3q7_ade/supp_d3q7_ade.tex
	modified:   latex/header.tex
	modified:   scripts/build.sh
	modified:   src/lbm.h
	modified:   src/tpde.hpp
	modified:   supp_d1q3_ade.pdf
	modified:   supp_d1q3_nse.pdf
	modified:   supp_d2q5_ade.pdf
	modified:   supp_d2q9_nse.pdf
	modified:   supp_d3q27_nse.pdf
	modified:   supp_d3q7_ade.pdf
parent 8073b588
......@@ -22,6 +22,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{MRT}
......@@ -30,6 +32,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{CLBM}
......@@ -38,6 +42,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\section{Comparison of SRT, MRT, and CLBM}
......
......@@ -22,9 +22,13 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsection{MRT}
......@@ -33,9 +37,13 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
......@@ -45,9 +53,13 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\section{Comparison of SRT, MRT, and CLBM}
......
......@@ -22,6 +22,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{MRT1}
......@@ -30,6 +32,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{MRT2}
......@@ -38,6 +42,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{CLBM1}
......@@ -46,6 +52,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{CLBM2}
......@@ -54,6 +62,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\section{Comparison of SRT, MRT, and CLBM}
......
......@@ -22,10 +22,18 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsection{MRT}
......@@ -34,10 +42,18 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsection{MRT2}
......@@ -46,10 +62,19 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsection{CLBM1}
......@@ -58,10 +83,18 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsection{CLBM2}
......@@ -70,10 +103,18 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsection{CuLBM1}
......@@ -82,10 +123,18 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsection{CuLBM2}
......@@ -94,10 +143,18 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\section{Comparison of SRT, MRT, CLBM, and CuLBM}
......
......@@ -22,12 +22,23 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsubsection{Conservation of momentum: $\rho \ow$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_03.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_03}
\subsection{MRT}
......@@ -36,12 +47,23 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsubsection{Conservation of momentum: $\rho \ow$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_03.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_03}
\subsection{MRT2}
......@@ -50,12 +72,23 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsubsection{Conservation of momentum: $\rho \ow$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_03.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_03}
\subsection{CLBM1}
......@@ -64,12 +97,23 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsubsection{Conservation of momentum: $\rho \ow$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_03.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_03}
\subsection{CLBM2}
......@@ -78,12 +122,23 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsubsection{Conservation of momentum: $\rho \ow$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_03.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_03}
\subsection{CuLBM1}
......@@ -92,12 +147,23 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsubsection{Conservation of momentum: $\rho \ow$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_03.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_03}
\subsection{CuLBM2}
......@@ -106,12 +172,23 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass: $\rho$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsubsection{Conservation of momentum: $\rho \ou$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_01.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_01}
\subsubsection{Conservation of momentum: $\rho \ov$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_02.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_02}
\subsubsection{Conservation of momentum: $\rho \ow$}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_03.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_03}
\section{Comparison of SRT, MRT, CLBM, and CuLBM}
......
......@@ -22,6 +22,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{MRT1}
......@@ -30,6 +32,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{MRT2}
......@@ -38,6 +42,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{CLBM1}
......@@ -46,6 +52,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\subsection{CLBM2}
......@@ -54,6 +62,8 @@
\input{\path/supp_\model_\colmod_defs}
\subsubsection{Conservation of mass equation}
\attachtxt{output_\model_\pde_\colmod_symbolic_pde_00.txt}
\input{output_\model_\pde_\colmod_symbolic/spatial_EPDE/pde_00}
\section{Comparison of SRT, MRT, and CLBM}
......
......@@ -14,8 +14,17 @@
\usepackage{subcaption}
\usepackage{algorithm}
\usepackage[left=2cm,right=2cm,a4paper]{geometry}
%\let\pdfxform\saveboxresource % fix for lualatex and attachfile
%\usepackage{embedfile}
\usepackage{attachfile}
\newcommand{\attachtxt}[1]{\attachfile[icon=Paperclip, mimetype=text/plain, description={#1}, author={LBMAT}]{#1} attached text file: \texttt{#1}}
\sloppy
\usepackage{breqn} % break equations automatically
%\breqnsetup{breakdepth={1},compact}
\usepackage{framed}
......
......@@ -44,10 +44,38 @@ for a in output*symbolic; do
if [ ! -f "${a}/spatial_EPDE/pde_00.tex" ]; then
../scripts/script_backsubst_searcher.py ${a}
"$CODEDIR/export_latex" ${a}
# check whether latex files exist
if [ ! -f "${a}/spatial_EPDE/pde_00.tex" ]; then
echo "${a}/spatial_EPDE/pde_00.tex not found"
exit;
fi
# check whether latex files exist
if [ ! -f "${a}/spatial_EPDE/pde_00.tex" ]; then
echo "${a}/spatial_EPDE/pde_00.tex not found"
exit;
fi
# for compatibility reasons
if [ ! -f "${a}/spatial_EPDE/pde_00.txt" ]; then
../scripts/script_backsubst_searcher.py ${a}
"$CODEDIR/export_latex" ${a}
fi
if [ ! -f "${a}/spatial_EPDE/pde_00.txt" ]; then
echo "${a}/spatial_EPDE/pde_00.txt not found"
exit;
fi
# check symbolic links
if [ ! -f "${a}_pde_00.txt" ]; then
ln "${a}/spatial_EPDE/pde_00.txt" "${a}_pde_00.txt" -s
fi
if [ ! -f "${a}_pde_01.txt" ]; then
if [ -f "${a}/spatial_EPDE/pde_01.txt" ]; then
ln "${a}/spatial_EPDE/pde_01.txt" "${a}_pde_01.txt" -s
fi
fi
if [ ! -f "${a}_pde_02.txt" ]; then
if [ -f "${a}/spatial_EPDE/pde_02.txt" ]; then
ln "${a}/spatial_EPDE/pde_02.txt" "${a}_pde_02.txt" -s
fi
fi
if [ ! -f "${a}_pde_03.txt" ]; then
if [ -f "${a}/spatial_EPDE/pde_03.txt" ]; then
ln "${a}/spatial_EPDE/pde_03.txt" "${a}_pde_03.txt" -s
fi
fi
done
......@@ -57,7 +85,8 @@ done
# call latex
buf_size=100000000
TEXCOMP="lualatex"
#TEXCOMP="lualatex"
TEXCOMP="pdflatex"
TEXFILE=supp_${MODEL}
${TEXCOMP} ${TEXFILE}
......
......@@ -18,7 +18,7 @@ struct LBM
matrix M, Mi, E[MQ], I;
symbol sym_o=symbol("omega","\\omega");
symbol sym_o=symbol("o","\\omega");
symbol sym_os[MQ];
symbol sym_cs=symbol("c_s","c_s");
symbol x=symbol("x","\\x"),y=symbol("y","\\y"),z=symbol("z","\\z"),t=symbol("t");
......@@ -141,7 +141,7 @@ struct LBM
char desc[100], desc_tex[100];
for (int q=0;q<MQ;q++)
{
sprintf(desc,"omega_%d",q+1);
sprintf(desc,"o_%d",q+1);
sprintf(desc_tex,"\\omega_{%d}",q+1);
sym_os[q]=symbol(desc,desc_tex);
}
......
......@@ -128,7 +128,9 @@ void TPDE::print(ostringstream &out)
str1 = ReplaceString( str1, "fu", "u");
str1 = ReplaceString( str1, "fv", "v");
str1 = ReplaceString( str1, "fw", "w");
str1 = ReplaceString( str1, "frho", "rho");
str1 = ReplaceString( str1, "frho", "frh");
str1 = ReplaceString( str1, "o", "omega");
str1 = ReplaceString( str1, "frh", "rho");
out << str1 << endl;
}
......@@ -408,22 +410,31 @@ void TPDE::save_raw(const char* dirname, const char* subdirname, const char*fmt,
sprintf(fn,"%s/%s/%s.txt",dirname,subdirname,name);
ostringstream sout_all;
sout_all << "Complete PDE that can be directly imported to Maxima:" << endl;
sout_all << raw() << endl << endl;
sout_all << "List of all coefficients of each partial derivatives:" << endl;
sout_all << sout.str() << endl;
ostringstream sout2;
sout2 << raw() << "=0" << endl;
// sout_all << "Complete PDE that can be directly imported to Maxima (for instance):" << endl;
// sout_all << raw() << endl << endl;
// sout_all << "List of all coefficients of each partial derivatives:" << endl;
// sout_all << sout.str() << endl;
string str1 = sout_all.str();
string str1 = sout2.str();
// str1 = ReplaceString( str1, "(t,x,y,z)", "");
// str1 = ReplaceString( str1, "(t,\\x,\\y,\\z)", "");
str1 = ReplaceString( str1, "frho", "rho");
str1 = ReplaceString( str1, "rho", "rh");
str1 = ReplaceString( str1, "o", "omega");
str1 = ReplaceString( str1, "fu", "u");
str1 = ReplaceString( str1, "fv", "v");
str1 = ReplaceString( str1, "fw", "w");
str1 = ReplaceString( str1, "frho", "rho");
str1 = ReplaceString( str1, "rh", "rho");
ofstream fout(fn);
fout << str1;
// fout << str1;
fout << "Complete expression of PDE that can be directly imported to Maxima:" << endl;
fout << str1 << endl << endl;
fout << "List of all coefficients of each partial derivatives:" << endl;
fout << sout.str() << endl;
fout.close();
}
......
No preview for this file type
No preview for this file type
No preview for this file type
No preview for this file type
No preview for this file type
No preview for this file type
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment