No part of the Journal may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior written permission from Summary Shifted linear systems are of the form (A − σ k I)x k = b, (1) where A ∈ C N ×N , b ∈ C N and {σ k } Nσ k=1 ∈ C is a sequence of numbers, called shifts. In order to so...

background: motorcycle accidents and sport accidents lead yearly to many head injuries like head fractures and concussion. so finding the most proper helmet for reducing the injuries to head can be very helpful for head protection in such cases. methods: after 3d modeling of the helmet and head and meshing the model, a compressive impacting load of 1.31mpa was exerted on head and the model was ...

The Lanczos method is an iterative procedure to compute an orthogonal basis for the Krylov subspace generated by a symmetric matrix A and a starting vector v. An interesting application of this method is the computation of the matrix exponential exp(−τA)v. This vector plays an important role in the solution of parabolic equations where A results from some form of discretization of an elliptic o...

To improve the computational efficiencies of real-space orbital-free density functional theory, this work develops a new single-grid solver by directly providing closed-form solution to inner iteration and using an improved bisection method accelerate line search process in outer iteration, extended multi-grid solver. Numerical examples show that while maintaining high level accuracy, proposed ...

Linearly constrained optimization problems with simple bounds are considered in the present work. First, a preconditioned spectral gradient method is defined for the case in which no simple bounds are present. This algorithm can be viewed as a quasiNewton method in which the approximate Hessians satisfy a weak secant equation. The spectral choice of steplength is embedded into the Hessian appro...

When rails are laid in a railway track, small gaps are left between the ends of adjoining rails to allow for the expansion of the rails when the temperature rises. Rail joints are provided to form a continuous rail path which are normally laid in standard lengths bolted together by fishplates. When two rail joints are exactly opposing each other, they are called square (side by side) joints...

In this paper, a new multi-parameter iterative algorithm is proposed to address the PageRank problem based on multi-splitting iteration method. The method solves two linear subsystems at each by splitting coefficient matrix, considering therefore inner and outer find approximate solutions of these subsystems. It can be shown that sequence generated finally converges vector when parameters satis...

For the solution of large sparse nonnegative constrained least squares (NNLS) problems, a new iterative method is proposed which uses the CGLS method for the inner iterations and the modulus iterative method for the outer iterations to solve the linear complementarity problem resulting from the Karush-Kuhn-Tucker condition of the NNLS problem. Theoretical convergence analysis including the opti...

In this paper, we have proposed a numerical method for singularly perturbed  fourth order ordinary differential equations of convection-diffusion type. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and  finite difference method. In order to get a numerical solution for the derivative of the solution, the given interval is divided  in...