|| || || || || ||

---

ComputeMatrix

Matrix::Mat* Impala::Application::Precompute::ComputeMatrix ( DistributedAccess *  da, Feature::FeatureTable *  devel, Feature::FeatureTable *  test, String  distanceFunction   ) Input one or two feature tables and out comes the kernel distance matrix. If there are two tables then the first is the devel set and the second is the test set. The testset always makes up the rows of the matrix while the devel set makes up the columns of the matrix. So variable horizontal (which stores the part of the develfeatures) loops over the columns. Please note that the function parameters row and column specify the place of this node in the partcount x partcount 'grid' of nodes. So therefore horizontal = GetPartial(f1, column, partcount); Note: : this function deletes f and f2 : f and f2 may point to the same table. Therefore care is taken not to delete it twice. Definition at line 149 of file mainPrecomputeKernelMatrix.cpp. References Impala::Core::Vector::Apply(), Impala::Core::Vector::Chi2Distance(), Impala::Core::Vector::Chi2DistanceSSE(), Impala::Core::Vector::DotProduct(), Impala::Core::Table::TableTem< Col1T, Col2T, Col3T, Col4T, Col5T, Col6T, Col7T, Col8T, Col9T >::GetColumn2(), GetPartial(), Impala::Core::Vector::HistogramIntersectionDouble(), ILOG_VAR, Impala::Core::Vector::VectorSet< ArrayT >::SetSize(), and Impala::Core::Table::Table::Size(). Referenced by Precompute(). { ILOG_VAR(Impala.Application.Precompute.ComputeMatrix); Feature::FeatureTable* horizontal; Feature::FeatureTable* vertical; if (da->GetTotalPartCount() == 1) { horizontal = devel; vertical = test; } else { int myPart = Link::Mpi::MyId(); int column = da->GetColumnPartOfPart(myPart); int row = da->GetRowPartOfPart(myPart); horizontal = GetPartial(devel, column, da->GetColumnPartCount()); vertical = GetPartial(test, row, da->GetRowPartCount()); delete devel; if (test != devel) delete test; } // fill a matrix with chi2 distance typedef Feature::FeatureTable::ColumnVectorSet VectorSet; VectorSet* vectorsH = horizontal->GetColumn2(); VectorSet* vectorsV = vertical->GetColumn2(); vectorsH->SetSize(horizontal->Size()); vectorsV->SetSize(vertical->Size()); Matrix::Mat* distanceMatrix = 0; if (distanceFunction == "dot") { distanceMatrix = Apply(&Vector::DotProduct, vectorsH, vectorsV); } else if (distanceFunction == "histogramintersection") { distanceMatrix = Apply(&Vector::HistogramIntersectionDouble, vectorsH, vectorsV); } else { #ifdef SSE_USED distanceMatrix = Apply(&Vector::Chi2DistanceSSE, vectorsH, vectorsV); #else distanceMatrix = Apply(&Vector::Chi2Distance, vectorsH, vectorsV); #endif } delete horizontal; if (horizontal != vertical) delete vertical; return distanceMatrix; } Here is the call graph for this function:

---