Skip to content
Snippets Groups Projects
Commit 93be8a6a authored by Konstantin Kutscher's avatar Konstantin Kutscher
Browse files

add tolerance

parent 6d25b921
No related branches found
No related tags found
No related merge requests found
...@@ -92,7 +92,7 @@ void CreateDemObjectsCoProcessor::createGeoObjects() ...@@ -92,7 +92,7 @@ void CreateDemObjectsCoProcessor::createGeoObjects()
for (int i = 0; i < size; i++) for (int i = 0; i < size; i++)
{ {
SPtr<GbSphere3D> sphere = std::dynamic_pointer_cast<GbSphere3D>(geoObjectPrototypeVector[i]); SPtr<GbSphere3D> sphere = std::dynamic_pointer_cast<GbSphere3D>(geoObjectPrototypeVector[i]);
if (demCoProcessor->isSpheresIntersection(sphere->getX1Centroid(), sphere->getX2Centroid(), sphere->getX3Centroid(), sphere->getRadius()*2.0*(1.0-toleranz))) if (demCoProcessor->isSpheresIntersection(sphere->getX1Centroid(), sphere->getX2Centroid(), sphere->getX3Centroid(), sphere->getRadius()*2.0*(1.0-tolerance)))
{ {
continue; continue;
} }
......
...@@ -32,6 +32,8 @@ public: ...@@ -32,6 +32,8 @@ public:
void addGeoObject(SPtr<GbObject3D> geoObjectPrototype, Vector3D initalVelocity); void addGeoObject(SPtr<GbObject3D> geoObjectPrototype, Vector3D initalVelocity);
void clearGeoObjects(); void clearGeoObjects();
void createGeoObjects(); void createGeoObjects();
double getToleranz() const { return tolerance; }
void setToleranz(double val) { tolerance = val; }
protected: protected:
private: private:
SPtr<Communicator> comm; SPtr<Communicator> comm;
...@@ -42,7 +44,7 @@ private: ...@@ -42,7 +44,7 @@ private:
SPtr<BCAdapter> velocityBcParticleAdapter; SPtr<BCAdapter> velocityBcParticleAdapter;
SPtr<Reconstructor> extrapolationReconstructor; SPtr<Reconstructor> extrapolationReconstructor;
int demCounter; int demCounter;
double toleranz; double tolerance;
#ifdef TIMING #ifdef TIMING
UbTimer timer; UbTimer timer;
#endif #endif
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment