Introduction
Basilar artery is an unpaired medium-sized artery formed by the confluence of right and left vertebral arteries at the pontomedullary junction and extends to the pontomesencephalic junction. It forms the spine of posterior cerebral circulation which is constituted by the vertebrobasilar system and its branches. Normal morphology of the basilar artery forms an essential component of cerebral circulation. The present study aims to measure the level of formation & termination, length, diameter, and angle of formation of the basilar artery. The data presented are relevant to understanding human variations and would be a good anatomical reference for clinicians, anatomists, and medical students. Caliber, length, and angle of bifurcation of the basilar artery help in assessing the feasibility and approach for various surgical procedures and predict cerebro-vascular diseases.
Materials & methods
96 adult human brain specimens were studied. (78 male, 18 female) (Age range: 19–80 y; Average age: 47.66 y). Measurements were taken using Vernier calliper. Data was analysed using Microsoft Excel and SPSS software version 17.
Results
The basilar artery was formed by the confluence of two vertebral arteries in all specimens extending from the pontomedullary junction to the pontomesencephalic junction in 2/3rd of the cases. The left vertebral artery was found to be dominant in 62.5% specimens. The basilar artery showed an average length of 3.1 cm (demonstrating positive correlation with age), average diameter of 3.6–3.9 mm at different levels, and average angle of formation as 65.38° in males and 62.22° in females. Fetal type posterior cerebral artery was noticed in 9.4% cases. 3.1% and 6.3% cases were seen on the right and left sides respectively. Basilar artery fenestration was noted in 2 percent specimens.
Conclusion
Basilar artery morphology was studied in 96 human adult cadavers. Basilar artery formation and termination was normal in more than 2/3rd cases. Variations were noted in its origin, vessel hypoplasia, presence of fenestrations, and fetal patterns. The data obtained from this study are relevant for anatomists, medical students, interventional radiologists, and neurosurgeons.
Rajagopalan Ramaswamy, Gunaseelan Mani .
Application of fixed point result to solve integral equation in the setting of graphical Branciari $ {\aleph } $-metric spaces. AIMS Mathematics, 2024, 9(11): 32945-32961.
doi: 10.3934/math.20241576
[2]
Muhammad Sarwar, Muhammad Fawad, Muhammad Rashid, Zoran D. Mitrović, Qian-Qian Zhang, Nabil Mlaiki .
Rational interpolative contractions with applications in extended $ b $-metric spaces. AIMS Mathematics, 2024, 9(6): 14043-14061.
doi: 10.3934/math.2024683
[3]
Gunaseelan Mani, Arul Joseph Gnanaprakasam, Absar Ul Haq, Imran Abbas Baloch, Choonkil Park .
On solution of Fredholm integral equations via fuzzy $ b $-metric spaces using triangular property. AIMS Mathematics, 2022, 7(6): 11102-11118.
doi: 10.3934/math.2022620
[4]
Aftab Hussain .
Fractional convex type contraction with solution of fractional differential equation. AIMS Mathematics, 2020, 5(5): 5364-5380.
doi: 10.3934/math.2020344
[5]
Hasanen A. Hammad, Hassen Aydi, Choonkil Park .
Fixed point approach for solving a system of Volterra integral equations and Lebesgue integral concept in F$ _{\text{CM}} $-spaces. AIMS Mathematics, 2022, 7(5): 9003-9022.
doi: 10.3934/math.2022501
[6]
Nizar Souayah, Nabil Mlaiki, Salma Haque, Doaa Rizk, Amani S. Baazeem, Wasfi Shatanawi .
A new type of three dimensional metric spaces with applications to fractional differential equations. AIMS Mathematics, 2022, 7(10): 17802-17814.
doi: 10.3934/math.2022980
[7]
Zhenhua Ma, Jamshaid Ahmad, Abdullah Eqal Al-Mazrooei .
Fixed point results for generalized contractions in controlled metric spaces with applications. AIMS Mathematics, 2023, 8(1): 529-549.
doi: 10.3934/math.2023025
[8]
Saif Ur Rehman, Iqra Shamas, Shamoona Jabeen, Hassen Aydi, Manuel De La Sen .
A novel approach of multi-valued contraction results on cone metric spaces with an application. AIMS Mathematics, 2023, 8(5): 12540-12558.
doi: 10.3934/math.2023630
[9]
Mian Bahadur Zada, Muhammad Sarwar, Reny George, Zoran D. Mitrović .
Darbo-Type $ \mathcal{Z}_{\rm{m}} $ and $ \mathcal{L}_{\rm{m}} $ contractions and its applications to Caputo fractional integro-differential equations. AIMS Mathematics, 2021, 6(6): 6340-6355.
doi: 10.3934/math.2021372
[10]
Hongyan Guan, Jinze Gou, Yan Hao .
On some weak contractive mappings of integral type and fixed point results in $ b $-metric spaces. AIMS Mathematics, 2024, 9(2): 4729-4748.
doi: 10.3934/math.2024228
Abstract
Introduction
Basilar artery is an unpaired medium-sized artery formed by the confluence of right and left vertebral arteries at the pontomedullary junction and extends to the pontomesencephalic junction. It forms the spine of posterior cerebral circulation which is constituted by the vertebrobasilar system and its branches. Normal morphology of the basilar artery forms an essential component of cerebral circulation. The present study aims to measure the level of formation & termination, length, diameter, and angle of formation of the basilar artery. The data presented are relevant to understanding human variations and would be a good anatomical reference for clinicians, anatomists, and medical students. Caliber, length, and angle of bifurcation of the basilar artery help in assessing the feasibility and approach for various surgical procedures and predict cerebro-vascular diseases.
Materials & methods
96 adult human brain specimens were studied. (78 male, 18 female) (Age range: 19–80 y; Average age: 47.66 y). Measurements were taken using Vernier calliper. Data was analysed using Microsoft Excel and SPSS software version 17.
Results
The basilar artery was formed by the confluence of two vertebral arteries in all specimens extending from the pontomedullary junction to the pontomesencephalic junction in 2/3rd of the cases. The left vertebral artery was found to be dominant in 62.5% specimens. The basilar artery showed an average length of 3.1 cm (demonstrating positive correlation with age), average diameter of 3.6–3.9 mm at different levels, and average angle of formation as 65.38° in males and 62.22° in females. Fetal type posterior cerebral artery was noticed in 9.4% cases. 3.1% and 6.3% cases were seen on the right and left sides respectively. Basilar artery fenestration was noted in 2 percent specimens.
Conclusion
Basilar artery morphology was studied in 96 human adult cadavers. Basilar artery formation and termination was normal in more than 2/3rd cases. Variations were noted in its origin, vessel hypoplasia, presence of fenestrations, and fetal patterns. The data obtained from this study are relevant for anatomists, medical students, interventional radiologists, and neurosurgeons.
1.
Introduction
Consider the following generalized absolute value equations (GAVE)
Ax−B|x|=b,
(1.1)
where A,B∈Rn×n are two given matrices and b∈Rn is a given vector, and |⋅| denotes the componentwise absolute value. Especially, if B=I, where I stands for the identity matrix, then the GAVE (1.1) reduces to the standard absolute value equations (AVE)
Ax−|x|=b.
(1.2)
The GAVE (1.1) and the AVE (1.1) arise in many scientific computing and engineering problems, including the linear programming problems, the linear complementarity problems (LCP), the bimatrix games, the quadratic programming and so on, see [1,2,3,4] for more details. Taking the well-known LCP as an example: for a given matrix M∈Rn×n and a given vector q∈Rn, find two vectors z,w∈Rn such that
z≥0,w=Mz+q≥0andzTw=0.
(1.3)
Here and thereafter, (⋅)T denotes the transpose of either a vector or a matrix. By simply letting z=|x|−x and w=|x|+x, then the LCP (1.3) can be equivalently transformed into the GAVE
(M+I)x−(M−I)|x|=qwithx=12((M−I)z+q).
(1.4)
In fact, the GAVE (1.4) can also be transformed into the LCP (1.3) [5,6]. For more details of the relation between the GAVE (1.1) and the LCP (1.3), please see [7,8,9].
In recent decades, much more attention has been paid to the GAVE (1.1) and the AVE (1.2). On one hand, some sufficient conditions have been studied to guarantee the existence and uniqueness of the solution of the GAVE (1.1) [10,11,12,13,14]. Rohn showed that the GAVE (1.1) is uniquely solvable for any b∈Rn if σmin(A)>σmax(|B|)[10]. Wu and Li extended the results of [10], and exhibited that the GAVE (1.1) is uniquely solvable for any b∈Rn if σmin(A)>σmax(B)[11,12]. In [13], Rohn et al. showed that the GAVE (1.1) is uniquely solvable for any b∈Rn if ρ(|A−1B|)<1. Here and in the sequel, σmin(⋅), σmax(⋅) and ρ(⋅) denote the minimal singular value, the maximal singular value and the spectral radius of the corresponding matrix, respectively. On the other hand, many efficient iteration methods, including the concave minimization algorithm [15,16], the sign accord algorithm [17], the optimization algorithm [18], the hybrid algorithm [19], the preconditioned AOR iterative method [20], the Picard-HSS iteration method [21], the Newton-type method [22,23,24] and so on, have been studied for solving the GAVE (1.1).
Due to the existence of the nonlinear term B|x|, the GAVE (1.1) can be regarded as a system of nonlinear equations
F(x)=0withF(x)=Ax−B|x|−b.
(1.5)
As a result, the well-known Newton iteration method
xk+1=xk−F′(xk)−1F(xk),k=0,1,2,⋯,
(1.6)
can be used provided that the Jacobian matrix F′(x) of F(x) exists and is invertible. However, the Newton iteration method (1.6) can not be used directly to solve the GAVE (1.1) since F(x)=Ax−B|x|−b is non-differentiable. For the special case, i.e., B=I, by considering F(x)=Ax−|x|−b as a piece-wise linear vector function, Mangasarian in [22] used the generalized Jacobian ∂|x| of |x| based on a subetaadient of its components and presented the following generalized Newton (GN) iteration method
x(k+1)=(A−D(x(k)))−1b,k=0,1,2,⋯
(1.7)
to get an approximate solution of the AVE (1.2), where D(x(k))=∂|x|=diag(sign(x(k))) and sign(x) stands for a vector with components equal to 1, 0, or −1 depending on whether the corresponding component of x is positive, zero or negative. Theoretical analysis showed that the GN iteration method (1.7) is globally linearly convergent under certain conditions [22]. Hu et al. extended the GN iteration scheme (1.7) to solve the GAVE (1.1) and proposed a weaker convergence condition [25]. For a general matrix B, the specific GN iteration scheme is
x(k+1)=(A−BD(x(k)))−1b,k=0,1,2,⋯.
(1.8)
Recently, convergence results of the GN iteration schemes (1.7) and (1.8) have been further discussed in [26,27,28]. From the GN iteration scheme (1.7) or (1.8), we can see that the coefficient matrix A−D(xk) or A−BD(xk) is changed at each iteration step. For large problems, it is very difficult especially if the coefficient matrix is ill-conditioned. In addition, if the generalized Jacobian matrix is singular, then the GN iteration method fails. To remedy these, a lot of effective improvements have been presented, such as a stable and quadratic locally convergent iteration scheme [28], the generalized Traub's method [29], the modified GN iteration method [24,30], the inexact semi-smooth Newton iteration method [31], a new two-step iterative method [32] and so on. All these improvements greatly accelerate the convergence rate of the GN iteration method. However, when the singular generalized Jacobian matrix happens, these newly developed iteration methods fail, too.
In this paper, by introducing a relaxation iteration parameter, we propose a relaxed generalized Newton (RGN) iteration method to solve the GAVE (1.1). In fact, the RGN iteration method is a generalization of the GN iteration method [22] and the Picard iteration method [13] studied recently. The advantage of the new RGN iteration method is twofold. By introducing suitable iteration parameter, it not only can avoid singularity of the generalized Jacobian matrix, but also improves the convergence rate. Theoretically, we prove that the RGN iteration method is well defined and globally linearly convergent under certain conditions. Moreover, a specific sufficient convergence condition is presented when the coefficient matrix A is symmetric positive definite. With two numerical examples, we show that the new proposed RGN iteration method is much more efficient than some existing Newton-type iteration methods.
The rest of this paper is organized as follows. In Section 2, the RGN iteration method is introduced to solve the GAVE (1.1). Convergence analyses are studied in detail in Section 3. In Section 4, two numerical examples from the LCP (1.3) are presented to demonstrate the effectiveness of our new method. Finally, we end this paper with some conclusions and outlook in Section 5.
2.
The relaxed generalized Newton iteration method
In this section, a new relaxed generalized Newton iteration method is introduced to solve the GAVE (1.1).
Let θ≥0 be a nonnegative real parameter. Based on the Newton iteration scheme (1.6) and the ideas studied [22,30], a new iteration scheme is introduced to solve the GAVE (1.1)
F(xk)+(∂F(xk)+(1−θ)BD(xk))(xk+1−xk)=0.
(2.1)
Substituting F(xk)=Axk−B|xk|−b (1.5) and the generalized Jacobian ∂F(xk)=A−BD(xk) into (2.1), we obtain
Axk−B|xk|−b+(A−θBD(xk))(xk+1−xk)=0.
(2.2)
Since D(x)=diag(sign(x)) is a diagonal matrix and satisfies D(xk)xk=|xk|[22]. Then the new iteration scheme (2.1) or (2.2) is simplified into the following final form
(A−θBD(xk))xk+1=b+(1−θ)B|xk|,k=0,1,2,⋯.
(2.3)
Here, the iteration parameter θ can be regarded as a role of relaxation, which can avoid the singularity problems and adjust the condition number of the coefficient matrix A−θBD(xk) so as to improve the convergence rate of the GN iteration method (1.8). So, we call the new iteration method (2.3) the relaxed generalized Newton (RGN) iteration method. In particular, if θ=1, then the RGN iteration method (2.3) reduces to the GN iteration method (1.8). If θ=0, then the RGN iteration method (2.3) becomes
Axk+1=b+B|xk|,k=0,1,2,⋯,
(2.4)
which is known as the Picard iteration method [7,13].
The RGN iteration method (2.3) is well defined provided that the coefficient matrix A−θBD(xk) is nonsingular at each iteration step. The following theorem presents a sufficient condition. To this end, we first define a set of matrices
since the diagonal matrix D(x)=diag(sign(x)) may change at each iteration step.
Theorem 2.1.Let A,B∈Rn×n, θ≥0 be a nonnegative real parameter and D∈Rn×n be any matrix in the set D (2.5). Let λmin(ATA) and λmax(BTB) be the smallest eigenvalue of ATA and the largest eigenvalue of BTB, respectively. If
λmin(ATA)λmax(BTB)>θ2,
(2.6)
then A−θBD is nonsingular and the RGN iteration method (2.3) is well defined.
Proof. We argue it by contradiction. If A−θBD is singular, then there exists a nonzero vector x such that
(A−θBD)x=0.
In addition, since D∈Rn×n is a diagonal matrix with each diagonal element being 1, −1 or 0, DTD is a diagonal matrix, too, and each diagonal element is either 1 or 0. Thus, it holds
Therefore, A−θBD is nonsingular and the RGN iteration method (2.3) is well defined provided that the condition (2.6) holds.
Remark 2.1.It should be noted that the condition given in Theorem 2.1 is a theoretical generalization of some recent works. In particular, if B=I and θ=1, then the condition (2.6) becomes λmin(ATA)>1, which means that all singular values of A exceed 1 and is in good agreement with [22, Lemma 2.1]. If only θ=1, then the condition (2.6) is the one given in [25, Theorem 3.1]. In addition, if θ=0, then the condition (2.6) is equivalent to show that the matrix A is nonsingular, which clearly shows that the Picard iteration method (2.4) is well defined.
To better implement the new RGN iteration method (2.3) in actual computations, we present an algorithmic version of the RGN iteration method (2.3) as follows. Here, ‖⋅‖2 denotes the Euclidean norm of either a vector or a matrix.
Algorithm 2.1.(The RGN iteration method)
1). Choose an arbitrary initial vector x0 and a nonnegative parameter θ. Given ε and set k=0;
2). If ‖Axk−B|xk|−b‖2≤ε‖b‖2, stop;
3). Compute D(xk)=diag(sign(xk));
4). Solve the following linear system to obtain xk+1
(A−θBD(xk))xk+1=b+(1−θ)B|xk|;
5). Set k=k+1. Go to Step 2.
3.
Convergence analysis
In this section, we will establish the convergence theory of the RGN iteration method (2.3) for solving the GAVE (1.1). Specially, two general convergence conditions of the RGN iteration method (2.3) will be presented firstly. Then, a sufficient convergence condition is proposed when the system matrix A is symmetric positive definite. In addition, as the special cases of our new RGN iteration method (2.3), the convergence conditions of the GN iteration method (1.8) and the Picard iteration method (2.4) can be acquired immediately.
3.1. General sufficient convergence conditions
In this subsection, we first study some sufficient convergence conditions only when the RGN iteration method (2.3) is well defined.
Theorem 3.1.Let A,B∈Rn×n, θ≥0 be a nonnegative real parameter and satisfy the condition (2.6). Let D∈Rn×n be any matrix in the set D (2.5). If
‖(A−θBD)−1‖2‖B‖2<11+θ,
(3.1)
then the RGN iteration method (2.3) converges linearly from any starting point to a solution x∗ of the GAVE (1.1).
Proof. Let x∗ be a solution of the GAVE (1.1), then it satisfies
where the inequality ‖|xk|−|x∗|‖2≤‖xk−x∗‖2 is used. From (3.4), we can see that the RGN iteration method (2.3) converges linearly from any starting point to a solution x∗ of the GAVE (1.1) provided that the condition (3.1) is satisfied.
Theorem 3.2.Under the conditions in Theorem 3.1. Further assume that A is nonsingular. If
‖A−1‖2‖B‖2<11+2θ,
(3.5)
then the RGN iteration method (2.3) converges linearly from any starting point to a solution x∗ of the GAVE (1.1).
Proof. According to Theorem 3.1, we only need to verify the condition (3.1). Under the condition (3.5), by the Banach perturbation lemma (see [33,Lemma 2.3.3]), we have
Therefore, the RGN iteration method (2.3) converges linearly from any starting point to a solution x∗ of the GAVE (1.1) if the condition (3.5) is satisfied.
As mentioned in Section 2, the well-known GN iteration method (1.8) and the Picard iteration method (2.4) are special cases of the new RGN iteration method (2.3) when the relaxation parameter θ is 1 and 0, respectively. By simply letting θ=1 and 0, we can obtain the following two corollaries, which describe the convergence conditions of the GN iteration method (1.8) and the Picard iteration method (2.4), respectively, for solving the GAVE (1.1).
Corollary 3.1.Let A,B∈Rn×n, and D∈Rn×n be any matrix in the set D (2.5). Assume that A−BD is nonsingular. If
‖(A−BD)−1‖2‖B‖2<12,
or if A is nonsingular and
‖A−1‖2‖B‖2<13,
then the GN iteration method (1.8) converges linearly from any starting point to a solution x∗ of the GAVE (1.1).
Corollary 3.2.Let A,B∈Rn×n. Assume that A is nonsingular. If
‖A−1‖2‖B‖2<1,
then the Picard iteration method (2.4) converges linearly from any starting point to a solution x∗ of the GAVE (1.1).
3.2. The case of symmetric positive definite
In this subsection, we turn to discuss the convergence conditions of the RGN iteration method (2.3) for solving the GAVE (1.1) when the system matrix A is symmetric positive definite.
Theorem 3.3.Let A∈Rn×n be a symmetric positive definite matrix, B∈Rn×n, D∈Rn×n be any matrix in the set D (2.5), θ be a positive constant and satisfy (2.6). Further assume that μmin is the smallest eigenvalue of the matrix A and ‖B‖2=τ. If
μmin>(1+2θ)τ,
(3.6)
then the RGN iteration method (2.3) converges linearly from any starting point to a solution x∗ of the GAVE (1.1).
Proof. Since the matrix A is symmetric positive definite, it is easy to check that
‖A−1‖2‖B‖2≤τμmin.
If μmin and τ further satisfy the condition (3.6), we have
‖A−1‖2‖B‖2≤τμmin<11+2θ.
Therefore, by Theorem 3.2, we obtain that the RGN iteration method (2.3) converges linearly from any starting point to a solution x∗ of the GAVE (1.1). This completes the proof.
In Theorem 3.3, by setting θ=1 and 0, we have the following two corollaries to guarantee the convergence of the GN iteration method (1.8) and the Picard iteration method (2.4) for solving the GAVE (1.1), respectively.
Corollary 3.3.Let A∈Rn×n be a symmetric positive definite matrix and μmin be its smallest eigenvalue. Let B∈Rn×n and ‖B‖2=τ. Let D∈Rn×n be any matrix in the set D (2.5). If
μmin>3τ,
then the GN iteration method (1.8) converges linearly from any starting point to a solution x∗ of the GAVE (1.1).
Corollary 3.4.Let A∈Rn×n be a symmetric positive definite matrix and μmin be its smallest eigenvalue. Let B∈Rn×n and ‖B‖2=τ. If
μmin>τ,
then the Picard iteration method (2.4) converges linearly from any starting point to a solution x∗ of the GAVE (1.1).
4.
Numerical experiments
In this section, we present some numerical examples arising from two types of LCP (1.3) to show the effectiveness of the RGN iteration method (2.6), and demonstrate the advantages of the new RGN iteration method over the well-known Lemke's method, the Picard iteration method (2.4), the GN iteration method (1.8) and the modified GN iteration (MGN) method [30] from aspects of the iteration counts (denoted by “IT”) and the elapsed CPU times (denoted by “CPU”). The first LCP comes from [4], which have been studied for standard test problem by many researchers. The second LCP arising from the practical traffic single bottleneck models [34], which are usually solved by the well-known Lemke's method. Note that the Lemke's method is one of the most efficient direct methods for solving the LCP (1.3).
In our experiments, the initial guess vector is the zero vector. All runs are terminated once
RES(x(k)):=‖Ax(k)−B|x(k)|−b‖2‖b‖2≤10−7
or if the prescribed iteration number kmax=5000 is exceeded. At each iteration step of the RGN iteration method, the Picard iteration method, the GN iteration method and the MGN iteration method, we need to solve a system of linear equations with the coefficient matrix A−θBD, A, A−BD and A+I−BD respectively. Here, these systems of linear equations are solved by the sparse LU factorization when these coefficient matrices are nonsymmetric and by the sparse Cholesky factorization when these coefficient matrices are symmetric positive definite. To efficiently implement the RGN iteration method, we need to choose the relaxation iteration parameter θ in advance. The convergence rates of all parameter dependent iteration methods heavily depend on the particular choice of the iteration parameter. The analytic determination of the value θ which results in the fastest convergence of the RGN iteration method appears to be quite a difficult problem. Here, the relaxation iteration parameter θ used in the new RGN iteration method is chosen to be the experimentally optimal one θexp, which leads to the smallest iteration step. In the following tables, “-” means that the corresponding iteration method does not converge to the approximate solution within kmax iteration steps or even diverges. All computations are run in MATLAB (version R2014a) in double precision, Intel(R) Core(TM) (i5-3337U CPU, 8G RAM) Windows 8 system. Here, we use the Matlab codes presented in https://ww2.mathworks.cn/matlabcentral/fileexchange/41485 to test the Lemke's method.
Example 4.1.([4]) Consider the LCP (1.3), in which M∈Rn×n is given by M=ˆM+μI∈Rn×n and q∈Rn is given by q=−Mz(∗), where
is the unique solution of the LCP (1.3). Here n=m2. From the discussion presented in Section 1, we can see that the LCP (1.3) can be equivalently expressed as the GAVE (1.1), where A=M+I and B=M−I. The exact solution of the GAVE (1.1) is
x(∗)=(−0.6,−0.6,−0.6,⋯,−0.6,⋯)T∈Rn.
For the first example, we take two cases of the parameter μ, i.e., μ=0 and μ=−1, in actual computations. For each parameter μ, four increasing sizes, i.e., m=30,60,90,120, are considered. The corresponding dimensions for each problem are n=900,3600,8100,14400, respectively. Note that for the case μ=0, both the matrices M and A are symmetric positive definite, while for the case μ=−1, the matrix M is symmetric indefinite and the matrix A is symmetric positive definite. In Tables 1 and 2, we list the numerical results of different methods for μ=0 and μ=−1, respectively.
Table 1.
Numerical results for Example 4.1 with μ=0.
From Tables 1 and 2, we can see that the GN iteration method and the new RGN iteration method perform much better than the other three computational methods in terms of both iteration steps and elapsed CPU times. For the case μ=0, the Lemke's method can not converge to a satisfactory solution for large problems. Although the Picard iteration method and the MGN iteration method converge, the numerical results show that the convergence rates of these two iteration methods are very slow. We have also noticed that the iteration steps of the Picard iteration step and the MGN iteration step are almost the same, but the elapsed CPU times show that the MGN iteration method costs much more expensive than the Picard iteration method. The reason is that the coefficient matrix of the MGN iteration method is changed at each iteration step. The best choice of the relaxation iteration parameter in the RGN iteration method is θexp=1, which means that the GN iteration method is the best one. For the case μ=−1, the Lemke's method computes the exact solution for all test problems, but the iteration steps and the elapsed CPU times show that this method is not competitive in actual computations. The Picard iteration method diverges. This is because the matrix M is indefinite and the convergence conditions in Corollary 3.2 and Corollary 3.4 can not be satisfied. Among the GN iteration method, the MGN iteration method and the RGN iteration method, we can see from the numerical results that the new RGN iteration method is the best one.
Example 4.2.([4]) Consider the LCP (1.3), in which M=ˆM++μI∈Rn×n and q=−Mz(∗). Different from Example 4.1, we assume that the matrix ˆM in the second example is nonsymmetric, i.e.,
Similar to Example 4.1, the second example can also be equivalently expressed as the GAVE (1.1). Note that Example 4.2 has the same exact solution as that of Example 4.1. For the second example, we still take two parameters for μ, i.e., μ=0 and μ=−1. For each parameter μ, we consider four increasing sizes, i.e., m=30,60,90,120. Thus, the total dimensions for each problem are n=900,3600,8100,14400, respectively. Different from Example 4.1, both the matrices M and A are nonsymmetric positive definite for the case μ=0. For the case μ=−1, the matrix M is nonsymmetric indefinite and the matrix A is nonsymmetric positive definite.
Tables 3 and 4 list the corresponding numerical results of different methods for μ=0 and μ=−1, respectively. These numerical results further confirm the observations obtained from Tables 1 and 2, i.e., the GN iteration method and the new RGN iteration method are superior to other three computational methods in terms of computing efficiency. For the case μ=0, the Lemke's method converge very slow for small problems and do not converge within the given maximum iteration step for large problems. The other four computational methods are convergent. However, the Picard iteration method and the MGN iteration method converge very slow. The GN iteration method and the new RGN iteration method have the same computational results and converge very fast. That means the GN iteration method is the best one. Most important, the iteration steps of both the GN iteration method and the new proposed RGN iteration method are constant as the problem sizes grow. For the case μ=−1, the Lemke's method performs much better than itself for the case μ=0. However, the computational results show that the Lemke's method is still not competitive in real applications. The Picard iteration method is still divergent. The reason is the same as that in Example 4.1. The GN iteration method, the MGN iteration method and the proposed RGN iteration method converge faster than the Lemke's method. From these numerical results, we see again that the RGN iteration method performs best among the three Newton-based iteration methods.
Table 3.
Numerical results for Example 4.2 with μ=0.
Example 4.3.([34]) The second example comes from the single bottleneck model with both the homogeneous commuters and the heterogeneous commuters. The dynamic equilibrium conditions for the single bottleneck model can be transformed into the LCP (1.3), in which the system matrix M and the vector q have the following specific structure
M=[0M1M2−MT3−MT1S000M1I0M3000]andq=[q1s1q2q3],
where the submatrices M1∈R(ΥG)×Υ, M2∈R(ΥG)×(ΥG), M3∈RG×(ΥG), S∈RΥ×Υ are
Here, τ∈T=[0,1,2,⋯,Υ] and g∈G=[1,2,⋯,G] are the indexes for the time interval and the user group, respectively. When G=1 and G>1, then the LCP (1.3) can be used to study the homogeneous case and the heterogeneous case, respectively. s denotes the bottleneck capacity with units given by number of vehicles per time interval and Ng denotes the number of individuals in group g. αg, βg and γg are the unit costs (or value) of the travel time, arriving early to work and arriving late to work in group g, respectively. τ∗g is the preferred arrival time in group g. 1 stands for a vector of all ones.
For the second example, the total dimension is n=2ΥG+Υ+G. It is proved in [34] that the system matrix M is a copositive matrix and the LCP (1.3) has a unique solution. In [34], the authors used the Lemke's method to solve such problems and simulate the single bottleneck model for both the homogeneous case (G=1) and the heterogenous case (G=3). In addition, they took the total demand N=25, the bottleneck capacity s=3 for time unit and the preferred arrival time τ∗=7. The time duration is 10 time units. For the homogeneous case, i.e., G=1, the unit costs are taken as α=2, β=1 and γ=4 for per time unit. For the heterogeneous case, i.e., G=3, the unit costs are taken as α:β:γrations=2:1:4 and τ∗g=6,7,8 for groups 1−3, respectively. For more detailed selection of these parameters, please see [34]. Here, the corresponding LCP is further equivalently transformed into the GAVE (1.1), where A=M+I and B=M−I. Then the Picard iteration method, the GN iteration method, the MGN iteration method and the new RGN iteration method are applied to solve the GAVE.
Numerical results of different computational methods are listed in Tables 5 and 6 for G=1 and G=3, respectively. From these two tables, we can see that the Lemke's method is successfully applied to solve the single bottleneck model, but the elapsed CPU times indicate that the Lemke's method costs very expensive. The GN iteration method fails to solve the test problems. This is because the singularity of the coefficient matrix A−BD(xk) occurs in implementing the GN iteration method. The Picard iteration method and the MGN iteration method can only be applied to some small problems. For large problems, these two iteration methods fail to converge within the prescribed largest iteration number. Our new RGN iteration method can be successfully applied to solve all the test problems and costs very cheap. Therefore, the new RGN iteration method is a powerful computational method to solve the GAVE (1.1).
Table 5.
Numerical results for Example 4.3 with homogeneous case (G=1).
In this paper, by introducing a relaxation iteration parameter, a new relaxed generalized Newton (RGN) iteration method is proposed to solve the generalized absolute value equations. We have proved that the RGN iteration method is well defined and converges globally under certain conditions. Two numerical examples, both arising from the well-known LCP problem, are used to illustrate the efficiency of the new computational method. Numerical results show that the RGN iteration method converges and has much better computing efficiency than some existing methods provided that suitable relaxation iteration parameters are chosen.
Just like most of the parameter-based iteration methods, the choice of the iteration parameter is an open and a challenging problem. Here, the RGN iteration method is proved to be only linearly convergent. In some recent works, the GN iteration method has been modified to be a globally and quadratically iteration method under very strong conditions. Therefore, how to improve the RGN iteration method needs further-in-depth studies. In addition, the generalized absolute value equations with general nonlinear term, which arise in nonlinear complementarity problems [35], implicit complementarity problems [36,37], quasi complementarity problems [38,39], are of great interesting. Future work should focus on estimating the quasi-optimal value of the relaxation iteration parameter, finding globally and quadratically convergent RGN iteration method, extending to solve more applications and so on.
Acknowledgments
The author are very much indebted to An Wang for writing the Matlab codes. This work is supported by the National Natural Science Foundation of China (Nos. 11771225, 61771265, 71771127), the Humanities and Social Science Foundation of the Ministry of Education of China (No. 18YJCZH274), the Science and Technology Project of Nantong City (No. JC2018142) and the ‘226’ Talent Scientific Research Project of Nantong City.
Conflict of interest
The authors declare there is no conflict of interest.
Acknowledgments
We thank the entire teaching and nonteaching faculty of Department of Forensic Medicine, Government Medical College Thiruvananthapuram for their kind support during the period of this study.
Conflicts of interest
The authors declare there is no conflicts of interest.
References
[1]
Standring S, Ellis H, Wigley C (2005) Gray's anatomy: The anatomical basis of clinical practice Edinburgh; New York: Elsevier Churchill Livingstone.
[2]
Wankhede HA, Hosmani PB, Nimje DA (2014) Morphological study of the basilar artery in adult human cadavers. Int J Anat Res 2: 497-502.
[3]
Efendic A, Isakovic E, Delic J, et al. (2014) Vascular geometry of vertebrobasilar tree with and without aneurysm. Med Glas (Zenica) 11: 252-257.
[4]
Ravensbergen J, Hillen B, Krijger JK, et al. (1996) The influence of the angle of confluence on the flow in a Vertebrobasilar junction model. J Biomech 29: 281-299. doi: 10.1016/0021-9290(95)00064-X
[5]
Tanaka M, Sakaguchi M, Miwa K, et al. (2013) Basilar artery diameter is an independent predictor of incident cardiovascular events. Arterioscler Thromb Vasc Biol 33: 2240-2244. doi: 10.1161/ATVBAHA.113.301467
[6]
Aoki J, Iguchi Y, Kimura K, et al. (2010) Diameter of the basilar artery may be associated with neurological deterioration in acute pontine infarction. Eur Neurol 63: 221-226. doi: 10.1159/000279619
[7]
Menshawi K, Gutierrez J, Mohr JP (2015) A functional perspective on embryology and anatomy of cerebral blood supply. J Stroke 17: 144-158. doi: 10.5853/jos.2015.17.2.144
[8]
van Raamt AF, van Laar PJ, Mali WP, et al. (2006) The fetal variant of the circle of Willis and its influence on the cerebral collateral circulation. Cerebrovasc Dis 22: 217-224. doi: 10.1159/000094007
[9]
Romanes GJ (1986) Cunningham's Manual of Practical Anatomy Oxford University Press.
[10]
Padmavathi G, Niranjana Murthy KV, Rajeshwari T (2011) Study of the variations in the origin and termination of basilar artery. Anatomica Karnataka 5: 54-59.
[11]
Yaşargil MG (2009) Microneurosurgery: Principles, applications, and training. Practical Handbook of Neurosurgery Vienna: Springer.
[12]
Fisher CM, Gore I, Okabe N, et al. (1965) Atherosclerosis of the carotid and vertebral arteries—extracranial and intracranial. J Neuropathol Exp Neurol 24: 455-476. doi: 10.1097/00005072-196507000-00007
[13]
Hong JM, Chung CS, Bang OY, et al. (2009) Vertebral artery dominance contributes to basilar artery curvature and peri-vertebrobasilar junctional infarcts. J Neurol Neurosurg Psychiatry 80: 1087-1092. doi: 10.1136/jnnp.2008.169805
[14]
Mamatha H, D'souza AS, Suhani S (2012) Human cadaveric study of the morphology of the basilar artery. Singapore Med J 53: 760-763.
[15]
Iqbal S (2013) Vertebro-basilar variants and their basic clinical implications. Int J Med Res Health Sci 2: 799-808.
[16]
Wankhede HA, Hosmani PB, Nimje DA (2014) Morphological study of the basilar artery in adult human cadavers. Int J Anat Res 2: 497-502.
[17]
Pai BS, Varma RG, Kulkarni RN, et al. (2007) Microsurgical anatomy of the posterior circulation. Neurology India 55: 31-41. doi: 10.4103/0028-3886.30424
[18]
Glagov S, Zarins C, Giddens DP, et al. (1988) Hemodynamics and atherosclerosis. Insights and perspectives gained from studies of human arteries. Arch Pathol Lab Med 112: 1018-1031.
[19]
Lee SH, Hur N, Jeong SK (2012) Geometric analysis and blood flow simulation of basilar artery. J Atheroscler Thromb 19: 397-405. doi: 10.5551/jat.10694
[20]
Passero S, Filosomi G (1998) Posterior circulation infarcts in patients with vertebrobasilar dolichoectasia. Stroke 29: 653-659. doi: 10.1161/01.STR.29.3.653
[21]
Sacks JG, Lindenberg R (1970) Dolicho-ectatic intracranial arteries. Symptomatology and pathogenesis of arterial elongation and distension. Johns Hopkins Med J 125: 95-105.
[22]
Goldstein SJ, Sacks JG, Lee C, et al. (1983) Computed tomographic findings in cerebral artery ectasia. AJNR 4: 501-504.
[23]
Saeki N, Rhoton AL (1977) Anatomy of the posterior circle of Willis. J Neurosurg 46: 564-578. doi: 10.3171/jns.1977.46.5.0563
[24]
Kamath SA (1979) A study of dimensions of the basilar artery in South Indian subjects. J Anat Soc India 28: 45-64.
[25]
Pico F, Labreuche J, Gourfinkel-An I, et al. (2006) Basilar artery diameter and 5-year mortality in patients with stroke. Stroke 37: 2342-2347. doi: 10.1161/01.STR.0000236058.57880.03
[26]
Black SP, Ansbacher LE (1984) Saccular aneurysm associated with segmental duplication of the basilar artery. A morphological study. J Neurosurg 6: 1005-1008. doi: 10.3171/jns.1984.61.6.1005
[27]
Finlay HM, Canham PB (1994) The layered fabric of cerebral artery fenestrations. Stroke 25: 1799-1806. doi: 10.1161/01.STR.25.9.1799
[28]
Stopford JS (1916) The arteries of the pons and medulla oblongata. J Anat Physiol 50: 131-164.
This article has been cited by:
1.
Wan-Chen Zhao, Xin-Hui Shao,
New matrix splitting iteration method for generalized absolute value equations,
2023,
8,
2473-6988,
10558,
10.3934/math.2023536
2.
Chen-Can Zhou, Qin-Qin Shen, Geng-Chen Yang, Quan Shi,
A general modulus-based matrix splitting method for quasi-complementarity problem,
2022,
7,
2473-6988,
10994,
10.3934/math.2022614
3.
Yiming Zhang, Dongmei Yu, Yifei Yuan,
On the Alternative SOR-like Iteration Method for Solving Absolute Value Equations,
2023,
15,
2073-8994,
589,
10.3390/sym15030589
4.
Luotao Wang, Kalidoss Rajakani,
Pattern Generation and Design of Floral Patterns Based on Newton’s Iterative Algorithm,
2022,
2022,
1530-8677,
1,
10.1155/2022/2116403
5.
Xin-Hui Shao, Wan-Chen Zhao,
Relaxed modified Newton-based iteration method for generalized absolute value equations,
2023,
8,
2473-6988,
4714,
10.3934/math.2023233
6.
Hongjun Sun, Tianyu Yang, Hongbing Ding, Jinxia Li, Wenqiang Zhang,
Online Measurement of Gas and Liquid Flow Rates in Wet Gas Using Vortex Flowmeter Coupled With Conductance Ring Sensor,
2022,
71,
0018-9456,
1,
10.1109/TIM.2021.3129222
7.
Yan-Xia Dai, Ren-Yi Yan, Ai-Li Yang,
Minimum Residual BAS Iteration Method for Solving the System of Absolute Value Equations,
2024,
2096-6385,
10.1007/s42967-024-00403-z
8.
Lu-Xin Wang, Yang Cao, Qin-Qin Shen,
Two Variants of Robust Two-Step Modulus-Based Matrix Splitting Iteration Methods for Mixed-Cell-Height Circuit Legalization Problem,
2024,
2096-6385,
10.1007/s42967-024-00400-2
9.
Lin Zheng, Yangxin Tang, Luigi Rarità,
A Two‐Step Matrix‐Splitting Iterative Method for Solving the Generalized Absolute Value Equation,
2024,
2024,
2314-4629,
10.1155/2024/8396895
10.
Rashid Ali, Fuad A. Awwad, Emad A. A. Ismail,
The development of new efficient iterative methods for the solution of absolute value equations,
2024,
9,
2473-6988,
22565,
10.3934/math.20241098
11.
Chen-Can Zhou, Yang Cao, Qin-Qin Shen, Quan Shi,
A modified Newton-based matrix splitting iteration method for generalized absolute value equations,
2024,
442,
03770427,
115747,
10.1016/j.cam.2023.115747
12.
Ximing Fang, Minhai Huang,
Convergence of the modified Newton-type iteration method for the generalized absolute value equation,
2025,
2193-5343,
10.1007/s40065-025-00500-8
Figure 1. A specimen showing normal variant of basilar artery and its branches at the base of brain. It is situated on the ventral aspect of pons in a shallow median groove called basilar sulcus. Abbreviations used: VA(R): Right vertebral artery ;VA(L): Left vertebral artery; BA: basilar artery; AICA: anterior inferior cerebellar artery; LA: labyrinthine artery; SCA: superior cerebellar artery; PCA: posterior cerebral artery; P.Comm.A: posterior communicating artery; ICA(R): right internal carotid artery; ICA(L): left internal carotid artery
Figure 2. Pie chart showing age distribution: Among 96 human brain specimens, 1% belonged to below 20 age group. 11% belonged to 21–30 age group. 19% specimens were from 31–40 age group. 29 % specimens were from 41–50 age group. 25% specimens belonged to 51–60 age group. In 61–70 age group, 10% specimens were present. Between 71–80 age group, 5% specimens were present
Figure 3. Box plot diagram representing angle of formation: The lower and upper end of the whisker represents minimum and maximum angle of formation respectively. Red line represents the median angle. Lower border of yellow box represents 25th percentile and upper border of the yellow box represents 75th percentile
Figure 4. Hypoplasia of right vertebral artery: right vertebral artery is hypoplastic compared to the left vertebral artery. Duplication of superior cerebellar artery is seen on right side. Abbreviations used: VA(L): left vertebral artery; BA: basilar artery; SCA: superior cerebellar artery; PCA: posterior cerebral artery
Figure 5. Hypoplastic left vertebral artery: The right vertebral artery is dominant compared to the left vertebral artery. A proximal basilar artery fenestration was noted in this specimen. Abbreviations used: BA: basilar artery; VA(R): right vertebral artery VA(L): left vertebral artery; AICA: anterior inferior cerebellar artery; SCA: superior cerebellar artery; PCA: posterior cerebral artery
Figure 6. A communicating branch between right and left vertebral arteries: A short communicating branch was observed between right and left vertebral arteries. Abbreviations used: BA: basilar artery; VA(R): right vertebral artery; VA(L): left vertebral artery
Figure 7. Scatter plot diagram showing the relationship between age and length of the basilar artery
Figure 8. Pie chart showing variation in the level of formation of the basilar artery (Abbreviation used: P-M junction: Ponto-medullary junction)
Figure 9. Pie chart showing variation in the level of termination of the basilar artery
Figure 10. Photograph showing the fetal type of posterior cerebral artery on the left side: Posterior cerebral artery arises from the internal carotid artery (fetal PCA) on the left side. On the right side, posterior cerebral artery shows normal origin from the basilar artery. Abbreviations used: BA: basilar artery; SCA: superior cerebellar artery; PCA: posterior cerebral artery; P.Comm.artery: posterior communicating artery; ICA(R): right internal carotid artery; ICA(L): left internal carotid artery