Researcher profile

Mário M. Graça

Mário M. Graça contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
6works
0followers
1topics
1close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2015arXiv

Higher-order root distillers

Recursive maps of high order of convergence $m$ (say $m=2^{10}$ or $m=2^{20}$) induce certain monotone step functions from which one can filter relevant information needed to globally separate and compute the real roots of a function on a given interval $[a,b]$. The process is here called a root distiller. A suitable root distiller has a powerful preconditioning effect enabling the computation, on the whole interval, of accurate roots of an high degree polynomial. Taking as model high-degree inexact Chebyshev polynomials and using the {\sl Mathematica} system, worked numerical examples are given detailing our distiller algorithm.

preprint2015arXiv

Maps for global separation of roots

Two simple predicates are adopted and certain real-valued piecewise continuous functions are constructed from them. This type of maps will be called quasi-step maps and aim to separate the fixed points of an iteration map in an interval. The main properties of these maps are studied. Several worked examples are given where appropriate quasi-step maps for Newton and Halley iteration maps illustrate the main features of quasi-step maps as tools for global separation of roots.

preprint2014arXiv

Recursive families of higher order iterative maps

To approximate a simple root of an equation we construct families of iterative maps of higher order of convergence. These maps are based on model functions which can be written as an inner product. The main family of maps discussed is defined recursively and is called {\it Newton-barycentric}. We illustrate the application of Newton-barycentric maps in two worked examples, one dealing with a typical least squares problem and the other showing how to locate simultaneously a great number of extrema of the Ackley's function.

preprint2012arXiv

A Simple Derivation of Newton-Cotes Formulas with Realistic Errors

In order to approximate the integral $I(f)=\int_a^b f(x) dx$, where $f$ is a sufficiently smooth function, models for quadrature rules are developed using a given {\it panel} of $n (n\geq 2)$ equally spaced points. These models arise from the undetermined coefficients method, using a Newton's basis for polynomials. Although part of the final product is algebraically equivalent to the well known closed Newton-Cotes rules, the algorithms obtained are not the classical ones. In the basic model the most simple quadrature rule $Q_n$ is adopted (the so-called left rectangle rule) and a correction $\tilde E_n$ is constructed, so that the final rule $S_n=Q_n+\tilde E_n$ is interpolatory. The correction $\tilde E_n$, depending on the divided differences of the data, might be considered a {\em realistic correction} for $Q_n$, in the sense that $\tilde E_n$ should be close to the magnitude of the true error of $Q_n$, having also the correct sign. The analysis of the theoretical error of the rule $S_n$ as well as some classical properties for divided differences suggest the inclusion of one or two new points in the given panel. When $n$ is even it is included one point and two points otherwise. In both cases this approach enables the computation of a {\em realistic error} $\bar E_{S_n}$ for the {\it extended or corrected} rule $S_n$. The respective output $(Q_n,\tilde E_n, S_n, \bar E_{S_n})$ contains reliable information on the quality of the approximations $Q_n$ and $S_n$, provided certain conditions involving ratios for the derivatives of the function $f$ are fulfilled. These simple rules are easily converted into {\it composite} ones. Numerical examples are presented showing that these quadrature rules are useful as a computational alternative to the classical Newton-Cotes formulas.

preprint2012arXiv

A unified framework for the computation of polynomial quadrature weights and errors

For the class of polynomial quadrature rules we show that conveniently chosen bases allow to compute both the weights and the theoretical error expression of a $n$-point rule via the undetermined coefficients method. As an illustration, the framework is applied to some classical rules such as Newton-Cotes, Adams-Bashforth, Adams-Moulton and Gaussian rules.

preprint2012arXiv

Quadrature as a least-squares and minimax problem

The vector of weights of an interpolatory quadrature rule with $n$ preassigned nodes is shown to be the least-squares solution $ω$ of an overdetermined linear system here called {\em the fundamental system} of the rule. It is established the relation between $ω$ and the minimax solution $\stackrel{\ast}{z}$ of the fundamental system, and shown the constancy of the $\infty$-norms of the respective residual vectors which are equal to the {\em principal moment} of the rule. Associated to $ω$ and $\stackrel{\ast}{z}$ we define several parameters, such as the angle of a rule, in order to assess the main properties of a rule or to compare distinct rules. These parameters are tested for some Newton-Cotes, Fejér, Clenshaw-Curtis and Gauss-Legendre rules.