Source author record

Zydrunas Gimbutas

Zydrunas Gimbutas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

7 published item(s)

preprint2019arXiv

A fast simple algorithm for computing the potential of charges on a line

We present a fast method for evaluating expressions of the form $$ u_j = \sum_{i = 1,i \not = j}^n \frac{α_i}{x_i - x_j}, \quad \text{for} \quad j = 1,\ldots,n, $$ where $α_i$ are real numbers, and $x_i$ are points in a compact interval of $\mathbb{R}$. This expression can be viewed as representing the electrostatic potential generated by charges on a line in $\mathbb{R}^3$. While fast algorithms for computing the electrostatic potential of general distributions of charges in $\mathbb{R}^3$ exist, in a number of situations in computational physics it is useful to have a simple and extremely fast method for evaluating the potential of charges on a line; we present such a method in this paper, and report numerical results for several examples.

preprint2019arXiv

Evaluation of Abramowitz functions in the right half of the complex plane

A numerical scheme is developed for the evaluation of Abramowitz functions $J_n$ in the right half of the complex plane. For $n=-1,\, \ldots,\, 2$, the scheme utilizes series expansions for $|z|<1$ and asymptotic expansions for $|z|>R$ with $R$ determined by the required precision, and modified Laurent series expansions which are precomputed via a least squares procedure to approximate $J_n$ accurately and efficiently on each sub-region in the intermediate region $1\le |z| \le R$. For $n>2$, $J_n$ is evaluated via a recurrence relation. The scheme achieves nearly machine precision for $n=-1, \ldots, 2$, with the cost about four times of evaluating a complex exponential per function evaluation.

preprint2014arXiv

A fast multipole method for the evaluation of elastostatic fields in a half-space with zero normal stress

In this paper, we present a fast multipole method (FMM) for the half-space Green's function in a homogeneous elastic half-space subject to zero normal stress, for which an explicit solution was given by Mindlin (1936). The image structure of this Green's function is unbounded, so that standard outgoing representations are not easily available. We introduce two such representations here, one involving an expansion in plane waves and one involving a modified multipole expansion. Both play a role in the FMM implementation.

preprint2014arXiv

Randomized methods for rank-deficient linear systems

We present a simple, accurate method for solving consistent, rank-deficient linear systems, with or without addi- tional rank-completing constraints. Such problems arise in a variety of applications, such as the computation of the eigenvectors of a matrix corresponding to a known eigenvalue. The method is based on elementary linear algebra combined with the ob- servation that if the matrix is rank-k deficient, then a random rank-k perturbation yields a nonsingular matrix with probability 1.

preprint2014arXiv

The Decoupled Potential Integral Equation for Time-Harmonic Electromagnetic Scattering

We present a new formulation for the problem of electromagnetic scattering from perfect electric conductors. While our representation for the electric and magnetic fields is based on the standard vector and scalar potentials ${\bf A},ϕ$ in the Lorenz gauge, we establish boundary conditions on the potentials themselves, rather than on the field quantities. This permits the development of a well-conditioned second kind Fredholm integral equation which has no spurious resonances, avoids low frequency breakdown, and is insensitive to the genus of the scatterer. The equations for the vector and scalar potentials are decoupled. That is, the unknown scalar potential defining the scattered field, $ϕ^{Sc}$, is determined entirely by the incident scalar potential $ϕ^{In}$. Likewise, the unknown vector potential defining the scattered field, ${\bf A}^{Sc}$, is determined entirely by the incident vector potential ${\bf A}^{In}$. This decoupled formulation is valid not only in the static limit but for arbitrary $ω\ge 0$.

preprint2012arXiv

A consistency condition for the vector potential in multiply-connected domains

A classical problem in electromagnetics concerns the representation of the electric and magnetic fields in the low-frequency or static regime, where topology plays a fundamental role. For multiply connected conductors, at zero frequency the standard boundary conditions on the tangential components of the magnetic field do not uniquely determine the vector potential. We describe a (gauge-invariant) consistency condition that overcomes this non-uniqueness and resolves a longstanding difficulty in inverting the magnetic field integral equation.

preprint2011arXiv

Fast multi-particle scattering: a hybrid solver for the Maxwell equations in microstructured materials

A variety of problems in device and materials design require the rapid forward modeling of Maxwell's equations in complex micro-structured materials. By combining high-order accurate integral equation methods with classical multiple scattering theory, we have created an effective simulation tool for materials consisting of an isotropic background in which are dispersed a large number of micro- or nano-scale metallic or dielectric inclusions.