Instance *CreateGraphRefInstance(Instance *instance1, ReferenceGraph *refGraph)
{
int i;
Instance *instance2 =
AllocateInstance(instance1->numVertices, instance1->numEdges);
for(i=0;i<instance1->numVertices;i++)
{
// assign the vertices to the full graph
instance2->vertices[i] = refGraph->vertices[instance1->vertices[i]].map;
}
// reorder indices based on the vertices in the full graph
SortIndices(instance2->vertices, 0, instance2->numVertices-1);
for(i=0;i<instance1->numEdges;i++)
{
// assign the vertices to the full graph
instance2->edges[i] = refGraph->edges[instance1->edges[i]].map;
}
// reorder indices based on the vertices in the full graph
SortIndices(instance2->edges, 0, instance2->numEdges-1);
return instance2;
}
void SortIndices(ULONG *A, long p, long r)
{
long q;
if (p < r)
{
q = Partition(A, p, r);
SortIndices(A, p, q-1);
SortIndices(A, q+1, r);
}
}
void CalculateRadii(vtkPolyData* mesh, const arma::sp_mat& edge_lengths, arma::vec& radii) {
const int npts = mesh->GetNumberOfPoints();
const int ncells = mesh->GetNumberOfCells();
radii.set_size(npts);
vtkSmartPointer<vtkIdList> cellIdList = vtkSmartPointer<vtkIdList>::New();
vtkSmartPointer<vtkIdList> pointIdList = vtkSmartPointer<vtkIdList>::New();
//for every face calculate gamma_ijk
arma::vec gamma_percell(ncells);
for (vtkIdType i = 0; i < ncells; i++) {
pointIdList->Reset();
mesh->GetCellPoints(cellIdList->GetId(i), pointIdList);
//!!! assert(pointIdList->GetNumberOfIds() == 3);
int ind_i = pointIdList->GetId(0);
int ind_j = pointIdList->GetId(1);
int ind_k = pointIdList->GetId(2);
SortIndices(i, ind_i, ind_j, ind_k);
const double l_ij = edge_lengths(ind_i, ind_j);
const double l_ki = edge_lengths(ind_k, ind_i);
const double l_jk = edge_lengths(ind_j, ind_k);
gamma_percell(i) = (l_ij + l_ki - l_jk) / 2;
}
//for every vertex get minimum gamma
for (vtkIdType i = 0; i < npts; i++) {
cellIdList->Reset();
mesh->GetPointCells(i, cellIdList);
radii(i) = DBL_MAX;
for (vtkIdType j = 0; j < cellIdList->GetNumberOfIds(); j++) {
if (gamma_percell(j) < radii(i)) radii(i) = gamma_percell(j); //get minimum of gammas from surrounding cells
}
}
}
void CalculateFaceAnglesE2(vtkPolyData* mesh, const arma::vec& edge_lengths, arma::vec& angles) {
const int ncells = mesh->GetNumberOfCells();
angles.set_size(ncells);
vtkSmartPointer<vtkIdList> pointIdList = vtkSmartPointer<vtkIdList>::New();
for (vtkIdType i = 0; i < ncells; i++) {
pointIdList->Reset();
mesh->GetCellPoints(i, pointIdList);
int ind_i = pointIdList->GetId(0);
int ind_j = pointIdList->GetId(1);
int ind_k = pointIdList->GetId(2);
SortIndices(i, ind_i, ind_j, ind_k);
const int l_ij = edge_lengths(ind_i, ind_j);
const int l_ki = edge_lengths(ind_k, ind_i);
const int l_kj = edge_lengths(ind_k, ind_j);
angles(i) = acos((l_ij * l_ij + l_ki * l_ki - l_kj * l_kj) / (2 * l_ki * l_ij));
}
}
std::vector<int> StreamPower::Indexx(std::vector<double>& arr)
{
return SortIndices(arr);
}
