static void show_cb(Widget w, XtPointer data, XtPointer cb)
{
if (!rens_initialized) {
char r[30];
sprintf(r, "widget=%ld", (long)drawa);
dp_init_ren("myrena", "gl", r, "");
fprintf(stderr, "%s: renderer A initialized\n",exename);
dp_print_ren("myrena");
sprintf(r, "widget=%ld", (long)drawb);
dp_init_ren("myrenb", "gl", r, "");
fprintf(stderr, "%s: renderer B initialized\n",exename);
dp_print_ren("myrenb");
rens_initialized= 1;
build_geom();
}
fprintf(stderr, "%s: about to render\n",exename);
rot_step= 0;
XtAddWorkProc(my_work_proc,NULL);
}
Graph::Graph(const int *cartan,
const std::vector<Word>& gens, //namesake
const std::vector<Word>& v_cogens,
const std::vector<Word>& e_gens,
const std::vector<Word>& f_gens,
const Vect& weights)
{
//define symmetry group relations
std::vector<Word> words = words_from_cartan(cartan);
{
const Logging::fake_ostream& os = logger.debug();
os << "relations =";
for (int w=0; w<6; ++w) {
Word& word = words[w];
os << "\n ";
for (unsigned i=0; i<word.size(); ++i) {
os << word[i];
}
}
os |0;
}
//check vertex stabilizer generators
{
const Logging::fake_ostream& os = logger.debug();
os << "v_cogens =";
for (unsigned w=0; w<v_cogens.size(); ++w) {
const Word& jenn = v_cogens[w]; //namesake
os << "\n ";
for (unsigned t=0; t<jenn.size(); ++t) {
int j = jenn[t];
os << j;
Assert (0<=j and j<4,
"generator out of range: letter w["
<< w << "][" << t << "] = " << j );
}
}
os |0;
}
//check edge generators
{
const Logging::fake_ostream& os = logger.debug();
os << "e_gens =";
for (unsigned w=0; w<e_gens.size(); ++w) {
const Word& edge = e_gens[w];
os << "\n ";
for (unsigned t=0; t<edge.size(); ++t) {
int j = edge[t];
os << j;
Assert (0<=j and j<4,
"generator out of range: letter w["
<< w << "][" << t << "] = " << j );
}
}
os |0;
}
//check face generators
{
const Logging::fake_ostream& os = logger.debug();
os << "f_gens =";
for (unsigned w=0; w<f_gens.size(); ++w) {
const Word& face = f_gens[w];
os << "\n ";
for (unsigned t=0; t<face.size(); ++t) {
int j = face[t];
os << j;
Assert (0<=j and j<4,
"generator out of range: letter w["
<< w << "][" << t << "] = " << j );
}
}
os |0;
}
//build symmetry group
Group group(words);
logger.debug() << "group.ord = " << group.ord |0;
//build subgroup
std::vector<int> subgroup; subgroup.push_back(0);
std::set<int> in_subgroup; in_subgroup.insert(0);
for (unsigned g=0; g<subgroup.size(); ++g) {
int g0 = subgroup[g];
for (unsigned j=0; j<gens.size(); ++j) {
int g1 = group.left(g0,gens[j]);
if (in_subgroup.find(g1) != in_subgroup.end()) continue;
subgroup.push_back(g1);
in_subgroup.insert(g1);
}
}
logger.debug() << "subgroup.ord = " << subgroup.size() |0;
//build cosets and count ord
std::map<int,int> coset; //maps group elements to cosets
ord = 0; //used as coset number
for (unsigned g=0; g<subgroup.size(); ++g) {
int g0 = subgroup[g];
if (coset.find(g0) != coset.end()) continue;
int c0 = ord++;
coset[g0] = c0;
std::vector<int> members(1, g0);
std::vector<int> others(0);
for (unsigned i=0; i<members.size(); ++i) {
int g1 = members[i];
for (unsigned w=0; w<v_cogens.size(); ++w) {
int g2 = group.left(g1, v_cogens[w]);
if (coset.find(g2) != coset.end()) continue;
coset[g2] = c0;
members.push_back(g2);
}
}
}
logger.info() << "cosets table built: " << " ord = " << ord |0;
//build edge lists
std::vector<std::set<int> > neigh(ord);
for (unsigned g=0; g<subgroup.size(); ++g) {
int g0 = subgroup[g];
int c0 = coset[g0];
for (unsigned w=0; w<e_gens.size(); ++w) {
int g1 = group.left(g0, e_gens[w]);
Assert (in_subgroup.find(g1) != in_subgroup.end(),
"edge leaves subgroup");
int c1 = coset[g1];
if (c0 != c1) neigh[c0].insert(c1);
}
}
// make symmetric
for (int c0=0; c0<ord; ++c0) {
const std::set<int>& n = neigh[c0];
for (std::set<int>::iterator c1=n.begin(); c1!=n.end(); ++c1) {
neigh[*c1].insert(c0);
}
}
// build edge table
adj.resize(ord);
for (int c=0; c<ord; ++c) {
adj[c].insert(adj[c].begin(), neigh[c].begin(), neigh[c].end());
}
neigh.clear();
deg = adj[0].size();
logger.info() << "edge table built: deg = " << deg |0;
//define faces
for (unsigned g=0; g<f_gens.size(); ++g) {
const Word& face = f_gens[g];
logger.debug() << "defining faces on " << face |0;
Logging::IndentBlock block;
//define basic face in group
Ring basic(1,0);
// g = 0;
int g0 = 0;
for (unsigned c=0; true; ++c) {
g0 = group.left(g0, face[c%face.size()]);
if (c >= face.size() and g0 == 0) break;
if (in_subgroup.find(g0) != in_subgroup.end() and g0 != basic.back()) {
basic.push_back(g0);
}
}
for (unsigned c=0; c<basic.size(); ++c) {
logger.debug() << " corner: " << basic[c] |0;
}
logger.debug() << "sides/face (free) = " << basic.size() |0;
//build orbit of basic face
std::vector<Ring> faces_g; faces_g.push_back(basic);
FaceRecognizer recognized; recognized(basic);
for (unsigned i=0; i<faces_g.size(); ++i) {
const Ring f = faces_g[i];
for (unsigned j=0; j<gens.size(); ++j) {
//right action of group on faces
Ring f_j(f.size());
for (unsigned c=0; c<f.size(); ++c) {
f_j[c] = group.right(f[c],gens[j]);
}
//add face
if (not recognized(f_j)) {
faces_g.push_back(f_j);
//logger.debug() << "new face: " << f_j |0;
} else {
//logger.debug() << "old face: " << f_j|0;
}
}
}
//hom face down to quotient graph
recognized.clear();
for (unsigned f=0; f<faces_g.size(); ++f) {
const Ring face_g = faces_g[f];
Ring face;
face.push_back(coset[face_g[0]]);
for (unsigned i=1; i<face_g.size(); ++i) {
int c = coset[face_g[i]];
if (c != face.back() and c != face[0]) {
face.push_back(c);
}
}
if (face.size() < 3) continue;
if (not recognized(face)) {
faces.push_back(face);
}
}
}
ord_f = faces.size();
logger.info() << "faces defined: order = " << ord_f |0;
//define vertex coset
std::vector<Word> vertex_coset;
for (unsigned g=0; g<subgroup.size(); ++g) {
int g0 = subgroup[g];
if (coset[g0]==0) vertex_coset.push_back(group.parse(g0));
}
//build geometry
std::vector<Mat> gen_reps(gens.size());
points.resize(ord);
build_geom(cartan, vertex_coset, gens, v_cogens, weights,
gen_reps, points[0]);
std::vector<int> pointed(ord,0);
pointed[0] = true;
logger.debug() << "geometry built" |0;
//build point sets
std::vector<int> reached(1,0);
std::set<int> is_reached;
is_reached.insert(0);
for (unsigned g=0; g<subgroup.size(); ++g) {
int g0 = reached[g];
for (unsigned j=0; j<gens.size(); ++j) {
int g1 = group.right(g0,gens[j]);
if (is_reached.find(g1) == is_reached.end()) {
if (not pointed[coset[g1]]) {
vect_mult(gen_reps[j], points[coset[g0]],
points[coset[g1]]);
pointed[coset[g1]] = true;
}
reached.push_back(g1);
is_reached.insert(g1);
}
}
}
logger.debug() << "point set built." |0;
//build face normals
normals.resize(ord_f);
for (int f=0; f<ord_f; ++f) {
Ring& face = faces[f];
Vect &a = points[face[0]];
Vect &b = points[face[1]];
Vect &c = points[face[2]];
Vect &n = normals[f];
cross4(a,b,c, n);
normalize(n);
/*
Assert1(fabs(inner(a,n)) < 1e-6,
"bad normal: <n,a> = " << fabs(inner(a,n)));
Assert1(fabs(inner(b,n)) < 1e-6,
"bad normal: <n,b> = " << fabs(inner(b,n)));
Assert1(fabs(inner(c,n)) < 1e-6,
"bad normal: <n,b> = " << fabs(inner(c,n)));
*/
}
logger.debug() << "face normals built." |0;
}
main(int argc, char *argv[])
{
int i;
exename= argv[0];
if (argc>1) { /* Check for command line arguments */
i= 1;
while ((--argc > 0) && (argv[i][0]=='-')) {
int c= argv[i][1];
switch (c) {
case 'a':
automanage= 1;
fprintf(stderr,"%s: Windows will be automanaged.\n",argv[0]);
break;
default:
fprintf(stderr,"%s: illegal option %c\n",argv[0],c);
fprintf(stderr,"usage: %s [-a]\n",argv[0]);
argc= 0;
break;
}
i++;
}
}
if (!automanage) {
Arg args[16];
Widget toplevel, main_window, menu_bar, menu_pane, button;
Widget draw_form;
XmString title_string;
int width=300, height=300;
i = 0;
toplevel = XtInitialize("p3d", "xpainter", NULL, 0, &i, NULL);
i= 0;
main_window= XmCreateMainWindow(toplevel, "main_window", args, i);
i= 0;
menu_bar= XmCreateMenuBar(main_window, "menu_bar", args, i);
XtManageChild(menu_bar);
i= 0;
menu_pane= XmCreatePulldownMenu(menu_bar, "menu_pane", args, i);
title_string= XmStringCreateSimple("File");
i= 0;
XtSetArg(args[i], XmNsubMenuId, menu_pane); i++;
XtSetArg(args[i], XmNlabelString, title_string); i++;
button= XmCreateCascadeButton(menu_bar,"file",args,i);
if (title_string) XmStringFree(title_string);
XtManageChild(button);
title_string= XmStringCreateSimple("Show");
i= 0;
XtSetArg(args[i], XmNlabelString, title_string); i++;
button= XmCreatePushButtonGadget(menu_pane,"show",args,i);
XtAddCallback(button,XmNactivateCallback,show_cb,NULL);
if (title_string) XmStringFree(title_string);
XtManageChild(button);
title_string= XmStringCreateSimple("Exit");
i= 0;
XtSetArg(args[i], XmNlabelString, title_string); i++;
button= XmCreatePushButtonGadget(menu_pane,"exit",args,i);
XtAddCallback(button,XmNactivateCallback,exit_cb,NULL);
if (title_string) XmStringFree(title_string);
XtManageChild(button);
i= 0;
draw_form= XmCreateForm(main_window,"draw_form", args, i);
XtManageChild(draw_form);
i = 0;
XtSetArg(args[i], XtNwidth, width); i++;
XtSetArg(args[i], XtNheight, height); i++;
XtSetArg(args[i], XmNtopAttachment, XmATTACH_FORM); i++;
XtSetArg(args[i], XmNbottomAttachment, XmATTACH_FORM); i++;
XtSetArg(args[i], XmNleftAttachment, XmATTACH_FORM); i++;
#ifdef USE_OPENGL
XtSetArg(args[i], GLwNrgba, TRUE); i++;
XtSetArg(args[i], GLwNdoublebuffer, TRUE); i++;
XtSetArg(args[i], GLwNdepthSize, 16); i++;
XtSetArg(args[i], GLwNallocateBackground, TRUE); i++;
XtSetArg(args[i], GLwNinstallBackground, TRUE); i++;
drawa= GLwCreateMDrawingArea(draw_form, "spina", args, i);
#else
XtSetArg(args[i], GlxNglxConfig, glxconfig); i++;
drawa= GlxCreateMDraw(draw_form, "spina", args, i);
#endif
XtManageChild(drawa);
i = 0;
XtSetArg(args[i], XtNwidth, width); i++;
XtSetArg(args[i], XtNheight, height); i++;
XtSetArg(args[i], XmNtopAttachment, XmATTACH_FORM); i++;
XtSetArg(args[i], XmNbottomAttachment, XmATTACH_FORM); i++;
XtSetArg(args[i], XmNleftAttachment, XmATTACH_WIDGET); i++;
XtSetArg(args[i], XmNleftWidget, drawa); i++;
XtSetArg(args[i], XmNrightAttachment, XmATTACH_FORM); i++;
#ifdef USE_OPENGL
XtSetArg(args[i], GLwNrgba, TRUE); i++;
XtSetArg(args[i], GLwNdoublebuffer, TRUE); i++;
XtSetArg(args[i], GLwNdepthSize, 16); i++;
XtSetArg(args[i], GLwNallocateBackground, TRUE); i++;
XtSetArg(args[i], GLwNinstallBackground, TRUE); i++;
drawb= GLwCreateMDrawingArea(draw_form, "spinb", args, i);
#else
XtSetArg(args[i], GlxNglxConfig, glxconfig); i++;
drawb= GlxCreateMDraw(draw_form, "spinb", args, i);
#endif
XtManageChild(drawb);
XtManageChild(main_window);
XtRealizeWidget(toplevel);
}
else { /* automanage renderers */
dp_init_ren("myrena", "gl","",
"title=\"Ren A\",size=\"256x256+200+200\"");
dp_init_ren("myrenb", "gl","",
"title=\"Ren B\",size=\"300x400+500+500\"");
fprintf(stderr, "%s: renderer initialized\n",exename);
dp_print_ren("myrena");
dp_print_ren("myrenb");
rens_initialized= 1;
build_geom();
}
if (!automanage) {
XtMainLoop();
}
else {
fprintf(stderr, "%s: about to render\n",exename);
for (i=0; i<180; i++) make_new_snaps();
finish_and_shutdown();
}
}
