Algorithms 70-71. Determination of a quadratic spline function with given values of the integrals in subintervals
Algorithms 73-76. Four algorithms for evaluation of improper integrals
Computer implementation of a coupled boundary and finite element methods for the steady exterior Oseen problem.
Fast and correctly rounded logarithms in double-precision
This article is a case study in the implementation of a portable, proven and efficient correctly rounded elementary function in double-precision. We describe the methodology used to achieve these goals in the crlibm library. There are two novel aspects to this approach. The first is the proof framework, and in general the techniques used to balance performance and provability. The second is the introduction of processor-specific optimization to get performance equivalent to the best current...
Is GPU the future of Scientific Computing ?
These past few years, new types of computational architectures based on graphics processors have emerged. These technologies provide important computational resources at low cost and low energy consumption. Lots of developments have been done around GPU and many tools and libraries are now available to implement efficiently softwares on those architectures.This article contains the two contributions of the mini-symposium about GPU organized by Loïc Gouarin (Laboratoire de Mathématiques d’Orsay),...
Numerical analysis of the MFS for certain harmonic problems
The Method of Fundamental Solutions (MFS) is a boundary-type meshless method for the solution of certain elliptic boundary value problems. In this work, we investigate the properties of the matrices that arise when the MFS is applied to the Dirichlet problem for Laplace’s equation in a disk. In particular, we study the behaviour of the eigenvalues of these matrices and the cases in which they vanish. Based on this, we propose a modified efficient numerical algorithm for the solution of the problem...
Numerical analysis of the MFS for certain harmonic problems
The Method of Fundamental Solutions (MFS) is a boundary-type meshless method for the solution of certain elliptic boundary value problems. In this work, we investigate the properties of the matrices that arise when the MFS is applied to the Dirichlet problem for Laplace's equation in a disk. In particular, we study the behaviour of the eigenvalues of these matrices and the cases in which they vanish. Based on this, we propose a modified efficient numerical algorithm for the solution of the problem...
Object oriented design philosophy for scientific computing
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
Object oriented design philosophy for scientific computing
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
Some new technics regarding the parallelisation of ZéBuLoN, an object oriented finite element code for structural mechanics
A finite element code, called ZéBuLoN was parallelised some years ago. This code is entirely written using an object oriented framework (C++ is the support language). The aim of this paper is to present some problems which arose during the parallelization, and some innovative solutions. Especially, a new concept of message passing is presented which allows to take into account SMP machines while still using the parallel virtual machine abstraction.
Some new technics regarding the parallelisation of ZéBuLoN, an object oriented finite element code for structural mechanics
A finite element code, called ZéBuLoN was parallelised some years ago. This code is entirely written using an object oriented framework (C++ is the support language). The aim of this paper is to present some problems which arose during the parallelization, and some innovative solutions. Especially, a new concept of message passing is presented which allows to take into account SMP machines while still using the parallel virtual machine abstraction.