• 제목/요약/키워드: computational mathematics

검색결과 3,169건 처리시간 0.026초

ANALYSIS OF THE UPPER BOUND ON THE COMPLEXITY OF LLL ALGORITHM

  • PARK, YUNJU;PARK, JAEHYUN
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • 제20권2호
    • /
    • pp.107-121
    • /
    • 2016
  • We analyze the complexity of the LLL algorithm, invented by Lenstra, Lenstra, and $Lov{\acute{a}}sz$ as a a well-known lattice reduction (LR) algorithm which is previously known as having the complexity of $O(N^4{\log}B)$ multiplications (or, $O(N^5({\log}B)^2)$ bit operations) for a lattice basis matrix $H({\in}{\mathbb{R}}^{M{\times}N})$ where B is the maximum value among the squared norm of columns of H. This implies that the complexity of the lattice reduction algorithm depends only on the matrix size and the lattice basis norm. However, the matrix structures (i.e., the correlation among the columns) of a given lattice matrix, which is usually measured by its condition number or determinant, can affect the computational complexity of the LR algorithm. In this paper, to see how the matrix structures can affect the LLL algorithm's complexity, we derive a more tight upper bound on the complexity of LLL algorithm in terms of the condition number and determinant of a given lattice matrix. We also analyze the complexities of the LLL updating/downdating schemes using the proposed upper bound.

튜링의 기계주의에 대한 괴델의 비평 (G$\ddot{o}$del's Critique of Turings Mechanism)

  • 현우식
    • 한국수학사학회지
    • /
    • 제17권4호
    • /
    • pp.27-36
    • /
    • 2004
  • 이 논문에서는 튜링의 기계주의에 대한 괴델의 비평을 다룬다. 여기에서 튜링의 기계주의란 튜링기계의 기호배열이 인간의 마음의 각 상태에 대응된다는 것을 의미한다. 첫째 부분에서는 계산으로서의 인지과정에 대한 튜링의 분석을 검토한다. 두 번째 부분에서는 튜링기계의 개념을 살펴보고, 세 번째 부분에서는 인지적 체계로서의 튜링기계가 갖는 계산적 한계를 설명한다. 네 번째 부분에서는 괴델이 튜링의 기계주의에 동의하지 않았음을 보이고, 마지막으로 오라클 튜링기계과 그 함의에 대하여 논의한다.

  • PDF

COMPACTNESS AND DIRICHLET'S PRINCIPLE

  • Seo, Jin Keun;Zorgati, Hamdi
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • 제18권2호
    • /
    • pp.193-207
    • /
    • 2014
  • In this paper we explore the emergence of the notion of compactness within its historical beginning through rigor versus intuition modes in the treatment of Dirichlet's principle. We emphasize on the intuition in Riemann's statement on the principle criticized by Weierstrass' requirement of rigor followed by Hilbert's restatement again criticized by Hadamard, which pushed the ascension of the notion of compactness in the analysis of PDEs. A brief overview of some techniques and problems involving compactness is presented illustrating the importance of this notion. Compactness is discussed here to raise educational issues regarding rigor vs intuition in mathematical studies. The concept of compactness advanced rapidly after Weierstrass's famous criticism of Riemann's use of the Dirichlet principle. The rigor of Weierstrass contributed to establishment of the concept of compactness, but such a focus on rigor blinded mathematicians to big pictures. Fortunately, Poincar$\acute{e}$ and Hilbert defended Riemann's use of the Dirichlet principle and found a balance between rigor and intuition. There is no theorem without rigor, but we should not be a slave of rigor. Rigor (highly detailed examination with toy models) and intuition (broader view with real models) are essentially complementary to each other.

A LOCALIZED GLOBAL DEFORMATION MODEL TO TRACK MYOCARDIAL MOTION USING ECHOCARDIOGRAPHY

  • Ahn, Chi Young
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • 제18권2호
    • /
    • pp.181-192
    • /
    • 2014
  • In this paper, we propose a robust real-time myocardial border tracking algorithm for echocardiography. Commonly, after an initial contour of LV border is traced at one or two frame from the entire cardiac cycle, LV contour tracking is performed over the remaining frames. Among a variety of tracking techniques, optical flow method is the most widely used for motion estimation of moving objects. However, when echocardiography data is heavily corrupted in some local regions, the errors bring the tracking point out of the endocardial border, resulting in distorted LV contours. This shape distortion often occurs in practice since the data acquisition is affected by ultrasound artifacts, dropout or shadowing phenomena of cardiac walls. The proposed method deals with this shape distortion problem and reflects the motion realistic LV shape by applying global deformation modeled as affine transform partitively to the contour. We partition the tracking points on the contour into a few groups and determine each affine transform governing the motion of the partitioned contour points. To compute the coefficients of each affine transform, we use the least squares method with equality constraints that are given by the relationship between the coefficients and a few contour points showing good tracking results. Many real experiments show that the proposed method supports better performance than existing methods.

A TWO-DIMENSIONAL FINITE VOLUME METHOD FOR TRANSIENT SIMULATION OF TIME- AND SCALE-DEPENDENT TRANSPORT IN HETEROGENEOUS AQUIFER SYSTEMS

  • Liu, F.;Turner, I.;Ahn, V.;Su, N.
    • Journal of applied mathematics & informatics
    • /
    • 제11권1_2호
    • /
    • pp.215-241
    • /
    • 2003
  • In this paper, solute transport in heterogeneous aquifers using a modified Fokker-Planck equation (MFPE) is investigated. This newly developed mathematical model is characterised with a time-, scale-dependent dispersivity. A two-dimensional finite volume quadrilateral mesh method (FVQMM) based on a quadrilateral background interpolation mesh is developed for analysing the model. The FVQMM transforms the coupled non-linear partial differential equations into a system of differential equations, which is solved using backward differentiation formulae of order one through five in order to advance the solution in time. Three examples are presented to demonstrate the model verification and utility. Henry's classic benchmark problem is used to show that the MFPE captures significant features of transport phenomena in heterogeneous porous media including enhanced transport of salt in the upper layer due to its parameters that represent the dependence of transport processes on scale and time. The time and scale effects are investigated. Numerical results are compared with published results on the some problems.

UNITARY INTERPOLATION FOR VECTORS IN TRIDIAGONAL ALGEBRAS

  • Jo, Young-Soo
    • Journal of applied mathematics & informatics
    • /
    • 제11권1_2호
    • /
    • pp.431-436
    • /
    • 2003
  • Given vectors x and y in a Hilbert space, an interpolating operator is a bounded operator T such that Tx = y. An interpolating operator for n vectors satisfies the equation $Tx_i\;:\;y_i,\;for\;i\;=\;1,\;2,\;{\cdots},\;n$. In this article, we obtained the following : $Let\;x\;=\;\{x_i\}\;and\;y=\{y_\}$ be two vectors in a separable complex Hilbert space H such that $x_i\;\neq\;0$ for all $i\;=\;1,\;2;\cdots$. Let L be a commutative subspace lattice on H. Then the following statements are equivalent. (1) $sup\;\{\frac{\$\mid${\sum_{k=1}}^l\;\alpha_{\kappa}E_{\kappa}y\$\mid$}{\$\mid${\sum_{k=1}}^l\;\alpha_{\kappa}E_{\kappa}x\$\mid$}\;:\;l\;\in\;\mathbb{N},\;\alpha_{\kappa}\;\in\;\mathbb{C}\;and\;E_{\kappa}\;\in\;L\}\;<\;\infty\;and\;$\mid$y_n\$\mid$x_n$\mid$^{-1}\;=\;1\;for\;all\;n\;=\;1,\;2,\;\cdots$. (2) There exists an operator A in AlgL such that Ax = y, A is a unitary operator and every E in L reduces, A, where AlgL is a tridiagonal algebra.

AN EFFICIENT AND STABLE ALGORITHM FOR NUMERICAL EVALUATION OF HANKEL TRANSFORMS

  • Singh, Om P.;Singh, Vineet K.;Pandey, Rajesh K.
    • Journal of applied mathematics & informatics
    • /
    • 제28권5_6호
    • /
    • pp.1055-1071
    • /
    • 2010
  • Recently, a number of algorithms have been proposed for numerical evaluation of Hankel transforms as these transforms arise naturally in many areas of science and technology. All these algorithms depend on separating the integrand $rf(r)J_{\upsilon}(pr)$ into two components; the slowly varying component rf(r) and the rapidly oscillating component $J_{\upsilon}(pr)$. Then the slowly varying component rf(r) is expanded either into a Fourier Bessel series or various wavelet series using different orthonormal bases like Haar wavelets, rationalized Haar wavelets, linear Legendre multiwavelets, Legendre wavelets and truncating the series at an optimal level; or approximating rf(r) by a quadratic over the subinterval using the Filon quadrature philosophy. The purpose of this communication is to take a different approach and replace rapidly oscillating component $J_{\upsilon}(pr)$ in the integrand by its Bernstein series approximation, thus avoiding the complexity of evaluating integrals involving Bessel functions. This leads to a very simple efficient and stable algorithm for numerical evaluation of Hankel transform.

A MULTI-SERVER RETRIAL QUEUEING MODEL WITH POISSON SIGNALS

  • CHAKRAVARTHY, SRINIVAS R.
    • Journal of applied mathematics & informatics
    • /
    • 제39권5_6호
    • /
    • pp.601-616
    • /
    • 2021
  • Retrial queueing models have been studied extensively in the literature. These have many practical applications, especially in service sectors. However, retrial queueing models have their own limitations. Typically, analyzing such models involve level-dependent quasi-birth-and-death processes, and hence some form of a truncation or an approximate method or simulation approach is needed to study in steady-state. Secondly, in general, the customers are not served on a first-come-first-served basis. The latter is the case when a new arrival may find a free server while prior arrivals are waiting in the retrial orbit due to the servers being busy during their arrivals. In this paper, we take a different approach to the study of multi-server retrial queues in which the signals are generated in such a way to provide a reasonably fair treatment to all the customers seeking service. Further, this approach makes the study to be level-independent quasi-birth-and-death process. This approach is different from any considered in the literature. Using matrix-analytic methods we analyze MAP/M/c-type retrial queueing models along with Poisson signals in steady-state. Illustrative numerical examples including a comparison with previously published retrial queues are presented and they show marked improvements in providing a quality of service to the customers.

THE WEIBULL MARSHALL-OLKIN LOMAX DISTRIBUTION WITH APPLICATIONS TO BLADDER AND HEAD CANCER DATA

  • KUMAR, DEVENDRA;KUMAR, MANEESH;ABD EL-BAR, AHMED M.T.;LIMA, MARIA DO CARMO S.
    • Journal of applied mathematics & informatics
    • /
    • 제39권5_6호
    • /
    • pp.785-804
    • /
    • 2021
  • The proposal of new families has been worked out by many authors over recent years. Many ways to generate new families have been developed as the methods of addition, linear combination, composition and, one of the newer, the T-X family of distributions. Using this latter method, Korkmaz et al. (2018) proposed a new class called Weibull Marshall-Olkin-G (WMO-G) family. In the present work, we propose a new distribution, based on the WMO-G family, using the Lomax distribution as baseline, called Weibull Marshall-Olkin Lomax (WMOL) distribution. The hazard rate function of this distribution can be increasing, decreasing, bathtub-shaped, decreasing-increasing-decreasing and unimodal. Some properties of the proposed model are developed. Besides that, we consider method of maximum likelihood for estimating the unknown parameters of the WMOL distribution. We provide a simulation study in order to verify the asymptotic properties of the maximum likelihood estimates. The applicability of the new distribution to modeling real life data is proved by two real data sets.

APPROXIMATED SEPARATION FORMULA FOR THE HELMHOLTZ EQUATION

  • Lee, Ju-Hyun;Jeong, Nayoung;Kang, Sungkwon
    • 호남수학학술지
    • /
    • 제41권2호
    • /
    • pp.403-420
    • /
    • 2019
  • The Helmholtz equation represents acoustic or electromagnetic scattering phenomena. The Method of Lines are known to have many advantages in simulation of forward and inverse scattering problems due to the usage of angle rays and Bessel functions. However, the method does not account for the jump phenomena on obstacle boundary and the approximation includes many high order Bessel functions. The high order Bessel functions have extreme blow-up or die-out features in resonance region obstacle boundary. Therefore, in particular, when we consider shape reconstruction problems, the method is suffered from severe instabilities due to the logical confliction and the severe singularities of high order Bessel functions. In this paper, two approximation formulas for the Helmholtz equation are introduced. The formulas are new and powerful. The derivation is based on Method of Lines, Huygen's principle, boundary jump relations, Addition Formula, and the orthogonality of the trigonometric functions. The formulas reduce the approximation dimension significantly so that only lower order Bessel functions are required. They overcome the severe instability near the obstacle boundary and reduce the computational time significantly. The convergence is exponential. The formulas adopt the scattering jump phenomena on the boundary, and separate the boundary information from the measured scattered fields. Thus, the sensitivities of the scattered fields caused by the boundary changes can be analyzed easily. Several numerical experiments are performed. The results show the superiority of the proposed formulas in accuracy, efficiency, and stability.