product_of_sines_field.h

Go to the documentation of this file.

#if !defined(KRATOS_ETHIR_FLOW_FIELD_H)
#define KRATOS_ETHIR_FLOW_FIELD_H
 
// /* External includes */
 
// System includes
 
// Project includes
#include "includes/variables.h"
#include "swimming_DEM_application.h"
 
/* System includes */
#include <limits>
#include <iostream>
#include <iomanip>
 
/* External includes */
#ifdef _OPENMP
#include <omp.h>
#endif
 
/* Project includes */
#include "includes/define.h"
#include "utilities/openmp_utils.h"
#include "real_functions.h"
#include "velocity_field.h"
 
 
namespace Kratos
{
class ProductOfSines : public VelocityField
{
public:
 
KRATOS_CLASS_POINTER_DEFINITION(ProductOfSines);
 
 
ProductOfSines():VelocityField()
{
    unsigned int number_of_threads = ParallelUtilities::GetNumThreads();
    ResizeVectorsForParallelism(number_of_threads);
    mOmega = Globals::Pi;
}
 
ProductOfSines(const double period)
{
    unsigned int number_of_threads = ParallelUtilities::GetNumThreads();
    ResizeVectorsForParallelism(number_of_threads);
 
    KRATOS_ERROR_IF(period <= 0.0) << "The period must be non-negative." << std::endl;
 
    mOmega = Globals::Pi / period;
}
 
 
 
virtual ~ProductOfSines(){}
 
 
//***************************************************************************************************************
//***************************************************************************************************************
void ResizeVectorsForParallelism(const int n_threads) override;
 
void UpdateCoordinates(const double time, const array_1d<double, 3>& coor, const int i_thread = 0) override;
 
void UpdateCoordinates(const double time, const DenseVector<double>& coor, const int i_thread = 0) override;
 
void LockCoordinates(const int i_thread = 0) override;
 
void UnlockCoordinates(const int i_thread = 0) override;
//***************************************************************************************************************
//***************************************************************************************************************
 
 
 
 
 
virtual std::string Info() const override
{
    return "";
}
 
 
virtual void PrintInfo(std::ostream& rOStream) const override
{
}
 
 
virtual void PrintData(std::ostream& rOStream) const override
{
}
 
 
 
 
protected:
 
 
 
 
 
 
// Values
 
double U0(const int i_thread = 0) override;
 
// First-order derivatives
 
double U0D0(const int i_thread = 0) override;
double U0D1(const int i_thread = 0) override;
double U0D2(const int i_thread = 0) override;
 
// Second-order derivatives
 
double U0D0D0(const int i_thread = 0) override;
double U0D0D1(const int i_thread = 0) override;
double U0D0D2(const int i_thread = 0) override;
double U0D1D1(const int i_thread = 0) override;
double U0D1D2(const int i_thread = 0) override;
double U0D2D2(const int i_thread = 0) override;
 
 
 
 
 
 
 
 
 
private:
 
 
double mOmega;
std::vector<int> mCoordinatesAreUpToDate;
std::vector<double> mSin0;
std::vector<double> mCos0;
std::vector<double> mSin1;
std::vector<double> mCos1;
std::vector<double> mSin2;
std::vector<double> mCos2;
 
 
 
 
 
 
 
 
 
ProductOfSines & operator=(ProductOfSines const& rOther);
 
 
 
}; // Class ProductOfSines
 
 
 
 
 
} // namespace Kratos.
 
#endif // KRATOS_ETHIR_FLOW_FIELD_H defined