Skip to content
Snippets Groups Projects
Commit 1d3be1e4 authored by Mohcine Chraibi's avatar Mohcine Chraibi
Browse files

change definition of method D velocity

not actually sure if this is ok.
But, the actual definition looks not correct.
parent d273b6d2
Branches
No related tags found
No related merge requests found
...@@ -366,7 +366,7 @@ double Method_D::GetVoronoiVelocity(const vector<polygon_2d>& polygon, const vec ...@@ -366,7 +366,7 @@ double Method_D::GetVoronoiVelocity(const vector<polygon_2d>& polygon, const vec
intersection(measureArea, polygon_iterator, v); intersection(measureArea, polygon_iterator, v);
if(!v.empty()) if(!v.empty())
{ {
meanV+=(Velocity[temp]*area(v[0])/area(measureArea)); meanV+=Velocity[temp]*area(v[0])/area(polygon_iterator);
if((area(v[0]) - area(polygon_iterator))>J_EPS) if((area(v[0]) - area(polygon_iterator))>J_EPS)
{ {
std::cout<<"this is a wrong result in calculating velocity\t"<<area(v[0])<<'\t'<<area(polygon_iterator)<< " (diff=" << area(v[0]) - area(polygon_iterator) << ")"<< std::endl; std::cout<<"this is a wrong result in calculating velocity\t"<<area(v[0])<<'\t'<<area(polygon_iterator)<< " (diff=" << area(v[0]) - area(polygon_iterator) << ")"<< std::endl;
...@@ -374,6 +374,8 @@ double Method_D::GetVoronoiVelocity(const vector<polygon_2d>& polygon, const vec ...@@ -374,6 +374,8 @@ double Method_D::GetVoronoiVelocity(const vector<polygon_2d>& polygon, const vec
} }
temp++; temp++;
} }
meanV = meanV/area(measureArea)*CMtoM*CMtoM;
return meanV; return meanV;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment