Search
Filters
total: 7785
-
Catalog
- Publications 6029 available results
- Journals 67 available results
- Conferences 7 available results
- People 67 available results
- Inventions 2 available results
- Projects 2 available results
- Laboratories 1 available results
- e-Learning Courses 86 available results
- Events 6 available results
- Open Research Data 1518 available results
displaying 1000 best results Help
Search results for: GENERALIZED EIGENVALUE PROBLEM, FEM, COMPLEX-VALUED SPARSE MATRIX PENCIL, GPU, MAXWELL’S EQUATIONS.
-
A GPU Solver for Sparse Generalized Eigenvalue Problems with Symmetric Complex-Valued Matrices Obtained Using Higher-Order FEM
PublicationThe paper discusses a fast implementation of the stabilized locally optimal block preconditioned conjugate gradient (sLOBPCG) method, using a hierarchical multilevel preconditioner to solve nonHermitian sparse generalized eigenvalue problems with large symmetric complex-valued matrices obtained using the higher-order finite-element method (FEM), applied to the analysis of a microwave resonator. The resonant frequencies of the low-order...
-
A memory efficient and fast sparse matrix vector product on a Gpu
PublicationThis paper proposes a new sparse matrix storage format which allows an efficient implementation of a sparse matrix vector product on a Fermi Graphics Processing Unit (GPU). Unlike previous formats it has both low memory footprint and good throughput. The new format, which we call Sliced ELLR-T has been designed specifically for accelerating the iterative solution of a large sparse and complex-valued system of linear equations arising...
-
Block Conjugate Gradient Method with Multilevel Preconditioning and GPU Acceleration for FEM Problems in Electromagnetics
PublicationIn this paper a GPU-accelerated block conjugate gradient solver with multilevel preconditioning is presented for solving large system of sparse equations with multiple right hand-sides (RHSs) which arise in the finite-element analysis of electromagnetic problems. We demonstrate that blocking reduces the time to solution significantly and allows for better utilization of the computing power of GPUs, especially when the system matrix...
-
Single and Dual-GPU Generalized Sparse Eigenvalue Solvers for Finding a Few Low-Order Resonances of a Microwave Cavity Using the Finite-Element Method
PublicationThis paper presents two fast generalized eigenvalue solvers for sparse symmetric matrices that arise when electromagnetic cavity resonances are investigated using the higher-order finite element method (FEM). To find a few loworder resonances, the locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm with null-space deflation is applied. The computations are expedited by using one or two graphical processing...
-
Tuning matrix-vector multiplication on GPU
PublicationA matrix times vector multiplication (matvec) is a cornerstone operation in iterative methods of solving large sparse systems of equations such as the conjugate gradients method (cg), the minimal residual method (minres), the generalized residual method (gmres) and exerts an influence on overall performance of those methods. An implementation of matvec is particularly demanding when one executes computations on a GPU (Graphics...
-
GPU-accelerated finite element method
PublicationIn this paper the results of the acceleration of computations involved in analysing electromagnetic problems by means of the finite element method (FEM), obtained with graphics processors (GPU), are presented. A 4.7-fold acceleration was achieved thanks to the massive parallelization of the most time-consuming steps of FEM, namely finite-element matrix-generation and the solution of a sparse system of linear equations with the...
-
Preconditioners with Low Memory Requirements for Higher-Order Finite-Element Method Applied to Solving Maxwell’s Equations on Multicore CPUs and GPUs
PublicationThis paper discusses two fast implementations of the conjugate gradient iterative method using a hierarchical multilevel preconditioner to solve the complex-valued, sparse systems obtained using the higher order finite-element method applied to the solution of the time-harmonic Maxwell equations. In the first implementation, denoted PCG-V, a classical V-cycle is applied and the system of equations on the lowest level is solved...
-
Tuning a Hybrid GPU-CPU V-Cycle Multilevel Preconditioner for Solving Large Real and Complex Systems of FEM Equations
PublicationThis letter presents techniques for tuning an accelerated preconditioned conjugate gradient solver with a multilevel preconditioner. The solver is optimized for a fast solution of sparse systems of equations arising in computational electromagnetics in a finite element method using higher-order elements. The goal of the tuning is to increase the throughput while at the same time reducing the memory requirements in order to allow...
-
A Task-Scheduling Approach for Efficient Sparse Symmetric Matrix-Vector Multiplication on a GPU
PublicationIn this paper, a task-scheduling approach to efficiently calculating sparse symmetric matrix-vector products and designed to run on Graphics Processing Units (GPUs) is presented. The main premise is that, for many sparse symmetric matrices occurring in common applications, it is possible to obtain significant reductions in memory usage and improvements in performance when the matrix is prepared in certain ways prior to computation....
-
Implementation of algebraic procedures on the GPU using CUDA architecture on the example of generalized eigenvalue problem
Publication -
Finite element matrix generation on a GPU
PublicationThis paper presents an efficient technique for fast generation of sparse systems of linear equations arising in computational electromagnetics in a finite element method using higher order elements. The proposed approach employs a graphics processing unit (GPU) for both numerical integration and matrix assembly. The performance results obtained on a test platform consisting of a Fermi GPU (1x Tesla C2075) and a CPU (2x twelve-core...
-
GPU-Accelerated LOBPCG Method with Inexact Null-Space Filtering for Solving Generalized Eigenvalue Problems in Computational Electromagnetics Analysis with Higher-Order FEM
PublicationThis paper presents a GPU-accelerated implementation of the Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) method with an inexact nullspace filtering approach to find eigenvalues in electromagnetics analysis with higherorder FEM. The performance of the proposed approach is verified using the Kepler (Tesla K40c) graphics accelerator, and is compared to the performance of the implementation based on functions from...
-
GENERAL DYNAMIC PROJECTING OF MAXWELL EQUATIONS
PublicationA complete – system of Maxwell equations is splitting into independent subsystems by means of a special dynamic projecting technique. The technique relies upon a direct link between field components that determine correspondent subspaces. The explicit form of links and corresponding subspace evolution equations are obtained in conditions of certain symmetry, it is illustrated by examples of spherical and quasi-one-dimensional waves.
-
Accuracy, Memory and Speed Strategies in GPU-based Finite-Element Matrix-Generation
PublicationThis paper presents strategies on how to optimize GPU-based finite-element matrix-generation that occurs in the finite-element method (FEM) using higher order curvilinear elements. The goal of the optimization is to increase the speed of evaluation and assembly of large finite-element matrices on a single GPU (Graphics Processing Unit) while maintaining the accuracy of numerical integration at the desired level. For this reason,...
-
Jacobi and gauss-seidel preconditioned complex conjugate gradient method with GPU acceleration for finite element method
PublicationIn this paper two implementations of iterative solvers for solving complex symmetric and sparse systems resulting from finite element method applied to wave equation are discussed. The problem under investigation is a dielectric resonator antenna (DRA) discretized by FEM with vector elements of the second order (LT/QN). The solvers use the preconditioned conjugate gradient (pcg) method implemented on Graphics Processing Unit (GPU)...
-
On tracking properties of real-valued generalized adaptive notch filters
PublicationGeneralized adaptive notch filters (GANFs) are used for identification/tracking of quasi-periodically varying dynamic systems and can be considered an extension, to the system case, of classical adaptive notch filters. The paper presents results of local performance analysis of a real-valued GANF algorithm, i.e., algorithm designed to track parameters of a real-valued system. This is an extension of the previous work which focused...
-
GPU-Accelerated Finite-Element Matrix Generation for Lossless, Lossy, and Tensor Media [EM Programmer's Notebook]
PublicationThis paper presents an optimization approach for limiting memory requirements and enhancing the performance of GPU-accelerated finite-element matrix generation applied in the implementation of the higher-order finite-element method (FEM). It emphasizes the details of the implementation of the matrix-generation algorithm for the simulation of electromagnetic wave propagation in lossless, lossy, and tensor media. Moreover, the impact...
-
Generalized adaptive notch smoothers for real-valued signals and systems
PublicationSystems with quasi-periodically varying coefficients can be tracked using the algorithms known as generalized adaptive notch filters (GANFs). GANF algorithms can be considered an extension, to the system case, of classical adaptive notch filters (ANFs). We show that estimation accuracy of the existing algorithms, as well as their robustness to the choice of design parameters, can be considerably improved by means of compensating...
-
Anisotropic Orlicz–Sobolev spaces of vector valued functions and Lagrange equations
PublicationIn this paper we study some properties of anisotropic Orlicz and Orlicz–Sobolev spaces of vector valued functions for a special class of G-functions. We introduce a variational setting for a class of Lagrangian Systems. We give conditions which ensure that the principal part of variational functional is finitely defined and continuously differentiable on Orlicz–Sobolev space.
-
Quasi-solutions for generalized second order differential equations with deviating arguments
PublicationThis paper deal with boundary value problems for generalized second order differential equations with deviating arguments. Existence of quasi-solutions and solutions are proved by monotone iterative method. Examples with numerical results are added.
-
Communication and Load Balancing Optimization for Finite Element Electromagnetic Simulations Using Multi-GPU Workstation
PublicationThis paper considers a method for accelerating finite-element simulations of electromagnetic problems on a workstation using graphics processing units (GPUs). The focus is on finite-element formulations using higher order elements and tetrahedral meshes that lead to sparse matrices too large to be dealt with on a typical workstation using direct methods. We discuss the problem of rapid matrix generation and assembly, as well as...
-
Multimodal Genetic Algorithm with Phase Analysis to Solve Complex Equations of Electromagnetic Analysis
PublicationIn this contribution, a new genetic-algorithm-based method of finding roots and poles of a complex function of a complex variable is presented. The algorithm employs the phase analysis of the function to explore the complex plane with the use of the genetic algorithm. Hence, the candidate regions of root and pole occurrences are selected and verified with the use of discrete Cauchy's argument principle. The algorithm is evaluated...
-
Inverse Nonlinear Eigenvalue Problem Framework for the Synthesis of Coupled-Resonator Filters With Nonresonant Nodes and Arbitrary Frequency-Variant Reactive Couplings
PublicationA novel, general circuit-level description of coupledresonator microwave filters is introduced in this article. Unlike well-established coupling-matrix models based on frequency-invariant couplings or linear frequency-variant couplings (LFVCs), a model with arbitrary reactive frequencyvariant coupling (AFVC) networks is proposed. The engineered formulation is more general than prior-art ones—with the only restriction that the coupling...
-
Estimation and tracking of complex-valued quasi-periodically varying systems
PublicationW artykule rozważany jest problem identyfikacji obiektów o parametrach zmieniających się w sposób pseudookresowy. Przedstawiono w nim algorytm oparty o metodę funkcji bazowych umożliwiający śledzenie takich obiektów oraz pokazano atrakcyjne z punktu widzenia złożoności obliczeń jego wersje zdekomponowane. Przydatność rozważanych algorytmów uzasadniono porównując je z rozwiązaniami innych autorów.
-
koło naukowe - Maxwell
e-Learning CoursesKoło naukowe Maxwell
-
DEM-FEM model of highly saturated soil motion due to seepage force
PublicationArtykuł poświęcono nowej propozycji rozwiązania problemu modelowania zjawiska transportu, które występuje wewnątrz wysoko-nawodnionego gruntu ziarnistego na skutek ciśnienia spływowego. Teoria oparta została na połączeniu metody elementów skończonych (FEM - Finite Element Method) - rozwiązania równań poro-elastycznych oraz metody poszczególnych elementów (DEM - Distinct Element Method), które dały podwójny model numeryczny. FEM...
-
Green's function for the wavized Maxwell fish-eye problem
PublicationRozpatrzono niezależne od czasu skalarne równanie falowe dla ośrodka typu ''rybie oko'' Maxwella w przestrzeni R^N (N >=2). Pokazano, że równanie to posiada unikalne własności transformacyjne względem inwersji w pewnej klasie hipersfer. Wykorzystano ten fakt do znalezienia zamkniętej postaci funkcji Greena, oraz uogólnionej funkcji Greena, dla wyjściowego równania.
-
GPU-Accelerated 3D Mesh Deformation for Optimization Based on the Finite Element Method
PublicationThis paper discusses a strategy for speeding up the mesh deformation process in the design-byoptimization of high-frequency components involving electromagnetic field simulations using the 3D finite element method (FEM). The mesh deformation is assumed to be described by a linear elasticity model of a rigid body; therefore, each time the shape of the device is changed, an auxiliary elasticity finite-element problem must be solved....
-
A Novel Coupling Matrix Synthesis Technique for Generalized Chebyshev Filters With Resonant Source–Load Connection
PublicationThis paper reports a novel synthesis method for microwave bandpass filters with resonant source–load connection. In effect, a network realizing N+1 transmission zeros (where N is the number of reflection zeros) is obtained. The method is based on a prototype transversal coupling matrix (N+2, N+2) with source and load connected by a resonant circuit formed by a capacitor in parallel with a frequency-invariant susceptance. To complement...
-
Efficiency of acoustic heating in the Maxwell fluid
PublicationThe nonlinear effects of sound in a fluid describing by the Maxwell model of the viscous stress tensor is the subject of investigation. Among other, viscoelastic biological media belong to this non-newtonian type of fluids. Generation of heating of the medium caused by nonlinear transfer of acoustic energy, is discussed in details. The governing equation of acoustic heating is derived by means of the special linear combination...
-
Efficiency of acoustic heating in the Maxwell fluid
PublicationThe nonlinear effects of sound in a fluid describing by the Maxwell model of the viscous stress tensor is the subject of investigation. Among other, viscoelastic biological media belong to this non-newtonian type of fluids. Generation of heating of the medium caused by nonlinear transfer of acoustic energy, is discussed in details. The governing equation of acoustic heating is derived by means of the special linear combination...
-
Solution of coupled integral equations for quantum scattering in the presence of complex potentials
PublicationIn this paper, we present a method to compute solutions of coupled integral equations for quantum scattering problems in the presence of a complex potential. We show how the elastic and absorption cross sections can be obtained from the numerical solution of these equations in the asymptotic region at large radial distances.
-
Modal FEM Analysis of Ferrite Resonant Structures
PublicationThe finite-element method (FEM) is applied for modal analysis of ferrite-loaded spherical resonators. To improve the efficiency of the numerical calculations, the body-of-revolution (BOR) technique is utilized. Due to the frequency-dependent ferrite permeability, FEM leads to a nonlinear eigenvalue problem that is challenging to solve. To this end, Beyn’s method is proposed. The effectiveness of the proposed approach is confirmed...
-
Performance evaluation of parallel background subtraction on GPU platforms
PublicationImplementation of the background subtraction algorithm on parallel GPUs is presented. The algorithm processes video streams and extracts foreground pixels. The work focuses on optimizing parallel algorithm implementation by taking into account specific features of the GPU architecture, such as memory access, data transfers and work group organization. The algorithm is implemented in both OpenCL and CUDA. Various optimizations of...
-
The maximum edge-disjoint paths problem in complete graphs
PublicationRozważono problem ścieżek krawędziowo rozłącznych w grafach pełnych. Zaproponowano wielomianowe algorytmy: 3.75-przybliżony (off-line) oraz 6.47-przybliżony (on-line), poprawiając tym samym wyniki wcześniej znane z literatury [P. Carmi, T. Erlebach, Y. Okamoto, Greedy edge-disjoint paths in complete graphs, in: Proc. 29th Workshop on Graph Theoretic Concepts in Computer Science, in: LNCS, vol. 2880, 2003, pp. 143-155]. Ponadto...
-
Main complications connected with detection, identification and determination of trace organic constituents in complex matrix samples
PublicationIt is well known that some problems with the determination of organic analytes at trace level can occur. This issue is connected with contamination during each stage of the analytical procedure from sampling to sample preparation up to chromatographic analysis, which often leads to false-positive or overestimated results. Another problem associated with determination of analytes occurs at trace- and ultra-trace level is a background...
-
Investigation of Vibro-Acoustical Properties Using Eigenvalue Problem Solution and Wavelets
Publication -
Comparison of Optimization Techniques for Coupling Matrix Synthesis Using Eigenvalue Based Approach
Publication -
Systems of differential equations with a singular matrix
PublicationPraca dotyczy układów równań różniczkowych z macierzą osobliwą. Szczególnymi przypadkami takich układów są równania różniczkowo-algebraiczne. Korzystając z odpowiedniej nierówności różniczkowej dla układów równań i metody quasi linearyzacji pokazano, iż odpowiednie monotoniczne ciągi są zbieżne do rozwiązania naszego problemu i jest to zbieżność kwadratowa. O funkcji f występującej po prawej stronie zagadnienia zakładano...
-
Multilevel model order reduction with generalized compression of boundaries for 3-d FEM electromagnetic analysis
PublicationThis paper presents a multilevel Model Order Reduction technique for a 3-D electromagnetic Finite Element Method analysis. The reduction process is carried out in a hierarchical way and involves several steps which are repeated at each level. This approach brings about versatility and allows one to efficiently analyze complex electromagnetic structures. In the proposed multilevel reduction the entire computational domain is covered...
-
Vortex flow caused by periodic and aperiodic sound in a relaxing maxwell fluid
PublicationThis paper concerns the description of vortex flow generated by periodic and aperiodic sound in relaxing Maxwell fluid. The analysis is based on governing equation of vorticity mode, which is a result of decomposition of the hydrodynamic equations for fluid flow with relaxation and thermal conductivity into acoustical and non-acoustical parts. The equation governing vorticity mode uses only instantaneous, not averaged over sound...
-
General quasilinearization method for systems of differential equations with a singular matrix
PublicationStosując metodę kwazilinearyzacji badano problem rozwiązań przybliżonych dla układów równań różniczkowych z macierzą osobliwą. Pokazano zbieżność kwadratową odpowiednich ciągów monotonicznych.
-
Complex Variables and Elliptic Equations
Journals -
Generalized quazilinearization for systems of degenerate singular perturbation problem
PublicationPraca dotyczy ogólnej metody kwazilinearyzacji dla układów równań różniczkowych z parametrami.
-
Sparse autoregressive modeling
PublicationIn the paper the comparison of the popular pitch determination (PD) algorithms for thepurpose of elimination of clicks from archive audio signals using sparse autoregressive (SAR)modeling is presented. The SAR signal representation has been widely used in code-excitedlinear prediction (CELP) systems. The appropriate construction of the SAR model is requiredto guarantee model stability. For this reason the signal representation...
-
On the Chaplyghin method for generalized solutions of partial differential functional equations
PublicationW pracy tej rozważa się zagadnienia początkowo-brzegowe dla semiliniowych równań różniczkowo-funkcyjnych cząstkowych pierwszego rzędu. Przy naturalnych założeniach o funkcjach danych, istnieje ciąg Czapłygina i jest on zbieżny do rozwiązania w sensie Caratheodory'ego wyjściowego zagadnienia. Podane jest oszacowanie błędu rozwiązania przybliżonego. Wykazuje się, iż metoda Czapłygina dla powyższego zagadnienia jest równoważna metodzie...
-
Differential equations with a singular matrix and convergence of approximate solutions
PublicationZastosowano metodę kwazilinearyzacji dla równań różniczkowych z macierzą osobliwą aby pokazać, że odpowiednie ciągi rozwiązań przybliżonych są kwadratowo zbieżne do rozwiązania naszego problemu.Metodę tę zastosowano do układów takich równań różniczkowych.
-
Degree of entaglement as a physically ill-posted problem: The case of entaglement with vacuum
PublicationAnalizujemy przypadek fotonu w superpozycji różnych modów i zadajemy pytanie o stopień ich splątania z próżnią. Problem okazuje się być źle postawiony, gdyż nie wiemy którą reprezentację algebry CCR wybrać dla kwantowania pola. Gdy dokonamy wyboru jednoznacznie możemy rozwiązać zagadnienie splątania. Tak więc trudność nie leży w matematyce lecz w fizyce problemu.
-
Macromodels for efficient FEM simulations of waveguides and resonators
PublicationThis paper introduces a novel technique for enhancing the efficiency of the finite element method (FEM) by incorporating special modules, called macromodels, into the standard eigenvalue formulation. The number of unknowns in the separated macromodel subdomain can be significantly reduced by orthogonal projection, using the efficient nodal order reduction algorithm. The idea of macromodels implementation is demonstrated on a simple...
-
Numerical Methods for Partial Differential Equations
e-Learning CoursesCourse description: This course focuses on modern numerical techniques for linear and nonlinear elliptic, parabolic and hyperbolic partial differential equations (PDEs), and integral equations fundamental to a large variety of applications in science and engineering. Topics include: formulations of problems in terms of initial and boundary value problems; finite difference and finite element discretizations; boundary element approach;...
-
Dynamic GPU power capping with online performance tracing for energy efficient GPU computing using DEPO tool
PublicationGPU accelerators have become essential to the recent advance in computational power of high- performance computing (HPC) systems. Current HPC systems’ reaching an approximately 20–30 mega-watt power demand has resulted in increasing CO2 emissions, energy costs and necessitate increasingly complex cooling systems. This is a very real challenge. To address this, new mechanisms of software power control could be employed. In this...
-
Generalized Euler method for first order partial differential functional equations
PublicationW pracy prezentowana jest nowa klasa metod numerycznych dla nieliniowych równań różniczkowo funkcyjnych pierwszego rzędu.Rozwiązania klasyczne zagadnień początkowo brzegowych przybliżane są w tej pracy przez rozwiązania odpowiedniego układu quasilininowego równań różnicowych. Podajemy kompletną analizę zbieżności metod i pokazujemy na przykładach, iż nowa metoda jest zauważalnie lepsza niż klasyczne schematy różnicowe. Dowód stabilności...
-
The generalized quasilinearization for integro-differential equations of Volterra type on time scales
PublicationBadano równania całkowo-różniczkowe on ''time scales'' i podano warunkidostateczne na zbieżność metody kwazilinearyzacji do jego rozwiązania. Podano warunki na to, aby zbieżność ta była kwadratową.
-
Generalized method of lines for nonlinear first order partial differential equations
PublicationKlasyczne rozwiązania zagadnień początkowych oraz początkowo brzegowych są przybliżane za pomocą rozwiązań równań różniczkowo różnicowych. Skonstruowana jest metoda prostych polegająca na dyskretyzacji wyjściowego równania względem zmiennych przestrzennych. Przedstawiony w pracy schemat bazuje na metodzie linearyzacji dla zagadnień nieliniowych. W pracy zastosowano metodę quasilinearyzacji polegającą na zamianie nieliniowego równania...
-
Generalized Euler method for nonlinear first order partial differential equations.
PublicationKlasyczne rozwiązania nieliniowych równań różniczkowych cząstkowych pierwszego rzędu są aproksymowane w tej pracy za pomocą rozwiązań quasiliniowych układów równań różnicowych. Podstawowa idea pracy jest oparta na teorii charakterystyk. Podane są warunki wystarczające dla zbieżności metody. Dowód stabilności schematu różnicowego wykorzystuje metodę porównawczą z nieliniowymi oszacowaniami typu Perrona dla danych funkcji.Podane...
-
Multimodal Particle Swarm Optimization with Phase Analysis to Solve Complex Equations of Electromagnetic Analysis
PublicationIn this paper, a new meta-heuristic method of finding roots and poles of a complex function of a complex variable is presented. The algorithm combines an efficient space exploration provided by the particle swarm optimization (PSO) and the classification of root and pole occurrences based on the phase analysis of the complex function. The method initially generates two uniformly distributed populations of particles on the complex...
-
Communication Model Order Reduction in Hybrid Methods Involving Generalized Impedance Matrix
PublicationA novel strategy for the efficient analysis of frequency-domain scattering electromagnetic problems in open and closed domains is presented. A fully automatic model-order reduction technique, called the enhanced reduced-basis method, is applied to increase the efficiency of the hybrid approach, which combines the finite-element and mode-matching methods. Numerical tests show that the proposed algorithm yields reliable and highly...
-
The generalized Banach match-box problem: Application in disc storage management
Publication -
Remarks on the convergence of an iterative method of solution of generalized least squares problem
PublicationW pracy przedstawiona jest metoda iteracyjna znajdowania regularyzowanego (w sensie Tichonowa) rozwiązania układu równań Ax=b z dowolną macierzą A. Dla danej liczby alfa i wektora g daje ona ciąg przybliżeń zbieżny do rozwiązania (w sensie najmniejszych kwadratów) tego układu. Rozwiązanie to minimalizuje odległość zbioru wszystkich rozwiązań średniokwadratowych układu Ax=b od wektora g. Podane zostało również oszacowanie szybkości...
-
Efficient and Systematic Solution of Real and Complex Eigenvalue Problems Employing Simplex Chain Vertices Searching Procedure
PublicationW artykule zaprezentowano nowatorską metodę, pozwalającą w bardzo efektywny sposób rozwiązywać rzeczywiste i zespolone zagadnienia własne (które przekształcić można w równania nieliniowe) często spotykane w elektrodynamice obliczeniowej. Idea metody oparta jest o śledzenie miejsc zerowych otrzymanego równania nieliniowego. W proponowanej technice zastosowano regularne wielościany, budując z nich łańcuch (ścieżkę) poszukiwanych...
-
GPU Acceleration of Multilevel Solvers for Analysis of Microwave Components With Finite Element Method
PublicationThe letter discusses a fast implementation of the conjugate gradient iterative method with ${rm E}$-field multilevel preconditioner applied to solving real symmetric and sparse systems obtained with vector finite element method. In order to accelerate computations, a graphics processing unit (GPU) was used and significant speed-up (2.61 fold) was achieved comparing to a central processing unit (CPU) based approach. These results...
-
Generalized Gradient Equivariant Multivalued Maps, Approximation and Degree
PublicationConsider the Euclidean space Rn with the orthogonal action of a compact Lie group G. We prove that a locally Lipschitz G-invariant mapping f from Rn to R can be uniformly approximated by G-invariant smooth mappings g in such a way that the gradient of g is a graph approximation of Clarke’s generalized gradient of f . This result enables a proper development of equivariant gradient degree theory for a class of set-valued gradient...
-
A Self-Adaptive Complex Root Tracing Algorithm for the Analysis of Propagation and Radiation Problem
PublicationAn improved complex root tracing algorithm for radiation and propagation issues is proposed. The approach is based on a self-adaptive discretization of Cauchy’s argument principle for a C × R space and requires a reduced number of function calls in comparison to other procedures presented in the literature. A few different examples concerning propagation and radiation problems have been considered to verify the validity and efficiency...
-
Generalized adaptive notch smoothing revisited
PublicationThe problem of identification of quasi-periodically varying dynamic systems is considered. This problem can be solved using generalized adaptive notch filtering (GANF) algorithms. It is shown that the accuracy of parameter estimates can be significantly increased if the results obtained from GANF are further processed using a cascade of appropriately designed filters. The resulting generalized adaptive notch smoothing (GANS) algorithm...
-
Parallelization of large vector similarity computations in a hybrid CPU+GPU environment
PublicationThe paper presents design, implementation and tuning of a hybrid parallel OpenMP+CUDA code for computation of similarity between pairs of a large number of multidimensional vectors. The problem has a wide range of applications, and consequently its optimization is of high importance, especially on currently widespread hybrid CPU+GPU systems targeted in the paper. The following are presented and tested for computation of all vector...
-
Equivariant degree of convex-valued applied to set-valued BVP
PublicationW pracy określono stopień ekwiwariantny dla odwzorowań o wartościach zwartych i wypukłych. Zastosowano go do dowodu twierdzenia o istnieniu wielunietrywialnych rozwiązań problemu brzegowego typu Dirichleta dla inkluzji różniczkowych z warunkami wzrostu typu Bernsteina-Nagumo, z dodatkowymi symetriami
-
Generalized solutions of mixed problems for first-order partial functional differential equations
PublicationW pracy udowodniono twierdzenie o istnieniu rozwiązań i ich ciągłej zależności od warunków początkowo brzegowych. Do przekształcenia zagadnienia oryginalnego w układ równań całkowo funkcyjnych typu Volterry użyto metody bicharakterystyk. Istnienie rozwiązań udowodnione jest metodą kolejnych przybliżeń, przy użyciu twierdzeń o nierównościach całkowych. Rozwiązania klasyczne układów równań całkowo funkcyjnych prowadzą do rozwiązań...
-
Efficient model order reduction for FEM analysis of waveguide structures and resonators
PublicationAn efficient model order reduction method for three-dimensional Finite Element Method (FEM) analysis of waveguide structures is proposed. The method is based on the Efficient Modal Order Reduction (ENOR) algorithm for creating macro-elements in cascaded subdomains. The resulting macro-elements are represented by very compact submatrices, leading to significant reduction of the overall number of unknowns. The efficiency of the model...
-
A more colorful hat problem
PublicationThe topic is the hat problem in which each of n players is randomly fitted with a blue or red hat. Then everybody can try to guess simultaneously his own hat color by looking at the hat colors of the other players. The team wins if at least one player guesses his hat color correctly, and no one guesses his hat color wrong; otherwise the team loses. The aim is to maximize the probability of winning. We consider a generalized hat...
-
Estimation of the amplitude of the signal for the active optical gesture sensor with sparse detectors
PublicationIn this paper we deal with the problem of precise gesture recognition for the active optical proximity sensor with sparse 8 photodiodes. We particularly focus on developing the method of estimating the real, usually not observable, maximum signal value representing maximum intensity of light reflected from an obstacle present in the front of the sensor. Different configurations of the fingers were used as an obstacle. The Monte Carlo...
-
Equivariant degree of convex-valued maps applied to set-valued BVP
Publication -
Modelling and simulation of GPU processing in the MERPSYS environment
PublicationIn this work, we evaluate an analytical GPU performance model based on Little's law, that expresses the kernel execution time in terms of latency bound, throughput bound, and achieved occupancy. We then combine it with the results of several research papers, introduce equations for data transfer time estimation, and finally incorporate it into the MERPSYS framework, which is a general-purpose simulator for parallel and distributed...
-
Existence of unbounded solutions to parabolic equations with functional dependence
PublicationThe Cauchy problem for nonlinear parabolic differential-functional equations is considered. Under natural generalized Lipschitz-type conditions with weights, the existence and uniqueness of unbounded solutions is obtained in three main cases: (i) the functional dependence u(·); (ii) the functional dependence u(·) and ∂xu(·); (iii) the functional dependence u(·)and the pointwise dependence ∂xu(t,x).
-
Thermal self-action effects for acoustic beams containing fronts in a Maxwell relaxing fluid
PublicationThis paper examines the thermal self-action of acoustic beams in a Maxwell relaxing fluid. This type of thermal self-action differs from that in a Newtonian fluid and behaves differently depending on a ratio of sound period and time of thermodynamic relaxation. The self-action which relates to sound beams containing shock fronts is also discussed. In addition, stationary and non-stationary types of self-action are considered.
-
Electromagnetic Simulations with 3D FEM and Intel Optane Persistent Memory
PublicationAbstract—Intel Optane persistent memory has the potential to induce a change in how high-performance calculations requiring a large system memory capacity are conducted. This article presents what this change may look like in the case of factorization of large sparse matrices describing electromagnetic problems arising in the 3D FEM analysis of passive highfrequency components. In numerical tests, the Intel oneAPI MKL PARDISO was...
-
Equations with Separated Variables on Time Scales
PublicationWe show that the well-known theory for classical ordinary differential equations with separated variables is not valid in case of equations on time scales. Namely, the uniqueness of solutions does not depend on the convergence of appropriate integrals.
-
New achievements in the field of extraction of trace analytes from samples characterized by complex composition of the matrix
PublicationWithout any doubt, the monitoring of compounds present in samples at trace or ultra-trace level usually requires a preliminary step of isolation and/or enrichment of analytes due to the fact that majority of analytical techniques are not sensitive enough for direct determination of trace compounds. On the other hand, sample preparation is considered as crucial part of whole analytical procedures, in particular in samples characterized...
-
Mathematical analysis of the lasing eigenvalue problem for the optical modes in a layered dielectric cavity with a quantum well and distributed Bragg reflectors
Publication -
Acceleration of the DGF-FDTD method on GPU using the CUDA technology
PublicationWe present a parallel implementation of the discrete Green's function formulation of the finite-difference time-domain (DGF-FDTD) method on a graphics processing unit (GPU). The compute unified device architecture (CUDA) parallel computing platform is applied in the developed implementation. For the sake of example, arrays of Yagi-Uda antennas were simulated with the use of DGF-FDTD on GPU. The efficiency of parallel computations...
-
Advanced Scalar-valued Intensity Measures for Residual Drift Prediction of SMRFs with Fluid Viscous Dampers
PublicationMaximum Residual Inter-story Drift Ratio (RIDRmax) plays an important role to specify the state of a structure after severe earthquake and the possibility of repairing the structure. Therefore, it is necessary to predict the RIDRmax of Steel Moment-Resisting Frames (SMRFs) with high reliability by employing powerful Intensity Measures (IMs). This study investigates the efficiency and sufficiency of scalar-valued IMs for predicting...
-
A New Type of Macro-Elements for Efficient Two-Dimensional FEM Analysis
PublicationThis letter deals with a model order reduction technique applicable for driven and eigenvalue problems solved using the finite element method (FEM). It allows one to efficiently compute electromagnetic parameters of structures comprising small features that require strong local mesh refinement. The subdomains of very fine mesh are separated from the global domain as so called macro-elements that undergo model reduction. The macro-elements...
-
Fractional differential equations with causal operators
PublicationWe study fractional differential equations with causal operators. The existence of solutions is obtained by applying the successive approximate method. Some applications are discussed including also the case when causal operator Q is a linear operator. Examples illustrate some results.
-
Generalized Chebyshev Bandpass Filters With Frequency-Dependent Couplings Based on Stubs
PublicationThis paper presents an accurate synthesis method for inline and cross-coupled generalized Chebyshev bandpass filters with frequency-dependent couplings implemented via open and short stubs. The technique involves the synthesis of a lumped-element prototy pe in the form of a coupling matrix with a frequency-dependent term and the conversion of this prototype to a distributed-element mode l composed of sections of TEM lines. This...
-
Performance of Vector-valued Intensity Measures for Estimating Residual Drift of Steel MRFs with Viscous Dampers
PublicationViscous Dampers (VDs) are widely used as passive energy dissipation system for improving seismic performance levels especially in retrofitting of buildings. Residual Inter-story Drift Ratio (R-IDR) is another important factor that specifies the condition of building after earthquake. The values of R-IDR illustrates the possibility of retrofitting and repairing of a building. Therefore, this study aims to explore the vector-valued...
-
On root finding algorithms for complex functions with branch cuts
PublicationA simple and versatile method is presented, which enhances the complex root finding process by eliminating branch cuts and branch points in the analyzed domain. For any complex function defined by a finite number of Riemann sheets, a pointwise product of all the surfaces can be obtained. Such single-valued function is free of discontinuity caused by branch cuts and branch points. The roots of the new function are the same as the...
-
Parabolic Equations with Functional Dependence
PublicationWe consider the Cauchy problem for nonlinear parabolic equations with functional dependence and prove theorems on the existence of solutions to parabolic differential-functional equations.
-
Single and Series of Multi-valued Decision Diagrams in Representation of Structure Function
PublicationStructure function, which defines dependency of performance of the system on performance of its components, is a key part of system description in reliability analysis. In this paper, we compare two approaches for representation of the structure function. The first one is based on use of a single Multi-valued Decision Diagram (MDD) and the second on use of a series of MDDs. The obtained results indicate that the series of MDDs...
-
Elimination of clicks from archive speech signals using sparse autoregressive modeling
PublicationThis paper presents a new approach to elimination of impulsivedisturbances from archive speech signals. The proposedsparse autoregressive (SAR) signal representation is given ina factorized form - the model is a cascade of the so-called formantfilter and pitch filter. Such a technique has been widelyused in code-excited linear prediction (CELP) systems, as itguarantees model stability. After detection of noise pulses usinglinear...
-
On the convergence of a finite-difference discretization à la Mickens of the generalized Burgers–Huxley equation
PublicationIn this note, we establish the property of convergence for a finite-difference discretization of a diffusive partial differential equation with generalized Burgers convective law and generalized Hodgkin–Huxley reaction. The numerical method was previously investigated in the literature and, amongst other features of interest, it is a fast and nonlinear technique that is capable of preserving positivity, boundedness and monotonicity....
-
Parallel implementation of the DGF-FDTD method on GPU Using the CUDA technology
PublicationThe discrete Green's function (DGF) formulation of the finite-difference time-domain method (FDTD) is accelerated on a graphics processing unit (GPU) by means of the Compute Unified Device Architecture (CUDA) technology. In the developed implementation of the DGF-FDTD method, a new analytic expression for dyadic DGF derived based on scalar DGF is employed in computations. The DGF-FDTD method on GPU returns solutions that are compatible...
-
A Cost-Effective Method for Reconstructing City-Building 3D Models from Sparse Lidar Point Clouds
PublicationThe recent popularization of airborne lidar scanners has provided a steady source of point cloud datasets containing the altitudes of bare earth surface and vegetation features as well as man-made structures. In contrast to terrestrial lidar, which produces dense point clouds of small areas, airborne laser sensors usually deliver sparse datasets that cover large municipalities. The latter are very useful in constructing digital...
-
Generalized solution of mixed problems for first order partial differential equations with state dependent delays
PublicationW pracy zostało udowodnione twierdzenie o istnieniu i jednoznaczności rozwiązań dla zagadnień początkowo-brzegowych z cząstkowym równaniem różniczkowo-funkcyjnym z opóźnionym argumentem zależnym od funkcji niewiadomej. Użyto metody bicharakterystyk. Jednoznaczność rozwiązań wykazano metodą porównawczą, istnienie - metodą ciągów przybliżeń.
-
Electromagnetic Control and Dynamics of Generalized Burgers’ Nanoliquid Flow Containing Motile Microorganisms with Cattaneo–Christov Relations: Galerkin Finite Element Mechanism
PublicationIn our research work, we have developed a model describing the characteristics of the bio-convection and moving microorganisms in the flows of a magnetized generalized Burgers’ nanoliquid with Fourier’s and Fick’s laws in a stretchable sheet. Considerations have been made to Cattaneo–Christov mass and heat diffusion theory. According to the Cattaneo–Christov relation, the Buongiorno phenomenon for the motion of a nanoliquid in...
-
Multimode systems of nonlinear equations: derivation, integrability, and numerical solutions
PublicationWe consider the propagation of electromagnetic pulses in isotropic media taking a third-order nonlinearityinto account. We develop a method for transforming Maxwell's equations based on a complete set ofprojection operators corresponding to wave-dispersion branches (in a waveguide or in matter) with thepropagation direction taken into account. The most important result of applying the method is a systemof equations describing the...
-
Generalized temperature dependence model for anammox process kinetics
PublicationTemperature is a key operational factor influencing the anammox process kinetics. In particular, at temperatures below 15 °C, the specific anammox activity (SAA) considerably decreases. This study aimed to describe the temperature dependence of the anammox process kinetics in the temperature range from 10 to 55 °C, including the specific characteristics of“cold anammox”. The commonly used Arrhenius and extended and modifiedRatkowsky...
-
Unifying Aspects of Generalized Calculus
PublicationNon-Newtonian calculus naturally unifies various ideas that have occurred over the years in the field of generalized thermostatistics, or in the borderland between classical and quantum information theory. The formalism, being very general, is as simple as the calculus we know from undergraduate courses of mathematics. Its theoretical potential is huge, and yet it remains unknown or unappreciated.
-
Macro-elements and Model Order Reduction for Efficient Three-Dimensional FEM Analysis
PublicationAn efficient model order reduction (MOR) methodology for three dimensional vector finite element method (FEM) is developed to accelerate simulations of the structures containing features that cause strong variations of mesh density. As the result of presented algorithm, FEM subsystems of equations corresponding to the selected refined region are converted into a very compact sets of linear equations, called macro-elements.Numerical...
-
Efficient parallel implementation of crowd simulation using a hybrid CPU+GPU high performance computing system
PublicationIn the paper we present a modern efficient parallel OpenMP+CUDA implementation of crowd simulation for hybrid CPU+GPU systems and demonstrate its higher performance over CPU-only and GPU-only implementations for several problem sizes including 10 000, 50 000, 100 000, 500 000 and 1 000 000 agents. We show how performance varies for various tile sizes and what CPU–GPU load balancing settings shall be preferred for various domain...
-
Justyna Signerska-Rynkowska dr inż.
PeopleI am currently an assistant professor (adjunct) at Gdansk University of Technology (Department of Differential Equations and Mathematics Applications). My scientific interests include dynamical systems theory, chaos theory and their applications to modeling of biological phenomena, especially to neurosciences. In June 2013 I completed PhD in Mathematics at the Institute of Mathematics of Polish Academy of Sciences (IMPAN) (thesis...
-
On the convergence of a nonlinear finite-difference discretization of the generalized Burgers–Fisher equation
PublicationIn this note, we establish analytically the convergence of a nonlinear finite-difference discretization of the generalized Burgers-Fisher equation. The existence and uniqueness of positive, bounded and monotone solutions for this scheme was recently established in [J. Diff. Eq. Appl. 19, 1907{1920 (2014)]. In the present work, we prove additionally that the method is convergent of order one in time, and of order two in space. Some...