Commit 5c917579 authored by Thomas Kummer's avatar Thomas Kummer
Browse files

Problem with bcBase idList solved

parent 53f4e355
......@@ -443,6 +443,8 @@ BCBase::BCBase ( const BCBase& bcBase ) :
{
M_bcVector = bcBase.M_bcVector->clone();
}
std::cout << 1 << std::endl;
std::cout << M_name << " " << M_idVector.size() << " " << bcBase.M_idVector.size() << std::endl;
// Important!!: The set member M_idSet is always empty at this point, it is just
// an auxiliary container used at the moment of the boundary update (see BCHandler::bcUpdate)
......
......@@ -159,8 +159,11 @@ BCHandler::addBC ( const bcName_Type& name,
const bcMode_Type& mode,
BCVectorBase& bcVector )
{
M_bcList.push_back ( BCBase ( name, flag, type, mode, bcVector ) );
std::sort ( M_bcList.begin(), M_bcList.end() );
std::cout << "bcHandler: " << name << std::endl;
M_bcList.push_back ( BCBase ( name, flag, type, mode, bcVector ) );
// std::sort ( M_bcList.begin(), M_bcList.end() );
std::cout << "bcHandler: " << name << std::endl;
}
void
......
......@@ -12,4 +12,4 @@ TRIBITS_ADD_EXECUTABLE(
TRIBITS_COPY_FILES_TO_BINARY_DIR(data_example_EMCoupling
SOURCE_FILES ParamList.xml dataP1BiVent dataP2BiVent dataP1SiVent inputfile singleVentricle.mesh FiberDirection.h5 SheetsDirection.h5
SOURCE_DIR ${CMAKE_CURRENT_SOURCE_DIR}
)
\ No newline at end of file
)
......@@ -365,10 +365,11 @@ int main (int argc, char** argv)
// Todo: Normal boundary condition!!
//solver.bcInterfacePtr() -> handler() -> addBC("LvPressure", LVFlag, Natural, Full, *pLvBCVectorPtr, 3); // BC for using function which keeps bc normal
solver.bcInterfacePtr() -> handler() -> addBC("LvPressure", LVFlag, Natural, Normal, *pLvBCVectorPtr, 3); // BC which is only at time zero normal
solver.bcInterfacePtr() -> handler() -> addBC("LvPressure", LVFlag, Natural, Normal, *pLvBCVectorPtr); // BC which is only at time zero normal
solver.bcInterfacePtr() -> handler() -> bcUpdate( *solver.structuralOperatorPtr() -> dispFESpacePtr() -> mesh(), solver.structuralOperatorPtr() -> dispFESpacePtr() -> feBd(), solver.structuralOperatorPtr() -> dispFESpacePtr() -> dof() );
solver.bcInterfacePtr() -> handler() -> showMe();
//********************************************//
// Volume integrators
//********************************************//
......
......@@ -579,7 +579,7 @@ EMSolver<Mesh, ElectroSolver>::setupMechanicalBC (std::string data_file_name,
M_bcInterfacePtr.reset (new bcInterface_Type() );
M_bcInterfacePtr->createHandler();
M_bcInterfacePtr->fillHandler ( data_file_name, "solid" );
M_bcInterfacePtr->handler()->bcUpdate ( *dFESpace->mesh(), dFESpace->feBd(), dFESpace->dof() );
// M_bcInterfacePtr->handler()->bcUpdate ( *dFESpace->mesh(), dFESpace->feBd(), dFESpace->dof() );
}
......
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