74,539 research outputs found
FitSKIRT: genetic algorithms to automatically fit dusty galaxies with a Monte Carlo radiative transfer code
We present FitSKIRT, a method to efficiently fit radiative transfer models to
UV/optical images of dusty galaxies. These images have the advantage that they
have better spatial resolution compared to FIR/submm data. FitSKIRT uses the
GAlib genetic algorithm library to optimize the output of the SKIRT Monte Carlo
radiative transfer code. Genetic algorithms prove to be a valuable tool in
handling the multi- dimensional search space as well as the noise induced by
the random nature of the Monte Carlo radiative transfer code. FitSKIRT is
tested on artificial images of a simulated edge-on spiral galaxy, where we
gradually increase the number of fitted parameters. We find that we can recover
all model parameters, even if all 11 model parameters are left unconstrained.
Finally, we apply the FitSKIRT code to a V-band image of the edge-on spiral
galaxy NGC4013. This galaxy has been modeled previously by other authors using
different combinations of radiative transfer codes and optimization methods.
Given the different models and techniques and the complexity and degeneracies
in the parameter space, we find reasonable agreement between the different
models. We conclude that the FitSKIRT method allows comparison between
different models and geometries in a quantitative manner and minimizes the need
of human intervention and biasing. The high level of automation makes it an
ideal tool to use on larger sets of observed data.Comment: 14 pages, 10 figures; accepted for publication in Astronomy and
Astrophysic
Development of Three-Dimensional Neoclassical Transport Simulation Code with High Performace Fortran on a Vector-Parallel Computer
A neoclassical transport simulation code (FORTEC-3D) applicable to threedimensional configurations has been developed using High Performance Fortran (HPF). Adoption of computing techniques for parallelization and a hybrid simulation model to the delta f Monte-Carlo method transport simulation, including non-local transport effects in three-dimensional configurations, makes it possible to simulate the dynamism of global, non-local transport phenomena with a self-consistent radial electric field within a reasonable computation time. In this paper, development of the transport code using HPF is reported. Optimization techniques in order to achieve both high vectorization and parallelization efficiency, adoption of a parallel random number generator, and also benchmark results, are shown
DMC and VMC Calculations of the Electric Dipole Moment and the Ground-State Total Energy of Hydrazine Molecule Using CASINO-Code
In this study quantum Monte Carlo (QMC) calculations of the electric dipole moment and ground-state total energy of hydrazine (N2H4) molecule using CASINO-code is presented. By employing the restricted Hartree-Fock (RHF) scheme, two QMC techniques were used in this work: variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC) techniques. The optimization of the Slater-Jastrow trial wave-function was done using variance-minimization scheme. The simulations require that the configurations must evolve on the time scale of the electronic motion, and after equilibration, the estimated effective time-step be obtained. In this study, the electric dipole moment of N2H4 molecule was calculated using only the DMC technique. The result obtained gives an electric dipole moment value of 2.0 D, which is in good agreement with the experimental value of 1.85 D. Similarly, the ground-state total energy of N2H4 molecule was calculated using both VMC and DMC methods. It was observed that the result obtained from the VMC technique agrees very-well with the best theoretical value while the DMC technique gave a ground-state total energy value lower than all other theoretical values in literature, suggesting that the DMC result –111.842774 ± 0.00394 a.u., should be the exact ground-state total energy of hydrazine molecule. However, the results from this study are found to be precisely approaching the required order of chemical accuracy. Keywords: Hydrazine (N2H4), QMC, VMC, DMC, CASINO-Code.
An Efficient Monte Carlo-based Probabilistic Time-Dependent Routing Calculation Targeting a Server-Side Car Navigation System
Incorporating speed probability distribution to the computation of the route
planning in car navigation systems guarantees more accurate and precise
responses. In this paper, we propose a novel approach for dynamically selecting
the number of samples used for the Monte Carlo simulation to solve the
Probabilistic Time-Dependent Routing (PTDR) problem, thus improving the
computation efficiency. The proposed method is used to determine in a proactive
manner the number of simulations to be done to extract the travel-time
estimation for each specific request while respecting an error threshold as
output quality level. The methodology requires a reduced effort on the
application development side. We adopted an aspect-oriented programming
language (LARA) together with a flexible dynamic autotuning library (mARGOt)
respectively to instrument the code and to take tuning decisions on the number
of samples improving the execution efficiency. Experimental results demonstrate
that the proposed adaptive approach saves a large fraction of simulations
(between 36% and 81%) with respect to a static approach while considering
different traffic situations, paths and error requirements. Given the
negligible runtime overhead of the proposed approach, it results in an
execution-time speedup between 1.5x and 5.1x. This speedup is reflected at
infrastructure-level in terms of a reduction of around 36% of the computing
resources needed to support the whole navigation pipeline
- …