spring_contact.cc File Reference

#include "generic.h"

Namespaces
	GlobalFct

Functions
void	GlobalFct::reset (Vector< double > &unknowns)
	Reset unknowns to some nontrivial garbage. More...
void	GlobalFct::get_residuals (const Vector< double > &param, const Vector< double > &unknowns, Vector< double > &residuals)
	Global residual fct. More...
void	GlobalFct::set_spurious_solution (Vector< double > &unknowns)
	Set unknowns to spurious solution. More...
void	GlobalFct::plot_it (const std::string filename)
	Plot "landscape" of residuals (only for 2D problems!) More...
int	main (int argc, char **argv)

Variables
double	GlobalFct::Force =0.0
	Force. More...
double	GlobalFct::Stiffness =1.0
	Spring stiffness. More...
double	GlobalFct::U_max =5.5
	Max. displacement. More...
unsigned	GlobalFct::N_primary =1

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

/////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 CommandLineArgs::specify_command_line_flag("--bilinear");
  
 CommandLineArgs::specify_command_line_flag("--completely_smooth");
 
 CommandLineArgs::specify_command_line_flag("--single_kink");
 
 CommandLineArgs::specify_command_line_flag("--old_version");
 
 CommandLineArgs::specify_command_line_flag("--kuhn_tucker");
 
 CommandLineArgs::specify_command_line_flag("--n_primary",
                                            &GlobalFct::N_primary);
 
 // Doc available command line flags
 CommandLineArgs::doc_available_flags();
 
 // Parse command line
 CommandLineArgs::parse_and_assign(); 
 
 // Doc what has actually been specified on the command line
 CommandLineArgs::doc_specified_flags();
 
 // Number of steps for load incrementation
 unsigned nstep=10;
 unsigned nstep_total=2*nstep;
 
 BlackBoxFDNewtonSolver::Max_iter=30;
 
 ofstream outfile2;
 outfile2.open("trace.dat");
 
 BlackBoxFDNewtonSolver::Doc_Progress=true;
 Vector<double> param;
 Vector<double> unknowns;
 GlobalFct::reset(unknowns);
 
 GlobalFct::Force=0.0;
 for (unsigned j=0;j<nstep_total;j++)
  {
   
   oomph_info << "\n\nSolving for force: " << GlobalFct::Force
              << std::endl << std::endl;
   
   char name[100];   
   sprintf(name,"landscape%i.dat",j);
   std::string filename(name);
   GlobalFct::plot_it(filename);
   try
    {
     BlackBoxFDNewtonSolver::black_box_fd_newton_solve(
      GlobalFct::get_residuals,param,unknowns);
     
     oomph_info << "Number of Newton iterations to convergence: "
                << BlackBoxFDNewtonSolver::N_iter_taken++ << std::endl;
     
     char name[100];   
     sprintf(name,"soln%i.dat",j);
     ofstream outfile;
     outfile.open(name);
     outfile <<  unknowns[0] << " " 
             <<  unknowns[1] << "\n";
     outfile.close();
     
     unsigned n=unknowns.size();
     outfile2 << GlobalFct::Force << " ";
     for (unsigned k=0;k<n;k++)
      {
       outfile2 << unknowns[k] << " ";
      }
     outfile2 << std::endl;
    }
   catch(OomphLibError) {}
   
   // Sweep up and then down
   if (j<nstep)
    {
     GlobalFct::Force+=1.0;
    }
   else
    {
     GlobalFct::Force-=1.0;
    }
  }
 
 outfile2.close();
 
 
 //----------------------------------------
 
 oomph_info << "\n\n\n\nTRY POSSIBLY SPURIOUS SOLUTION\n"
            <<         "==============================\n\n";
 
 // Try spurious solution
 GlobalFct::set_spurious_solution(unknowns);
 
 BlackBoxFDNewtonSolver::black_box_fd_newton_solve(
  GlobalFct::get_residuals,param,unknowns);
 
 oomph_info << "Number of Newton iterations to convergence: "
            << BlackBoxFDNewtonSolver::N_iter_taken++ << std::endl;
 
 oomph_info << "\nPossibly spurious solution: "
            << unknowns[0] << " " <<  unknowns[1] << "\n\n";
 
 
 
} // end_of_main

References oomph::BlackBoxFDNewtonSolver::black_box_fd_newton_solve(), oomph::CommandLineArgs::doc_available_flags(), oomph::BlackBoxFDNewtonSolver::Doc_Progress, oomph::CommandLineArgs::doc_specified_flags(), MergeRestartFiles::filename, GlobalFct::Force, GlobalFct::get_residuals(), j, k, oomph::BlackBoxFDNewtonSolver::Max_iter, n, oomph::BlackBoxFDNewtonSolver::N_iter_taken, GlobalFct::N_primary, plotDoE::name, oomph::oomph_info, oomph::CommandLineArgs::parse_and_assign(), GlobalFct::plot_it(), GlobalFct::reset(), GlobalFct::set_spurious_solution(), Flag_definition::setup(), oomph::CommandLineArgs::specify_command_line_flag(), and oomph::Global_string_for_annotation::string().