Bevel

Edge Bevel: Requirements

If you have an edge that is part of exactly two faces, an edge bevel should replace that edge with two parallel edges, each moved some distance along their respective faces. Usually the distance moved along each face is the same.

Bevel Edge

Bevel Amount

In the Bevel Edge figure the red edge is beveled.
Four different ways of measuring the bevel amount are shown by colored lines on the diagram on the left:

Offset (magenta line) - the distance each line moves on its respective face

Width (green line) - the distance across the newly made face between the two lines

Depth (yellow line) - the distance from the old edge to the new face, along the line that bisects the angle between the two faces

Percentage (percentage of pink line) - the distance moved along each face is a percentage of the edge along which the new edge endpoint slides

Note that the Percentage case is strange in several ways: because it moves each of the 4 endpoints of the two new lines a distance that depends on the length of an adjacent edge, the resulting new edges may not be parallel, and the distance that the two lines move along their respective faces may not be the same. Some users have asked for this as an option, however, in order to create special modeling effects (there's an airplane wing modeling tutorial that uses it, for example).

When beveling a single edge, the first three methods are equivalent in that they all produce the same final models, just with different actual values of bevel amount. For a given final result, the amounts are related to each other by simple trigonometric functions of the angle between the faces. The current Bevel in Blender uses the Offset method for measuring bevel amount, by default.

However when beveling more than one edge at once, on a model that isn't completely symmetrical and with 90 degree angles, the different measurement methods yield different results.

Bevel Offset Methods

In the Bevel Offset Methods figure, you can see the effect of using the same offset but different offset methods on faces that meet in obtuse or acute angles. The same colors as in the Bevel Edge figure show where the specified distance (the same here for all methods except percentage) is measured. For the percentage case, the size and shape of the offset depends on how long the face is in the direction away from the beveled edge.

From the above diagram, it looks like the Width method gives more even looking results.
However, when you bevel all edges around a face, and all around a vertex, the Offset method has advantages.

Bevel Offset Methods on Face

Notice that in Bevel Offset Methods On Face figure, the Offset method results in a new inside face that has a similar shape to the original face. This is because it insets each edge the same amount. Whereas the Width method has a different shape from the original face because each edge moves in a different amount, depending on the angle between the top face and the face on the other side of the beveled edge.

Probably artists would like the option of either Offset or Width. The argument for Depth is more theoretical: it corresponds to the physical act of planing of the edge and going to the same planing depth on each edge. That might make sense to some artists, though in practice it doesn't seem very useful. The Percentage method, as mentioned above, seems to have specialized uses that some artists want sometimes.

Non-manifold Edges

What to do about beveling edges that are non-manifold - that is, they are not connected to exactly two faces? For example, an edge that is the outer boundary of a plane, or an edge where three or more planes meet in a kind of fan. While one can imagine things to do in these cases, for now the decision in Blender is to regard Edge Bevel as only working on manifold edges. Non-manifold edges will be ignored when beveling.

Beveled Edge Ends

It is not always clear how to finish off the ends of beveled edges.

One case where it is pretty clear is when the vertex at the end of the edge is on an outside boundary, so the only two faces involved are the faces adjacent to the edge.

Bevel Ends - ends on boundary

At a vertex where three edges meet, with three faces between them, and only one edge is beveled, most people would expect the result on the left in the next figure. But the result on the right is another option.

Bevel Ends - one of three edges beveled at each end

In the left result, the faces between the non-beveled edges are carved away by the end of the beveled edge.
In the right result, those faces are instead split by the ends of the beveled edge.
Let us call the left method Vertex Deletion and the right method Vertex Preservation, where this refers to whether the vertex at the end of
the original beveled edge is removed from the final result or not.
One could make these two choices an option, but in situations like this Vertex Perservation doesn't seem very useful,
so for now Blender will always use Vertex Deletion in this one-of-three-edges-beveled case.

A vertex where only two edges meet and only one of the edges is beveled presents another problem.

Bevel Ends - one of two edges beveled

In the figure above, (a) is the original, showing a split edge in a plane, and (b) - (e) show four choices of how to end the bevel.
Of these choices, only (b) and (e) make a new bevel face that is a quad, something that seems desirable.
Also, choices (c), (d), and (e) will make non-planar ngons if the upper and lower faces in (a) are not coplanar.
Blender currently always uses choice (b) for situations like this.

Even with choice (b) there is another subchoice to be made: what size angle should we make at that triangle tip
(or, equivalently, how far along the unbeveled edge should we slide to get the position of the new end of that edge)?
The current choice is to slide it by the offset (the Offset distance, see above).

In the nonplanar case, this may not be exactly what is desired:

Bevel Ends - one of two edges beveled, nonplanar

An artist might like the option to specify the slide distance as zero in this case, to make a nice vertical transition between
the beveled and unbeveled edges.

Now consider the case of more than edges into a vertex, all beveled.

Bevel Ends - six of six edges beveled

In this figure, where all six side edges of an irregular six-sided cone are beveled, it is pretty clear what to do.
On the left side you see the original mesh, with dotted yellow lines showing the construction of the bevel.
Whichever way the bevel amount is specified by the user, that can be translated into the Offset amount for each edge.
Even if two adjacent bevel edges have different offsets, as they will in general when anything but the Width method has been used,
there is always an intersection point on the common face.
The old vertex gets replaced by a six-sided face that is formed by joining those intersection points.
One problem with this is that that six-sided face may not be planar - and is not, in this case:

Detail of previous bevel

The only way to avoid this would be to move the intersection points to a common plane.
Doing so would distort the desired widths, and in a case like this, would lead to uneven widths.
Failing that, the polygon could be triangulated if it is nonplanar.
One interesting way to triangulate would be to put a vertex at the center and join all the others to it.
Of course, in many common cases, only three edges meet at a vertex and the triangle that replaces the vertex cannot be non-planar.
Another case that offers no difficulties is when all the beveled edges are coplanar.

The harder case is when only some of the edges into a vertex are beveled.
Consider this case, where in figure (a), the orange edges are beveled and the black ones are not.

Bevel Ends - Not all Edges Beveled

In figure (a), the yellow dotted lines show where the beveled edge sides should go in order to make all bevels have the same width.
As you can see, they do not meet at a common point on the unbeveled edge between two beveled edges.
Since we need to keep that non-beveled edge, what can we do? Some choices:

Allow Jogs: Let the two beveled edges on either side meet the non-beveled edges at different points. In the diagram (a), this would mean making the magenta and aqua vertices and joining them. Very ugly, and won't consider this option further.

Preserve Widths: Keep the specified widths correct by changing the angle of the non-beveled edge. In the diagram (a), this would mean making the yellow vertex as the meeting point. The result is diagram (b). This could really distort the gross outline of the model but is not so bad if all of these edges are coplanar.

Preserve Angles: Keep the angles of the non-beveled edges as they are by changing the widths of one of the bevels. In the diagram (a), this would mean choosing either the aqua point or the magenta point as the meeting point (or even some other point, such as the average of those two). The result in diagram (c) is when the aqua point is chosen. Note that the bevel amount for the horizontal beveled edge has now been adjusted and is different on each side of it: the top half is as it was before but the bottom half is now a lot narrower. One could also think of a variant where both sides are changed to the same value, but then that would force a cascade of changes all around the vertex, and would lead to further inconsistencies, so let us not pursue this variant.

Blender used to use the Preserve Widths method.
But this lead to unexpected behavior, reported as bugs by some artists.
The issue is that you might expect that when there are unbeveled edges involved,
they stay as they were and the new bevel edges slide along them (like in a loop slide operation):

Bevel showing (b) Preserve Widths and (c) Preserve Angles

You can see in the above that using Preserve Widths, shown in (b), looks pretty strange,
even though the widths of the angled bevels (as measured perpendicularly across them)
are even.
Using Preserve Angles, shown in (c), is more likely what artists would like to see.

Another problem with Preserve Widths is that it is sometimes impossible to make a meeting point
for the offset edges that gives the desired width on both edges.
In the Bevel Ends - Not all Edges Beveled diagram above, the yellow point that is the
intersection of the two offset lines looks like a point that can be used to preserve widths.
However, that point only exists if, as in that diagram, everything is in the same plane.
If the two beveled edges are in different planes, then the lines, in general will not meet:

Bevel, Preserve Widths, on Bent Plane

In the above diagram, the orange lines show the needed offset lines to preserve widths.
Those lines have been artificially extended in the three images below, shown from different viewpoints,
to show that they do not meet at an intersection point.
In such a case, one has to choose what to do, and falling back on Preserve Angles is probably the
best of a number of not-great choices.

So now the method used by Blender is Preserve Angles, always, though the hook remains
in the code to use Preserve Widths (only useful in planar cases), should we decide to offer
that as a user option.

Using Preserve Angles has problems of its own:
as soon as you start changing widths, you would like to change the widths at the other end of
the edge to match, else you get an uneven effect like this:

Uneven Bevel If Change Widths

In the above diagram, the mostly-vertical beveled edge has had its width changed at the center vertex
in order to preserve angles, but the width at the bottom of that edge remains what the user specified.
Clearly we should change the width at the bottom here to match the change made at the center.
But notice that this means we cannot just process each vertex in isolation (as we are able to do
if we use Preserve Widths, which makes no width changes).
We could try to just change the width at the other end whenever we find we have to adjust a width,
but what if that vertex has already been processed? E.g., in the above diagram, suppose we process
all the vertices on the boundary before the middle one -- then we get the picture shown.
We either need a more intelligent order of processing vertices, or we need a different algorithm
that does a global width consistency pass before finally processing all the vertices.

Another question that arises when only some edges into a vertex are beveled is illustrated
by the following figure:

Bevel Choices: One Of Six Edges

In the above figure, (a) shows a single edge out of six beveled.
Whenever there are at least two unbeveled edges between beveled ones entering a vertex,
we have a choice between these two options:

Full Vertex Bevel: as shown in (b), make new vertices sliding along the unbeveled edges so as to make a vertex polygon to finish off the beveled edge(s)

Partial Vertex Bevel: as shown in (c), and from another angle in detail in (d), leave the original vertex intact and only affect faces and edges adjacent to beveled edges. Sometimes the reflex face in (d) can look strange.

When you use Full Vertex Bevel, there is the question of how much to slide the new vertices along unbeveled edges.
Blender currently uses the amount that the previously encountered beveled edge slid its vertex, which makes for an even-looking polygon.

In general, there may be more than one group of beveled edges that need to get a vertex polygon,
with each group separated by two or more unbeveled edges. E.g.:

Bevel Choices: Three Of Ten Edges, in Two Groups

For simplicity of code, and because the reflex face faces can look strange, Blender
currently always uses the Full Vertex Bevel method. It is possible that artists
might like Partial Vertex Bevel as an option.

Multiple Segments

So far the discussion has been about bevels where each edge to be beveled turns into a single new face.
Blender also allows for bevels with multiple segments, where a segment is a new face.
For example, here is an edge that has been beveled into 3 segments:

Beveled Edge With 3 Segments

Profile

The multiple segments should have parallel edges.
Viewed end-on, the multisegment beveled edge has a profile, which is typically formed
by evenly spacing the segments+1 vertices needed for the profile along a smooth curve:

Profile of Beveled Edge With 3 Segments

What smooth curve should Blender use?
In the above diagram, where there is a 90 degree angle between the faces adjacent to
the beveled edge and the offsets along those faces are equal, the natural curve
used is the quarter-circle arc with radius equal to the offset.

Artists have asked for more options for the profile, for instance to have concave profiles
or have a profile that follows the original faces (but with multiple segments).
One way to do this would be to use the parameterized family of Superellipses.
Another option would be to allow an arbitrary Curve object to be specified as the profile.
The 2.70 release of Blender implements the superellipse idea.
There is now a profile parameter, which is really just 1/4 of the superellipse exponent:

Effect of Profile Parameter

Yet another thing to consider is the effect of running the subdivision algorithm on the result of a multisegment bevel.
If you want the limit surface (as the number of subdivisions goes to infinity) to be, say, a circular arc, then the vertices
should not be placed directly on a circular arc, but rather in different, further out, positions that could be calculated mathematically.

When the adjacent faces do not form a 90 degree angle, or the offsets along each face are different,
then we need to distort the quarter-circle arc somehow.

Transform Non-Square Profile to Square Profile

The way to do this is shown in the above diagram.
On the left, the edge to be beveled is at C, viewed end-on, and the desired offset distances along the adjacent faces (in red)
are offset1 and offset2. Letter A is placed so as to make ABCD a parallelogram.
Now we can make a transform matrix that maps ABCD onto A'B'C'D', a square with sides of length 1.
Finally, we can place vertices evenly on the quarter-circle arc (in blue) in the unit square,
and then use the inverse transform matrix to map them back onto the
distorted quarter-circle arc in the original figure on the left.

Multiple Segment Corners

When all edges into a vertex are beveled into multiple segments, we need a way of making a vertex mesh
which matches the profiles coming in and looks smooth.

Here are some choices for how to do that at a three-edge vertex.

Patterns for 3-Edge, 4-segment Vertex Mesh

For reference, let us give some names to the patterns in the above diagram:

Poly - shown in (a), make a polygon out of the ends of the beveled edges; in general this polygon will be non-planar

LatLng - shown in (b), the pattern is like the latitude and longitude lines on a globe

None - shown in (c), there is no separate mesh for the vertex but rather the ends of the segments are joined up with their adjacent edge segments; this requires cutting the ends of the beveled edges in two halves

Adj - shown in (d), the pattern shares quads between adjacent edges of the vertex polygon in a way that will be explained more below.

The Poly pattern is the one currently used for one-segment bevels. It is a good pattern if the beveled edges
are all in one plane, so that the vertex polygon will be planar, but it does not look good in the non-planar, multi-segment case.
Artists probably would prefer that the mesh around the vertex somehow continues the silhouette of the profiles of
the beveled edges, when viewed from different angles.
For a quarter-circle arc profile, this would lead to the natural choice of a one-eighth-sphere octant for the case where
three beveled edges meet at 90 degree lines (as in the corner of a cube).
Patterns (b) and (d) attempt to do that, putting a number of points on the sphere octant and joining them up into smaller polygons.
The LatLng pattern is easy to compute, but it has several disadvantages: it is asymmetrical, requiring the choice of
one vertex as the pole; it is not all quads, pinching into triangles at the pole; and it is hard to see how to generalize
to cases with more than 3 beveled edges into a vertex.
The Adj pattern has none of these disadvantages, and so Blender currently uses that in preference to LatLng.
The None pattern gives a kind of sharp point that may be desired in special cases, but Blender currently
only uses it in some cases where only two edges into a vertex are beveled (the code calls this a "weld").

As you can see, the pattern is different for odd and even nseg (number of segments).
For odd nsegs, there is a center polygon with the same number of sides as the number of beveled edges.
For even nsegs, there is a single vertex in the center.
The rule to create the pattern is:

Form the boundary of the vertex mesh as in the Poly pattern (more on this below); the end of each beveled edge will have nseg segments and thus nseg+1 vertices across it (including the two ends). Call the group of vertices and edges that are on the profile of one beveled edge a boundary arc.

Let ns2 be nseg /2, rounding down if nseg is odd. Then each boundary arc has ns2 rings of quads going from the previous boundary arc to the succeeding boundary arc. The tops and bottoms of the rings go from the last half of vertices of the preceding boundary arc to the first half of vertices of the succeeding boundary arc. The sides of the ring quads are the tops and bottoms for the ring quads of the preceding and succeeding boundary arcs.

If nseg is odd, form a center polygon from the tops of the center quads in the innermost rings.

The above rules give the topology of the Adj pattern, but not how to set the coordinates of the vertices of the mesh.
I have not been able to find a formula for placing the vertices on an octant of a sphere;
it would be great if someone could find one.
The current method used by Blender is described in the implementation section below.
It can almost certainly be improved upon.

To make the boundary arcs for the Poly pattern and the outer rim of the Adj pattern,
see the following diagram:

Construction of Boundary Arcs

In the diagram, the bevel edge is part of the red line, and the offset lines are shown in dotted teal.
Points A and B have been set as described in the Bevel Edge Ends section for the case of one segment.
We need to find the position of point C.
The diagram shows a natural choice, which is the point C on the beveled edge that is closest to line AB -
that is, it makes a right angle with each of lines AB and CD.
This forms a triangle ABC. Note that in general the plane containing ABC may be tilted, and not
the plane that has the beveled edge as normal.
Now that we have triangle ABC, we can form the parallelogram that extends it and transform
to a unit square, as shown in the diagram Transform Non-Square Profile to Square Profile above.
The quarter-circle arc can be divided into nseg equal pieces and used to place the internal vertices
of the boundary arc in the unit square, and then transformed back into 3D space in their original positions.

The plane ABC is called the "profile" plane because it is the plane in which the multisegment profile
(cross section end) is inscribed.

There are cases where the above method for defining the profile plane and center arc point C does
not yield what artists expect. One example is when the edges to be beveled are as shown in (A) here:

(A) Original (B) Tilted profiles (C) Non-tilted profiles

The image in (B) shows the result of using the above algorithm to set the profile planes.
Notice that the profile at the bottom lifts off the ground plane slightly,
and the profile at the mid-level is also tilted out of the plane that contains all of
the unbeveled horizontal edges at mid-level.
The current Blender code recognizes these as special cases --- that is, when there
are two unbeveled edges that either cross or terminate the beveled edge ---
and adjusts the bevel plane to that of those unbeveled edges, as shown in (C).

Edge Bevel: Implementation

Data Structure for Vertex at End of Beveled Edge

Consider a vertex that has a mixture of beveled and unbeveled edges leading into it.

Vertex With Edges a, b, and e Beveled

On the left, the vertex in the middle is a vertex involved in a bevel,
and is represented by a BevVert struct, which contains
an array of EdgeHalf structs.
Each EdgeHalf represents, as its name implies, a half of an edge --
the half that is attached to this BevVert.
So in the example shown, the array has EdgeHalfs for edges a, b, c, d, e, and f.
In general, there is no intrinsic order to edges around a vertex:
think about a bunch of wire edges sticking out at all angles in space from a vertex.
But a very common case is that if you look down from the vertex normal,
all of the edges are connected to exactly two faces and you can find a counterclockwise
ordering of edges that is defined by going to the next edge by going across a common face.
The EdgeHalf array for a BevVert is arranged in that order
(starting from some arbitrary edge).
In the above diagram, if the vertex normal is coming out of the page, then a, b, c, d, e, f
would be a possible ordering in the array.
If there are gaps between some edges (no faces at all), the code tries to do the best it can
to order the subsequences that can be so ordered.

What happens after you bevel edges a, b, and e is that we create new vertices 0, 1, 2, and 3
as shown on the right.
These new vertices are called Boundary Vertices and are represented by a BoundVert struct.
The BoundVerts are linked together in a circular list, in counterclockwise order
(from the vertex normal side, the same as the edge order).
More than one EdgeHalf can be attached to the same BoundVert,
and an EdgeHalf that is for a beveled edge will be attached to two BoundVerts.
The following diagram shows some of the pointers that link the EdgeHalfs and
BoundVerts.

BoundVerts and EdgeHalfs for previous diagram

An EdgeHalf has a left side and a right side,
which are on the left and right sides of the edge when you look along the edge towards the BevVert
(from the vertex normal side).
The vleft and vright pointers of an EdgeHalf
point to the BoundVerts on its left and right side.
When the edge isn't beveled, they will both point to the same BoundVert,
but when it is beveled there will be two different BoundVerts.
For example, in the above diagram, beveled EdgeHalf b has
vleft = 0 and vright = 1.
Unbeveled edge d has vleft = vright = 2.
EdgeHalfs also have fprev and fnext
pointers which point to the faces between the edge and the previous and next edges, respectively,
in counterclockwise order. Those pointers may be NULL if there is no face between the edges.

A BoundVert has EdgeHalf pointers efirst and
elast, which point to the first and last EdgeHalfs in counterclockwise
order attached to the BoundVert.
If it is attached to the left side of a beveled edge then there will also be a ebev pointer that
points to that EdgeHalf.
For example, in the above diagram, BoundVert 2 has efirst = d,
elast = ebev = 3.

The vertex mesh for a BevVert is kept in a VMesh struct.
The main component of that is an array of NewVert structs, each of which
holds a 3D coordinate and, eventually, the BMesh vertex instantiated for it.
The NewVerts are indexed in the array according to the mesh pattern.
For the Adj pattern, they are indexed by (i, j, k) where i is the index of the BoundVert,
j is the index of the ring, and k is the index of the segment along a boundary arc between the
given BoundVert and its counterclockwise successor.
The vertices in the pattern have two or more indices according to this scheme.
Of course there should only be one BMVert allocated for a vertex,
so some of the NewVerts will share a BMesh vertex.

VMesh Adj Pattern indexing (even nseg)

VMesh Adj Pattern indexing (odd nseg)

The above two diagrams show the indexing patterns for an even and odd number of segments,
and three vertices in order, A, B, and C.
In the diagrams, it is assumed that the predecessor boundary arc to A and the successor boundary arc to C
are not for beveled edges.
When a boundary arc is not for a beveled edge then only segment indices 0 and nseg are used.

Update: Now even boundary arcs for non-beveled edges are divided into nseg segments.
This yields a generally more pleasing vertex mesh, and is also what is needed for multisegment vertex bevel.
The downside is that now the adjacent old faces will be higher-degree ngons when rebuild; but they should be flat so maybe this isn't such a big downside.
Leaving the above diagrams there in case some future developer wants to try going back to the mixed mode of boundary arcs.

Calculating Vertex Positions for Adj Vertex Mesh Pattern

The original method for calculating the vertex positions for the Adj vertex mesh pattern
was:

Calculate the position of the boundary arc vertices using the triangle formed by the two end vertices and the closest point on the beveled edge to the line between them, as described above in the Multiple Segment Corners section.

For the inner rings of a given boundary arc, use the same method: find the endpoints on the adjacent boundary arcs and then form the triangle with those two points and the closest point on the beveled edge, and use the map-to-unit-square method to map a quarter-circle arc into position.

Now there will be several choices for each vertex, because they can be indexed in more than one way (see the previous section), due to being on a ring for more than one boundary arc. Solve this dilemma by averaging the coordinate positions of all aliased indices for a vertex.

While this produced results that looked OK in many cases, there were also a lot of cases where it didn't work well. For instance, for some corner angles, as the number of segments increased, the contour lines looked uneven and could even cross.
Also, this method did not seem to converge to a solution where the vertex mesh smoothly transitioned into the adjacent sides (in other words, the normals of the planes at the edges of the vertex mesh didn't converge to the normals of the adjacent sides).

For these reasons, the literature was consulted for a better solution.
Here are papers about two similar approaches:

Both use the 'Adj' pattern, and use subdivision schemes with modification in order to match the boundary curves and the tangencies there.
The Hwang/Chuang paper uses a quadratric subdivision (like Doo-Sabin), while the Levin paper uses cubic subdivision (like Catmull-Clark).

After implementing both, I ended up choosing to use the Levin method.

A key part of the Levin method is to do a Catmull-Clark subdivision step.
He recommends a variant that uses weights calculated by Sabin.
As a concrete example of our indexing scheme, here is a picture of indexing of the unit cube with BoundVerts at (1,0,0), (0,1,0), and (0,0,1) and quarter-circle arcs on the faces of orthogonal planes through the origin:

Cube corner with one step of subdivision

The mesh on the left in the above diagram shows how the initial construction works:
each boundary arc is divided in two, and one central vertex is made to join to the middle vertices.
The coordinates of the middle vertex are up to ``the designer.
We move it along the direction from the original vertex to the centroid of the boundary vertices.
The amount to move depends on the profile parameter.

For the general formulas, here is one step of a an internal quad turning into four quads with one step:

General interior step of subdivision

Let M be the mesh on the left and M’ be the mesh on the right.
M' is formed by adding new edges-vertices (crosses) and a new face-vertex (circle).
Also, the coordinates of the old vertices(vertex-vertices) are changed.
The Catmull-Clark rules give the following. Assume we make all the face vertices first, then all the edge ones, then all the vertex ones.

New face vertex M’(2j+1, 2k+1):

f(j,k)

= avg(M(j,k) + M(j,k+1) + M(j+1,k) + M(j+1,k+1))

Note that f(j,k) is stored in M’(2j+1, 2k+1), which will be used below.

There is also the new center vertex which is like the New vertex vertex formula but needs to gather new edge vertices and new face vertices all around it, and combine it with α, β, γ that depend on the number of sides of the hole.

Finally, the border and corner vertices need special adjustment.
I implemented the formulas in the paper and they didn't seem to work.
They seemed too extreme and not convergent to the desired boundary curves.
Furthermore, reading Levin's thesis, which describes the same process, it appears the formulas
are different.
I probably am misunderstanding, but it is also possible there is an error in one or the other.

Instead, what I did is resample the boundary curves at even parameter spacing and set them to that.

Now there are two problems with the above method that still need to be resolved:

As you iterate, you always double the number of segments. So you will only get powers of 2 for the number of segments. We need to have an arbitrary number of segments.

The spacing between the lines does not come out even (except where forced by resampling at the boundary).

Two solve these problems, the adj_vmesh() subdivides until the power of 2 exceeds the desired number of segments, and then uses interp_vmesh() to sample the required number of segments at even intervals along two dimensions.

Even Sampling of Superellipses

The introduction of profile control via the superellipse parameter brought a new problem to the implementation: how to sample the profile so that you get evenly spaced chords.

With a circular arc for the profile, it is easy, using the trigonometric parameterization of the unit circle:

P(Θ) = (cos Θ, sin Θ)

gives the P = (x, y) coordinates for points on the arc as a function of the angle Θ.
This has the nice property that if you divide a quarter circle angle into n equal pieces, and feed the cumulative angles into the above formula, then the arc lengths (and chord lengths) will all be equal.
This is what an artist would expect along a boundary curve.

Unfortunately, the same doesn't hold for the trigonometric parameterization of a superellipse:

P(Θ) = (cos2/r Θ, sin2/r Θ)

when sampled at even intervals of Θ, gives an uneven spacing of arc length or chord length, when r is not equal to 2.
The following picture shows this:

Superellipses with Trigonometric Parameterization

On the left is a superellipse with parameter r=0.3, and on the right is one with parameter r=3.0. The dots show the points along the arc using even spacing of Θ. For parameters < 2, the dots bunch together near the ends; for parameters > 2, the dots bunch together in the middle.

There are two other ways of parameterizing superellipses. The angle-center parameterization is:

P(Θ) = (R cos Θ, R sin Θ) where R = (cosr Θ + sinr Θ)-1/r

And there's the explicit equation for y in terms of x:

y = (1 - xr)1/r

I'll spare you the graphs, but they show that neither of these parameterizations yield an even spacing of the dots for an even spacing of the parameter.
In fact, a search of the literature failed to find a closed form parameterization of a superellipse that would give an even spacing.

Ideally, we would like, for a given nseg and superellipse parameter r, to find a way of getting the ith (x, y) point along the quarter arc in the first quadrant, where i varies from Θ to nseg, such that the chord lengths between successive points are all equal.

The code does this by precomputing a series of Θ values such that evaluating at each of them in turn yields equal-length chords. (Actually, it precalculates a series of a parameter u that is a scaled version of Θ that avoids dealing with π so much.)
It currently implements this as follows:

The function superellipse_co(u, r, r_co) evaluates the trigonometric parameterization of the superellipse with parameter r at u = Θ/4 (so as u goes from 0 to 2, the arc traverses the first quadrant).

The function find_superellipse_chord_u(u0, d2goal, r) finds a value of u > u0 so that the square of the chord length from the point at u0 to the point at u equals d2goal. It works by binary search on suitable u, until the chord length squareds match to suitable tolerance.

The function find_even_superellipse_params(n, r, r_params) finds and returns an array of u parameters to give n even chords in the first quadrant. It does this by binary searching for the chord length squared that makes when repeated n times goes from the beginning to the end of the quadrant. The find_superellipse_chord_u function is used as a helper in this binary search.

For a given invocation of Bevel, which has a particular nseg, we need to precalculate the superellipse u's for that nseg.
Also, because we need to do superellipse sampling during the Catmull-Clark subdivision steps, we also need to precalculate the u's for all n that are powers of 2 up to a value that equals or exceeds nseg.
While this may not be exactly right, we just do this for the highest power of 2 needed, and then when we need the intervals for a power of 2 less than that, we suitably subsample the values for the highest power of two.

Consistent Widths for Even Bevels

To show the problem of choosing consistent edge widths, consider the following diagram:

Mesh With Widths Needing Matching

In the diagram above, the black dots and lines are vertices and edges of a mesh,
and the dotted red lines show where we would like to move the offset edges if beveling all of the edges
of the two larger triangles but not beveling the edges joining to the center dots.
Notice that some of the red lines happen to meet on the spoke edges, or at least come close,
but in general they do not, and if we want to preserve angles, we will have to adjust some bevel widths.
The letters in the diagrams label widths (between a black line and a dotted red line).
A beveled edge might be forced to have four different widths: on each side and each end.
For example, e, f, g, and h widths might all have to be different in the general case.
(There are also widths on the other sides of the outer beveled edges too, not shown in this diagram.)

We can make a graph showing which widths are dependent on other widths:

Width Dependencies for Mesh With Widths Needing Matching

In the above diagram, the red lines connect widths that are part of the same bevel edge side.
The magenta lines connect widths that are part of the same bevel edge end (on opposite sides of the original edge).
And the blue dotted lines connect widths that are on beveled edges separated by one unbeveled edge.
We would like the widths connected by red lines to be the same, so that you don't get the tapered-edge effect on beveled edges.
We would like the widths connected by magenta lines to be the same so that the new bevel face evenly straddles the original edge.
We are forced to make the edges connected by dotted blue lines to be related according to a function of the angles involved
(really: so that the meet point is somewhere sliding along the unbeveled intermediate edge).
You can see that because there are cycles in the graph, it may not be possible to satisfy all the constraints;
there may be situations in which we will have to accept some tapered beveled edges.
But by using the graph to adjust the edges in a certain order, we can minimize the number of places where that needs to happen,
and in many cases eliminate the need for that to happen.
Note: the condition that widths connected by magenta lines should be equal is probably not very important.
What matters more is that the final new bevel face looks even - the memory of how it straddled the original face is probably gone.
An alternate desired constraint might be that the sum of the widths connected by magenta lines should be as close as possible to
the width originally specified by the user.

What can we do when there are conflicts? Consider this diagram:

Two beveled edges with conflicting widths

The desired offsets for the two beveled edges are w1 and w2, which would like us to put the vertex at A and B respectively. We are going to have to pick some C at or between A and B at which to compromise; this will lead to new offsets w1’ and w2’. We would like to minimize either |w1 - w1’| + |w2 - w2’|, or (w1 - w1’)2 + (w2 - w2’)2. Let δ = the distance between A and B. We want to choose a λ between 0 and 1 and place C a distance λδ from A along the line AB. So to minimize the sum of absolute values, we minimize:

The diagram also illustrates that there are going to be some special cases. If both θs are zero (so all three lines are the same), then can only have all widths be zero. If either θ is straight or reflex then there is no way for any offset line to meet the unbeveled edge; in that case we need to use the ‘offset-in-two-planes’ method to set coordinate C (and that method might fail because the lines don’t meet, and if so, we can’t fall back to the ‘on-edge’ method! need to just choose something like the average of the closes points on each line).