void PLUDecomposite(PMatrix *P) {
for (int k = 1; k <= n()-1; ++k) {
int m = k;
for (int i = k+1; i <= n(); ++i) {
if (fabs(M(i, k)) > fabs(M(m, k))) {
m = i;
}
}
if (fabs(M(m, k) + 0.0) < 1e-15) {
_reg = false;
return;
}
if (m != k) {
swapRows(m, k);
P -> swapRows(m, k);
}
for (int i = k+1; i <= n(); ++i) {
setM(i, k, M(i, k) / M(k, k));
for (int j = k+1; j <= n(); ++j) {
setM(i, j, M(i, j) - M(i, k) * M (k, j));
}
}
}
_reg = fabs(M(n(), n())) >= 1e-15;
}
void swapRows(int i1, int i2) {
for (int i = 1; i <= n(); ++i) {
double tmp = M(i1, i);
setM(i1, i, M(i2, i));
setM(i2, i, tmp);
}
}
bool LL1Parser::constructLL1Table()
{
m.clear();
g.computeFollow();
vector<string> v = g.v;
vector<string> t = g.t;
vector<Production> p = g.p;
bool succeed = true;
m.assign(v.size(), vector<int>(t.size() + 1, -1));
// row: variables, column: terminals + $
// elements: production indices
for(size_t i = 0; i < p.size(); ++i)
{
// A -> alpha
string A = p[i].left;
int AIndex = g.getVariableIndex(A);
vector<string> firstAlpha = g.computeFirst(p[i].right);
for(size_t j = 0; j < firstAlpha.size(); ++j)
{
string a = firstAlpha[j];
if(a != "")
{
if(!setM(AIndex, g.getTerminalIndex(a), i))
{
succeed = false;
}
}
}
if(in(string(""), firstAlpha))
{
vector<string> followA = g.follow[g.getFollowIndex(A)].second;
for(size_t j = 0; j < followA.size(); ++j)
{
if(followA[j] != "$")
{
if(!setM(AIndex, g.getTerminalIndex(followA[j]), i))
{
succeed = false;
}
}
else
{
if(!setM(AIndex, t.size(), i))
{
succeed = false;
}
}
}
}
}
return succeed;
}
Mahalanobis::Mahalanobis(arma::mat M) {
setM(M / M(0,0));
}
static void
create_markerss (void)
{
GtkWidget *window = NULL;
GtkWidget *box1;
GtkWidget *box2;
GtkWidget *close_button;
GtkWidget *box;
GtkWidget *table;
GtkWidget *label;
GtkWidget *separator;
gfloat *X;
gfloat *Y;
gint i;
GtkDataboxGraph *graph;
GdkColor color;
window = gtk_window_new (GTK_WINDOW_TOPLEVEL);
gtk_widget_set_size_request (window, 600, 600);
g_signal_connect (GTK_OBJECT (window), "destroy",
G_CALLBACK (gtk_main_quit), NULL);
gtk_window_set_title (GTK_WINDOW (window), "GtkDatabox: Markers");
gtk_container_set_border_width (GTK_CONTAINER (window), 0);
box1 = gtk_vbox_new (FALSE, 0);
gtk_container_add (GTK_CONTAINER (window), box1);
label =
gtk_label_new
("Some markerss (green) are shown below.\nThe text corresponds to the position of the text, relative to the markers.\n\nThere are also some texts (yellow).\n\nData points are red\n");
gtk_box_pack_start (GTK_BOX (box1), label, FALSE, FALSE, 0);
separator = gtk_hseparator_new ();
gtk_box_pack_start (GTK_BOX (box1), separator, FALSE, FALSE, 0);
gtk_databox_create_box_with_scrollbars_and_rulers (&box, &table,
TRUE, TRUE, TRUE, TRUE);
gtk_box_pack_start (GTK_BOX (box1), table, TRUE, TRUE, 0);
color.red = 16383;
color.green = 16383;
color.blue = 16383;
gtk_widget_modify_bg (box, GTK_STATE_NORMAL, &color);
X = g_new0 (gfloat, 20 + 4 + 5);
Y = g_new0 (gfloat, 20 + 4 + 5);
for (i = 0; i < 20 + 4 + 5; ++i)
{
X[i] = i % 4;
Y[i] = i / 4;
}
X[20] = -0.5;
Y[20] = -0.5;
X[21] = -0.5;
Y[21] = 4.5;
X[22] = 3.5;
Y[22] = -0.5;
X[23] = 3.5;
Y[23] = 4.5;
X[24] = -1.;
Y[24] = -1.;
X[25] = -1.;
Y[25] = +5.;
X[26] = +4.;
Y[26] = -1.;
X[27] = +4.;
Y[27] = +5.;
X[28] = +1.5;
Y[28] = +2.;
color.red = 65535;
color.green = 0;
color.blue = 0;
graph = gtk_databox_points_new (20 + 4 + 5, X, Y, &color, 3);
gtk_databox_graph_add (GTK_DATABOX (box), graph);
color.red = 0;
color.green = 65535;
color.blue = 0;
graph = gtk_databox_markers_new (20, X, Y, &color, 15,
GTK_DATABOX_MARKERS_TRIANGLE);
gtk_databox_graph_add (GTK_DATABOX (box), graph);
setM (graph, 0, GTK_DATABOX_MARKERS_W, GTK_DATABOX_MARKERS_TEXT_W, "West", FALSE);
setM (graph, 1, GTK_DATABOX_MARKERS_N, GTK_DATABOX_MARKERS_TEXT_N, "North", FALSE);
setM (graph, 2, GTK_DATABOX_MARKERS_S, GTK_DATABOX_MARKERS_TEXT_S, "South", FALSE);
setM (graph, 3, GTK_DATABOX_MARKERS_E, GTK_DATABOX_MARKERS_TEXT_E, "East", FALSE);
setM (graph, 4, GTK_DATABOX_MARKERS_W, GTK_DATABOX_MARKERS_TEXT_W, "West", TRUE);
setM (graph, 5, GTK_DATABOX_MARKERS_N, GTK_DATABOX_MARKERS_TEXT_N, "North", TRUE);
setM (graph, 6, GTK_DATABOX_MARKERS_S, GTK_DATABOX_MARKERS_TEXT_S, "South", TRUE);
setM (graph, 7, GTK_DATABOX_MARKERS_E, GTK_DATABOX_MARKERS_TEXT_E, "East", TRUE);
setM (graph, 8, GTK_DATABOX_MARKERS_E, GTK_DATABOX_MARKERS_TEXT_W, "West", TRUE);
setM (graph, 9, GTK_DATABOX_MARKERS_S, GTK_DATABOX_MARKERS_TEXT_N, "North", TRUE);
setM (graph, 10, GTK_DATABOX_MARKERS_N, GTK_DATABOX_MARKERS_TEXT_S, "South", TRUE);
setM (graph, 11, GTK_DATABOX_MARKERS_W, GTK_DATABOX_MARKERS_TEXT_E, "East", TRUE);
setM (graph, 12, GTK_DATABOX_MARKERS_W, GTK_DATABOX_MARKERS_TEXT_S, "South", TRUE);
setM (graph, 13, GTK_DATABOX_MARKERS_N, GTK_DATABOX_MARKERS_TEXT_W, "West", TRUE);
setM (graph, 14, GTK_DATABOX_MARKERS_S, GTK_DATABOX_MARKERS_TEXT_E, "East", TRUE);
setM (graph, 15, GTK_DATABOX_MARKERS_E, GTK_DATABOX_MARKERS_TEXT_N, "North", TRUE);
setM (graph, 16, GTK_DATABOX_MARKERS_W, GTK_DATABOX_MARKERS_TEXT_SW, "SW", TRUE);
setM (graph, 17, GTK_DATABOX_MARKERS_N, GTK_DATABOX_MARKERS_TEXT_NW, "NW", TRUE);
setM (graph, 18, GTK_DATABOX_MARKERS_S, GTK_DATABOX_MARKERS_TEXT_SE, "SE", TRUE);
setM (graph, 19, GTK_DATABOX_MARKERS_E, GTK_DATABOX_MARKERS_TEXT_NE, "NE", TRUE);
graph = gtk_databox_markers_new (4, X + 20, Y + 20, &color, 1,
GTK_DATABOX_MARKERS_DASHED_LINE);
gtk_databox_graph_add (GTK_DATABOX (box), graph);
setL (graph, 0, GTK_DATABOX_MARKERS_N, GTK_DATABOX_MARKERS_TEXT_W, "West", TRUE);
setL (graph, 1, GTK_DATABOX_MARKERS_E, GTK_DATABOX_MARKERS_TEXT_S, "South", TRUE);
setL (graph, 2, GTK_DATABOX_MARKERS_W, GTK_DATABOX_MARKERS_TEXT_N, "North", TRUE);
setL (graph, 3, GTK_DATABOX_MARKERS_S, GTK_DATABOX_MARKERS_TEXT_E, "East", TRUE);
color.red = 65535;
color.green = 65535;
color.blue = 0;
graph = gtk_databox_markers_new (5, X + 20 + 4, Y + 20 + 4, &color, 1,
GTK_DATABOX_MARKERS_NONE);
gtk_databox_graph_add (GTK_DATABOX (box), graph);
setT (graph, 0, GTK_DATABOX_MARKERS_TEXT_N, "North", FALSE);
setT (graph, 1, GTK_DATABOX_MARKERS_TEXT_E, "East", FALSE);
setT (graph, 2, GTK_DATABOX_MARKERS_TEXT_S, "South", FALSE);
setT (graph, 3, GTK_DATABOX_MARKERS_TEXT_W, "West", FALSE);
setT (graph, 4, GTK_DATABOX_MARKERS_TEXT_CENTER, "Centered", TRUE);
gtk_databox_auto_rescale (GTK_DATABOX (box), 0.05);
separator = gtk_hseparator_new ();
gtk_box_pack_start (GTK_BOX (box1), separator, FALSE, TRUE, 0);
box2 = gtk_vbox_new (FALSE, 10);
gtk_container_set_border_width (GTK_CONTAINER (box2), 10);
gtk_box_pack_end (GTK_BOX (box1), box2, FALSE, TRUE, 0);
close_button = gtk_button_new_with_label ("close");
g_signal_connect_swapped (GTK_OBJECT (close_button), "clicked",
G_CALLBACK (gtk_main_quit), GTK_OBJECT (box));
gtk_box_pack_start (GTK_BOX (box2), close_button, TRUE, TRUE, 0);
gtk_widget_set_can_default(close_button, GTK_CAN_DEFAULT);
gtk_widget_grab_default (close_button);
gtk_widget_show_all (window);
}