Enrico Bertolazzi

Enrico Bertolazzi, Marco Frego

An effective solution to the problem of Hermite $G^1$ interpolation with a clothoid curve is provided. At the beginning the problem is naturally formulated as a system of nonlinear equations with multiple solutions that is generally difficult to solve numerically. All the solutions of this nonlinear system are reduced to the computation of the zeros of a single nonlinear equation. A simple strategy, together with the use of a good and simple guess function, permits to solve the single nonlinear equation with a few iterations of the Newton--Raphson method. The computation of the clothoid curve requires the computation of Fresnel and Fresnel related integrals. Such integrals need asymptotic expansions near critical values to avoid loss of precision. This is necessary when, for example, the solution of interpolation problem is close to a straight line or an arc of circle. Moreover, some special recurrences are deduced for the efficient computation of asymptotic expansion. The reduction of the problem to a single nonlinear function in one variable and the use of asymptotic expansions make the solution algorithm fast and robust.

Enrico Bertolazzi, Marco Frego

The basic module for the solution of the minimum time optimal control of a car-like vehicle is herein presented. The vehicle is subject to the effect of laminar (linear) and aerodynamic (quadratic) drag, taking into account the asymmetric bounded longitudinal accelerations. This module is studied and designed to be fast and robust in sight to be the fundamental building block of a more extended optimal control problem that considers a given clothoid as the trajectory and the presence of a constraint on the lateral acceleration of the vehicle. The nonlinear dynamics and the different possible boundary conditions yield different analytical solutions of the differential equations, hence they by themselves a particular attention. The study of the numeric stability of the computation for limit values of the parameters is essential as showed in the numerical tests.

Enrico Bertolazzi, Marco Frego

A new robust algorithm for the numerical computation of biarcs, i.e. $G^1$ curves composed of two arcs of circle, is presented. Many algorithms exist but are based on geometric constructions, which must consider many geometrical configurations. The proposed algorithm uses an algebraic construction which is reduced to the solution of a single $2$ by $2$ linear system. Singular angles configurations are treated smoothly by using the pseudoinverse matrix when solving the linear system. The proposed algorithm is compared with the Matlab's routine \texttt{rscvn} that solves geometrically the same problem. Numerical experiments show that Matlab's routine sometimes fails near singular configurations and does not select the correct solution for large angles, whereas the proposed algorithm always returns the correct solution. The proposed solution smoothly depends on the geometrical parameters so that it can be easily included in more complex algorithms like splines of biarcs or least squares data fitting.

Ana Alonso Rodríguez, Enrico Bertolazzi, Alberto Valli

New variational formulations are devised for the curl--div system, and the corresponding finite element approximations are shown to converge. Curl--free and divergence--free finite elements are employed for discretizing the problem.

Enrico Bertolazzi, Marco Frego

A new polynomial preconditioner for symmetric complex linear systems based on Hermitian and skew-Hermitian splitting (HSS) for complex symmetric linear systems is herein presented. It applies to Conjugate Orthogonal Conjugate Gradient (COCG) or Conjugate Orthogonal Conjugate Residual (COCR) iterative solvers and does not require any estimation of the spectrum of the coefficient matrix. An upper bound of the condition number of the preconditioned linear system is provided. Moreover, to reduce the computational cost, an inexact variant based on incomplete Cholesky decomposition or orthogonal polynomials is proposed. Numerical results show that the present preconditioner and its inexact variant are efficient and robust solvers for this class of linear systems. A stability analysis of the method completes the description of the preconditioner.