5 #include "Epetra_ConfigDefs.h" 8 #include "Epetra_MpiComm.h" 10 #include "Epetra_SerialComm.h" 13 #include "Teuchos_CommandLineProcessor.hpp" 14 #include "Teuchos_StandardCatchMacros.hpp" 15 #include "Teuchos_ParameterList.hpp" 16 #include "Teuchos_XMLParameterListCoreHelpers.hpp" 20 int main(
int argc,
char *argv[]){
22 std::string xmlInFileName =
"";
23 std::string extraXmlFile =
"";
24 std::string xmlOutFileName =
"paramList.out";
26 Teuchos::CommandLineProcessor clp(
false);
27 clp.setOption(
"xml-in-file",&xmlInFileName,
"The XML file to read into a parameter list");
28 clp.setDocString(
"TO DO.");
30 Teuchos::CommandLineProcessor::EParseCommandLineReturn
31 parse_return = clp.parse(argc,argv);
32 if( parse_return != Teuchos::CommandLineProcessor::PARSE_SUCCESSFUL ) {
33 std::cout <<
"\nEnd Result: TEST FAILED" << std::endl;
38 MPI_Init(&argc, &argv);
39 Epetra_MpiComm Comm(MPI_COMM_WORLD);
41 Epetra_SerialComm Comm;
44 Teuchos::RCP<Teuchos::ParameterList> paramList = Teuchos::rcp(
new Teuchos::ParameterList);
45 if(xmlInFileName.length()) {
46 Teuchos::updateParametersFromXmlFile(xmlInFileName, inoutArg(*paramList));
50 paramList->print(std::cout,2,
true,
true);
54 Teuchos::RCP<newtonRaphson> Newton = Teuchos::rcp(
new newtonRaphson(*my_interface,*paramList));
56 Newton->Initialization();
57 Newton->setParameters(*paramList);
58 int error = Newton->Solve_with_Aztec(
true);
60 std::string path =
"/home/s/staber/Trilinos_results/arteries/gmrf_neumann/";
61 std::string filename1 = path +
"disp_mean_model.mtx";
62 Newton->print_newton_solution(filename1);
63 std::string filename2 = path +
"stress_mean_model";
64 my_interface->compute_center_cauchy_stress(*Newton->x,filename2);
int main(int argc, char *argv[])