Homec4science
Diffusion ConSol (master) Browse

ConSol/ConSol/ConSol/Mesh/BC/MixedBC.cpp
36dd85413a37master

Last Change
Raw File

MixedBC.cpp
View Options

//
// MixedBC.cpp
// ConSol
//
// Created by Fabian Moenkeberg on 12.07.17.
// Copyright © 2017 EPFL. All rights reserved.
//
#include "MixedBC.hpp"
#include <assert.h>
void MixedBC::updateBoundary(std::vector<std::vector<double> > &U, MeshBase *mesh){
switch (mesh->getNdims()) {
case 1:
for (int k = 0; k < sysDim; k++){
BC[k].updateBoundary(U, mesh);
}
break;
default:
assert(0 && "MixedBC: higher dimensions are not implemented, yet.");
break;
}
}
void MixedBcComponent::updateBoundary(std::vector<std::vector<double> > &U, MeshBase *mesh){
switch (mesh->getNdims()) {
case 1:{
std::vector<int> internal = mesh->getInternal();
int indMax = *std::max_element(internal.begin(),internal.end());
int indMin = *std::min_element(internal.begin(),internal.end());
std::vector<std::vector<int>> gcMap1(2,std::vector<int>(indMin));
std::vector<std::vector<int>> gcMap2(2,std::vector<int>(indMin));
for (int i = 0; i < indMin; i++){
gcMap1[0][i] = i;
gcMap1[1][i] = 2*indMin - i;
gcMap2[0][i] = indMax + 1 + i;
gcMap2[1][i] = indMax - 1 - i;
}
switch (type[0]) {
case NEUMANN:{
for (int i = 0; i < indMin; i++){
U[gcMap1[0][i]][indComp] = U[gcMap1[1][i]][indComp];
// U[gcMap1[0][i]][indComp] = U[indMin][indComp] + 2.0*(U[gcMap1[1][i]][indComp] - U[indMin][indComp]);
}
break;
}
case DIRICHLET:{
for (int i = 0; i < indMin; i++){
U[gcMap1[0][i]][indComp] = -U[gcMap1[1][i]][indComp] + 2*bdryValues[0];
}
break;
}
default:
assert(0 && "MixedBcComponent: Bdry condition not implemented.");
break;
}
switch (type[1]) {
case NEUMANN:{
for (int i = 0; i < indMin; i++){
U[gcMap2[0][i]][indComp] = U[gcMap2[1][i]][indComp];
// U[gcMap2[0][i]][indComp] = U[indMax][indComp] + 2.0*(U[gcMap2[1][i]][indComp]- U[indMax][indComp]);
}
break;
}
case DIRICHLET:{
for (int i = 0; i < indMin; i++){
U[gcMap2[0][i]][indComp] = -U[gcMap2[1][i]][indComp] + 2*bdryValues[1];
}
break;
}
default:
assert(0 && "MixedBcComponent: Bdry condition not implemented.");
break;
}
break;
}
default:
assert( 0 && "MixedBcComponent: higher dimensions not implemented, yet.");
break;
}
}
c4science · Help