Subspace gradient descent method for linear tensor equations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Iannacito, Martina, Piccinini, Lorenzo, Simoncini, Valeria
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866915816589492224
author Iannacito, Martina
Piccinini, Lorenzo
Simoncini, Valeria
author_facet Iannacito, Martina
Piccinini, Lorenzo
Simoncini, Valeria
contents The numerical solution of algebraic tensor equations is a largely open and challenging task. Assuming that the operator is symmetric and positive definite, we propose two new gradient-descent type methods for tensor equations that generalize the recently proposed Subspace Conjugate Gradient (SS-CG), D. Palitta et al, SIAM J. Matrix Analysis and Appl (2025). As our interest is mainly in a modest number of tensor modes, the Tucker format is used to efficiently represent low-rank tensors. Moreover, mixed-precision strategies are employed in certain subtasks to improve the memory usage, and different preconditioners are applied to enhance convergence. The potential of our strategies is illustrated by experimental results on tensor-oriented discretizations of three-dimensional partial differential equations with separable coefficients. Comparisons with the state-of-the-art Alternating Minimal Energy (AMEn) algorithm confirm the competitiveness of the proposed strategies.
format Preprint
id arxiv_https___arxiv_org_abs_2602_21974
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Subspace gradient descent method for linear tensor equations
Iannacito, Martina
Piccinini, Lorenzo
Simoncini, Valeria
Numerical Analysis
65F10, 65F99
G.1.3
The numerical solution of algebraic tensor equations is a largely open and challenging task. Assuming that the operator is symmetric and positive definite, we propose two new gradient-descent type methods for tensor equations that generalize the recently proposed Subspace Conjugate Gradient (SS-CG), D. Palitta et al, SIAM J. Matrix Analysis and Appl (2025). As our interest is mainly in a modest number of tensor modes, the Tucker format is used to efficiently represent low-rank tensors. Moreover, mixed-precision strategies are employed in certain subtasks to improve the memory usage, and different preconditioners are applied to enhance convergence. The potential of our strategies is illustrated by experimental results on tensor-oriented discretizations of three-dimensional partial differential equations with separable coefficients. Comparisons with the state-of-the-art Alternating Minimal Energy (AMEn) algorithm confirm the competitiveness of the proposed strategies.
title Subspace gradient descent method for linear tensor equations
topic Numerical Analysis
65F10, 65F99
G.1.3
url https://arxiv.org/abs/2602.21974