NoReconstr_2D.cpp
a701db2d9ca1master

NoReconstr_2D.cpp
View Options

			//
			// Created by Fabian Moenkeberg on 30.07.18.
			//

			#include "NoReconstr_2D.hpp"

			void NoReconstr_2D::reconstruct(std::vector<std::vector<std::vector<std::vector<double>>>> &UReconstr, std::vector<std::vector<double>> &U,
			MeshBase *mesh, double t, std::vector<std::vector<double>> &Source){


			int dimSyst = mesh->getNelem();
			int nCells = mesh->getNcells();
			UReconstr.resize(2,std::vector<std::vector<std::vector<double>>>(1,std::vector<std::vector<double>>(mesh->getNedges(),std::vector<double>(dimSyst))));
			std::vector<double> source_loc;

			std::vector<int> c2e_loc(4);
			int nNodes;
			for (int ii = 0; ii < nCells; ii++){
			c2e_loc = mesh->getC2E_internal(ii);
			nNodes = c2e_loc.size();
			int internal_i = mesh->getInternal(ii);

			for (int j = 0; j < nNodes; j++){

			if (mesh->getRight(c2e_loc[j])== internal_i){
			UReconstr[0][c2e_loc[j]][0] = U[internal_i];
			}else{
			UReconstr[1][c2e_loc[j]][0] = U[internal_i];
			}
			}
			if (source->getIsUsed()) {
			std::vector<std::vector<double>> xy = mesh->getXCoordOfCell(internal_i);
			source_loc = source->source(U[internal_i], t, std::vector<double>{(xy[0][0] + xy[1][0] + xy[2][0]) / 3,
			(xy[0][1] + xy[1][1] + xy[2][1]) / 3});

			for (int k = 0; k < dimSyst; k++) {
			Source[internal_i][k] = source_loc[k];
			}
			}
			}
			mesh->updateEdges(UReconstr, t);
			}