Catalog footprint

What is connected

49works
34topics
4close 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

49 published item(s)

preprint2026arXiv

Ada-MK: Adaptive MegaKernel Optimization via Automated DAG-based Search for LLM Inference

When large language models (LLMs) serve real-time inference in commercial online advertising systems, end-to-end latency must be strictly bounded to the millisecond range. Yet every token generated during the decode phase triggers thousands of kernel launches, and kernel launch overhead alone can account for 14.6% of end-to-end inference time. MegaKernel eliminates launch overhead and inter-operator HBM round-trips by fusing multiple operators into a single persistent kernel. However, existing MegaKernel implementations face a fundamental tension between portability and efficiency on resource-constrained GPUs such as NVIDIA Ada: hand-tuned solutions are tightly coupled to specific architectures and lack portability, while auto-compiled approaches introduce runtime dynamic scheduling whose branch penalties are unacceptable in latency-critical settings. We observe that under a fixed deployment configuration, the optimal execution path of a MegaKernel is uniquely determined, and runtime dynamic decision-making can be entirely hoisted to compile time. Building on this insight, we propose Ada-MK: (1) a three-dimensional shared-memory constraint model combined with K-dimension splitting that reduces peak shared memory usage by 50%; (2) MLIR-based fine-grained DAG offline search that solidifies the optimal execution path, completely eliminating runtime branching; and (3) a heterogeneous hybrid inference engine that embeds MegaKernel as a plugin into TensorRT-LLM, combining high-throughput Prefill with low-latency Decode. On an NVIDIA L20, Ada-MK improves single-batch throughput by up to 23.6% over vanilla TensorRT-LLM and 50.2% over vLLM, achieving positive gains across all tested scenarios--the first industrial deployment of MegaKernel in a commercial online advertising system.

preprint2022arXiv

A Large-scale Comprehensive Dataset and Copy-overlap Aware Evaluation Protocol for Segment-level Video Copy Detection

In this paper, we introduce VCSL (Video Copy Segment Localization), a new comprehensive segment-level annotated video copy dataset. Compared with existing copy detection datasets restricted by either video-level annotation or small-scale, VCSL not only has two orders of magnitude more segment-level labelled data, with 160k realistic video copy pairs containing more than 280k localized copied segment pairs, but also covers a variety of video categories and a wide range of video duration. All the copied segments inside each collected video pair are manually extracted and accompanied by precisely annotated starting and ending timestamps. Alongside the dataset, we also propose a novel evaluation protocol that better measures the prediction accuracy of copy overlapping segments between a video pair and shows improved adaptability in different scenarios. By benchmarking several baseline and state-of-the-art segment-level video copy detection methods with the proposed dataset and evaluation metric, we provide a comprehensive analysis that uncovers the strengths and weaknesses of current approaches, hoping to open up promising directions for future works. The VCSL dataset, metric and benchmark codes are all publicly available at https://github.com/alipay/VCSL.

preprint2022arXiv

CATNet: Cross-event Attention-based Time-aware Network for Medical Event Prediction

Medical event prediction (MEP) is a fundamental task in the medical domain, which needs to predict medical events, including medications, diagnosis codes, laboratory tests, procedures, outcomes, and so on, according to historical medical records. The task is challenging as medical data is a type of complex time series data with heterogeneous and temporal irregular characteristics. Many machine learning methods that consider the two characteristics have been proposed for medical event prediction. However, most of them consider the two characteristics separately and ignore the correlations among different types of medical events, especially relations between historical medical events and target medical events. In this paper, we propose a novel neural network based on attention mechanism, called cross-event attention-based time-aware network (CATNet), for medical event prediction. It is a time-aware, event-aware and task-adaptive method with the following advantages: 1) modeling heterogeneous information and temporal information in a unified way and considering temporal irregular characteristics locally and globally respectively, 2) taking full advantage of correlations among different types of events via cross-event attention. Experiments on two public datasets (MIMIC-III and eICU) show CATNet can be adaptive with different MEP tasks and outperforms other state-of-the-art methods on various MEP tasks. The source code of CATNet will be released after this manuscript is accepted.

preprint2022arXiv

Cauchy, normal and correlations versus heavy tails

A surprising result of Pillai and Meng (2016) showed that a transformation $\sum_{j=1}^n w_j X_j/Y_j$ of two iid centered normal random vectors, $(X_1,\ldots, X_n)$ and $(Y_1,\ldots, Y_n)$, $n>1$, for any weights $0\leq w_j\leq 1$, $ j=1,\ldots, n$, $\sum_{j=1}^n w_j=1$, has a Cauchy distribution regardless of any correlations within the normal vectors. The correlations appear to lose out in the competition with the heavy tails. To clarify how extensive this phenomenon is, we analyze two other transformations of two iid centered normal random vectors. These transformations are similar in spirit to the transformation considered by Pillai and Meng (2016). One transformation involves absolute values: $\sum_{j=1}^n w_j X_j/|Y_j|$. The second involves randomly stopped Brownian motions: $\sum_{j=1}^n w_j X_j\bigl(Y_j^{-2}\bigr)$, where $\bigl\{\bigl( X_1(t),\ldots, X_n(t)\bigr), \, t\geq 0\bigr\},\ n>1,$ is a Brownian motion with positive variances; $(Y_1,\ldots, Y_n)$ is a centered normal random vector with the same law as $( X_1(1),\ldots, X_n(1))$ and independent of it; and $X(Y^{-2})$ is the value of the Brownian motion $X(t)$ evaluated at the random time $t=Y^{-2}$. All three transformations result in a Cauchy distribution if the covariance matrix of the normal components is diagonal, or if all the correlations implied by the covariance matrix equal 1. However, while the transformation Pillai and Meng (2016) considered produces a Cauchy distribution regardless of the normal covariance matrix. the transformations we consider here do not always produce a Cauchy distribution. The correlations between jointly normal random variables are not always overwhelmed by the heaviness of the marginal tails. The mysteries of the connections between normal and Cauchy laws remain to be understood.

preprint2022arXiv

DeepShovel: An Online Collaborative Platform for Data Extraction in Geoscience Literature with AI Assistance

Geoscientists, as well as researchers in many fields, need to read a huge amount of literature to locate, extract, and aggregate relevant results and data to enable future research or to build a scientific database, but there is no existing system to support this use case well. In this paper, based on the findings of a formative study about how geoscientists collaboratively annotate literature and extract and aggregate data, we proposed DeepShovel, a publicly-available AI-assisted data extraction system to support their needs. DeepShovel leverages the state-of-the-art neural network models to support researcher(s) easily and accurately annotate papers (in the PDF format) and extract data from tables, figures, maps, etc. in a human-AI collaboration manner. A follow-up user evaluation with 14 researchers suggested DeepShovel improved users' efficiency of data extraction for building scientific databases, and encouraged teams to form a larger scale but more tightly-coupled collaboration.

preprint2022arXiv

Diverse Preference Augmentation with Multiple Domains for Cold-start Recommendations

Cold-start issues have been more and more challenging for providing accurate recommendations with the fast increase of users and items. Most existing approaches attempt to solve the intractable problems via content-aware recommendations based on auxiliary information and/or cross-domain recommendations with transfer learning. Their performances are often constrained by the extremely sparse user-item interactions, unavailable side information, or very limited domain-shared users. Recently, meta-learners with meta-augmentation by adding noises to labels have been proven to be effective to avoid overfitting and shown good performance on new tasks. Motivated by the idea of meta-augmentation, in this paper, by treating a user's preference over items as a task, we propose a so-called Diverse Preference Augmentation framework with multiple source domains based on meta-learning (referred to as MetaDPA) to i) generate diverse ratings in a new domain of interest (known as target domain) to handle overfitting on the case of sparse interactions, and to ii) learn a preference model in the target domain via a meta-learning scheme to alleviate cold-start issues. Specifically, we first conduct multi-source domain adaptation by dual conditional variational autoencoders and impose a Multi-domain InfoMax (MDI) constraint on the latent representations to learn domain-shared and domain-specific preference properties. To avoid overfitting, we add a Mutually-Exclusive (ME) constraint on the output of decoders to generate diverse ratings given content data. Finally, these generated diverse ratings and the original ratings are introduced into the meta-training procedure to learn a preference meta-learner, which produces good generalization ability on cold-start recommendation tasks. Experiments on real-world datasets show our proposed MetaDPA clearly outperforms the current state-of-the-art baselines.

preprint2022arXiv

Empirical Bayes Multistage Testing for Large-Scale Experiments

Modern application of A/B tests is challenging due to its large scale in various dimensions, which demands flexibility to deal with multiple testing sequentially. The state-of-the-art practice first reduces the observed data stream to always-valid p-values, and then chooses a cut-off as in conventional multiple testing schemes. Here we propose an alternative method called AMSET (adaptive multistage empirical Bayes test) by incorporating historical data in decision-making to achieve efficiency gains while retaining marginal false discovery rate (mFDR) control that is immune to peeking. We also show that a fully data-driven estimation in AMSET performs robustly to various simulation and real data settings at a large mobile app social network company.

preprint2022arXiv

foREST: A Tree-based Approach for Fuzzing RESTful APIs

Representational state transfer (REST) is a widely employed architecture by web applications and cloud. Users can invoke such services according to the specification of their application interfaces, namely RESTful APIs. Existing approaches for fuzzing RESTful APIs are generally based on classic API-dependency graphs. However, such dependencies are inefficient for REST services due to the explosion of dependencies among APIs. In this paper, we propose a novel tree-based approach that can better capture the essential dependencies and largely improve the efficiency of RESTful API fuzzing. In particular, the hierarchical information of the endpoints across multiple APIs enables us to construct an API tree, and the relationships of tree nodes can indicate the priority of resource dependencies, \textit{e.g.,} it's more likely that a node depends on its parent node rather than its offspring or siblings. In the evaluation part, we first confirm that such a tree-based approach is more efficient than traditional graph-based approaches. We then apply our tool to fuzz two real-world RESTful services and compare the performance with two state-of-the-art tools, EvoMaster and RESTler. Our results show that foREST can improve the code coverage in all experiments, ranging from 11.5\% to 82.5\%. Besides, our tool finds 11 new bugs previously unknown.

preprint2022arXiv

Iterative Adaptively Regularized LASSO-ADMM Algorithm for CFAR Estimation of Sparse Signals: IAR-LASSO-ADMM-CFAR Algorithm

The least-absolute shrinkage and selection operator (LASSO) is a regularization technique for estimating sparse signals of interest emerging in various applications and can be efficiently solved via the alternating direction method of multipliers (ADMM), which will be termed as LASSO-ADMM algorithm. The choice of the regularization parameter has significant impact on the performance of LASSO-ADMM algorithm. However, the optimization for the regularization parameter in the existing LASSO-ADMM algorithms has not been solved yet. In order to optimize this regularization parameter, we propose an efficient iterative adaptively regularized LASSO-ADMM (IAR-LASSO-ADMM) algorithm by iteratively updating the regularization parameter in the LASSO-ADMM algorithm. Moreover, a method is designed to iteratively update the regularization parameter by adding an outer iteration to the LASSO-ADMM algorithm. Specifically, at each outer iteration the zero support of the estimate obtained by the inner LASSO-ADMM algorithm is utilized to estimate the noise variance, and the noise variance is utilized to update the threshold according to a pre-defined const false alarm rate (CFAR). Then, the resulting threshold is utilized to update both the non-zero support of the estimate and the regularization parameter, and proceed to the next inner iteration. In addition, a suitable stopping criterion is designed to terminate the outer iteration process to obtain the final non-zero support of the estimate of the sparse measurement signals. The resulting algorithm is termed as IAR-LASSO-ADMM-CFAR algorithm. Finally, simulation results have been presented to show that the proposed IAR-LASSO-ADMM-CFAR algorithm outperforms the conventional LASSO-ADMM algorithm and other existing algorithms in terms of reconstruction accuracy, and its sparsity order estimate is more accurate than the existing algorithms.

preprint2022arXiv

On structure theorems and non-saturated examples

For any minimal system $(X,T)$ and $d\geq 1$ there is an associated minimal system $(N_{d}(X), \mathcal{G}_{d}(T))$, where $\mathcal{G}_{d}(T)$ is the group generated by $T\times\cdots\times T$ and $T\times T^2\times\cdots\times T^{d}$ and $N_{d}(X)$ is the orbit closure of the diagonal under $\mathcal{G}_{d}(T)$. It is known that the maximal $d$-step pro-nilfactor of $N_d(X)$ is $N_d(X_d)$, where $X_d$ is the maximal $d$-step pro-nilfactor of $X$. In this paper, we further study the structure of $N_d(X)$. We show that the maximal distal factor of $N_d(X)$ is $N_d(X_{dis})$ with $X_{dis}$ being the maximal distal factor of $X$, and prove that as minimal systems $(N_{d}(X), \mathcal{G}_{d}(T))$ has the same structure theorem as $(X,T)$. In addition, a non-saturated metric example $(X,T)$ is constructed, which is not $T\times T^2$-saturated and is a Toeplitz minimal system.

preprint2022arXiv

Resonant tunneling in disordered borophene nanoribbons with line defects

Very recently, borophene has been attracting extensive and ongoing interest as the new wonder material with structural polymorphism and superior attributes, showing that the structural imperfection of line defects (LDs) occurs widely at the interface between $ν_{1/5}$ ($χ_3$) and $ν_{1/6}$ ($β_{12}$) boron sheets. Motivated by these experiments, here we present a theoretical study of electron transport through two-terminal disordered borophene nanoribbons (BNRs) with random distribution of LDs. Our results indicate that LDs could strongly affect the electron transport properties of BNRs. In the absence of LDs, both $ν_{1/5}$ and $ν_{1/6}$ BNRs exhibit metallic behavior, in agreement with experiments. While in the presence of LDs, the overall electron transport ability is dramatically decreased, but some resonant peaks of conductance quantum can be found in the transmission spectrum of any disordered BNR with arbitrary arrangement of LDs. These disordered BNRs exhibit metal-insulator transition by varying nanoribbon width with tunable transmission gap in the insulating regime. Furthermore, the bond currents present fringe patterns and two evolution phenomena of resonant peaks are revealed for disordered BNRs with different widths. These results may help for understanding structure-property relationships and designing LD-based nanodevices.

preprint2022arXiv

STELLA: Sparse Taint Analysis for Enclave Leakage Detection

Intel SGX (Software Guard Extension) is a promising TEE (trusted execution environment) technique that can protect programs running in user space from being maliciously accessed by the host operating system. Although it provides hardware access control and memory encryption, the actual effectiveness also depends on the quality of the software. In particular, improper implementation of a code snippet running inside the enclave may still leak private data due to the invalid use of pointers. This paper serves as a first attempt to study the privacy leakage issues of enclave code and proposes a novel static sparse taint analysis approach to detect them. We first summarize five common patterns of leakage code. Based on these patterns, our approach performs forward analysis to recognize all taint sinks and then employs a backward approach to detect leakages. Finally, we have conducted experiments with several open-source enclave programs and found 78 vulnerabilities previously unknown in 13 projects.

preprint2022arXiv

The China Trade Shock and the ESG Performances of US firms

How does import competition from China affect engagement on ESG initiatives by US corporates? On the one hand, reduced profitability due to import competition and lagging ESG performance of Chinese exporters can disincentivize US firms to put more resources to ESG initiatives. On the other hand, the shift from labor-intensive production to capital/technology-intensive production along with offshoring may improve the US company's ESG performance. Moreover, US companies have incentives to actively pursue more ESG engagement to differentiate from Chinese imports. Exploiting a trade policy in which US congress granted China the Permanent Normal Trade Relations and the resulting change in expected tariff rates on Chinese imports, we find that greater import competition from China leads to an increase in the US company's ESG performance. The improvement primarily stems from "doing more positives" and from more involvement on environmental initiatives. Indirect and direct evidence shows that the improvement is not driven by the change in production process or offshoring, but is consistent with product differentiation. Our results suggest that the trade shock from China has significant impact on the US company's ESG performance.

preprint2022arXiv

The structure of pointwise recurrent expansive homeomorphisms

Let $X$ be a compact metric space and let $f:X\rightarrow X$ be a homeomorphism on $X$. We show that if $f$ is both pointwise recurrent and expansive, then the dynamical system $(X, f)$ is topologically conjugate to a subshift of some symbolic system. Moreover, if $f$ is pointwise positively recurrent, then the subshift is semisimple; a counterexample is given to show the necessity of positive recurrence to ensure the semisimilicity.

preprint2022arXiv

The structures of higher rank lattice actions on dendrites

Let $Γ$ be a higher rank lattice acting on a nondegenerate dendrite $X$ with no infinite order points. We show that there exists a nondegenerate subdendrite $Y$ which is $Γ$-invariant and satisfies the following items: (1) There is an inverse system of finite actions $\{(Y_i, Γ):i=1,2,3,\cdots\}$ with monotone bonding maps $ϕ_i: Y_{i+1}\rightarrow Y_i$ and with each $Y_i$ being a dendrite, such that $(Y, Γ|Y)$ is topologically conjugate to the inverse limit $(\underset{\longleftarrow}{\lim}(Y_i, Γ), Γ)$. (2) The first point map $r:X\rightarrow Y$ is a factor map from $(X, Γ)$ to $(Y, Γ|Y)$; if $x\in X\setminus Y$, then $r(x)$ is an end point of $Y$ with infinite orbit; for each $y\in Y$, $r^{-1}(y)$ is contractible, that is there is a sequence $g_i\in Γ$ with ${\rm diam}(g_ir^{-1}(y))\rightarrow 0$.

preprint2021arXiv

DROID: Minimizing the Reality Gap using Single-Shot Human Demonstration

Reinforcement learning (RL) has demonstrated great success in the past several years. However, most of the scenarios focus on simulated environments. One of the main challenges of transferring the policy learned in a simulated environment to real world, is the discrepancy between the dynamics of the two environments. In prior works, Domain Randomization (DR) has been used to address the reality gap for both robotic locomotion and manipulation tasks. In this paper, we propose Domain Randomization Optimization IDentification (DROID), a novel framework to exploit single-shot human demonstration for identifying the simulator's distribution of dynamics parameters, and apply it to training a policy on a door opening task. Our results show that the proposed framework can identify the difference in dynamics between the simulated and the real worlds, and thus improve policy transfer by optimizing the simulator's randomization ranges. We further illustrate that based on these same identified parameters, our method can generalize the learned policy to different but related tasks.

preprint2021arXiv

Explore missing flow dynamics by physics-informed deep learning: the parameterised governing systems

Gaining and understanding the flow dynamics have much importance in a wide range of disciplines, e.g. astrophysics, geophysics, biology, mechanical engineering and biomedical engineering. As a reliable way in practice, especially for turbulent flows, regional flow information such as velocity and its statistics, can be measured experimentally. Due to the poor fidelity or other technical limitations, some information may not be resolved in a region of interest. On the other hand, detailed flow features are described by the governing equations, e.g. the Navier-Stokes equations for viscous fluid, and can be resolved numerically, which is heavily dependent on the capability of either computing resources or modelling. Alternatively, we address this problem by employing the physics-informed deep learning, and treat the governing equations as a parameterised constraint to recover the missing flow dynamics. We demonstrate that with limited data, no matter from experiment or others, the flow dynamics in the region where the required data is missing or not measured, can be reconstructed with the parameterised governing equations. Meanwhile, a richer dataset, with spatial distribution of the control parameter (e.g. eddy viscosity of turbulence modellings), can be obtained. The method provided in this paper may shed light on data-driven scale-adaptive turbulent structure recovering and understanding of complex fluid physics, and can be extended to other parameterised governing systems beyond fluid mechanics.

preprint2021arXiv

Memory-Safety Challenge Considered Solved? An In-Depth Study with All Rust CVEs

Rust is an emerging programing language that aims at preventing memory-safety bugs without sacrificing much efficiency. The claimed property is very attractive to developers, and many projects start using the language. However, can Rust achieve the memory-safety promise? This paper studies the question by surveying 186 real-world bug reports collected from several origins which contain all existing Rust CVEs (common vulnerability and exposures) of memory-safety issues by 2020-12-31. We manually analyze each bug and extract their culprit patterns. Our analysis result shows that Rust can keep its promise that all memory-safety bugs require unsafe code, and many memory-safety bugs in our dataset are mild soundness issues that only leave a possibility to write memory-safety bugs without unsafe code. Furthermore, we summarize three typical categories of memory-safety bugs, including automatic memory reclaim, unsound function, and unsound generic or trait. While automatic memory claim bugs are related to the side effect of Rust newly-adopted ownership-based resource management scheme, unsound function reveals the essential challenge of Rust development for avoiding unsound code, and unsound generic or trait intensifies the risk of introducing unsoundness. Based on these findings, we propose two promising directions towards improving the security of Rust development, including several best practices of using specific APIs and methods to detect particular bugs involving unsafe code. Our work intends to raise more discussions regarding the memory-safety issues of Rust and facilitate the maturity of the language.

preprint2021arXiv

Randomization Inference for Composite Experiments with Spillovers and Peer Effects

Group-formation experiments, in which experimental units are randomly assigned to groups, are a powerful tool for studying peer effects in the social sciences. Existing design and analysis approaches allow researchers to draw inference from such experiments without relying on parametric assumptions. In practice, however, group-formation experiments are often coupled with a second, external intervention, that is not accounted for by standard nonparametric approaches. This note shows how to construct Fisherian randomization tests and Neymanian asymptotic confidence intervals for such composite experiments, including in settings where the second intervention exhibits spillovers. We also propose an approach for designing optimal composite experiments.

preprint2021arXiv

Sensitive group actions on regular curves of almost $\leq n$ order

Let $X$ be a regular curve and $n$ be a positive integer such that for every nonempty open set $U\subset X$, there is a nonempty connected open set $V\subset U$ with the cardinality $|\partial_X(V)|\leq n$. We show that if $X$ admits a sensitive action of a group $G$, then $G$ contains a free subsemigroup and the action has positive geometric entropy. As a corollary, $X$ admits no sensitive nilpotent group action.

preprint2020arXiv

Deep Reinforcement Learning in Fluid Mechanics: a promising method for both Active Flow Control and Shape Optimization

In recent years, Artificial Neural Networks (ANNs) and Deep Learning have become increasingly popular across a wide range of scientific and technical fields, including Fluid Mechanics. While it will take time to fully grasp the potentialities as well as the limitations of these methods, evidence is starting to accumulate that point to their potential in helping solve problems for which no theoretically optimal solution method is known. This is particularly true in Fluid Mechanics, where problems involving optimal control and optimal design are involved. Indeed, such problems are famously difficult to solve effectively with traditional methods due to the combination of non linearity, non convexity, and high dimensionality they involve. By contrast, Deep Reinforcement Learning (DRL), a method of optimization based on teaching empirical strategies to an ANN through trial and error, is well adapted to solving such problems. In this short review, we offer an insight into the current state of the art of the use of DRL within fluid mechanics, focusing on control and optimal design problems.

preprint2020arXiv

Invariant Radon measures and minimal sets for subgroups of $\text{Homeo}_+(\mathbb{R})$

Let $G$ be a subgroup of $\text{Homeo}_+(\mathbb{R})$ without crossed elements. We show the equivalence among three items: (1) existence of $G$-invariant Radon measures on $\mathbb R$; (2) existence of minimal closed subsets of $\mathbb R$; (3) nonexistence of infinite towers covering the whole line. For a nilpotent subgroup $G$ of $\text{Homeo}_+(\mathbb{R})$, we show that $G$ always has an invariant Radon measure and a minimal closed set if every element of $G$ is $C^{1+α} (α>0$); a counterexample of $C^1$ commutative subgroup of $\text{Homeo}_+(\mathbb{R})$ is constructed.

preprint2020arXiv

Topological conjugation classes of tightly transitive subgroups of $\text{Homeo}_{+}(\mathbb{S}^1)$

Let $\text{Homeo}_{+}(\mathbb{S}^1)$ denote the group of orientation preserving homeomorphisms of the circle $\mathbb{S}^1$. A subgroup $G$ of $\text{Homeo}_{+}(\mathbb{S}^1)$ is tightly transitive if it is topologically transitive and no subgroup $H$ of $G$ with $[G: H]=\infty$ has this property; is almost minimal if it has at most countably many nontransitive points. In the paper, we determine all the topological conjugation classes of tightly transitive and almost minimal subgroups of $\text{Homeo}_{+}(\mathbb{S}^1)$ which are isomorphic to $\mathbb{Z}^n$ for any integer $n\geq 2$.

preprint2019arXiv

Nektar++: enhancing the capability and application of high-fidelity spectral/$hp$ element methods

Nektar++ is an open-source framework that provides a flexible, high-performance and scalable platform for the development of solvers for partial differential equations using the high-order spectral/$hp$ element method. In particular, Nektar++ aims to overcome the complex implementation challenges that are often associated with high-order methods, thereby allowing them to be more readily used in a wide range of application areas. In this paper, we present the algorithmic, implementation and application developments associated with our Nektar++ version 5.0 release. We describe some of the key software and performance developments, including our strategies on parallel I/O, on in situ processing, the use of collective operations for exploiting current and emerging hardware, and interfaces to enable multi-solver coupling. Furthermore, we provide details on a newly developed Python interface that enables a more rapid introduction for new users unfamiliar with spectral/$hp$ element methods, C++ and/or Nektar++. This release also incorporates a number of numerical method developments - in particular: the method of moving frames, which provides an additional approach for the simulation of equations on embedded curvilinear manifolds and domains; a means of handling spatially variable polynomial order; and a novel technique for quasi-3D simulations to permit spatially-varying perturbations to the geometry in the homogeneous direction. Finally, we demonstrate the new application-level features provided in this release, namely: a facility for generating high-order curvilinear meshes called NekMesh; a novel new AcousticSolver for aeroacoustic problems; our development of a 'thick' strip model for the modelling of fluid-structure interaction problems in the context of vortex-induced vibrations. We conclude by commenting some directions for future code development and expansion.

preprint2019arXiv

Oil drop deposition on solid surfaces in mixed polymer-surfactant solutions in relation to hair- and skin-care applications

The deposition of oil drops on solid substrates from mixed solutions of surfactants and cationic polymer is investigated. The used anionic surfactants are sodium laurylethersulfate (SLES) and sulfonated methyl esters (SME); the zwitterionic surfactant is cocamidopropyl betaine (CAPB). A new method, called the pressed drop method (PDM), was proposed to study the drop adhesion to substrates of different hydrophobicity. The PDM allows one to detect the presence or absence of drop adhesion at different degrees of dilution of the initial solution and, thus, to determine the threshold concentration of drop adhesion. The results show that the increase of the fraction of CAPB in the mixture with the anionic surfactant suppresses the oil-drop deposition; SME provides easier drop adhesion than SLES; the addition of NaCl enhances, whereas coco-fatty-acid-monoethanolamide (CMEA) suppresses the drop deposition; no drop adhesion is observed in the absence of polymer. The drop-to-substrate adhesion is interpreted in terms of the acting surface forces: polymer bridging attraction; hydrophobic attraction between segments of adsorbed polymer brushes and electrostatic forces. From viewpoint of applications, the PDM experiments enable one to compare the performance of various components in personal care formulations and to optimize their composition with respect to the oil-drop deposition.

preprint2016arXiv

Dramatic increase of the onset critical temperature and critical field of elemental Sn in the form of thin nanowires

Sn is a well-known classical superconductor on the border between type I and type II with critical temperature of 3.722K and critical field of 0.031T. We show by means of specific heat and electric magneto-transport data that its critical parameters can be dramatically increased if it is brought in the form of loosely bound bundles of thin nanowires. The specific heat displays a pronounced double phase transition at 3.7K and 5.5K, which we attribute to the inner 'bulk' contribution of the nanowires and to the surface contribution, respectively. The latter is visible only because of the large volume fraction of the surface layer in relation to their bulk inner volume. The upper transition coincides with the onset of the resistive transition, while zero resistance is gradually approached below the lower transition. The large coherence length of 230nm at 0K likely actuates a Josephson coupling between adjacent neighboring nanowires and thus suppresses the effect of 1D phase fluctuations along the nanowires, and stabilizes 3D phase coherence throughout the entire network with zero resistance. A magnetic field of more than 3T is required to restore the normal state, which means that the critical field is enhanced by about two orders of magnitude with respect to Sn in its bulk form.

preprint2016arXiv

Influence of localised smooth steps on the instability of a boundary layer

We consider a smooth forward facing step defined by the Gauss error function of height 4-30\% and four times the width of the local boundary layer thickness $δ_{99}$. The boundary layer flow over a smooth forward-facing stepped plate is studied with particular emphasis on stabilisation and destabilisation of the Tollmien-Schlichting (TS) waves and subsequently on transition. The interaction between TS waves at a range of frequencies and a base flow over a single/two forward facing smooth steps is conducted by linear analysis. The results indicate that for a high frequency TS wave, the amplitude of the TS wave is attenuated in the unstable regime of the neutral stability curve corresponding to a flat plate boundary layer. Furthermore, it is observed that two smooth forward facing steps lead to a more acute reduction of the amplitude of the TS wave. When the height of a step is increased to more than 20\% of the local boundary layer thickness for a fixed width parameter, the TS wave is amplified and thereby a destabilisation effect is introduced. Therefore, stabilisation or destabilisation effect of a smooth step is typically dependent on its shape parameters. To validate the results of the linear stability analysis, where a high-frequency TS wave is damped by the forward facing smooth steps direct numerical simulation (DNS) is performed. The results of the DNS correlate favorably with the linear analysis and show that for the investigated high frequency TS wave, the K-type transition process is altered whereas the onset of the H-type transition is postponed. The results of the DNS suggest that for a high-frequency perturbation $\mathcal{F}=150$ and in the absence of other external perturbations, two forward facing steps of height 5\% and 12\% of the boundary layer thickness delayed H-type transition scenario and completely suppresses it for the K-type transition.

preprint2016arXiv

Integrated Energy Exchange Scheduling for Multi-microgrid System with Electric Vehicles

Electric vehicles (EVs) can be considered as flexible mobile battery storages in microgrids. For multiple microgrids in an area, coordinated scheduling on charging and discharging are required to avoid power exchange spikes between the multimicrogrid system and the main grid. In this paper, a two-stage integrated energy exchange scheduling strategy for multimicrogrid system is presented, which considers EVs as storage devices. Then several dual variables, which are representative of the marginal cost of proper constraints, are utilized to form an updated price, thereby being a modification on the original electricity price. With this updated price signal, a price-based decentralized scheduling strategy is presented for the Microgrid Central Controller (MGCC). Simulation results show that the two-stage scheduling strategy reduces the electricity cost and avoids frequent transitions between battery charging/discharging states. With the proposed decentralized scheduling strategy, each microgrid only needs to solve its local problem and limits the total power exchange within the safe range.

preprint2016arXiv

IntelliAd Understanding In-APP Ad Costs From Users Perspective

Ads are an important revenue source for mobile app development, especially for free apps, whose expense can be compensated by ad revenue. The ad benefits also carry with costs. For example, too many ads can interfere the user experience, leading to less user retention and reduced earnings ultimately. In the paper, we aim at understanding the ad costs from users perspective. We utilize app reviews, which are widely recognized as expressions of user perceptions, to identify the ad costs concerned by users. Four types of ad costs, i.e., number of ads, memory/CPU overhead, traffic usage, and bettery consumption, have been discovered from user reviews. To verify whether different ad integration schemes generate different ad costs, we first obtain the commonly used ad schemes from 104 popular apps, and then design a framework named IntelliAd to automatically measure the ad costs of each scheme. To demonstrate whether these costs indeed influence users reactions, we finally observe the correlations between the measured ad costs and the user perceptions. We discover that the costs related to memory/CPU overhead and battery consumption are more concerned by users, while the traffic usage is less concerned by users in spite of its obvious variations among different schemes in the experiments. Our experimental results provide the developers with suggestions on better incorporating ads into apps and, meanwhile, ensuring the user experience.

preprint2016arXiv

Minimax Design and Order Estimation of FIR Filters for Extending the Bandwidth of ADCs

The bandwidth of the sampling systems, especially for time-interleaved analog-to-digital converters, needs to be extended along with the rapid increase of the sampling rate. A digitally assisted technique becomes a feasible approach to extend the analog bandwidth, as it is impractical to implement the extension in analog circuits. This paper derives accurate order estimation formulas for the bandwidth extension filter, which is designed in the minimax sense with the ripple constraints as the design criteria. The derived filter order estimation is significant in evaluating the computational complexity from the viewpoint of the top-level system design. Moreover, with the proposed order estimates, one can conveniently obtain the minimal order that satisfies the given ripple constraints, which contributes to reducing the design time. Both the performance of the extension filter and its order estimation are illustrated and demonstrated through simulation examples.

preprint2016arXiv

On closedness under convolution roots related to an infinitely divisible distribution in the distribution class L(γ)

We consider questions related to the well-known conjecture due to Embrechts and Goldie on the closedness of different classes of heavy- and light-tailed distributions with respect to convolution roots. We show that the class L(γ)\cap OS is not closed under convolution roots related to an infinitely divisible distribution for any γ\ge0, i.e. we provide examples of infinitely divisible distributions belonging to this class such that the corresponding Levy spectral distribution does not. We also prove a similar statement for the class (L(γ)\cap OS)\ S(γ). In order to facilitate our analysis, we explore the structural properties of some of the classes of distributions, and study some properties of the well-known transformation from a heavy-tailed distribution to a light-tailed one.

preprint2016arXiv

The local asymptotic estimation for the supremum of a random walk with generalized strong subexponential summands

In this paper, the local asymptotic estimation for the supremum of a random walk and its applications are presented. The summands of the random walk have common long-tailed and generalized strong subexponential distribution. This distribution class and the corresponding generalized local subexponential distribution class are two new distribution classes with some good properties. Further, some long-tailed distributions with intuitive and concrete forms are found, which show that the intersection of the two above-mentioned distribution classes with long-tailed distribution class properly contain the strong subexponential distribution class and the locally subexponential distribution class, respectively.

preprint2016arXiv

The product of dependent random variables with applications to a discrete-time risk model

Let $X$ be a real valued random variable with an unbounded distribution $F$ and let $Y$ be a nonnegative valued random variable with a unbounded distribution $G$, which satisfy that \begin{eqnarray*} P(X>x|Y=y)\sim h(y)P(X>x) \end{eqnarray*} holds uniformly for $y\geq0$ as $x\to \infty$. Under the condition that $\overline{G}(bx)=o(\overline H(x))$ holds for all constant $b>0$, we proved that $F\in\mathcal{L}(γ)$ for some $γ\geq 0$ implied $H\in \mathcal{L}(γ/β_G)$ and that $F\in\mathcal{S}(γ)$ for some $γ\geq 0$ implied $H\in \mathcal{S}(γ/β_G)$, where $H$ is the distribution of the product $XY$, and $β_G$ is the right endpoint of $G$, that is, $β_G=\sup\{y:~G(y)<1\}\in (0,\infty],$ and when $β_G=\infty$, $γ/β_G$ is understood as 0. Furthermore, in a discrete-time risk model in which the net insurance loss and the stochastic discount factor are equipped with a dependence structure, a general asymptotic formula for the finite-time ruin probability is obtained when the net insurance loss has a subexponential tail.

preprint2016arXiv

Wireless Fractal Cellular Networks

With the seamless coverage of wireless cellular networks in modern society, it is interesting to consider the shape of wireless cellular coverage. Is the shape a regular hexagon, an irregular polygon, or another complex geometrical shape? Based on fractal theory, the statistical characteristic of the wireless cellular coverage boundary is determined by the measured wireless cellular data collected from Shanghai, China. The measured results indicate that the wireless cellular coverage boundary presents an extremely irregular geometrical shape, which is also called a statistical fractal shape. Moreover, the statistical fractal characteristics of the wireless cellular coverage boundary have been validated by values of the Hurst parameter estimated in angular scales. The statistical fractal characteristics of the wireless cellular coverage boundary can be used to evaluate and design the handoff scheme of mobile user terminals in wireless cellular networks.

preprint2015arXiv

Convolution and convolution-root properties of long-tailed distributions

We obtain a number of new general properties, related to the closedness of the class of long-tailed distributions under convolutions, that are of interest themselves and may be applied in many models that deal with "plus" and/or "max" operations on heavy-tailed random variables. We analyse the closedness property under convolution roots for these distributions. Namely, we introduce two classes of heavy-tailed distributions that are not long-tailed and study their properties. These examples help to provide further insights and, in particular, to show that the properties to be both long-tailed and so-called "generalised subexponential" are not preserved under the convolution roots. This leads to a negative answer to a conjecture of Embrechts and Goldie [10, 12] for the class of long-tailed and generalised subexponential distributions. In particular, our examples show that the following is possible: an infinitely divisible distribution belongs to both classes, while its Levy measure is neither long-tailed nor generalised subexponential.

preprint2015arXiv

N-Version Obfuscation: Impeding Software Tampering Replication with Program Diversity

Tamper-resistance is a fundamental software security research area. Many approaches have been proposed to thwart specific procedures of tampering, e.g., obfuscation and self-checksumming. However, to our best knowledge, none of them can achieve theoretically tamper-resistance. Our idea is to impede the replication of tampering via program diversification, and thus increasing the complexity to break the whole software system. To this end, we propose to deliver same featured, but functionally nonequivalent software copies to different machines. We formally define the problem as N-version obfuscation, and provide a viable means to solve the problem. Our evaluation result shows that the time required for breaking a software system is linearly increased with the number of software versions, which is O(n) complexity.

preprint2015arXiv

On the structure of a class of distributions obeying the principle of a single big jump

In this paper, we present several heavy-tailed distributions belonging to the new class J of distributions obeying the principle of a single big jump introduced by Beck et al. [1]. We describe the structure of this class from different angles. First, we show that heavy-tailed distributions in the class J are automatically strongly heavy-tailed and thus have tails which are not too irregular. Second, we show that such distributions are not necessarily weakly tail equivalent to a subexponential distribution. We also show that the class of heavy-tailed distributions in J which are neither long-tailed nor dominatedly-varying-tailed is not only non-empty but even quite rich in the sense that it has a nonempty intersection with several other well-established classes. In addition, the integrated tail distribution of some particular of these distributions shows that the Pakes-Veraverbeke-Embrechts Theorem for the class J in [1] does not hold trivially.

preprint2015arXiv

PersisDroid: Android Performance Diagnosis via Anatomizing Asynchronous Executions

Android applications (apps) grow dramatically in recent years. Apps are user interface (UI) centric typically. Rapid UI responsiveness is key consideration to app developers. However, we still lack a handy tool for profiling app performance so as to diagnose performance problems. This paper presents PersisDroid, a tool specifically designed for this task. The key notion of PersisDroid is that the UI-triggered asynchronous executions also contribute to the UI performance, and hence its performance should be properly captured to facilitate performance diagnosis. However, Android allows tremendous ways to start the asynchronous executions, posing a great challenge to profiling such execution. This paper finds that they can be grouped into six categories. As a result, they can be tracked and profiled according to the specifics of each category with a dynamic instrumentation approach carefully tailored for Android. PersisDroid can then properly profile the asynchronous executions in task granularity, which equips it with low-overhead and high compatibility merits. Most importantly, the profiling data can greatly help the developers in detecting and locating performance anomalies. We code and open-source release PersisDroid. The tool is applied in diagnosing 20 open-source apps, and we find 11 of them contain potential performance problems, which shows its effectiveness in performance diagnosis for Android apps.

preprint2014arXiv

On a transformation between distributions obeying the principle of a single big jump

Beck et al. (2013) introduced a new distribution class J which contains many heavy-tailed and light-tailed distributions obeying the principle of a single big jump. Using a simple transformation which maps heavy-tailed distributions to light-tailed ones, we find some light-tailed distributions, which belong to the class J but do not belong to the convolution equivalent distribution class and which are not even weakly tail equivalent to any convolution equivalent distribution. This fact helps to understand the structure of the light-tailed distributions in the class J and leads to a negative answer to an open question raised by the above paper.

preprint2012arXiv

Analysis of the absorbing layers for the weakly-compressible lattice Boltzmann schemes

It has been demonstrated that Lattice Boltzmann schemes (LBSs) are very efficient for Computational AeroAcoustics (CAA). In order to handle the issue of absorbing acoustic boundary conditions for LBS, three kinds of damping terms are proposed and added into the right hand sides of the governing equations of LBS. From the classical theory, these terms play an important role to absorb and minimize the acoustic wave reflections from computational boundaries. Meanwhile, the corresponding macroscopic equations with the damping terms are recovered for analyzing the macroscopic behaviors of the these damping terms and determining the critical absorbing strength. Further, in order to detect the dissipation and dispersion behaviors, the linearized LBS with the damping terms is derived and analyzed. The dispersive and dissipative properties are explored in the wave-number spaces via the Von Neumann analysis. The related damping strength critical values and the optimal absorbing term are addressed. Finally, some benchmark problems are implemented to assess the theoretical results.

preprint2012arXiv

Brownian Model Theory of Nonequilibrium Liquid Structure and Hydrodynamics of Strong Electrolyte Solutions in an External Field

In this paper, on the basis of the Onsager--Wilson theory of strong binary electrolyte solutions we completely work out the solutions of the governing equations (Onsager-Fuoss equations and Poisson equations) for nonequilibrium pair correlation functions and ionic potentials and the solutions for the Stokes equation for the velocity and pressure in the case of strong binary electrolyte solutions under the influence of an external electric field of arbitrary strength. The solutions are calculated in the configuration space as functions of coordinates and reduced field strength. Thus the axial and transversal components of the velocity and the accompanying nonequilibrium pressure are explicitly obtained and computed for all values of position coordinates. Computation of velocity profiles makes it possible to visualize the movement and distortion of ion atmosphere under the influence of an external electric field.

preprint2012arXiv

Wien Effect on Ionic Conductance of Binary Strong Electrolyte Solutions in a High External Electric Field

In the preceding paper, the exact solution of Stokes equation was obtained for a binary strong electrolyte solution in an external electric field. In the present paper, the solution is applied to calculate the Wien effect on deviation from the Coulombic law of conduction in high fields. One of the important aims of the present line of work was in removing or avoiding the divergence difficulty in calculating the electrophoretic and relaxation time coefficients. The present work achieves that aim by calculating on the basis of computing the axial velocity profiles the position of the center of the ion atmosphere as a function of the reduced field strength and therewith computing the electrophoretic and relaxation time coefficients for the migrating spherical ion atmosphere at each value of the reduced field strength. With the electrophoretic and relaxation time coefficients thus calculated along the trajectory of the center of the ion atmosphere with regard to the reduced field strength, the equivalent ionic conductance of a magnesium sulfate solution is calculated and compared with experimental data in very good accuracy. It is also compared with the prediction by Wilson's theory and is found very much improved over the latter. Moreover, the present method is divergence-free unlike Wilson's theory. It therefore not only sheds a new light on the nonequilibrium liquid structure of ionic solutions and the motion of ion atmosphere in an electric field, but also shows how the Onsager--Wilson theory of conductance in electrolyte solutions should be applied to account for the Wien effect on ionic conductance.

preprint2011arXiv

A Lifting Relation from Macroscopic Variables to Mesoscopic Variables in Lattice Boltzmann Method: Derivation, Numerical Assessments and Coupling Computations Validation

In this paper, analytic relations between the macroscopic variables and the mesoscopic variables are derived for lattice Boltzmann methods (LBM). The analytic relations are achieved by two different methods for the exchange from velocity fields of finite-type methods to the single particle distribution functions of LBM. The numerical errors of reconstructing the single particle distribution functions and the non-equilibrium distribution function by macroscopic fields are investigated. Results show that their accuracy is better than the existing ones. The proposed reconstruction operator has been used to implement the coupling computations of LBM and macro-numerical methods of FVM. The lid-driven cavity flow is chosen to carry out the coupling computations based on the numerical strategies of domain decomposition methods (DDM). The numerical results show that the proposed lifting relations are accurate and robust.

preprint2011arXiv

Optimal low-dispersion low-dissipation LBM schemes for computational aeroacoustics

Lattice Boltmzmann Methods (LBM) have been proved to be very effective methods for computational aeroacoustics (CAA), which have been used to capture the dynamics of weak acoustic fluctuations. In this paper, we propose a strategy to reduce the dispersive and disspative errors of the two-dimensional (2D) multi-relaxation-time lattice Boltzmann method (MRT-LBM). By presenting an effective algorithm, we obtain a uniform form of the linearized Navier-Stokes equations corresponding to the MRT-LBM in wave-number space. Using the matrix perturbation theory and the equivalent modified equation approach for finite difference methods, we propose a class of minimization problems to optimize the free-parameters in the MRT-LBM. We obtain this way a dispersion-relation-preserving LBM (DRP-LBM) to circumvent the minimized dispersion error of the MRT-LBM. The dissipation relation precision is also improved.And the stability of the MRT-LBM with the small bulk viscosity is guaranteed. Von Neuman analysis of the linearized MRT-LBM is performed to validate the optimized dispersion/dissipation relations considering monochromatic wave solutions. Meanwhile, dispersion and dissipation errors of the optimized MRT-LBM are quantitatively compared with the original MRT-LBM . Finally, some numerical simulations are carried out to assess the new optimized MRT-LBM schemes.

preprint2011arXiv

Sensitivity analysis and determination of free relaxation parameters for the weakly-compressible MRT-LBM schemes

It is well-known that there exist several free relaxation parameters in the MRT-LBM. Although these parameters have been tuned via linear analysis, the sensitivity analysis of these parameters and other related parameters are still not sufficient for detecting the behaviors of the dispersion and dissipation relations of the MRT-LBM. Previous researches have shown that the bulk dissipation in the MRT-LBM induces a significant over-damping of acoustic disturbances. This indicates that MRT-LBM cannot be used to obtain the correct behavior of pressure fluctuations because of the fixed bulk relaxation parameter. In order to cure this problem, an effective algorithm has been proposed for recovering the linearized Navier-Stokes equations from the linearized MRT-LBM. The recovered L-NSE appear as in matrix form with arbitrary order of the truncation errors with respect to $δt$. Then, in wave-number space, the first/second-order sensitivity analyses of matrix eigenvalues are used to address the sensitivity of the wavenumber magnitudes to the dispersion-dissipation relations. By the first-order sensitivity analysis, the numerical behaviors of the group velocity of the MRT-LBM are first obtained. Afterwards, the distribution sensitivities of the matrix eigenvalues corresponding to the linearized form of the MRT-LBM are investigated in the complex plane. Based on the sensitivity analysis and the recovered L-NSE, we propose some simplified optimization strategies to determine the free relaxation parameters in the MRT-LBM. Meanwhile, the dispersion and dissipation relations of the optimal MRT-LBM are quantitatively compared with the exact dispersion and dissipation relations. At last, some numerical validations on classical acoustic benchmark problems are shown to assess the new optimal MRT-LBM.

preprint2002arXiv

On the "Mandelbrot set" for a pair of linear maps and complex Bernoulli convolutions

We consider the "Mandelbrot set" $M$ for pairs of complex linear maps, introduced by Barnsley and Harrington in 1985 and studied by Bousch, Bandt and others. It is defined as the set of parameters $λ$ in the unit disk such that the attractor $A_λ$ of the IFS $\{λz-1, λz+1\}$ is connected. We show that a non-trivial portion of $M$ near the imaginary axis is contained in the closure of its interior (it is conjectured that all non-real points of $M$ are in the closure of the set of interior points of $M$). Next we turn to the attractors $A_λ$ themselves and to natural measures $ν_λ$ supported on them. These measures are the complex analogs of much-studied infinite Bernoulli convolutions. Extending the results of Erdös and Garsia, we demonstrate how certain classes of complex algebraic integers give rise to singular and absolutely continuous measures $ν_λ$. Next we investigate the Hausdorff dimension and measure of $A_λ$, for $λ$ in the set $M$, for Lebesgue-a.e. $λ$. We also obtain partial results on the absolute continuity of $ν_λ$ for a.e. $λ$ of modulus greater than $\sqrt{1/2}$.