A delta-convex (d.c.) function is one which can be written as the difference of two convex functions.

The space of d.c. functions includes all C^{2} functions, and is interesting because it allows many notions from differential geometry to be generalized to objects with "corners". For example, the "tangent cone" of directional derivatives exists at every point, and is actually a tangent plane except on a set of measure 0 (the corners). This means that the gradient exists except on a set of measure 0, and indeed so does the Hessian.

In general, there is no canonical way to decompose a d.c. function into a convex and concave part. However, there are some notable exceptions (from here on, "decompose" means "decompose in a canonical way"):

f : ℝ → ℝ is d.c. iff it is the integral of a function of locally bounded variation. Since such functions can be decomposed into monotone nondecreasing and nonincreasing parts, f itself can be decomposed into the integrals of said parts.

In the particular case where f : ℝ → ℝ is C

^{2}, we simply decompose f'' into non-negative and non-positive parts, and then integrate twice.A (homogeneous) quadratic function f : ℝ

^{n}→ ℝ is completely characterized by its Hessian; f is convex (resp. concave) iff its Hessian is positive (resp. negative) semi-definite. The inner product on ℝ^{n}decomposes any symmetric linear map into positive and negative semi-definite parts, so this gives a decomposition of f into convex and concave parts.

The last two examples suggest a common theme for C^{2} functions: if we can decompose the Hessian into positive and negative semi-definite parts, then we can decompose the function itself into convex and concave parts.

The decomposition of symmetric linear maps is associated with an inner product on the space of said maps, defined as

(A,B) = tr(AB).

Call the set of positive semi-definite maps P (so the negative semi-definite maps are -P). Then the inner product has the property that (A,B) ≥ 0 for any A, B ∈ P with equality iff all nonzero eigenspaces of A are orthogonal to those of B.

So, saying that any symmetric map C can be decomposed into positive and negative semi-definite parts is equivalent to saying that C can be written uniquely as A - B where A, B ∈ P are orthogonal. We can find A and B by noting that A + B ∈ (P + C) ∩ (P - C), a convex set, and it has minimal norm among elements of this set. (Technically we only need this norm, and not an inner product, but the properties of the inner product feel very desirable.)

If f = g - h where g, h are convex C^{2} functions whose Hessians are pointwise orthogonal, then g and h are certainly unique; and it would be nice to think of this decomposition as canonical. One way to make this rigorous might be to define an inner product on C^{2} functions, and find g + h by analogy with the above.

My initial instinct was to integrate the pointwise inner product of the Hessians. However, this generally doesn't give a finite value, even in one dimension. Since d.c. functions can be characterized as having "gradients of locally bounded variation (LBV)" - in the sense that the restriction of such a function to any line has LBV derivative - I thought of limiting the inner product to those functions whose gradients were globally BV. Such functions look a lot like tangent cones, from a great enough distance, and feel like they should be important in their own right. However, even for these functions the inner product defined above is finite only on finite-radius tubular neighbourhoods of curves whose derivatives are themselves BV.

Now I'm looking at an even more specific space of functions, those which are actually (everywhere) equal to their tangent cone at the origin. I'm convinced that these functions are key to understanding the functions of globally BV gradient, which in turn might eventually lead to insight into general d.c. functions. But so far, I've not got anywhere.

My question, then: does anyone else have any insight into the problem of canonically decomposing d.c. functions? Maybe even know some results? Am I wasting my time completely - should I move on to trying to generalize some basic results about manifolds and the calculus of variations, which was my original goal?

Any input will be greatly appreciated.