Source author record

François Cuvelier

François Cuvelier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Mathematical Software Numerical Analysis Applications math.AP math.PR

Catalog footprint

What is connected

4works

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.

preprint2015arXiv

An efficient way to assemble finite element matrices in vector languages

Efficient Matlab codes in 2D and 3D have been proposed recently to assemble finite element matrices. In this paper we present simple, compact and efficient vectorized algorithms, which are variants of these codes, in arbitrary dimension, without the use of any lower level language. They can be easily implemented in many vector languages (e.g. Matlab, Octave, Python, Scilab, R, Julia, C++ with STL,...). The principle of these techniques is general, we present it for the assembly of several finite element matrices in arbitrary dimension, in the P1 finite element case. We also provide an extension of the algorithms to the case of a system of PDE's. Then we give an extension to piecewise polynomials of higher order. We compare numerically the performance of these algorithms in Matlab, Octave and Python, with that in FreeFEM++ and in a compiled language such as C. Examples show that, unlike what is commonly believed, the performance is not radically worse than that of C : in the best/worst cases, selected vector languages are respectively 2.3/3.5 and 2.9/4.1 times slower than C in the scalar and vector cases. We also present numerical results which illustrate the computational costs of these algorithms compared to standard algorithms and to other recent ones.

preprint2014arXiv

Generalized Kirchhoff approximation for Helmholtz equation

We give integral formulas to approximate solutions of Dirichlet and Neumann problems for Helmholtz equation at high frequencies. These approximations are valid in the complementary of a union of convex compact obstacles. The first step of the iterative procedure is the classical Kirchhoff approximation. Convergence is proved by comparison with the geometrical optics asymptotics. The method is shown to be numerically stable.

preprint2013arXiv

An efficient way to perform the assembly of finite element matrices in Matlab and Octave

We describe different optimization techniques to perform the assembly of finite element matrices in Matlab and Octave, from the standard approach to recent vectorized ones, without any low level language used. We finally obtain a simple and efficient vectorized algorithm able to compete in performance with dedicated software such as FreeFEM++. The principle of this assembly algorithm is general, we present it for different matrices in the P1 finite elements case and in linear elasticity. We present numerical results which illustrate the computational costs of the different approaches

preprint2010arXiv

A probabilistic algorithm approximating solutions of a singular PDE of porous media type

The object of this paper is a one-dimensional generalized porous media equation (PDE) with possibly discontinuous coefficient $β$, which is well-posed as an evolution problem in $L^1(\mathbb{R})$. In some recent papers of Blanchard et alia and Barbu et alia, the solution was represented by the solution of a non-linear stochastic differential equation in law if the initial condition is a bounded integrable function. We first extend this result, at least when $β$ is continuous and the initial condition is only integrable with some supplementary technical assumption. The main purpose of the article consists in introducing and implementing a stochastic particle algorithm to approach the solution to (PDE) which also fits in the case when $β$ is possibly irregular, to predict some long-time behavior of the solution and in comparing with some recent numerical deterministic techniques.