
This paper presents our work on designing a parallel platform for largescale
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 distributedmemory parallel systems, where MPI (Message Passing Interface)
is employed for communications among computation nodes. It provides structured
grid due to its simplicity and cellcentered 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. Commonlyused 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 generalpurpose 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.

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 GK\"ahlerlike or K\"ahlerlike, 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 nonK\"ahler CalabiYau 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.

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 finegrained management. In this paper, for
balancing the supplydemand and reducing the loadshedding 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 loadshedding 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 NS3Matlab
cosimulator. The numerical results demonstrate the feasibility and
effectiveness of our solution.

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.

Exotic nonAbelian quasiparticles are believed to occur in certain fractional
quantum Hall (FQH) states when effective threebody correlations form between
spinpolarized electrons in the first excited Landau level. Inspired by recent
observations of exotic physics from Floquet engineering, we investigate
periodic driving of anisotropic twobody interactions as an alternative route
for realizing robust nonAbelian multicomponent FQH states. We develop an
analytic formalism to describe this Floquet FQH protocol, which is distinct
from previous proposals that modulate singlebody hoppings for bandstructure
engineering. Our Floquet mechanism is shown to lead to highlytunable
threebody 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
nonAbelian multicomponent FQH phases. Finally, we propose a realistic
implementation of our Floquet mechanism in optically dressed ultracold polar
molecules with modulated Rabi frequencies.

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, 3DPhysNet, which is able to
predict how a threedimensional 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 realtime applications from partial views. Experimental results show the
viability and the generalisation properties of the proposed architecture.

Modelling the physical properties of everyday objects is a fundamental
prerequisite for autonomous robots. We present a novel generative adversarial
network (DefoNet), able to predict body deformations under external forces
from a single RGBD 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. DefoNet inherits the generalisation properties of GANs. This means
that the network is able to reconstruct the whole 3D 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 realtime 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 RGBD camera can use the network to predict
terrain deformations under different payload configurations and use this to
avoid unsafe areas.

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 finegrained 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 realtime
interaction of fog computing. The proposed protocol takes advantage of
attributebased 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.

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 GALFAHI 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).

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 positionpositionvelocity 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 nonperiodic 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.

A study of highorder solitons in three nonlocal nonlinear Schr\"{o}dinger
equations is presented, which includes the \PTsymmetric, reversetime, and
reversespacetime nonlocal nonlinear Schr\"{o}dinger equations. General
highorder solitons in three different equations are derived from the same
RiemannHilbert solutions of the AKNS hierarchy, except for the difference in
the corresponding symmetry relations on the "perturbed" scattering data.
Dynamics of general highorder solitons in these equations is further analyzed.
It is shown that the highorder fundamentalsoliton 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 highorder multisolitons could have more complicated wave
structures and behave very differently from highorder fundamental solitons.
More interesting is the highorder hybridpattern solitons, which are derived
from combination of different size of block matrix in the RiemannHilbert
solutions and thus they can describe a nonlinear interaction between several
types of solitons.

In this paper, we propose a novel approach, 3DRecGAN++, 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 3DRecGAN++ 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 finegrained 3D structures of
objects in highdimensional voxel space. Extensive experiments on large
synthetic datasets and realworld Kinect datasets show that the proposed
3DRecGAN++ significantly outperforms the state of the art in single view 3D
object reconstruction, and is able to reconstruct unseen types of objects.

In this work, we aim at building a bridge from poor behavioral data to an
effective, quickresponse, 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
multidimensional lowquality data, e.g., offline checkins 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 realworld datasets: Foursquare and Yelp. In the widelyused
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 realtime online identity
authentication can be improved via modeling users' composite behavioral
patterns.

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
realtime taxibookingthat there is a welldefined 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 sublinearly. 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.

In this paper, we propose a noninvasive imaging method of microwave near
field using a diamond containing nitrogenvacancy 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 nitrogenvacancy 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 highresolution near
field imaging. The proposed method is expected to have applications in
monolithicmicrowaveintegrated circuit chip local diagnosis, antenna
characterization, and field mode imaging of microwave cavities and waveguides.

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 transienttime 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 twoside
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 mixedinteger 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 timecoupled constraints
and the longtail problem, we also propose a lowcomplexity nearoptimal
algorithm based on BD (LNBD). The results verify the effectivity of the
proposed methods and algorithms.

For the fast rotating quasitwodimensional 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 twobody 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.

In this work, we study the dynamics of rogue waves in the partially
$\cal{PT}$symmetric nonlocal DaveyStewartson(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 DSI equation, the solutions are obtained and
expressed in term of determinants. For the partially $\cal{PT}$symmetric DSII
equation, the solutions are represented as quasiGram 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
finitetime 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 antidark rational travelling waves
generates the novel hybridpattern waves. However, no highorder rogue waves
are found in this partially $\cal{PT}$symmetric nonlocal DS systems. Instead,
it can produce some highorder travelling waves from the highorder rational
solutions.

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 allelectric ships.
Firstly, the multizone 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.

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 quarterpinched sectional
curvature.

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 42\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 nonK\"ahler compact Bismut flat surfaces are exactly
those isosceles Hopf surfaces by \cite{WYZ}.

We analytically derived the effective twobody interaction for a finite
thickness quantum Hall system with a harmonic perpendicular confinement and an
inplane 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 inplane 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.

In this paper, we propose a novel 3DRecGAN 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 3DRecGAN 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
finegrained 3D structures of objects in highdimensional voxel space.
Extensive experiments on large synthetic datasets show that the proposed
3DRecGAN 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/3DRecGAN.

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 nonlinear
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 realtime basis, for which the involved vehicles are required to briefly
adjust their dynamics. This leads to a selforganised 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.

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 spacefilling curve
method. A spacefilling curve is a continuous curve and defines a map between a
onedimensional space and a multidimensional space. A Hilbert spacefilling
curve is one special spacefilling 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 multidimensional space are also close to
each other in a one dimensional space. This property can model communications
in gridbased parallel applications. The idea of the Hilbert spacefilling
curve method is to map a computational domain into a onedimensional space,
partition the onedimensional 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 highdimensional coordinates to a scalar value
and an efficient onedimensional partitioning module that divides a
onedimensional space and makes sure that all intervals have a similar
workload.
The Hilbert spacefilling curve method is compared with ParMETIS, a famous
graph partitioning package. The results show that our Hilbert spacefilling
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.