1
Fork 0
turbonss/cxxmph/mphtable.cc

190 lines
6.5 KiB
C++
Raw Normal View History

#include <limits>
2010-10-28 03:17:09 +03:00
#include <iostream>
using std::cerr;
using std::endl;
2010-06-28 22:01:18 +03:00
2010-10-05 17:51:17 +03:00
#include "mphtable.h"
using std::vector;
2010-10-28 03:17:09 +03:00
namespace {
2010-10-28 03:17:09 +03:00
static const cmph_uint8 kUnassigned = 3;
// table used for looking up the number of assigned vertices to a 8-bit integer
static cmph_uint8 kBdzLookupTable[] =
{
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
3, 3, 3, 2, 3, 3, 3, 2, 3, 3, 3, 2, 2, 2, 2, 1,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
3, 3, 3, 2, 3, 3, 3, 2, 3, 3, 3, 2, 2, 2, 2, 1,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
4, 4, 4, 3, 4, 4, 4, 3, 4, 4, 4, 3, 3, 3, 3, 2,
3, 3, 3, 2, 3, 3, 3, 2, 3, 3, 3, 2, 2, 2, 2, 1,
3, 3, 3, 2, 3, 3, 3, 2, 3, 3, 3, 2, 2, 2, 2, 1,
3, 3, 3, 2, 3, 3, 3, 2, 3, 3, 3, 2, 2, 2, 2, 1,
3, 3, 3, 2, 3, 3, 3, 2, 3, 3, 3, 2, 2, 2, 2, 1,
2, 2, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 1, 1, 1, 0
};
2010-10-05 17:51:17 +03:00
2010-10-28 03:17:09 +03:00
} // anonymous namespace
namespace cxxmph {
2010-11-05 08:40:15 +02:00
const cmph_uint8 MPHTable::valuemask[] = { 0xfc, 0xf3, 0xcf, 0x3f};
2010-10-29 09:26:37 +03:00
void MPHTable::clear() {
// TODO(davi) impolement me
}
2010-10-28 03:17:09 +03:00
bool MPHTable::GenerateQueue(
TriGraph* graph, vector<cmph_uint32>* queue_output) {
2010-10-05 17:51:17 +03:00
cmph_uint32 queue_head = 0, queue_tail = 0;
2010-10-28 03:17:09 +03:00
cmph_uint32 nedges = m_;
cmph_uint32 nvertices = n_;
2010-10-05 17:51:17 +03:00
// Relies on vector<bool> using 1 bit per element
2010-10-28 05:45:43 +03:00
vector<bool> marked_edge(nedges + 1, false);
2010-10-28 03:17:09 +03:00
vector<cmph_uint32> queue(nvertices, 0);
for (cmph_uint32 i = 0; i < nedges; ++i) {
const TriGraph::Edge& e = graph->edges()[i];
if (graph->vertex_degree()[e[0]] == 1 ||
graph->vertex_degree()[e[1]] == 1 ||
graph->vertex_degree()[e[2]] == 1) {
2010-10-05 17:51:17 +03:00
if (!marked_edge[i]) {
queue[queue_head++] = i;
2010-10-05 17:51:17 +03:00
marked_edge[i] = true;
}
}
}
2010-10-28 05:45:43 +03:00
for (unsigned int i = 0; i < marked_edge.size(); ++i) {
cerr << "vertex with degree " << static_cast<cmph_uint32>(graph->vertex_degree()[i]) << " marked " << marked_edge[i] << endl;
}
for (unsigned int i = 0; i < queue.size(); ++i) {
cerr << "vertex " << i << " queued at " << queue[i] << endl;
}
// At this point queue head is the number of edges touching at least one
// vertex of degree 1.
cerr << "Queue head " << queue_head << " Queue tail " << queue_tail << endl;
2010-11-09 02:02:18 +02:00
graph->DebugGraph();
2010-10-05 17:51:17 +03:00
while (queue_tail != queue_head) {
cmph_uint32 current_edge = queue[queue_tail++];
2010-10-05 17:51:17 +03:00
graph->RemoveEdge(current_edge);
const TriGraph::Edge& e = graph->edges()[current_edge];
2010-10-05 17:51:17 +03:00
for (int i = 0; i < 3; ++i) {
cmph_uint32 v = e[i];
if (graph->vertex_degree()[v] == 1) {
cmph_uint32 first_edge = graph->first_edge()[v];
if (!marked_edge[first_edge]) {
2010-10-05 17:51:17 +03:00
queue[queue_head++] = first_edge;
marked_edge[first_edge] = true;
}
}
}
}
2010-10-28 05:45:43 +03:00
for (unsigned int i = 0; i < queue.size(); ++i) {
cerr << "vertex " << i << " queued at " << queue[i] << endl;
}
int cycles = queue_head - nedges;
if (cycles == 0) queue.swap(*queue_output);
return cycles == 0;
2010-06-28 22:01:18 +03:00
}
2010-10-28 03:17:09 +03:00
void MPHTable::Assigning(
const vector<TriGraph::Edge>& edges, const vector<cmph_uint32>& queue) {
cmph_uint32 current_edge = 0;
2010-11-09 06:29:39 +02:00
vector<bool> marked_vertices(n_ + 1);
2010-10-28 03:17:09 +03:00
// Initialize vector of half nibbles with all bits set.
2010-10-28 05:45:43 +03:00
cmph_uint32 sizeg = static_cast<cmph_uint32>(ceil(n_/4.0));
vector<cmph_uint8>(sizeg, std::numeric_limits<cmph_uint8>::max()).swap(g_);
2010-11-09 06:29:39 +02:00
cmph_uint32 nedges = m_; // for legibility
for (int i = nedges - 1; i + 1 >= 1; --i) {
current_edge = queue[i];
2010-10-28 05:45:43 +03:00
cerr << "Current edge " << current_edge << " at queue pos " << i << endl;
const TriGraph::Edge& e = edges[current_edge];
2010-10-28 05:45:43 +03:00
cerr << "B: " << e[0] << " " << e[1] << " " << e[2] << " -> "
<< get_2bit_value(g_, e[0]) << " "
<< get_2bit_value(g_, e[1]) << " "
<< get_2bit_value(g_, e[2]) << " " << endl;
if (!marked_vertices[e[0]]) {
if (!marked_vertices[e[1]]) {
2010-10-28 03:17:09 +03:00
set_2bit_value(&g_, e[1], kUnassigned);
marked_vertices[e[1]] = true;
}
if (!marked_vertices[e[2]]) {
2010-10-28 03:17:09 +03:00
set_2bit_value(&g_, e[2], kUnassigned);
2010-11-09 06:29:39 +02:00
assert(marked_vertices.size() > e[2]);
marked_vertices[e[2]] = true;
}
2010-10-28 03:17:09 +03:00
set_2bit_value(&g_, e[0], (6 - (get_2bit_value(g_, e[1]) + get_2bit_value(g_, e[2]))) % 3);
marked_vertices[e[0]] = true;
2010-10-28 03:17:09 +03:00
} else if (!marked_vertices[e[1]]) {
if (!marked_vertices[e[2]]) {
set_2bit_value(&g_, e[2], kUnassigned);
marked_vertices[e[2]] = true;
}
2010-10-28 03:17:09 +03:00
set_2bit_value(&g_, e[1], (7 - (get_2bit_value(g_, e[0]) + get_2bit_value(g_, e[2]))) % 3);
marked_vertices[e[1]] = true;
} else {
2010-10-28 03:17:09 +03:00
set_2bit_value(&g_, e[2], (8 - (get_2bit_value(g_, e[0]) + get_2bit_value(g_, e[1]))) % 3);
marked_vertices[e[2]] = true;
2010-10-05 17:51:17 +03:00
}
2010-10-28 05:45:43 +03:00
cerr << "A: " << e[0] << " " << e[1] << " " << e[2] << " -> "
<< get_2bit_value(g_, e[0]) << " "
<< get_2bit_value(g_, e[1]) << " "
<< get_2bit_value(g_, e[2]) << " " << endl;
2010-10-05 17:51:17 +03:00
}
}
2010-10-28 03:17:09 +03:00
void MPHTable::Ranking() {
cmph_uint32 nbytes_total = static_cast<cmph_uint32>(ceil(n_ / 4.0));
cmph_uint32 size = k_ >> 2U;
2010-10-28 03:17:09 +03:00
cmph_uint32 ranktablesize = static_cast<cmph_uint32>(
ceil(n_ / static_cast<double>(k_)));
// TODO(davi) Change swap of member classes for resize + memset to avoid
// fragmentation
vector<cmph_uint32> (ranktablesize).swap(ranktable_);;
cmph_uint32 offset = 0;
cmph_uint32 count = 0;
2010-10-28 05:45:43 +03:00
cmph_uint32 i = 1;
while (1) {
2010-10-28 03:17:09 +03:00
if (i == ranktable_.size()) break;
cmph_uint32 nbytes = size < nbytes_total ? size : nbytes_total;
2010-10-28 03:17:09 +03:00
for (cmph_uint32 j = 0; j < nbytes; ++j) count += kBdzLookupTable[g_[offset + j]];
ranktable_[i] = count;
offset += nbytes;
nbytes_total -= size;
++i;
}
2010-10-05 17:51:17 +03:00
}
2010-10-28 03:17:09 +03:00
cmph_uint32 MPHTable::Rank(cmph_uint32 vertex) const {
cmph_uint32 index = vertex >> b_;
cmph_uint32 base_rank = ranktable_[index];
2010-10-28 03:17:09 +03:00
cmph_uint32 beg_idx_v = index << b_;
2010-10-28 05:45:43 +03:00
cmph_uint32 beg_idx_b = beg_idx_v >> 2;
cmph_uint32 end_idx_b = vertex >> 2;
while (beg_idx_b < end_idx_b) base_rank += kBdzLookupTable[g_[beg_idx_b++]];
beg_idx_v = beg_idx_b << 2;
2010-10-28 05:45:43 +03:00
cerr << "beg_idx_v: " << beg_idx_v << endl;
cerr << "base rank: " << base_rank << endl;
cerr << "G: ";
for (unsigned int i = 0; i < n_; ++i) {
cerr << get_2bit_value(g_, i) << " ";
}
while (beg_idx_v < vertex) {
2010-10-28 05:45:43 +03:00
cerr << get_2bit_value(g_, beg_idx_v) << " ";
if (get_2bit_value(g_, beg_idx_v) != kUnassigned) ++base_rank;
++beg_idx_v;
}
2010-10-28 05:45:43 +03:00
cerr << "Base rank: " << base_rank << endl;
return base_rank;
2010-10-05 17:51:17 +03:00
}
2010-06-28 22:01:18 +03:00
} // namespace cxxmph