adams_bashforth_strategy.h

Go to the documentation of this file.

//
// Author:
// Guillermo Casas gcasas@cimne.upc.edu
//
#if !defined(KRATOS_ADAMS_BASHFORTH_STRATEGY)
#define KRATOS_ADAMS_BASHFORTH_STRATEGY
 
// /* External includes */
 
// System includes
 
// Project includes
#include "utilities/timer.h"
#include "custom_elements/Particle_Contact_Element.h"
#include "includes/variables.h"
#include "includes/deprecated_variables.h"
 
/* System includes */
#include <limits>
#include <iostream>
#include <iomanip>
#include <time.h>
 
/* External includes */
#ifdef _OPENMP
#include <omp.h>
#endif
 
#define CUSTOMTIMER 0  // ACTIVATES AND DISABLES ::TIMER:::::
 
/* Project includes */
#include "includes/define.h"
#include "utilities/openmp_utils.h"
#include "includes/model_part.h"
#include "custom_strategies/strategies/explicit_solver_strategy.h"
 
#ifdef USING_CGAL
#include <CGAL/spatial_sort.h>
#endif
 
/* Timer defines */
#ifdef CUSTOMTIMER
#define KRATOS_TIMER_START(t) Timer::Start(t);
#define KRATOS_TIMER_STOP(t) Timer::Stop(t);
#else
#define KRATOS_TIMER_START(t)
#define KRATOS_TIMER_STOP(t)
#endif
 
namespace Kratos {
 
    class KRATOS_API(SWIMMING_DEM_APPLICATION) AdamsBashforthStrategy: public ExplicitSolverStrategy {
    public:
 
        typedef ExplicitSolverStrategy BaseType;
        typedef BaseType::NodesArrayType NodesArrayType;
        typedef BaseType::ElementsArrayType ElementsArrayType;
        typedef BaseType::ElementsIterator ElementsIterator;
        typedef BaseType::ConditionsArrayType ConditionsArrayType;
        typedef GlobalPointersVector<Element> ParticleWeakVectorType;
        typedef GlobalPointersVector<Element >::iterator ParticleWeakIteratorType;
        typedef ParticleWeakVectorType::ptr_iterator ParticleWeakIteratorType_ptr;
 
        using BaseType::mpInlet_model_part;
        using BaseType::mpCluster_model_part;
        using BaseType::mpContact_model_part;
        using BaseType::GetModelPart;
        using BaseType::GetFemModelPart;
        using BaseType::mNumberOfThreads;
        using BaseType::mListOfSphericParticles;
        using BaseType::mListOfGhostSphericParticles;
        using BaseType::SearchNeighbours;
        using BaseType::SetSearchRadiiOnAllParticles;
        KRATOS_CLASS_POINTER_DEFINITION(AdamsBashforthStrategy);
 
        AdamsBashforthStrategy(){}
 
        AdamsBashforthStrategy(ExplicitSolverSettings& settings,
                const double max_delta_time,
                const int n_step_search,
                const double safety_factor,
                const int delta_option,
                ParticleCreatorDestructor::Pointer p_creator_destructor,
                DEM_FEM_Search::Pointer p_dem_fem_search,
                SpatialSearch::Pointer pSpSearch,
                Parameters strategy_parameters):
                ExplicitSolverStrategy(settings,
                                       max_delta_time,
                                       n_step_search,
                                       safety_factor,
                                       delta_option,
                                       p_creator_destructor,
                                       p_dem_fem_search,
                                       pSpSearch,
                                       strategy_parameters)
        {
            mFirstStep = true;
            ExplicitSolverStrategy::GetParticleCreatorDestructor() = p_creator_destructor;
        }
        virtual ~AdamsBashforthStrategy() {}
 
        double SolveSolutionStep() override;
 
    protected:
 
        void ReconstructForces(ModelPart& r_model_part);
 
    private:
        bool mFirstStep;
    }; // Class AdamsBashforthStrategy
} // namespace Kratos.
 
#endif // KRATOS_ADAMS_BASHFORTH_STRATEGY  defined