Set algloop to high values instead when otherwise errors were thrown.

Author: qichenghua

SimulationRuntime/cpp/Solver/Nox/Nox.cpp

@@ -115,14 +115,16 @@ void Nox::solve()
 	//_solverParametersPtr->sublist("Direction").set("Method", "Steepest Descent");
 
 	//resetting Line search method to default (Full Step, ie. Standard Newton with lambda=1)
-	_solverParametersPtr->sublist("Line Search").set("Method", "Backtrack");
+	_solverParametersPtr->sublist("Line Search").set("Method", "Full Step");
+	//_solverParametersPtr->sublist("Line Search").sublist("Full Step").set("Full Step", 0.5);
+	//_solverParametersPtr->sublist("Line Search").set("Method", "Polynomial");
 	// Set the level of output
     if (_generateoutput){
 		_solverParametersPtr->sublist("Printing").set("Output Information", NOX::Utils::Error + NOX::Utils::Warning + NOX::Utils::OuterIteration + NOX::Utils::Details + NOX::Utils::Debug); //(there are also more options, but error and outer iteration are the ones that I commonly use.
 	}else{
 		_solverParametersPtr->sublist("Printing").set("Output Information", NOX::Utils::Error);
 	}
-	_solverParametersPtr->sublist("Printing").set("Output Information", NOX::Utils::Error + NOX::Utils::Warning + NOX::Utils::OuterIteration + NOX::Utils::InnerIteration + NOX::Utils::Details + NOX::Utils::Debug); //(there are also more options, but error and outer iteration are the ones that I commonly use.
+	// _solverParametersPtr->sublist("Printing").set("Output Information", NOX::Utils::Error + NOX::Utils::Warning + NOX::Utils::OuterIteration + NOX::Utils::InnerIteration + NOX::Utils::Details + NOX::Utils::Debug); //(there are also more options, but error and outer iteration are the ones that I commonly use.
 
 
 	if (_generateoutput) std::cout << "creating noxLapackInterface" << std::endl;
@@ -330,7 +332,13 @@ void Nox::solve()
 		int numberofdifferentnormtests=5;//this variable is used twice.
 
 		_algLoop->setReal(_y);
-		_algLoop->evaluate();
+		try{
+			_algLoop->evaluate();
+		}
+		catch(const std::exception &ex)
+		{
+			if (_generateoutput) std::cout << "algloop evaluation after solve failed with error message:" << std::endl << ex.what() << std::endl << "Trying to continue without. This should hopefully lead to statusTest::Failed." << std::endl;
+		}
 
 		if (_generateoutput) {
 			std::cout << "solutionvector=(";
@@ -684,7 +692,12 @@ void Nox::LocaHomotopySolve(int numberofhomotopytries)
 	}
 
 	_algLoop->setReal(_y);
-	_algLoop->evaluate();
+	try{
+		_algLoop->evaluate();
+	}catch(const std::exception &ex)
+	{
+		if (_generateoutput) std::cout << "algloop evaluation after solve failed with error message:" << std::endl << ex.what() << std::endl << "Trying to continue without. This should hopefully lead to statusTest::Failed." << std::endl;
+	}
 
 	if (_generateoutput) {
 		std::cout << "solutionvector=(";

SimulationRuntime/cpp/Solver/Nox/NoxLapackInterface.cpp

@@ -187,18 +187,38 @@ bool NoxLapackInterface::computeActualF(NOX::LAPACK::Vector& f, const NOX::LAPAC
 		}
 	}
 
-	_algLoop->setReal(xp);
-	_algLoop->evaluate();
-	_algLoop->getRHS(rhs);
-
 	if (_generateoutput) {
 		std::cout << "we are at position x=(";
 		for (int i=0;i<_dimSys;i++){
 			std::cout << std::setprecision (std::numeric_limits<double>::digits10 + 8) << x(i) << " ";
 		}
 		std::cout << ")" << std::endl;
 		std::cout << std::endl;
+	}
+
+	_algLoop->setReal(xp);
+	_algLoop->getRHS(rhs);
+	try{
+		_algLoop->evaluate();
+		_algLoop->getRHS(rhs);
+	}catch(const std::exception &ex)
+	{
+		if (_generateoutput) std::cout << "calculating right hand side failed with error message:" << std::endl << ex.what() << std::endl;
+		//the following should be done when some to be implemented flag like "continue if function evaluation fails" is activated.
+		if (_generateoutput) std::cout << "setting high values into right hand side:" << std::endl << "(";
+		for(int i=0;i<_dimSys;i++){
+			if (rhs[i]==0.0){
+				rhs[i]=1000000.0;
+			}else{
+				rhs[i]=1000000.0*rhs[i];
+			}
+			if (_generateoutput) std::cout << rhs[i] << " ";
+		}
+		if (_generateoutput) std::cout << ")" << std::endl;
+	}
 
+
+	if (_generateoutput) {
 		std::cout << "the right hand side is given by (";
 		for (int i=0;i<_dimSys;i++){
 			std::cout << rhs[i] << " ";