scholarly journals Matrix decomposition algorithms for arbitrary order C0 tensor product finite element systems

2015 ◽  
Vol 275 ◽  
pp. 162-182 ◽  
Author(s):  
Kui Du ◽  
Graeme Fairweather ◽  
Weiwei Sun
2008 ◽  
Vol 52 (1) ◽  
pp. 1-23 ◽  
Author(s):  
Bernard Bialecki ◽  
Graeme Fairweather ◽  
David B. Knudson ◽  
D. Abram Lipman ◽  
Que N. Nguyen ◽  
...  

CALCOLO ◽  
2021 ◽  
Vol 58 (2) ◽  
Author(s):  
Francesca Bonizzoni ◽  
Guido Kanschat

AbstractA finite element cochain complex on Cartesian meshes of any dimension based on the $$H^1$$ H 1 -inner product is introduced. It yields $$H^1$$ H 1 -conforming finite element spaces with exterior derivatives in $$H^1$$ H 1 . We use a tensor product construction to obtain $$L^2$$ L 2 -stable projectors into these spaces which commute with the exterior derivative. The finite element complex is generalized to a family of arbitrary order.


2018 ◽  
Vol 8 (1) ◽  
pp. 210-217
Author(s):  
Marek Stabrowski

Abstract This paper presents numerical experiments with assorted versions of parallel LU matrix decomposition algorithms (Gauss and Crout algorithm). The tests have been carried out on the hardware platform with fourcore Skylake processor featuring hyperthreading technology doubling virtually core number. Parallelization algorithms have been implemented with the aid of classic POSIX threads library. Experiments have shown that basic 4-thread acceleration of all parallel implementations is almost equal to the number of threads/processors. Both algorithms are worth considering in real-world applications (Florida University collection). Gauss algorithm is a better performer, with respect to timing, in the case of matrices with lower density of nonzeros, as opposed to higher density matrices. The latter are processed more efficiently with the aid of Crout algorithm implementation.


2012 ◽  
Vol 45 (2) ◽  
pp. 299-308 ◽  
Author(s):  
R. Oishi-Tomiyasu ◽  
M. Yonemura ◽  
T. Morishima ◽  
A. Hoshikawa ◽  
S. Torii ◽  
...  

Z-Rietveldis a program suite for Rietveld analysis and the Pawley method; it was developed for analyses of powder diffraction data in the Materials and Life Science Facility of the Japan Proton Accelerator Research Complex. Improvements have been made to the nonlinear least-squares algorithms ofZ-Rietveldso that it can deal with singular matrices and intensity non-negativity constraints. Owing to these improvements,Z-Rietveldsuccessfully executes the Pawley method without requiring any constraints on the integrated intensities, even in the case of severely or exactly overlapping peaks. In this paper, details of these improvements are presented and their advantages discussed. A new approach to estimate the number of independent reflections contained in a powder pattern is introduced, and the concept of good reflections proposed by Sivia [J. Appl. Cryst.(2000),33, 1295–1301] is shown to be explained by the presence of intensity non-negativity constraints, not the intensity linear constraints.


Sign in / Sign up

Export Citation Format

Share Document