MixedBC_2D.cpp
f9c1f927e1e5master

MixedBC_2D.cpp
View Options


			// Created by Fabian Moenkeberg on 09.08.18.
			//

			#include "MixedBC_2D.hpp"

			void MixedBC_2D::initialize(MeshBase *mesh, Configuration config) {}

			void MixedBC_2D::updateBoundary(std::vector<std::vector<double> > &U, MeshBase *mesh, double t){
			// std::cout << "testB0\n";
			// for (int k = 0; k < sysDim; k++){
			// std::cout << "testB"<< k<<"\n";
			// BC[k].updateBoundary(U, mesh);
			// }
			int nCellsTot = mesh->getNcellsTot();
			// std::cout << "testB1\n";
			for (int i = 0; i < nCellsTot; i++){
			// std::cout << "testBi\n";
			U[i] = U[mesh->getBCMap(i)];
			}
			}

			void MixedBC_2D::updateEdges(std::vector<std::vector<std::vector<std::vector<double>>>> &UReconstr, MeshBase *mesh, double t){
			std::vector<std::vector<int>> edgesAtBC = mesh->getEdgesAtBC();
			// std::cout <<"test\n";
			int nQuad = UReconstr[0][0].size();
			int nEdgesBC = edgesAtBC.size();
			for (int k = 0; k < nQuad; k++){
			for(int i = 0; i < nEdgesBC; i++) {
			if (edgesAtBC[i][2] == 0) { //means right side of edge is out
			UReconstr[0][edgesAtBC[i][0]][k] = UReconstr[1][edgesAtBC[i][0]][k];
			}else{
			UReconstr[1][edgesAtBC[i][0]][k] = UReconstr[0][edgesAtBC[i][0]][k];
			}
			}
			}
			}

			void MixedBcComponent_2D::updateBoundary(std::vector<std::vector<double> > &U, MeshBase *mesh, double t){
			// std::cout << "testB1a\n";
			int nCellsTot = mesh->getNcellsTot();
			// std::cout << "testB1\n";
			for (int i = 0; i < nCellsTot; i++){
			// std::cout << "testBi\n";
			U[i] = U[mesh->getBCMap(i)];
			}
			}