Search

In this paper, we introduce novel fast matrix inversion algorithms that leverage triangular decomposition and recurrent formalism, incorporating Strassen's fast matrix multiplication. Our research places particular emphasis on triangular matrices, where we propose a novel computational approach based on combinatorial techniques for finding the inverse of a general non-singular triangular matrix. Unlike iterative methods, our combinatorial approach for (block) triangular-type matrices enables direct computation of the matrix inverse through a nonlinear combination of carefully selected combinatorial entries from the initial matrix. This unique characteristic makes our proposed method fully parallelizable, offering significant potential for efficient implementation on parallel computing architectures. Our approach demonstrates intriguing features that allow the derivation of recurrent relations for constructing the matrix inverse. By combining the (block) combinatorial approach, with a recursive triangular split method for inverting triangular matrices, we develop potentially competitive algorithms that strike a balance between efficiency and accuracy. We provide rigorous mathematical proofs of the newly presented method. Additionally, we conduct extensive numerical tests to showcase its applicability and efficiency. The comprehensive evaluation and experimental results presented in this paper confirm the practical utility of our proposed algorithms, demonstrating their superiority over classical approaches in terms of computational efficiency.

Derivatives of the displacement tensor with respect to the independent model parameters of the subsurface, also called Fr\'echet derivatives (or sensitivity kernels), are a key ingredient for seismic full-waveform inversion with a local-search optimization algorithm. They provide a quantitative measure of the expected changes in the seismograms due to perturbations of the subsurface model parameters for a given survey geometry. Since 2.5-D wavefield modeling involves a real point source in a 2-D geological model with 3D (spherical) wave properties, it yields synthetic data much closer to the actual practical field data than the commonly used 2-D wave simulation does, which uses an unrealistic line source in which the waves spread cylindrically. Based on our recently developed general 2.5-D wavefield modeling scheme, we apply the perturbation method to obtain explicit analytic expressions for the derivatives of the displacement tensor for 2.5-D/2-D frequency-domain seismic full-waveform inversion in general viscoelastic anisotropic media. We then demonstrate the numerical calculations of all these derivatives in two common cases: (i) viscoelastic isotropic and (ii) viscoelastic tilted transversely isotropic (TTI) solids. Examples of the differing sensitivity patterns for the various derivatives are investigated and compared for four different homogeneous models involving 2-D and 2.5-D modeling., Comment: 50 pages, 15 figures, no published before

This paper presents a new algorithm for the parallel in time (PiT) numerical simulation of time dependent partial/ordinary differential equations. We propose a reliable alternative to the well know parareal in time algorithm, by formulating the parallel in time problem algebraically and solve it using an adapted Bi-Conjugate gradient stabilized method. The proposed Parallel in time Stable Bi-Conjugate algorithm (PiTSBiCG) has a great potential in stabilizing the parallel resolution for a variety of problems. In this work, we describe the mathematical approach to the new algorithm and provide numerical evidences that show its superiority to the standard parareal method.

This paper presents a general description of a parameter estimation inverse problem for systems governed by nonlinear differential equations. The inverse problem is presented using optimal control tools with state constraints, where the minimization process is based on a first-order optimization technique such as adaptive monotony-backtracking steepest descent technique and nonlinear conjugate gradient methods satisfying strong Wolfe conditions. Global convergence theory of both methods is rigorously established where new linear convergence rates have been reported. Indeed, for the nonlinear non-convex optimization we show that under the Lipschitz-continuous condition of the gradient of the objective function we have a linear convergence rate toward a stationary point. Furthermore, nonlinear conjugate gradient method has also been shown to be linearly convergent toward stationary points where the second derivative of the objective function is bounded. The convergence analysis in this work has been established in a general nonlinear non-convex optimization under constraints framework where the considered time-dependent model could whether be a system of coupled ordinary differential equations or partial differential equations. Numerical evidence on a selection of popular nonlinear models is presented to support the theoretical results. Nonlinear Conjugate gradient methods, Nonlinear Optimal control and Convergence analysis and Dynamical systems and Parameter estimation and Inverse problem

We present a regularization strategy that leads to well-conditioned boundary integral equation formulations of Helmholtz equations with impedance boundary conditions in two-dimensional Lipschitz domains. We consider both the case of classical impedance boundary conditions, as well as the case of transmission impedance conditions wherein the impedances are certain coercive operators. The latter type of problems is instrumental in the speed up of the convergence of Domain Decomposition Methods for Helmholtz problems. Our regularized formulations use as unknowns the Dirichlet traces of the solution on the boundary of the domain. Taking advantage of the increased regularity of the unknowns in our formulations, we show through a variety of numerical results that a graded-mesh based Nystr\"om discretization of these regularized formulations leads to efficient and accurate solutions of interior and exterior Helmholtz problems with impedance boundary conditions., Comment: arXiv admin note: text overlap with arXiv:1509.04415

We present a time-parallelization method that enables to accelerate the computation of quantum optimal control algorithms. We show that this approach is approximately fully efficient when based on a gradient method as optimization solver: the computational time is approximately divided by the number of available processors. The control of spin systems, molecular orientation and Bose-Einstein condensates are used as illustrative examples to highlight the wide range of application of this numerical scheme., Comment: Journal paper submitted to Physical Review A

In this paper we present an advanced numerical method to simulate a real life challenging industrial problem that consists of the non-destructive testing in steam generators. We develop a finite element technique that handles the big data numerical set of systems arising when a discretization of the eddy-current equation in three dimensional space is made. Using a high performance technique, our method becomes fully efficient. We provide numerical simulations using the software Freefem++ which has a powerful tool to handle finite element method and parallel computing. We show that our technique speeds up the simulation with a good efficiency factor., Comment: 8 pages towcolomnw, 11 figure

This work is motivated by the monitoring of conductive clogging deposits in steam generator at the level of support plates. One would like to use monoaxial coils measurements to obtain estimates on the clogging volume. We propose a 3D shape optimization technique based on simplified parametrization of the geometry adapted to the measurement nature and resolution. The direct problem is modeled by the eddy current approximation of time-harmonic Maxwell's equations in the low frequency regime. A potential formulation is adopted in order to easily handle the complex topology of the industrial problem setting. We first characterize the shape derivatives of the deposit impedance signal using an adjoint field technique. For the inversion procedure, the direct and adjoint problems have to be solved for each coil vertical position which is excessively time and memory consuming. To overcome this difficulty, we propose and discuss a steepest descent method based on a fixed and invariant triangulation. Numerical experiments are presented to illustrate the convergence and the efficiency of the method.

We consider the inverse problem of estimating the shape profile of an unknown deposit from a set of eddy current impedance measurements. The measurements are acquired with an axial probe, which is modeled by a set of coils that generate a magnetic field inside the tube. For the direct problem, we validate the method that takes into account the tube support plates, highly conductive part, by a surface impedance condition. For the inverse problem, finite element and shape sensitivity analysis related to the eddy current problem are provided in order to determine the explicit formula of the gradient of a least square misfit functional. A geometrical-parametric shape inversion algorithm based on cylindrical coordinates is designed to improve the robustness and the quality of the reconstruction. Several numerical results are given in the experimental part. Numerical experiments on synthetic deposits, nearby or far away from the tube, with different shapes are considered in the axisymmetric configuration., Comment: 34

In this paper we present a new steepest-descent type algorithm for convex optimization problems. Our algorithm pieces the unknown into sub-blocs of unknowns and considers a partial optimization over each sub-bloc. In quadratic optimization, our method involves Newton technique to compute the step-lengths for the sub-blocs resulting descent directions. Our optimization method is fully parallel and easily implementable, we first presents it in a general linear algebra setting, then we highlight its applicability to a parabolic optimal control problem, where we consider the blocs of unknowns with respect to the time dependency of the control variable. The parallel tasks, in the last problem, turn "on" the control during a specific time-window and turn it "off" elsewhere. We show that our algorithm significantly improves the computational time compared with recognized methods. Convergence analysis of the new optimal control algorithm is provided for an arbitrary choice of partition. Numerical experiments are presented to illustrate the efficiency and the rapid convergence of the method., Comment: 28 pages

We present a parareal in time algorithm for the simulation of neutron diffusion transient model. The method is made efficient by means of a coarse solver defined with large time steps and steady control rods model. Using finite element for the space discretization, our implementation provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch-Maurer-Werner (LMW) benchmark [1].

In this paper, we present a method that enables solving in parallel the Euler-Lagrange system associated with the optimal control of a parabolic equation. Our approach is based on an iterative update of a sequence of intermediate targets that gives rise to independent sub-problems that can be solved in parallel. This method can be coupled with the parareal in time algorithm. Numerical experiments show the efficiency of our method., Comment: 14 pages