Source author record

Pedro M. Lima

Pedro M. Lima appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Biological Physics cond-mat.dis-nn math.DS nlin.PS Numerical Analysis

Catalog footprint

What is connected

5works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Numerical investigation of a neural field model including dendritic processing

We consider a simple neural field model in which the state variable is dendritic voltage, and in which somas form a continuous one-dimensional layer. This neural field model with dendritic processing is formulated as an integro-differential equation. We introduce a computational method for approximating solutions to this nonlocal model, and use it to perform numerical simulations for neuro-biologically realistic choices of anatomical connectivity and nonlinear firing rate function. For the time discretisation we adopt an Implicit-Explicit (IMEX) scheme; the space discretisation is based on a finite-difference scheme to approximate the diffusion term and uses the trapezoidal rule to approximate integrals describing the nonlocal interactions in the model. We prove that the scheme is of first-order in time and second order in space, and can be efficiently implemented if the factorisation of a small, banded matrix is precomputed. By way of validation we compare the outputs of a numerical realisation to theoretical predictions for the onset of a Turing pattern, and to the speed and shape of a travelling front for a specific choice of Heaviside firing rate. We find that theory and numerical simulations are in excellent agreement.

preprint2020arXiv

Numerical solution of a class of third-kind Volterra integral equations using Jacobi wavelets

We propose a spectral collocation method, based on the generalized Jacobi wavelets along with the Gauss-Jacobi quadrature formula, for solving a class of third-kind Volterra integral equations. To do this, the interval of integration is first transformed into the interval [-1,1], by considering a suitable change of variable. Then, by introducing special Jacobi parameters, the integral part is approximated using the Gauss-Jacobi quadrature rule. An approximation of the unknown function is considered in terms of Jacobi wavelets functions with unknown coefficients, which must be determined. By substituting this approximation into the equation, and collocating the resulting equation at a set of collocation points, a system of linear algebraic equations is obtained. Then, we suggest a method to determine the number of basis functions necessary to attain a certain precision. Finally, some examples are included to illustrate the applicability, efficiency, and accuracy of the new scheme.

preprint2015arXiv

Numerical Simulations in Two-Dimensional Neural Fields

In the present paper we are concerned with a numerical algorithm for the approximation of the two-dimensional neural field equation with delay. We consider three numerical examples that have been analysed before by other authors and are directly connected with real world applications. The main purposes are 1) to test the performance of the mentioned algorithm, by comparing the numerical results with those obtained by other authors; 2) to analyse with more detail the properties of the solutions and take conclusions about their physical meaning.

preprint2015arXiv

Numerical Solution of the Neural Field Equation in the Two-dimensional Case

We are concerned with the numerical solution of a class integro-differential equations, known as Neural Field Equations, which describe the large-scale dynamics of spatially structured networks of neurons. These equations have many applications in Neuroscience and Robotics. We describe a numerical method for the approximation of solutions in the two-dimensional case, including a space-dependent delay in the integrand function. Compared with known algorithms for this type of equation we propose a scheme with higher accuracy in the time discretisation. Since computational efficiency is a key issue in this type of calculations, we use a new method for reducing the complexity of the algorithm. The convergence issues are discussed in detail and a number of numerical examples is presented, which illustrate the performance of the method.

preprint2014arXiv

Root Finding by High Order Iterative Methods Based on Quadratures

We discuss a recursive family of iterative methods for the numerical approximation of roots of nonlinear functions in one variable. These methods are based on Newton-Cotes closed quadrature rules. We prove that when a quadrature rule with $n+1$ nodes is used the resulting iterative method has convergence order at least $n+2$, starting with the case $n=0$ (which corresponds to the Newton's method).