Commit be73053c authored by thomaskummer's avatar thomaskummer
Browse files

intrepid lifev update

parent e62851cd
......@@ -519,7 +519,7 @@ int main (int argc, char** argv)
};
solver.fullMeshPtr()->showMe();
if ( 0 == comm->MyPID() ) solver.fullMeshPtr()->showMe();
Vector3D center1, center2;
Real radius1 = 2;
......@@ -726,7 +726,7 @@ int main (int argc, char** argv)
//std::cout << "\nJFe = " << JFe;
//std::cout << "\nJCirc = " << JCirc;
//std::cout << "\nJR = " << JR;
std::cout << "\n==============================================================="; }
std::cout << "\n===============================================================\n"; }
};
......
......@@ -612,8 +612,9 @@ EMSolver<Mesh, ElectroSolver>::setupMechanicalSolver ( GetPot& dataFile)
{
if (M_commPtr -> MyPID() == 0)
{
std::cout << "EMS - setting up mechanical solver\n";
std::cout << "\nEMSolver: setupMechanicalSolver ... ";
}
boost::shared_ptr<StructuralConstitutiveLawData> dataStructure (new StructuralConstitutiveLawData( ) );
dataStructure->setup (dataFile);
......@@ -642,6 +643,12 @@ EMSolver<Mesh, ElectroSolver>::setupMechanicalSolver ( GetPot& dataFile)
M_wteTotal.setup(dataStructure, dFESpace, dETFESpace, M_commPtr, 0, M_EMStructuralOperatorPtr->EMMaterial());
// M_wtePassive.setup(dataStructure, dFESpace, dETFESpace, M_commPtr, 0, "passive");
// M_wteActive.setup(dataStructure, dFESpace, dETFESpace, M_commPtr, 0, "active");
if (M_commPtr -> MyPID() == 0)
{
std::cout << "done";
}
}
/////////////////////
......
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