46 |
|
#include "math/SquareMatrix3.hpp" |
47 |
|
#include "math/DynamicRectMatrix.hpp" |
48 |
|
#include "primitives/Molecule.hpp" |
49 |
< |
#include "applications/hydrodynamics/MoleculeShape.hpp" |
49 |
> |
#include "applications/hydrodynamics/StuntDoubleShape.hpp" |
50 |
|
#include "utils/any.hpp" |
51 |
|
namespace oopse { |
52 |
|
struct HydrodynamicProps { |
66 |
|
|
67 |
|
class HydrodynamicsModel { |
68 |
|
public: |
69 |
< |
HydrodynamicsModel(Molecule* mol, const DynamicProperty& extraParams) : mol_(mol) {} |
70 |
< |
bool calcHydrodyanmicsProps(double eta); |
69 |
> |
HydrodynamicsModel(StuntDouble* sd, const DynamicProperty& extraParams); |
70 |
> |
bool calcHydrodyanmicsProps(); |
71 |
|
|
72 |
|
Vector3d getDiffCenter(); |
73 |
|
Mat3x3d getTransDiff(); |
76 |
|
void writeBeads(std::ostream& os); |
77 |
|
void writeDiffCenterAndDiffTensor(std::ostream& os); |
78 |
|
protected: |
79 |
< |
Molecule* mol_; |
79 |
> |
StuntDouble* sd_; |
80 |
|
private: |
81 |
|
virtual bool createBeads(std::vector<BeadParam>& beads) = 0; |
82 |
|
|
83 |
|
|
84 |
|
HydrodynamicProps props_; |
85 |
|
std::vector<BeadParam> beads_; |
86 |
+ |
double viscosity_; |
87 |
+ |
double temperature_; |
88 |
|
|
89 |
|
}; |
90 |
|
|