161 lines
3.0 KiB
C++
Executable File
161 lines
3.0 KiB
C++
Executable File
/*------------------------------- phasicFlow ---------------------------------
|
|
O C enter of
|
|
O O E ngineering and
|
|
O O M ultiscale modeling of
|
|
OOOOOOO F luid flow
|
|
------------------------------------------------------------------------------
|
|
Copyright (C): www.cemf.ir
|
|
email: hamid.r.norouzi AT gmail.com
|
|
------------------------------------------------------------------------------
|
|
Licence:
|
|
This file is part of phasicFlow code. It is a free software for simulating
|
|
granular and multiphase flows. You can redistribute it and/or modify it under
|
|
the terms of GNU General Public License v3 or any other later versions.
|
|
|
|
phasicFlow is distributed to help others in their research in the field of
|
|
granular and multiphase flows, but WITHOUT ANY WARRANTY; without even the
|
|
implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
|
|
|
|
-----------------------------------------------------------------------------*/
|
|
|
|
#ifndef __positionParticles_hpp__
|
|
#define __positionParticles_hpp__
|
|
|
|
#include "virtualConstructor.hpp"
|
|
#include "Vectors.hpp"
|
|
#include "dictionary.hpp"
|
|
|
|
namespace pFlow
|
|
{
|
|
|
|
class regionBase
|
|
{
|
|
public:
|
|
|
|
regionBase() = default;
|
|
|
|
regionBase(const regionBase&) = default;
|
|
|
|
regionBase& operator =(const regionBase&) = default;
|
|
|
|
virtual ~regionBase() = default;
|
|
|
|
virtual bool isInside(const realx3 point)const = 0;
|
|
|
|
virtual realx3 minPoint()const =0;
|
|
|
|
virtual realx3 maxPoint()const =0;
|
|
|
|
virtual word name()const =0;
|
|
|
|
};
|
|
|
|
template<typename T>
|
|
class region
|
|
:
|
|
public regionBase
|
|
{
|
|
protected:
|
|
|
|
T region_;
|
|
|
|
public:
|
|
|
|
region(const T& rgn)
|
|
:
|
|
region_(rgn)
|
|
{}
|
|
|
|
region(const dictionary& dict)
|
|
:
|
|
region_(dict)
|
|
{}
|
|
|
|
region(const region&) = default;
|
|
|
|
region& operator =(const region&) = default;
|
|
|
|
virtual ~region()=default;
|
|
|
|
bool isInside(const realx3 point) const override
|
|
{
|
|
return region_.isInside(point);
|
|
}
|
|
|
|
realx3 minPoint()const override
|
|
{
|
|
return region_.minPoint();
|
|
}
|
|
|
|
realx3 maxPoint()const override
|
|
{
|
|
return region_.maxPoint();
|
|
}
|
|
|
|
word name()const override
|
|
{
|
|
return region_.typeName();
|
|
}
|
|
|
|
};
|
|
|
|
class positionParticles
|
|
{
|
|
protected:
|
|
|
|
uniquePtr<regionBase> region_ = nullptr;
|
|
|
|
size_t maxNumberOfParticles_ = 10000;
|
|
|
|
Logical mortonSorting_;
|
|
|
|
static const size_t numReports_ = 40;
|
|
|
|
realx3Vector sortByMortonCode(realx3Vector& position)const;
|
|
|
|
public:
|
|
|
|
// - type Info
|
|
TypeInfo("positionParticles");
|
|
|
|
positionParticles(const dictionary& dict);
|
|
|
|
create_vCtor
|
|
(
|
|
positionParticles,
|
|
dictionary,
|
|
(const dictionary& dict),
|
|
(dict)
|
|
);
|
|
|
|
virtual ~positionParticles() = default;
|
|
|
|
//// - Methods
|
|
|
|
virtual label numPoints()const = 0;
|
|
|
|
virtual label size()const = 0;
|
|
|
|
virtual real maxDiameter() const = 0;
|
|
|
|
// - const access to position
|
|
virtual const realx3Vector& position()const = 0;
|
|
|
|
|
|
// - access to position
|
|
virtual realx3Vector& position() = 0;
|
|
|
|
virtual realx3Vector getFinalPosition();
|
|
|
|
static
|
|
uniquePtr<positionParticles> create(const dictionary & dict);
|
|
|
|
};
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
#endif // __positionParticles_hpp__
|