AxisymmEuler.cpp
View Options

			//
			// Created by Fabian Moenkeberg on 2020-01-29.
			//

			#include "AxisymmEuler.hpp"

			std::vector<double> AxisymmEuler::source(std::vector<std::vector<double>> &U, MeshBase *mesh, double t, int iCell){
			std::vector<double> ret(U.size(),0);

			long double rrho, rp, rm1, rm2, rm3, rE, rm, uD;

			rrho = U[iCell][0];
			rm1 = U[iCell][1];
			rm2 = U[iCell][2];
			rm3 = U[iCell][4];
			rE = U[iCell][3];
			rp = (gamma - 1)(rE - 0.5(rm1rm1 + rm2rm2)/rrho);

			std::vector<double> r = mesh->MidPoint(iCell);
			ret[2] = rp/r[2];
			return ret;
			}

ConSol/ConSol/ConSol/Source/AxisymmEuler.cpp70d9a38942f6master