1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
// equi2.cmd -- to color edges subject to equiangulation
// Programmer: Ken Brakke, brakke@susqu.edu, http://www.susqu.edu/brakke
// useage: showequi
showequi := {
local ecount, al,bl,cl,dl,el,nn;
ecount := 0;
foreach edge ee do
{ al := ee.length;
nn := 1;
foreach ee.facet ff do
{ foreach ff.edge fe do
{ if ( fe.id != ee.id ) then
{ if ( nn == 1 ) then bl := fe.length
else if ( nn == 2 ) then cl := fe.length
else if ( nn == 3 ) then dl := fe.length
else if ( nn == 4 ) then el := fe.length;
nn := nn + 1;
}
}
};
if ( nn == 5 ) then
if (bl^2+cl^2-al^2)/bl/cl + (dl^2+el^2-al^2)/dl/el < 0 then
{ set ee color green;
ecount := 1 + ecount
}
};
printf "Found %g edges.\n",ecount;
show_expr edges where 1
}
|
