BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Nick Vannieuwenhoven

Nick Vannieuwenhoven contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.AGmath.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2020arXiv

All secant varieties of the Chow variety are nondefective for cubics and quaternary forms

The Chow rank of a form is the length of its smallest decomposition into a sum of products of linear forms. For a generic form, this corresponds to finding the smallest secant variety of the Chow variety which fills the ambient space. We determine the Chow rank of generic cubics and quaternary forms by proving nondefectivity of all involved secant varieties. The main new ingredient in our proof is the generalization of a technique by [Brambilla and Ottaviani, On the Alexander--Hirschowitz theorem, J. Pure Appl. Algebra, 2008] that consists of employing Terracini's lemma and Newton's backward difference formula to compute the dimensions of secant varieties of arbitrary projective varieties. Via this inductive construction, the proof of nondefectivity ultimately reduces to proving a number of base cases. These are settled via a computer-assisted proof because of the large dimensions of the spaces involved. The largest base case required in our proof consisted of computing the dimension of a vector space constructed from the $400$th secant variety of a degree-$82$ Chow variety embedded in $\mathbb{P}^{98769}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The condition number of Riemannian approximation problems

We consider the local sensitivity of least-squares formulations of inverse problems. The sets of inputs and outputs of these problems are assumed to have the structures of Riemannian manifolds. The problems we consider include the approximation problem of finding the nearest point on a Riemannian embedded submanifold from a given point in the ambient space. We characterize the first-order sensitivity, i.e., condition number, of local minimizers and critical points to arbitrary perturbations of the input of the least-squares problem. This condition number involves the Weingarten map of the input manifold, which measures the amount by which the input manifold curves in its ambient space. We validate our main results through experiments with the $n$-camera triangulation problem in computer vision.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the average condition number of tensor rank decompositions

We compute the expected value of powers of the geometric condition number of random tensor rank decompositions. It is shown in particular that the expected value of the condition number of $n_1\times n_2 \times 2$ tensors with a random rank-$r$ decomposition, given by factor matrices with independent and identically distributed standard normal entries, is infinite. This entails that it is expected and probable that such a rank-$r$ decomposition is sensitive to perturbations of the tensor. Moreover, it provides concrete further evidence that tensor decomposition can be a challenging problem, also from the numerical point of view. On the other hand, we provide strong theoretical and empirical evidence that tensors of size $n_1~\times~n_2~\times~n_3$ with all $n_1,n_2,n_3 \ge 3$ have a finite average condition number. This suggests there exists a gap in the expected sensitivity of tensors between those of format $n_1\times n_2 \times 2$ and other order-3 tensors. For establishing these results, we show that a natural weighted distance from a tensor rank decomposition to the locus of ill-posed decompositions with an infinite geometric condition number is bounded from below by the inverse of this condition number. That is, we prove one inequality towards a so-called condition number theorem for the tensor rank decomposition.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A Riemannian Trust Region Method for the Canonical Tensor Rank Approximation Problem

The canonical tensor rank approximation problem (TAP) consists of approximating a real-valued tensor by one of low canonical rank, which is a challenging non-linear, non-convex, constrained optimization problem, where the constraint set forms a non-smooth semi-algebraic set. We introduce a Riemannian Gauss-Newton method with trust region for solving small-scale, dense TAPs. The novelty of our approach is threefold. First, we parametrize the constraint set as the Cartesian product of Segre manifolds, hereby formulating the TAP as a Riemannian optimization problem, and we argue why this parametrization is among the theoretically best possible. Second, an original ST-HOSVD-based retraction operator is proposed. Third, we introduce a hot restart mechanism that efficiently detects when the optimization process is tending to an ill-conditioned tensor rank decomposition and which often yields a quick escape path from such spurious decompositions. Numerical experiments show improvements of up to three orders of magnitude in terms of the expected time to compute a successful solution over existing state-of-the-art methods.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Pencil-based algorithms for tensor rank decomposition are not stable

We prove the existence of an open set of $n_1\times n_2 \times n_3$ tensors of rank $r$ on which a popular and efficient class of algorithms for computing tensor rank decompositions based on a reduction to a linear matrix pencil, typically followed by a generalized eigendecomposition, is arbitrarily numerically forward unstable. Our analysis shows that this problem is caused by the fact that the condition number of the tensor rank decomposition can be much larger for $n_1 \times n_2 \times 2$ tensors than for the $n_1\times n_2 \times n_3$ input tensor. Moreover, we present a lower bound for the limiting distribution of the condition number of random tensor rank decompositions of third-order tensors. The numerical experiments illustrate that for random tensor rank decompositions one should anticipate a loss of precision of a few digits.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The condition number of join decompositions

The join set of a finite collection of smooth embedded submanifolds of a mutual vector space is defined as their Minkowski sum. Join decompositions generalize some ubiquitous decompositions in multilinear algebra, namely tensor rank, Waring, partially symmetric rank and block term decompositions. This paper examines the numerical sensitivity of join decompositions to perturbations; specifically, we consider the condition number for general join decompositions. It is characterized as a distance to a set of ill-posed points in a supplementary product of Grassmannians. We prove that this condition number can be computed efficiently as the smallest singular value of an auxiliary matrix. For some special join sets, we characterized the behavior of sequences in the join set converging to the latter's boundary points. Finally, we specialize our discussion to the tensor rank and Waring decompositions and provide several numerical experiments confirming the key results.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Identifiability beyond Kruskal's bound for symmetric tensors of degree 4

We show how methods of algebraic geometry can produce criteria for the identifiability of specific tensors that reach beyond the range of applicability of the celebrated Kruskal criterion. More specifically, we deal with the symmetric identifiability of symmetric tensors in Sym$^4(\mathbb{C}^{n+1})$, i.e., quartic hypersurfaces in a projective space $\mathbb{P}^n$, that have a decomposition in 2n+1 summands of rank 1. This is the first case where the reshaped Kruskal criterion no longer applies. We present an effective algorithm, based on efficient linear algebra computations, that checks if the given decomposition is minimal and unique. The criterion is based on the application of advanced geometric tools, like Castelnuovo's lemma for the existence of rational normal curves passing through a finite set of points, and the Cayley-Bacharach condition on the postulation of finite sets. In order to apply these tools to our situation, we prove a reformulation of these results, hereby extending classical results such as Castelnuovo's lemma and the analysis of Geramita, Kreuzer, and Robbiano, "Cayley-Bacharach schemes and their canonical modules", Trans. Amer. Math. Soc. 339:443-452, 1993.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A condition number for the tensor rank decomposition

The tensor rank decomposition problem consists of recovering the unique set of parameters representing a robustly identifiable low-rank tensor when the coordinate representation of the tensor is presented as input. A condition number for this problem measuring the sensitivity of the parameters to an infinitesimal change to the tensor is introduced and analyzed. It is demonstrated that the absolute condition number coincides with the inverse of the least singular value of Terracini's matrix. Several basic properties of this condition number are investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effective criteria for specific identifiability of tensors and forms

In applications where the tensor rank decomposition arises, one often relies on its identifiability properties for interpreting the individual rank-$1$ terms appearing in the decomposition. Several criteria for identifiability have been proposed in the literature, however few results exist on how frequently they are satisfied. We propose to call a criterion effective if it is satisfied on a dense, open subset of the smallest semi-algebraic set enclosing the set of rank-$r$ tensors. We analyze the effectiveness of Kruskal's criterion when it is combined with reshaping. It is proved that this criterion is effective for both real and complex tensors in its entire range of applicability, which is usually much smaller than the smallest typical rank. Our proof explains when reshaping-based algorithms for computing tensor rank decompositions may be expected to recover the decomposition. Specializing the analysis to symmetric tensors or forms reveals that the reshaped Kruskal criterion may even be effective up to the smallest typical rank for some third, fourth and sixth order symmetric tensors of small dimension as well as for binary forms of degree at least three. We extended this result to $4 \times 4 \times 4 \times 4$ symmetric tensors by analyzing the Hilbert function, resulting in a criterion for symmetric identifiability that is effective up to symmetric rank $8$, which is optimal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On generic identifiability of symmetric tensors of subgeneric rank

We prove that the general symmetric tensor in $S^d {\mathbb C}^{n+1}$ of rank r is identifiable, provided that r is smaller than the generic rank. That is, its Waring decomposition as a sum of r powers of linear forms is unique. Only three exceptional cases arise, all of which were known in the literature. Our original contribution regards the case of cubics ($d=3$), while for $d\ge 4$ we rely on known results on weak defectivity by Ballico, Ciliberto, Chiantini, and Mella.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Most secant varieties of tangential varieties to Veronese varieties are nondefective

We prove a conjecture stated by Catalisano, Geramita, and Gimigliano in 2002, which claims that the secant varieties of tangential varieties to the $d$th Veronese embedding of the projective $n$-space $\mathbb{P}^n$ have the expected dimension, modulo a few well-known exceptions. As Bernardi, Catalisano, Gimigliano, and Idá demonstrated that the proof of this conjecture may be reduced to the case of cubics, i.e., $d=3$, the main contribution of this work is the resolution of this base case. The proposed proof proceeds by induction on the dimension $n$ of the projective space via a specialization argument. This reduces the proof to a large number of initial cases for the induction, which were settled using a computer-assisted proof. The individual base cases were computationally challenging problems. Indeed, the largest base case required us to deal with the tangential variety to the third Veronese embedding of $\mathbb{P}^{79}$ in $\mathbb{P}^{88559}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An algorithm for generic and low-rank specific identifiability of complex tensors

We propose a new sufficient condition for verifying whether generic rank-r complex tensors of arbitrary order admit a unique decomposition as a linear combination of rank-1 tensors. A practical algorithm is proposed for verifying this condition, with which it was established that in all spaces of dimension less than 15000, with a few known exceptions, listed in the paper, generic identifiability holds for ranks up to one less than the generic rank of the space. This is the largest possible rank value for which generic identifiability can hold, except for spaces with a perfect shape. The algorithm can also verify the identifiability of a given specific rank-r decomposition, provided that it can be shown to correspond to a nonsingular point of the r-th order secant variety. For sufficiently small rank, which nevertheless improves upon the known bounds for specific identifiability, some local equations of this variety are known, allowing us to verify this property. As a particular example of our approach, we prove the identifiability of a specific 5x5x5 tensor of rank 7, which cannot be handled by the conditions recently provided in [I. Domanov and L. De Lathauwer, On the Uniqueness of the Canonical Polyadic Decomposition of third-order tensors--Part II: Uniqueness of the overall decomposition, SIAM J. Matrix Anal. Appl. 34(3), 2013]. Finally, we also present a surprising new class of weakly-defective Segre varieties that nevertheless turns out to admit a generically unique decomposition.

0 citations0 reviewsNo code linkedOpen access