FreeFallRestartUnitTest.cpp File Reference

#include "DPMBase.h"
#include "Particles/BaseParticle.h"
#include "Walls/InfiniteWall.h"
#include <iostream>
#include <Species/LinearViscoelasticSpecies.h>
#include <Logger.h>

Go to the source code of this file.

Classes
class	FreeFall
	This code tests: 1) Restarting 2) Saving arcoss multiple files 3) and accepting command line argument. More...

Functions
void	runFreeFall (int argc, char *argv[])
int	main (int argc, char *argv[])

Function Documentation

int main	(	int argc	,
		char *	argv[]
	)

Definition at line 87 of file FreeFallRestartUnitTest.cpp.

References BaseHandler< T >::begin(), BaseHandler< T >::end(), FATAL, INFO, logger, DPMBase::particleHandler, DPMBase::readRestartFile(), and runFreeFall().

{
    
    //On first pass i.e. code being called as a selftest, it enters here and call the code again 3 times with different arguments.
    if (argc==1)
        {
            logger(INFO,"Case 1: not restarted");
            if (system("./FreeFallRestartUnitTest -tmax 0.4 -name FreeFallRestartUnitTest_NoRestart"))
                       logger(FATAL, "code did not run");

            logger(INFO,"Case 2: restarted at t=0.2");
            //restarted at t=0.2

            if (system("./FreeFallRestartUnitTest -tmax 0.2 -name FreeFallRestartUnitTest_Restarted"))
                       logger(FATAL, "code did not run");
        
            if (system("./FreeFallRestartUnitTest -r FreeFallRestartUnitTest_Restarted -tmax 0.4"))
                       logger(FATAL,"code did not run");

        
            logger(INFO,"Case 3: restarted at t=0.2; using separate data files");


            if (system("./FreeFallRestartUnitTest -tmax 0.2 -name FreeFallRestartUnitTest_SeparateFiles -fileTypeData 2"))
                       logger(FATAL,"code did not run");
            if (system("./FreeFallRestartUnitTest -r FreeFallRestartUnitTest_SeparateFiles -tmax 0.4 "))
                       logger(FATAL,"code did not run");
    } else
    {
        runFreeFall(argc, argv);
        return 0;
    }
    
    //To compare the code using gnuplot follow the instructions below.
    //gnuplot> plot 'free_fall_restart/free_fall_restart.ene' u 1:3 w l, 'free_fall_restart/free_fall_no_restart.ene' u 1:3 w l
    //../sc/fpdiff.py ./free_fall_restart.fstat ./free_fall_no_restart.fstat
      
    //Final stage now we check what we get.
    logger(INFO,"Finished running, now comparing");
    
    FreeFall FreeFallProblemReload;
    FreeFall FreeFallProblemReloadRestart;
    FreeFall FreeFallProblemReloadSplitFiles;
    
    FreeFallProblemReload.readRestartFile("FreeFallRestartUnitTest_NoRestart.restart");
    FreeFallProblemReloadRestart.readRestartFile("FreeFallRestartUnitTest_Restarted.restart");
    FreeFallProblemReloadSplitFiles.readRestartFile("FreeFallRestartUnitTest_SeparateFiles.restart");

    auto FreeFallProblemReloadRestartIt = FreeFallProblemReloadRestart.particleHandler.begin();
    auto FreeFallProblemReloadSplitFilesIt = FreeFallProblemReloadSplitFiles.particleHandler.begin();
   
    for (auto FreeFallProblemReloadIt = FreeFallProblemReload.particleHandler.begin(); FreeFallProblemReloadIt != FreeFallProblemReload.particleHandler.end(); ++FreeFallProblemReloadIt)
    {
      //  if (!(*FreeFallProblemReloadIt))->getTimeMax()
        
        if (!(*FreeFallProblemReloadIt)->getPosition().isEqualTo((*FreeFallProblemReloadRestartIt)->getPosition(),1e-6))
        {
            logger(FATAL,"Particles is not in the same place after restart. Before it was % and now it is %.",(*FreeFallProblemReloadIt)->getPosition(),(*FreeFallProblemReloadRestartIt)->getPosition());
        }
        if (!(*FreeFallProblemReloadIt)->getPosition().isEqualTo((*FreeFallProblemReloadSplitFilesIt)->getPosition(), 1e-10))
        {
             logger(FATAL,"Particles velocities are not the same place. Before it was % and now it is %.",(*FreeFallProblemReloadIt)->getVelocity(),(*FreeFallProblemReloadRestartIt)->getVelocity());
        }
        ++FreeFallProblemReloadRestartIt;
        ++FreeFallProblemReloadSplitFilesIt;
    }
    
    return 0;
}

void runFreeFall	(	int	argc,
		char *	argv[]
	)

Start off my solving the default problem

Definition at line 72 of file FreeFallRestartUnitTest.cpp.

References BaseHandler< T >::copyAndAddObject(), Files::dataFile, File::getFstream(), Files::setName(), Files::setSaveCount(), DPMBase::setTimeStep(), DPMBase::setXMax(), DPMBase::setYMax(), DPMBase::solve(), and DPMBase::speciesHandler.

Referenced by main().

{
    FreeFall FreeFallProblem;
    auto species = FreeFallProblem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    species->setDensity(2000);
    FreeFallProblem.setName("FreeFallRestartUnitTest");
    FreeFallProblem.setSaveCount(10000);
    FreeFallProblem.setTimeStep(1e-6);
    FreeFallProblem.setYMax(0.1);
    FreeFallProblem.setXMax(0.01);
    FreeFallProblem.dataFile.getFstream().precision(10);
    FreeFallProblem.solve(argc,argv);
}