• ### Development of A Scalable Platform for Large-scale Reservoir Simulations on Parallel computers(1602.05901)

Sept. 1, 2018 cs.CE
This paper presents our work on designing a parallel platform for large-scale reservoir simulations. Detailed components, such as grid and linear solver, and data structures are introduced, which can serve as a guide to parallel reservoir simulations and other parallel applications. The main objective of platform is to support implementation of various parallel reservoir simulators on distributed-memory parallel systems, where MPI (Message Passing Interface) is employed for communications among computation nodes. It provides structured grid due to its simplicity and cell-centered data is applied for each cell. The platform has a distributed matrix and vector module and a map module. The matrix and vector module is the base of our parallel linear systems. The map connects grid and linear system modules, which defines various mappings between grid and linear systems. Commonly-used Krylov subspace linear solvers are implemented, including the restarted GMRES method and the BiCGSTAB method. It also has an interface to a parallel algebraic multigrid solver, BoomerAMG from HYPRE. Parallel general-purpose preconditioners and special preconditioners for reservoir simulations are also developed. Various data structures are designed, such as grid, cell, data, linear solver and preconditioner, and some key default parameters are presented in this paper. The numerical experiments show that our platform has excellent scalability and it can simulate giant reservoir models with hundreds of millions of grid cells using thousands of CPU cores.
• ### On Curvature Tensors of Hermitian Manifolds(1602.01189)

July 29, 2018 math.DG
In this article, we examine the behavior of the Riemannian and Hermitian curvature tensors of a Hermitian metric, when one of the curvature tensors obeys all the symmetry conditions of the curvature tensor of a K\"ahler metric. We will call such metrics G-K\"ahler-like or K\"ahler-like, for lack of better terminologies. Such metrics are always balanced when the manifold is compact, so in a way they are more special than balanced metrics, which drew a lot of attention in the study of non-K\"ahler Calabi-Yau manifolds. In particular we derive various formulas on the difference between the Riemannian and Hermitian curvature tensors in terms of the torsion of the Hermitian connection. We believe that these formulas could lead to further applications in the study of Hermitian geometry with curvature assumptions.
• ### Distributed Load Shedding for Microgrid with Compensation Support via Wireless Network(1701.02195)

July 23, 2018 cs.SY
Due to the limited generation and finite inertia, microgrid suffers from the large frequency and voltage deviation which can lead to system collapse. Thus, reliable load shedding to keep frequency stable is required. Wireless network, benefiting from the high flexibility and low deployment cost, is considered as a promising technology for fine-grained management. In this paper, for balancing the supply-demand and reducing the load-shedding amount, a distributed load shedding solution via wireless network is proposed. Firstly, active power coordination of different priority loads is formulated as an optimisation problem. To solve it, a distributed load shedding algorithm based on subgradient method (DLSS) is developed for gradually shedding loads. Using this method, power compensation can be utilised and has more time to lower the power deficit so as to reduce the load-shedding amount. Secondly, to increase the response rate and enhance the reliability of our method, a multicast metropolis schedule based on TDMA (MMST) is developed. In this protocol, time slots are dedicatedly allocated and a checking and retransmission mechanism is utilised. Finally, the proposed solution is evaluated by NS3-Matlab co-simulator. The numerical results demonstrate the feasibility and effectiveness of our solution.
• ### Graph Bayesian Optimization: Algorithms, Evaluations and Applications(1805.01157)

May 9, 2018 cs.AI, cs.LG, stat.ML
Network structure optimization is a fundamental task in complex network analysis. However, almost all the research on Bayesian optimization is aimed at optimizing the objective functions with vectorial inputs. In this work, we first present a flexible framework, denoted graph Bayesian optimization, to handle arbitrary graphs in the Bayesian optimization community. By combining the proposed framework with graph kernels, it can take full advantage of implicit graph structural features to supplement explicit features guessed according to the experience, such as tags of nodes and any attributes of graphs. The proposed framework can identify which features are more important during the optimization process. We apply the framework to solve four problems including two evaluations and two applications to demonstrate its efficacy and potential applications.
• ### Floquet Mechanism for Non-Abelian Fractional Quantum Hall States(1805.01446)

Exotic non-Abelian quasiparticles are believed to occur in certain fractional quantum Hall (FQH) states when effective three-body correlations form between spin-polarized electrons in the first excited Landau level. Inspired by recent observations of exotic physics from Floquet engineering, we investigate periodic driving of anisotropic two-body interactions as an alternative route for realizing robust non-Abelian multicomponent FQH states. We develop an analytic formalism to describe this Floquet FQH protocol, which is distinct from previous proposals that modulate single-body hoppings for bandstructure engineering. Our Floquet mechanism is shown to lead to highly-tunable three-body interactions that can be repulsive as well as attractive. We systematically analyze the resulting interactions with generalized pseudopotentials, and numerically demonstrate that they support a variety of non-Abelian multicomponent FQH phases. Finally, we propose a realistic implementation of our Floquet mechanism in optically dressed ultracold polar molecules with modulated Rabi frequencies.
• ### 3D-PhysNet: Learning the Intuitive Physics of Non-Rigid Object Deformations(1805.00328)

April 25, 2018 cs.CV
The ability to interact and understand the environment is a fundamental prerequisite for a wide range of applications from robotics to augmented reality. In particular, predicting how deformable objects will react to applied forces in real time is a significant challenge. This is further confounded by the fact that shape information about encountered objects in the real world is often impaired by occlusions, noise and missing regions e.g. a robot manipulating an object will only be able to observe a partial view of the entire solid. In this work we present a framework, 3D-PhysNet, which is able to predict how a three-dimensional solid will deform under an applied force using intuitive physics modelling. In particular, we propose a new method to encode the physical properties of the material and the applied force, enabling generalisation over materials. The key is to combine deep variational autoencoders with adversarial training, conditioned on the applied force and the material properties. We further propose a cascaded architecture that takes a single 2.5D depth view of the object and predicts its deformation. Training data is provided by a physics simulator. The network is fast enough to be used in real-time applications from partial views. Experimental results show the viability and the generalisation properties of the proposed architecture.
• ### Defo-Net: Learning Body Deformation using Generative Adversarial Networks(1804.05928)

April 16, 2018 cs.RO
Modelling the physical properties of everyday objects is a fundamental prerequisite for autonomous robots. We present a novel generative adversarial network (Defo-Net), able to predict body deformations under external forces from a single RGB-D image. The network is based on an invertible conditional Generative Adversarial Network (IcGAN) and is trained on a collection of different objects of interest generated by a physical finite element model simulator. Defo-Net inherits the generalisation properties of GANs. This means that the network is able to reconstruct the whole 3-D appearance of the object given a single depth view of the object and to generalise to unseen object configurations. Contrary to traditional finite element methods, our approach is fast enough to be used in real-time applications. We apply the network to the problem of safe and fast navigation of mobile robots carrying payloads over different obstacles and floor materials. Experimental results in real scenarios show how a robot equipped with an RGB-D camera can use the network to predict terrain deformations under different payload configurations and use this to avoid unsafe areas.
• ### Assured Data Deletion with Fine-grained Access Control for Fog-based Industrial Applications(1804.02834)

April 9, 2018 cs.CR
The advances of cloud computing, fog computing and Internet of Things (IoT) make the industries more prosperous than ever. A wide range of industrial systems such as transportation systems and manufacturing systems have been developed by integrating cloud computing, fog computing and IoT successfully. Security and privacy issues are a major concern that hinders the wide adoptions of these novel techniques. In this paper, we focus on assured data deletion, an issue which is important but received less attention in academia and industry. We firstly propose a framework to integrate the cloud, the fog and the things together to manage the stored data from industries or individuals. We then focus on secure data deletion in this framework by proposing an assured data deletion scheme which fulfills fine-grained access control over sensitive data and verifiable data deletion. Only the data owners and the fog devices are involved when deleting a data key and validating the data deletion, which makes the protocol practical due to the features of low latency and real-time interaction of fog computing. The proposed protocol takes advantage of attribute-based encryption and is provably secure under the standard model. The theoretical analysis shows the good performance and functionality requirements while the implementation results demonstrate the feasibility of our proposal.
• ### Dispersion of velocity gradients: Mapping magnetization with the Velocity Gradient Technique(1802.02984)

April 9, 2018 astro-ph.GA, astro-ph.IM
Recent developments of the Velocity Gradient Technique (VGT) show that the velocity gradients provide a reliable tracing of magnetic field direction in turbulent plasmas. In this paper, we explore the ability of velocity gradients to measure the magnetization of interstellar medium. We demonstrate that the distribution of velocity gradient orientations provides a reliable estimation of the magnetization of the media. In particular, we determine the relation between Alfvenic Mach number $M_A$ and properties of the velocity gradient distribution, namely, with the dispersion of velocity gradient orientation as well as with the peak to base ratio of the amplitudes. We apply our technique for a selected GALFA-HI region and find the results consistent with the expected behavior of $M_A$. Using 3D MHD simulations we successfully compare the results with our new measure of magnetization that is based on the dispersion of starlight polarization. We demonstrate that, combined with the velocity dispersion along the line of sight direction, our technique is capable to deliver the magnetic field strength. The new technique opens a way to measure magnetization using other gradient measures such as synchrotron intensity gradients (SIGs) and synchrotron polarization gradients (SPGs).
• ### Statistical tracing of magnetic fields: comparing and improving the techniques(1804.02732)

Magnetohydrodynamic(MHD) turbulence displays anisotropic velocity and density features which reflect the direction of the magnetic field. This anisotropy has led to the development of a number of statistical techniques for studying magnetic fields in the interstellar medium. In this paper, we review and compare three techniques for determining magnetic field strength and morphology that use radio position-position-velocity data: the correlation function anisotropy (CFA), Principal Component Analysis of Anisotropies (PCAA), and the more recent Velocity Gradient Technique (VGT). We compare these three techniques and suggest improvements to the CFA and PCAA techniques to increase their accuracy and versatility. In particular, we suggest and successfully implement a much faster way of calculating non-periodic correlation functions for the CFA. We discuss possible improvements to the current implementation of the PCAA. We show the advantages of the VGT in terms of magnetic field tracing and stress the complementary nature with the other two techniques.
• ### Dynamics of high-order solitons in the nonlocal nonlinear Schr\"{o}dinger equations(1802.05498)

Feb. 15, 2018 nlin.PS
A study of high-order solitons in three nonlocal nonlinear Schr\"{o}dinger equations is presented, which includes the \PT-symmetric, reverse-time, and reverse-space-time nonlocal nonlinear Schr\"{o}dinger equations. General high-order solitons in three different equations are derived from the same Riemann-Hilbert solutions of the AKNS hierarchy, except for the difference in the corresponding symmetry relations on the "perturbed" scattering data. Dynamics of general high-order solitons in these equations is further analyzed. It is shown that the high-order fundamental-soliton is always moving on several different trajectories in nearly equal velocities, and they can be nonsingular or repeatedly collapsing, depending on the choices of the parameters. It is also shown that high-order multi-solitons could have more complicated wave structures and behave very differently from high-order fundamental solitons. More interesting is the high-order hybrid-pattern solitons, which are derived from combination of different size of block matrix in the Riemann-Hilbert solutions and thus they can describe a nonlinear interaction between several types of solitons.
• ### 3D Object Dense Reconstruction from a Single Depth View(1802.00411)

Feb. 1, 2018 cs.AI, cs.CV, cs.LG, cs.RO
In this paper, we propose a novel approach, 3D-RecGAN++, which reconstructs the complete 3D structure of a given object from a single arbitrary depth view using generative adversarial networks. Unlike existing work which typically requires multiple views of the same object or class labels to recover the full 3D geometry, the proposed 3D-RecGAN++ only takes the voxel grid representation of a depth view of the object as input, and is able to generate the complete 3D occupancy grid with a high resolution of 256^3 by recovering the occluded/missing regions. The key idea is to combine the generative capabilities of autoencoders and the conditional Generative Adversarial Networks (GAN) framework, to infer accurate and fine-grained 3D structures of objects in high-dimensional voxel space. Extensive experiments on large synthetic datasets and real-world Kinect datasets show that the proposed 3D-RecGAN++ significantly outperforms the state of the art in single view 3D object reconstruction, and is able to reconstruct unseen types of objects.
• ### Composite Behavioral Modeling for Identity Theft Detection in Online Social Networks(1801.06825)

Jan. 21, 2018 cs.CR, cs.SI
In this work, we aim at building a bridge from poor behavioral data to an effective, quick-response, and robust behavior model for online identity theft detection. We concentrate on this issue in online social networks (OSNs) where users usually have composite behavioral records, consisting of multi-dimensional low-quality data, e.g., offline check-ins and online user generated content (UGC). As an insightful result, we find that there is a complementary effect among different dimensions of records for modeling users' behavioral patterns. To deeply exploit such a complementary effect, we propose a joint model to capture both online and offline features of a user's composite behavior. We evaluate the proposed joint model by comparing with some typical models on two real-world datasets: Foursquare and Yelp. In the widely-used setting of theft simulation (simulating thefts via behavioral replacement), the experimental results show that our model outperforms the existing ones, with the AUC values $0.956$ in Foursquare and $0.947$ in Yelp, respectively. Particularly, the recall (True Positive Rate) can reach up to $65.3\%$ in Foursquare and $72.2\%$ in Yelp with the corresponding disturbance rate (False Positive Rate) below $1\%$. It is worth mentioning that these performances can be achieved by examining only one composite behavior (visiting a place and posting a tip online simultaneously) per authentication, which guarantees the low response latency of our method. This study would give the cybersecurity community new insights into whether and how a real-time online identity authentication can be improved via modeling users' composite behavioral patterns.
• ### Phase Transition in Taxi Dynamics and Impact of Ridesharing(1801.00462)

Jan. 16, 2018 physics.soc-ph, nlin.AO
We develop a numerical model using both artificial and empirical inputs to analyze taxi dynamics in an urban setting. More specifically, we quantify how the supply and demand for taxi services, the underlying road network, and the public acceptance of taxi ridesharing (TRS) affect the optimal number of taxis for a particular city, as well as commuters' average waiting time and trip time. Results reveal certain universal features of the taxi dynamics with real-time taxi-booking---that there is a well-defined transition between the oversaturated phase when demand exceeds supply, and the undersaturated phase when supply exceeds demand. The boundary between the two phases gives the optimal number of taxis a city should accommodate, given the specific demand, road network and commuter habits. Adding or removing taxis may affect commuter experience very differently in the two phases revealed. In the oversaturated phase the average waiting time is affected exponentially, while in the undersaturated phase it is affected sub-linearly. We analyze various factors that can shift the phase boundary, and show that an increased level of acceptance for TRS universally shifts the phase boundary by reducing the number of taxis needed. We discuss some of the useful insights on the benefits and costs of TRS, especially how under certain situations TRS will not only have economic benefits for commuters, but can also save the overall travel time for the shared parties, by significantly reducing the time commuters spend on waiting for taxis. Our simulations also suggest that simple artificial taxi systems can capture most of the universal features of the taxi dynamics. The relevance of the assumptions and the overall methodology are also illustrated using the empirical road network and taxi demand in Singapore.
• ### Non-Invasive Imaging Method of Microwave Near Field Based on Solid State Quantum Sensing(1801.01706)

In this paper, we propose a non-invasive imaging method of microwave near field using a diamond containing nitrogen-vacancy centers. We applied synchronous pulsed sequence combined with charge coupled device camera to measure the amplitude of the microwave magnetic field. A full reconstruction formulation of the local field vector, including the amplitude and phase, is developed by measuring both left and right circular polarizations along the four nitrogen-vacancy axes. Compared to the raster scanning approach, the two dimensional imaging method is promising for application to circuit failure analysis. A diamond film with micrometer thinness enables high-resolution near field imaging. The proposed method is expected to have applications in monolithic-microwave-integrated circuit chip local diagnosis, antenna characterization, and field mode imaging of microwave cavities and waveguides.
• ### Optimal Power Management for Failure Mode of AC/DC Microgrids in All-Electric Ships(1712.02552)

Jan. 2, 2018 cs.SY
Optimal power management of shipboard power system for failure mode (OPMSF) is a significant and challenging problem considering the safety of system and person. Many existing works focused on the transient-time recovery without consideration of the operating cost and the voyage plan. In this paper, we formulate the OPMSF problem considering the long- time scheduling and the faults at bus and generator. For reducing the fault effects, we adopt two-side adjustment methods including the load shedding and the reconfiguration. To address this non- convex problem, we transform the travel equality constraint into an inequality constraint and obtain a convex problem. Then, considering the infeasibility scenario affected by faults, a further relaxation is adopted to obtain a new problem with feasibility guaranteed. Furthermore, we derive a sufficient condition to guarantee that the new problem has the same optimal solution as the original one. Because of the mixed-integer nonlinear feature, we develop an optimal algorithm based on Benders decomposition (BD) to solve the new one. Due to the slow convergence caused by the time-coupled constraints and the long-tail problem, we also propose a low-complexity near-optimal algorithm based on BD (LNBD). The results verify the effectivity of the proposed methods and algorithms.
• ### Phase transition and intrinsic metric of the dipolar fermions in quantum Hall regime(1710.05597)

For the fast rotating quasi-two-dimensional dipolar fermions in the quantum Hall regime, the interaction between two dipoles breaks the rotational symmetry when the dipole moment has component in the the plane via being tuned by an external field. For the anisotropic two-body interaction, we expand it in a generalized pseudopotentials (PPs). With assuming that all the dipoles are polarized in the same direction, we perform the numerical diagonalization for finite size systems on a torus. We find that the most stable fractional quantum Hall (FQH) states in the lowest Landau level (LLL) and the first Landau level (1LL) are {\nu} = 1/3 and {\nu} = 2 + 1/5 Laughlin state respectively in the isotropic case. While rotating the dipolar angle, these FQH states reveal a robustness and finally enter into a molecule phase in which all the particles are attracted and form a bound state. The anisotropy and the phase transition are studied by the intrinsic metric, the wave function overlap and the nematic order parameter.
• ### Dynamics of Rogue Waves in the Partially PT-symmetric Nonlocal Davey-Stewartson Systems(1710.07061)

Oct. 19, 2017 math-ph, math.MP, nlin.SI
In this work, we study the dynamics of rogue waves in the partially $\cal{PT}$-symmetric nonlocal Davey-Stewartson(DS) systems. Using the Darboux transformation method, general rogue waves in the partially $\cal{PT}$-symmetric nonlocal DS equations are derived. For the partially $\cal{PT}$-symmetric nonlocal DS-I equation, the solutions are obtained and expressed in term of determinants. For the partially $\cal{PT}$-symmetric DS-II equation, the solutions are represented as quasi-Gram determinants. It is shown that the fundamental rogue waves in these two systems are rational solutions which arises from a constant background at $t\rightarrow -\infty$, and develops finite-time singularity on an entire hyperbola in the spatial plane at the critical time. It is also shown that the interaction of several fundamental rogue waves is described by the multi rogue waves. And the interaction of fundamental rogue waves with dark and anti-dark rational travelling waves generates the novel hybrid-pattern waves. However, no high-order rogue waves are found in this partially $\cal{PT}$-symmetric nonlocal DS systems. Instead, it can produce some high-order travelling waves from the high-order rational solutions.
• ### Hybrid Optimization Method for Reconfiguration of AC/DC Microgrids in All-Electric Ships(1709.05505)

Sept. 19, 2017 cs.SY
Since the limited power capacity, finite inertia, and dynamic loads make the shipboard power system (SPS) vulnerable, the automatic reconfiguration for failure recovery in SPS is an extremely significant but still challenging problem. It is not only required to operate accurately and optimally, but also to satisfy operating constraints. In this paper, we consider the reconfiguration optimization for hybrid AC/DC microgrids in all-electric ships. Firstly, the multi-zone medium voltage DC (MVDC) SPS model is presented. In this model, the DC power flow for reconfiguration and a generalized AC/DC converter are modeled for accurate reconfiguration. Secondly, since this problem is mixed integer nonlinear programming (MINLP), a hybrid method based on Newton Raphson and Biogeography based Optimization (NRBBO) is designed according to the characteristics of system, loads, and faults. This method facilitates to maximize the weighted load restoration while satisfying operating constraints. Finally, the simulation results demonstrate this method has advantages in terms of power restoration and convergence speed.
• ### A note on the almost one half holomorphic pinching(1709.02527)

Sept. 8, 2017 math.DG
Motivated by a previous work of Zheng and the second named author, we study pinching constants of compact K\"ahler manifolds with positive holomorphic sectional curvature. In particular we prove a gap theorem following the work of Petersen and Tao on Riemannian manifolds with almost quarter-pinched sectional curvature.
• ### On compact Hermitian manifolds with flat Gauduchon connections(1709.02530)

Sept. 8, 2017 math.DG
Given a Hermitian manifold $(M^n,g)$, the Gauduchon connections are the one parameter family of Hermitian connections joining the Chern connection and the Bismut connection. We will call $\nabla^s = (1-\frac{s}{2})\nabla^c + \frac{s}{2}\nabla^b$ the $s$-Gauduchon connection of $M$, where $\nabla^c$ and $\nabla^b$ are respectively the Chern and Bismut connections. It is natural to ask when a compact Hermitian manifold could admit a flat $s$-Gauduchon connection. This is related to a question asked by Yau \cite{Yau}. The cases with $s=0$ (a flat Chern connection) or $s=2$ (a flat Bismut connection) are classified respectively by Boothby \cite{Boothby} in the 1950s or by Q. Wang and the authors recently \cite{WYZ}. In this article, we observe that if either $s\geq 4+2\sqrt{3} \approx 7.46$ or $s\leq 4-2\sqrt{3}\approx 0.54$ and $s\neq 0$, then $g$ is K\"ahler. We also show that, when $n=2$, $g$ is always K\"ahler unless $s=2$. Note that non-K\"ahler compact Bismut flat surfaces are exactly those isosceles Hopf surfaces by \cite{WYZ}.
• ### Anisotropic Pseudopotential Characterization of Quantum Hall Systems under Tilted Magnetic Field(1709.00286)

We analytically derived the effective two-body interaction for a finite thickness quantum Hall system with a harmonic perpendicular confinement and an in-plane magnetic field. The anisotropic effective interaction in the lowest Landau level (LLL) and first Landau level (1LL) are expanded in the basis of the generalized pseudopotentials (PPs), and we analyze how the coefficients of some prominent isotropic and anisotropic PPs depend on the thickness of the sample and the strength of the in-plane magnetic field. We also investigate the stability of the topological quantum Hall states, especially the Laughlin state and its emergent guiding center metric, which we can now compute analytically. An interesting reorientation of the anisotropy direction of the Laughlin state in the 1LL is revealed, and we also discuss various possible experimental ramifications for this quantum Hall system with broken rotational symmetry.
• ### 3D Object Reconstruction from a Single Depth View with Adversarial Learning(1708.07969)

Aug. 26, 2017 cs.AI, cs.CV, cs.LG, cs.RO
In this paper, we propose a novel 3D-RecGAN approach, which reconstructs the complete 3D structure of a given object from a single arbitrary depth view using generative adversarial networks. Unlike the existing work which typically requires multiple views of the same object or class labels to recover the full 3D geometry, the proposed 3D-RecGAN only takes the voxel grid representation of a depth view of the object as input, and is able to generate the complete 3D occupancy grid by filling in the occluded/missing regions. The key idea is to combine the generative capabilities of autoencoders and the conditional Generative Adversarial Networks (GAN) framework, to infer accurate and fine-grained 3D structures of objects in high-dimensional voxel space. Extensive experiments on large synthetic datasets show that the proposed 3D-RecGAN significantly outperforms the state of the art in single view 3D object reconstruction, and is able to reconstruct unseen types of objects. Our code and data are available at: https://github.com/Yang7879/3D-RecGAN.
• ### Efficient Intersection Control for Minimally Guided Vehicles: A Self-Organised and Decentralized Approach(1708.04553)

Aug. 15, 2017 nlin.AO
An important question for the practical applicability of the highly efficient traffic intersection control is about the minimal level of intelligence the vehicles need to have so as to move beyond the traffic light control. We propose an efficient intersection traffic control scheme without the traffic lights, that only requires a majority of vehicles on the road to be equipped with a simple driver assistance system. The algorithm of our scheme is completely decentralized, and takes into full account the non-linear interaction between the vehicles at high density. For vehicles approaching the intersection in different directions, our algorithm imposes simple interactions between vehicles around the intersection, by defining specific conditions on the real-time basis, for which the involved vehicles are required to briefly adjust their dynamics. This leads to a self-organised traffic flow that is safe, robust, and efficient. We also take into account of the driver comfort level and study its effect on the control efficiency. The scheme has low technological barrier, minimal impact on the conventional driving behaviour, and can coexist with the traffic light control. It also has the advantages of being easily scalable, and fully compatible with both the conventional road systems as well as the futuristic scenario in which driverless vehicles dominate the road. The mathematical formulation of our scheme permits large scale realistic numerical simulations of busy intersections, allowing a more complete evaluation of the control performance, instead of just the collision avoidance at the intersection.
• ### Load Balancing using Hilbert Space-filling Curves for Parallel Reservoir Simulations(1708.01365)

Aug. 4, 2017 cs.DC
The goal of load balancing (grid partitioning) is to minimize overall computations and communications, and to make sure that all processors have a similar workload. Geometric methods divide a grid by using a location of a cell while topological methods work with connectivity of cells, which is generally described as a graph. This paper introduces a Hilbert space-filling curve method. A space-filling curve is a continuous curve and defines a map between a one-dimensional space and a multi-dimensional space. A Hilbert space-filling curve is one special space-filling curve discovered by Hilbert and has many useful characteristics, such as good locality, which means that two objects that are close to each other in a multi-dimensional space are also close to each other in a one dimensional space. This property can model communications in grid-based parallel applications. The idea of the Hilbert space-filling curve method is to map a computational domain into a one-dimensional space, partition the one-dimensional space to certain intervals, and assign all cells in a same interval to a MPI. To implement a load balancing method, a mapping kernel is required to convert high-dimensional coordinates to a scalar value and an efficient one-dimensional partitioning module that divides a one-dimensional space and makes sure that all intervals have a similar workload. The Hilbert space-filling curve method is compared with ParMETIS, a famous graph partitioning package. The results show that our Hilbert space-filling curve method has good partition quality. It has been applied to grids with billions of cells, and linear scalability has been obtained on IBM Blue Gene/Q.