Source author record

Yukihiro Komura

Yukihiro Komura 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

9works
2topics
2close 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

9 published item(s)

preprint2016arXiv

A generalized GPU-based connected component labeling algorithm

We propose a generalized GPU-based connected component labeling (CCL) algorithm that can be applied to both various lattices and to non-lattice environments in a uniform fashion. We extend our recent GPU-based CCL algorithm without the use of conventional iteration to the generalized method. As an application of this algorithm, we deal with the bond percolation problem. We investigate bond percolation on the honeycomb and triangle lattices to confirm the correctness of this algorithm. Moreover, we deal with bond percolation on the Bethe lattice as a substitute for a network structure, and demonstrate the performance of this algorithm on those lattices.

preprint2015arXiv

Probing phase transition order of $q$-state Potts models using Wang-Landau Algorithm

Phase transitions are ubiquitous phenomena, exemplified by the melting of ice and spontaneous magnetization of magnetic material. In general, a phase transition is associated with a symmetry breaking of a system; occurs due to the competition between coupling interaction and external fields such as thermal energy. If the phase transition occurs with no latent heat, the system experiences continuous transition, also known as second order phase transition. The ferromagnetic $q$-state Potts model with $r$ extra invisible states, introduced by Tamura, Tanaka, and Kawashima [Prog. Theor. Phys. 124, 381 (2010)], is studied by using the Wang-Landau method. The density of states difference (DOSD), $\ln g(E +ΔE) - \ln g(E)$, is used to investigate the order of the phase transition and examine the critical value of $r$ changing the second to the first order transition.

preprint2014arXiv

CUDA programs for GPU computing of Swendsen-Wang multi-cluster spin flip algorithm: 2D and 3D Ising, Potts, and XY models

We present sample CUDA programs for the GPU computing of the Swendsen-Wang multi-cluster spin flip algorithm. We deal with the classical spin models; the Ising model, the $q$-state Potts model, and the classical XY model. As for the lattice, both the 2D (square) lattice and the 3D (simple cubic) lattice are treated. We already reported the idea of the GPU implementation for 2D models [Comput. Phys. Commun. 183 (2012) 1155-1161]. We here explain the details of sample programs, and discuss the performance of the present GPU implementation for the 3D Ising and XY models. We also show the calculated results of the moment ratio for these models, and discuss phase transitions.

preprint2012arXiv

Difference of energy density of states in the Wang-Landau algorithm

Paying attention to the difference of density of states, Δln g(E) = ln g(E+ΔE) - ln g(E), we study the convergence of the Wang-Landau method. We show that this quantity is a good estimator to discuss the errors of convergence, and refer to the $1/t$ algorithm. We also examine the behavior of the 1st-order transition with this difference of density of states in connection with Maxwell's equal area rule. A general procedure to judge the order of transition is given.

preprint2012arXiv

GPU-based single-cluster algorithm for the simulation of the Ising model

We present the GPU calculation with the common unified device architecture (CUDA) for the Wolff single-cluster algorithm of the Ising model. Proposing an algorithm for a quasi-block synchronization, we realize the Wolff single-cluster Monte Carlo simulation with CUDA. We perform parallel computations for the newly added spins in the growing cluster. As a result, the GPU calculation speed for the two-dimensional Ising model at the critical temperature with the linear size L=4096 is 5.60 times as fast as the calculation speed on a current CPU core. For the three-dimensional Ising model with the linear size L=256, the GPU calculation speed is 7.90 times as fast as the CPU calculation speed. The idea of quasi-block synchronization can be used not only in the cluster algorithm but also in many fields where the synchronization of all threads is required.

preprint2012arXiv

GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional classical spin systems

We present the GPU calculation with the common unified device architecture (CUDA) for the Swendsen-Wang multi-cluster algorithm of two-dimensional classical spin systems. We adjust the two connected component labeling algorithms recently proposed with CUDA for the assignment of the cluster in the Swendsen-Wang algorithm. Starting with the q-state Potts model, we extend our implementation to the system of vector spins, the q-state clock model, with the idea of embedded cluster. We test the performance, and the calculation time on GTX580 is obtained as 2.51 nano sec per a spin flip for the q=2 Potts model (Ising model) and 2.42 nano sec per a spin flip for the q=6 clock model with the linear size L=4096 at the critical temperature, respectively. The computational speed for the q=2 Potts model on GTX580 is 12.4 times as fast as the calculation speed on a current CPU core. That for the q=6 clock model on GTX580 is 35.6 times as fast as the calculation speed on a current CPU core.

preprint2012arXiv

Large-scale Monte Carlo simulation of two-dimensional classical XY model using multiple GPUs

We study the two-dimensional classical XY model by the large-scale Monte Carlo simulation of the Swendsen-Wang multi-cluster algorithm using multiple GPUs on the open science supercomputer TSUBAME 2.0. Simulating systems up to the linear system size L=65536, we investigate the Kosterlitz-Thouless (KT) transition. Using the generalized version of the probability-changing cluster algorithm based on the helicity modulus, we locate the KT transition temperature in a self-adapted way. The obtained inverse KT temperature β_{KT} is 1.11996(6). We estimate the exponent to specify the multiplicative logarithmic correction, -2r, and precisely reproduce the theoretical prediction -2r=1/8.

preprint2012arXiv

Multi-GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional q-state Potts model

We present the multiple GPU computing with the common unified device architecture (CUDA) for the Swendsen-Wang multi-cluster algorithm of two-dimensional (2D) q-state Potts model. Extending our algorithm for single GPU computing [Comp. Phys. Comm. 183 (2012) 1155], we realize the GPU computation of the Swendsen-Wang multi-cluster algorithm for multiple GPUs. We implement our code on the large-scale open science supercomputer TSUBAME 2.0, and test the performance and the scalability of the simulation of the 2D Potts model. The performance on Tesla M2050 using 256 GPUs is obtained as 37.3 spin flips per a nano second for the q=2 Potts model (Ising model) at the critical temperature with the linear system size L=65536.

preprint2010arXiv

Phase transition of two-dimensional generalized XY model

We study the two-dimensional generalized XY model that depends on an integer $q$ by the Monte Carlo method. This model was recently proposed by Romano and Zagrebnov. We find a single Kosterlitz-Thouless (KT) transition for all values of $q$, in contrast with the previous speculation that there may be two transitions, one a regular KT transition and another a first-order transition at a higher temperature. We show the universality of the KT transitions by comparing the universal finite-size scaling behaviors at different values of $q$ without assuming a specific universal form in terms of the KT transition temperature $T_{\rm KT}$.