Search results for "physics.comp-ph"
showing 10 items of 115 documents
Unsupervised Anomaly and Change Detection With Multivariate Gaussianization
2022
Anomaly detection (AD) is a field of intense research in remote sensing (RS) image processing. Identifying low probability events in RS images is a challenging problem given the high dimensionality of the data, especially when no (or little) information about the anomaly is available a priori. While a plenty of methods are available, the vast majority of them do not scale well to large datasets and require the choice of some (very often critical) hyperparameters. Therefore, unsupervised and computationally efficient detection methods become strictly necessary, especially now with the data deluge problem. In this article, we propose an unsupervised method for detecting anomalies and changes …
Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model
2010
A Modern Graphics Processing unit (GPU) is able to perform massively parallel scientific computations at low cost. We extend our implementation of the checkerboard algorithm for the two-dimensional Ising model [T. Preis et al., Journal of Chemical Physics 228 (2009) 4468–4477] in order to overcome the memory limitations of a single GPU which enables us to simulate significantly larger systems. Using multi-spin coding techniques, we are able to accelerate simulations on a single GPU by factors up to 35 compared to an optimized single Central Processor Unit (CPU) core implementation which employs multi-spin coding. By combining the Compute Unified Device Architecture (CUDA) with the Message P…
Local Granger causality
2021
Granger causality is a statistical notion of causal influence based on prediction via vector autoregression. For Gaussian variables it is equivalent to transfer entropy, an information-theoretic measure of time-directed information transfer between jointly dependent processes. We exploit such equivalence and calculate exactly the 'local Granger causality', i.e. the profile of the information transfer at each discrete time point in Gaussian processes; in this frame Granger causality is the average of its local version. Our approach offers a robust and computationally fast method to follow the information transfer along the time history of linear stochastic processes, as well as of nonlinear …
Real-time computation of parameter fitting and image reconstruction using graphical processing units
2016
Abstract In recent years graphical processing units (GPUs) have become a powerful tool in scientific computing. Their potential to speed up highly parallel applications brings the power of high performance computing to a wider range of users. However, programming these devices and integrating their use in existing applications is still a challenging task. In this paper we examined the potential of GPUs for two different applications. The first application, created at Paul Scherrer Institut (PSI), is used for parameter fitting during data analysis of μ SR (muon spin rotation, relaxation and resonance) experiments. The second application, developed at ETH, is used for PET (Positron Emission T…
Point Group Symmetry Analysis of the Electronic Structure of Bare and Protected Metal Nanocrystals
2018
The electronic structures of a variety of experimentally identified gold and silver nanoclusters from 20 to 246 atoms, either unprotected or protected by several types of ligands, are characterized by using point group specific symmetry analysis. The delocalized electron states around the HOMO-LUMO energy gap, originating from the metal s-electrons in the cluster core, show symmetry characteristics according to the point group that describes best the atomic arrangement of the core. This indicates strong effects of the lattice structure and overall shape of the metal core to the electronic structure, which cannot be captured by the conventional analysis based on identification of spherical a…
Boundary-integral method for poloidal axisymmetric AC magnetic fields
2006
This paper presents a boundary-integral equation (BIE) method for the calculation of poloidal axisymmetric magnetic fields applicable in a wide range of ac frequencies. The method is based on the vector potential formulation and it uses the Green's functions of Laplace and Helmholtz equations for the exterior and interior of conductors, respectively. The work is particularly focused on a calculation of axisymmetric Green's function for the Helmholtz equation which is both simpler and more accurate compared to previous approaches. Three different approaches are used for calculation of the Green's function depending on the parameter range. For low and high dimensionless ac frequencies we use …
Robustness of Coherence: An Operational and Observable Measure of Quantum Coherence
2016
Quantifying coherence is an essential endeavour for both quantum foundations and quantum technologies. Here the robustness of coherence is defined and proven a full monotone in the context of the recently introduced resource theories of quantum coherence. The measure is shown to be observable, as it can be recast as the expectation value of a coherence witness operator for any quantum state. The robustness of coherence is evaluated analytically on relevant classes of states, and an efficient semidefinite program that computes it on general states is given. An operational interpretation is finally provided: the robustness of coherence quantifies the advantage enabled by a quantum state in a …
A numerical method for imaging of biological microstructures by VHF waves
2013
Imaging techniques give a fundamental support to medical diagnostics during the pathology discovery as well as for the characterization of biological structures. The imaging methods involve electromagnetic waves in a frequency range that spans from some Hz to GHz and over. Most of these methods involve scanning of wide human body areas even if only small areas need to be analyzed. In this paper, a numerical method to evaluate the shape of micro-structures for application in the medical field, with a very low invasiveness for the human body, is proposed. A flexible thin-wire antenna radiates the VHF waves and then, by measuring the spatial magnetic field distribution it is possible to recons…
Efficient and portable acceleration of quantum chemical many-body methods in mixed floating point precision using OpenACC compiler directives
2016
It is demonstrated how the non-proprietary OpenACC standard of compiler directives may be used to compactly and efficiently accelerate the rate-determining steps of two of the most routinely applied many-body methods of electronic structure theory, namely the second-order M{\o}ller-Plesset (MP2) model in its resolution-of-the-identity (RI) approximated form and the (T) triples correction to the coupled cluster singles and doubles model (CCSD(T)). By means of compute directives as well as the use of optimized device math libraries, the operations involved in the energy kernels have been ported to graphics processing unit (GPU) accelerators, and the associated data transfers correspondingly o…
Quantum dynamics by the constrained adiabatic trajectory method
2011
We develop the constrained adiabatic trajectory method (CATM) which allows one to solve the time-dependent Schr\"odinger equation constraining the dynamics to a single Floquet eigenstate, as if it were adiabatic. This constrained Floquet state (CFS) is determined from the Hamiltonian modified by an artificial time-dependent absorbing potential whose forms are derived according to the initial conditions. The main advantage of this technique for practical implementation is that the CFS is easy to determine even for large systems since its corresponding eigenvalue is well isolated from the others through its imaginary part. The properties and limitations of the CATM are explored through simple…