Coadjoint Formalism : Nonorthogonal Basis Problems

Using nonorthogonal bases in spectral methods demands considerable effort, because applying the Gram-Schmidt process is a fundamental condition for calculations. However, operational matrices numerical methods are being used in an increasing way and extensions for nonorthogonal bases appear, requiring simplified procedures. Here, extending previous work, an efficient tensorial method is presented, in order to simplify the calculations related to the use of nonorthogonal bases in spectral numerical problems. The method is called coadjoint formalism and is based on bracket Dirac’s formulation of quantum mechanics. Some examples are presented, showing how simple it is to use the method.


Introduction
Spectral numerical methods are increasingly used to solve differential equations, even when fractional derivatives appear.Concerning these methods, they frequently use orthogonal functions [1], in order to ease the calculations and to preserve the basis elements when dimensional expansions are necessary.This is an important advantage when operational matrices are used, because if it is necessary to extend an orthogonal basis of functions from dimensions  to  + 1, only the new element must be calculated with all the former elements preserved [2,3].If the basis functions are not orthogonal, this assumption is not true.
However, there are some situations that, even for numerical methods with operational matrices [4], solution of differential equations demands new basis of functions [5,6], with some of them being nonorthogonal as, for instance, Bernstein's polynomials, mainly in the cases where boundary or initial conditions must be considered [7][8][9].Some iterative methods, based on Krylov formulation, for instance, present the use of such bases and are implemented by using the Gram-Schmidt process [8,10,11] that, in most cases, is difficult to be operationalized.
Trying to simplify this procedure, an alternative operational method is presented with the first ideas developed in [12] that is described by using tensorial language and called coadjoint formalism [13], allowing the direct use of nonorthogonal bases without any kind of previous conditioning process.
In a nutshell, coadjoint formalism adapts Dirac's bracket notation [14] to spectral methods considering finite dimension complex vectorial spaces.The methodology is applied to nonorthogonal bases in finite dimension function spaces with a generalized tensorial approach [15][16][17] that simplifies the operational conditions.
In the next section, the theoretical fundamentals of the coadjoint formalism are presented, trying to connect it with the well-known Dirac bracket, followed by a section where two examples show how simple it is to apply the method developed here, with a conclusion section closing the work.

Coadjoint Formalism: Theoretical Foundations
This section presents the concepts and definitions used to build the coadjoint formalism.The bases to be considered are finite sets of complex functions of a real variable; that is, the bases elements {  },  = 1, .(i) A finite dimensional Hilbert space H.
(ii) The finite dimension space with dimension , generated by the series expansion with  terms for the considered function, denoted by O  and called order space.
It is assumed that kets are represented by column vectors and vice versa.It is the same for bras, represented by line vectors, understood as covectors, that is, dual space elements.Consequently, an equivalence relation can be established [14]: The kets can belong either to the function space, | ⟩ ∞ ∈ H or to the order space, | ⟩ 0  ∈ O  .The same is valid for bras.
Covariant and contravariant distinction are assumed [15], and the traditional notation of differential geometry and tensorial calculus is used: Quantities described by kets are in the original space; quantities described by bras are in the dual space, in the traditional way of linear algebra [14].Primitive space is related to the space where the quantity is defined; that is, if a bra defines a quantity, the dual space assumes the primitive space condition of the quantity.
Given a covariant basis composed of |  ⟩ ∈ V, a vector  belonging to the order space, O  (V), can be described by a ket in two different ways: , with [] = {|  ⟩}.The spaces V where the spectral methods are generally applied are H, C and C  .Given a reference basis [] from V  , for a basis composed of |  ⟩ ∈ V  and described in terms of the basis [], where V can be either the Hilbert space H or any other complex vectorial space, two representations are possible: ], with the inferior dot representing the covariant nature of the basis and M  represents the set of the -order square matrices.
In this kind of disposition, the bases are described by their coordinates in the reference basis [], but the reference basis can be omitted and the invariant representation is assumed.Consequently, there are two isomorphic representations, given by ] .
Given a vector |⟩ with its components expressed in a generical basis |  ⟩, the following operations are defined: (iv) Duality:

Inner Products. The inner products (IP
are defined in several different ways, depending on the space and the representation, as it is shown below.

Order Space
(1) Hermitian product is as follows: (2) Dual product is as follows: (i) coordinate representation: Functions Space.Considering the space of the functions with domain Hybrid Spaces.The mixed inner product ⟨ | ⟩ : (3)

Duality Relations. It is possible to obtain the dual basis
{|  ⟩} from an original basis {|  ⟩} by using the duality relations, given by It is possible to express the same vector in several ways, considering the formalism described here, considering four distinct bases.Here, these representations are called connected representations and are shown in Table 1.
Considering that the covariant representation in the original space is considered to be natural, the dot below the Table 2: Transformation relations between bases.

Basis Notation Transformation relations Original
basis symbol can be omitted and the transformation relations between the bases are shown in Table 2.
The covariant and contravariant components are obtained by applying (1) contravariant:

Series Expansions of Functions. Considering a basis
{|  ⟩ ∞ } from O  (H) and elements   : R → H, the infinite series expansion of a function  is given by with   being the expansion coefficients in that basis.Considering the mixed inner product definition, this expression can be modified as with the covector ⟨| ∈ O ⋆  (C) called coefficient covector.In the subspace O  (H), projectors can be defined by with the second one called adjoint projector.
The eigenprojectors in the O  (H) space are defined as the projectors in the -dimensional proper space; that is, P = P  and P † = P †  .In a way analogous to that followed in differential geometry, fundamental metric tensor of a basis, {|  ⟩ ∞ } from O  (H), is defined as with its  eciprocal metric tensor given by Consequently, the IPs between two functions ,  represented by series are given by with f and ĝ, respectively, representing the  and  expansion coefficients, in the basis {|⟩}.
The metric tensorial fundamental operator and its reciprocal are defined as Therefore, meaning that the reciprocal tensor matrix is the inverse of the metrical tensor matrix.Calling the first one  ⋆ : 2.6.Finite Expansions.The -order finite expansion of a function  can be expressed by the action of the reciprocal tensorial metric operator over the function; that is, Under these conditions, the coefficient covector ⟨| ∈ O ⋆  (C) can be expressed as Then, the expansion is given by the mixed IP: ] being the vector associated with the second representation chosen basis.Equation (18) can be expressed by using matrices, giving a useful computational expression, in order to find the coefficient covector.
Considering the coefficient covector it is possible to write As a consequence, the expanded function becomes

Application Examples
In this section, two examples are developed showing the expansion of functions using nonorthogonal bases, one using the canonical basis and the other using a basis of nonorthogonal complex functions.

Canonical Basis.
Here, the more common case of nonorthogonality is developed, considering the canonical basis of functions defined in the real interval, defined by the polynomials: The metric tensor can be calculated as Therefore, the analytical expression for the metric tensor elements in the domain  1 can be calculated, resulting in with Following the calculations for  = 5, the metric matrix is given by allowing, by using expression (17), the calculation of the matrix representing its reciprocal tensor: By using the methodology developed in the former section and applying the transformation relation (13), the reciprocal basis is composed of polynomials that, for  = 5, are listed below: ( It is important to notice that, if the basis is orthogonal, to obtain its reciprocal is simply scale changing.Considering a nonorthogonal basis, the procedure is described here and changing order implies a whole recalculation of the basis elements. In order to illustrate the ideas in a particular case, the expansion of the function () = 2 cos 1 − 2 2 cos  is considered.
The development for  = 5 gives According to ( 22), the matrix expression is or In these expressions, the coefficient covector is generating the expansion for the function Figure 1(a) shows the almost perfect superposition of the function () and its expansion  5 ().Figure 1(b) shows the local error for the expansion.

Defining the basis
the metric tensor is expressed by

Mathematical Problems in Engineering
For  = 5 and, for the reciprocal metrics, Expanding the given function for  = 5, where (42) Figure 2 shows the superposition of the function () and its expansion  5 ().It can be noticed that the five-term approximation is not good.From the explicit expression (26), the metrical matrix for  = 10 can be determined: with the reciprocal matrix given by [ 194.628   Consequently, the ten-term expansion is given by where   = ∑ where  ⋆ [] is given by the new matrix  or, analogously, calculating the coefficients by the equation Consequently, the new approximation is given by

Conclusion
Considering that using nonorthogonal bases is being increased in spectral numerical methods requiring Gram-Schmidt procedures that are generally difficult from the operational point of view, this paper presents a simpler method for expanding functions, based on bracket formalism and called coadjoint method.
The mathematical ideas of the coadjoint method were presented and the examples have shown its practicability.
It can be added that coadjoint method is an efficient and concise tool for nonorthogonal bases with low computational costs.As it is not necessary to have orthogonal bases, more general function can be used in the numerical methods, increasing the quality of the whole process.

𝑓 10 (Figure 3
Figure 3(a) shows the almost perfect superposition of the function () and the real part of its expansion  10 ().Figure 3(b) shows the local rest for the expansion.

Figure 3 (
Figure 3(a) shows the almost perfect superposition of the function () and the real part of its expansion  10 ().Figure 3(b) shows the local rest for the expansion.