TEST(Math, VectorEqual) {
Vector<3> vec1(0);
Vector<3> vec2(0);
ASSERT_EQ(vec1, vec2);
Vector<4> vec3(1);
Vector<4> vec4(1);
ASSERT_EQ(vec3, vec4);
Vector<2> vec5(0);
Vector<2> vec6(1);
ASSERT_NE(vec5, vec6);
}
//----------------------------------------------------------------------------------
//! Constructor
//----------------------------------------------------------------------------------
SyncFields::SyncFields ()
: BaseOp(0, 0)
{
ML_TRACE_IN("SyncFields::SyncFields ()");
handleNotificationOff();
FieldContainer *fieldC = getFieldContainer();
const char* syncModes[] = {"Float","MultiFloat","Vector","MultiVector"};
f_Mode = fieldC->addEnum("mode",syncModes,4);
f_Mode->setEnumValue( FLOAT );
f_On = fieldC->addBool("on");
f_On->setBoolValue( true );
f_Precision = fieldC->addDouble("precision");
f_Precision->setDoubleValue( 1e-3f );
for (unsigned int iFloat=0; iFloat < NFIELDS; ++iFloat ){
std::stringstream ss;
ss << iFloat;
std::string currentName = "float" + ss.str();
m_FloatsFields.push_back(fieldC->addDouble(currentName.c_str()) );
m_FloatsFields[iFloat]->setDoubleValue( 0.0f );
}
for (unsigned int iVector =0; iVector < NFIELDS; ++iVector ){
std::stringstream ss;
ss << iVector;
std::string currentName = "vector" + ss.str();
m_Vec6Fields.push_back( fieldC->addVec6f(currentName.c_str()) );
m_Vec6Fields[iVector]->setVec6fValue( vec6() );
}
m_Updating = false;
handleNotificationOn();
}
void CSODistance::_process()
{
ML_TRACE_IN("void CSODistance::_process()");
_outputXMarkerList->clearList();
// Delete all created CurveData objects
while ( !_outputCurveList->getCurveList().empty() ) {
ML_DELETE( _outputCurveList->getCurveList().back() );
_outputCurveList->getCurveList().pop_back();
}
if (_csoList0 != NULL) {
int nCSOs = _csoList0->numCSO();
double minDistance = 0.0;
double maxDistance = 0.0;
double avgDistance = 0.0;
double stdDevDistance = 0.0;
std::stringstream distances;
distances << _tableHeader << std::endl;
switch ( _modeFld->getEnumValue() ){
case FIRST2 :
{
if ((nCSOs >= 2) && (_csoList0->getCSOAt(0)->getIsFinished()) && (_csoList0->getCSOAt(1)->getIsFinished())) {
std::vector<vec3>pointSet1;
std::vector<vec3>pointSet2;
_csoList0->getCSOAt(0)->fillPathPointCoordinatesFlattened(pointSet1);
_csoList0->getCSOAt(1)->fillPathPointCoordinatesFlattened(pointSet2);
Vector3 minPoint1,minPoint2,maxPoint1,maxPoint2;
MinimalDistancePointClouds* pointSetsMinDist = NULL;
ML_CHECK_NEW(pointSetsMinDist, MinimalDistancePointClouds);
pointSetsMinDist->setPointSets(pointSet1, pointSet2);
pointSetsMinDist->setNumEntries(200);
pointSetsMinDist->computeDistance(minPoint1,minPoint2);
minDistance = sqrt((minPoint1[0]-minPoint2[0])*(minPoint1[0]-minPoint2[0]) +
(minPoint1[1]-minPoint2[1])*(minPoint1[1]-minPoint2[1]) +
(minPoint1[2]-minPoint2[2])*(minPoint1[2]-minPoint2[2]));
_minimumDistanceFld->setFloatValue( static_cast<float>(minDistance));
_outputXMarkerList->appendItem(XMarker( vec6(minPoint1[0],minPoint1[1],minPoint1[2],0.5f,0.0f,0.0f),
vec3(minPoint2[0]-minPoint1[0],minPoint2[1]-minPoint1[1],minPoint2[2]-minPoint1[2]),
0,""));
distances << _csoList0->getCSOAt(0)->getId() << ","
<< _csoList0->getCSOAt(1)->getId() << ","
<< minDistance << ","
<< maxDistance << ","
<< avgDistance << ","
<< stdDevDistance
<< std::endl;
} else {
_minimumDistancePoint1Fld->setVec3fValue(vec3(0.0,0.0,0.0));
_minimumDistancePoint2Fld->setVec3fValue(vec3(0.0,0.0,0.0));
_minimumDistanceFld->setFloatValue(0.0f);
_distancesFld->setStringValue( _tableHeader );
}
break;
}
case INPLANE:
case INDEX:
{
CurveData *outputCurve = new CurveData;
double *yValues = new double[ nCSOs ];
double minDist = ML_DOUBLE_MAX;
vec3 point1;
vec3 point2;
double averageMinDistance = 0.0;
double averageMeanDistance = 0.0;
double averageMaxDistance = 0.0;
for ( int iCSO = 0; iCSO<nCSOs; ++iCSO ){
CSO* currentCSO = _csoList0->getCSOAt( iCSO );
CSO* matchingCSO = _findMatchingCSO(iCSO);
if (!matchingCSO) {continue;}
std::vector<vec3>pointSet1;
std::vector<vec3>pointSet2;
currentCSO->fillPathPointCoordinatesFlattened(pointSet1);
matchingCSO->fillPathPointCoordinatesFlattened(pointSet2);
Vector3 minPoint1,minPoint2,maxPoint1,maxPoint2;
_getDistances( pointSet1, pointSet2,
minDistance,maxDistance,avgDistance,stdDevDistance,
minPoint1,minPoint2,maxPoint1,maxPoint2);
averageMinDistance += minDistance;
averageMeanDistance += avgDistance;
averageMaxDistance += maxDistance;
distances << currentCSO->getId() << ","
<< matchingCSO->getId() << ","
<< minDistance << ","
<< maxDistance << ","
<< avgDistance << ","
<< stdDevDistance
<< std::endl;
if ( minDistance < minDist ){
minDist = minDistance;
point1 = minPoint1;
point2 = minPoint2;
}
_outputXMarkerList->appendItem(XMarker( vec6(minPoint1[0],minPoint1[1],minPoint1[2],0.5f,0.0f,0.0f),
vec3(minPoint2[0]-minPoint1[0],minPoint2[1]-minPoint1[1],minPoint2[2]-minPoint1[2]),
currentCSO->getId(),""));
switch ( _curveStatistic->getEnumValue() ){
case MIN:
yValues[ iCSO ] = minDistance;
break;
case MAX:
yValues[ iCSO ] = maxDistance;
break;
case MEAN:
yValues[ iCSO ] = avgDistance;
break;
case STDEV:
yValues[ iCSO ] = stdDevDistance;
break;
default:
break;
}
} // iCSO
averageMinDistance /= (nCSOs != 0 ? nCSOs : 1.0);
averageMeanDistance /= (nCSOs != 0 ? nCSOs : 1.0);
averageMaxDistance /= (nCSOs != 0 ? nCSOs : 1.0);
outputCurve->setY( nCSOs, yValues);
delete[] yValues;
_outputCurveList->getCurveList().push_back( outputCurve );
_distancesFld->setStringValue( distances.str() );
_minimumDistancePoint1Fld->setVec3fValue(point1);
_minimumDistancePoint2Fld->setVec3fValue(point2);
_minimumDistanceFld->setFloatValue( static_cast<float>(minDistance) );
_AverageMinimumDistanceFld->setDoubleValue( averageMinDistance );
_AverageMeanDistanceFld->setDoubleValue( averageMeanDistance );
_AverageMaxDistanceFld->setDoubleValue( averageMaxDistance );
break;
}
default:
break;
}
} else {
_minimumDistancePoint1Fld->setVec3fValue(vec3(0.0,0.0,0.0));
_minimumDistancePoint2Fld->setVec3fValue(vec3(0.0,0.0,0.0));
_minimumDistanceFld->setFloatValue(0.0f);
_distancesFld->setStringValue( _tableHeader );
}
_outputXMarkerListFld->setBaseValue(_outputXMarkerList);
_outputCurveListFld->touch();
}
