Source author record

Xiaobing Feng

Xiaobing Feng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
18topics
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

40 published item(s)

preprint2026arXiv

Beyond Prompts: Space-Time Decoupling Control-Plane Jailbreaks in LLM Structured Output

Content Warning: This paper may contain unsafe or harmful content generated by LLMs that may be offensive to readers. Large Language Models (LLMs) are extensively used as tooling platforms through structured output APIs to ensure syntax compliance so that robust integration with existing software, like agent systems, can be achieved. However, the feature enabling the functionality of grammar-guided structured output presents significant security vulnerabilities. In this work, we reveal a critical control-plane attack surface orthogonal to traditional data-plane vulnerabilities. We introduce Constrained Decoding Attack (CDA), a novel jailbreak class that weaponizes structured output constraints to bypass both external auditing and internal safety alignment. Unlike prior attacks focused on input prompt designs, CDA operates by embedding malicious intent in schema-level grammar rules (control-plane) while maintaining benign surface prompts (data-plane). We instantiate this with two proof-of-concept attacks: EnumAttack, which embeds malicious content in enum fields; and the more evasive DictAttack, which decouples the malicious payload across a benign prompt and a dictionary-based grammar. Our evaluation spans a broad spectrum of 13 proprietary/open-weight models. In particular, DictAttack achieves 94.3--99.5% ASR across five benchmarks on gpt-5, gemini-2.5-pro, deepseek-r1, and gpt-oss-120b. Furthermore, we demonstrate the significant challenge in defending against these threats: while basic grammar auditing mitigates EnumAttack, the more sophisticated DictAttack maintains a 75.8% ASR even against multiple state-of-the-art jailbreak guardrails. This exposes a critical "semantic gap" in current safety architectures and underscores the urgent need for cross-plane defenses that can bridge the data and control planes to secure the LLM generation pipeline.

preprint2026arXiv

The New Compiler Stack: A Survey on the Synergy of LLMs and Compilers

This survey has provided a systematic overview of the emerging field of LLM-enabled compilation by addressing several key research questions. We first answered how LLMs are being integrated by proposing a comprehensive, multi-dimensional taxonomy that categorizes works based on their Design Philosophy (Selector, Translator, Generator), LLM Methodology, their operational Level of Code Abstraction, and the specific Task Type they address. In answering what advancements these approaches offer, we identified three primary benefits: the democratization of compiler development, the discovery of novel optimization strategies, and the broadening of the compiler's traditional scope. Finally, in addressing the field's challenges and opportunities, we highlighted the critical hurdles of ensuring correctness and achieving scalability, while identifying the development of hybrid systems as the most promising path forward. By providing these answers, this survey serves as a foundational roadmap for researchers and practitioners, charting the course for a new generation of LLM-powered, intelligent, adaptive and synergistic compilation tools.

preprint2022arXiv

A narrow-stencil framework for convergent numerical approximations of fully nonlinear second order PDEs

This paper develops a unified general framework for designing convergent finite difference and discontinuous Galerkin methods for approximating viscosity and regular solutions of fully nonlinear second order PDEs. Unlike the well-known monotone (finite difference) framework, the proposed new framework allows for the use of narrow stencils and unstructured grids which makes it possible to construct high order methods. The general framework is based on the concepts of consistency and g-monotonicity which are both defined in terms of various numerical derivative operators. Specific methods that satisfy the framework are constructed using numerical moments. Admissibility, stability, and convergence properties are proved, and numerical experiments are provided along with some computer implementation details.

preprint2022arXiv

A practical framework for multi-domain speech recognition and an instance sampling method to neural language modeling

Automatic speech recognition (ASR) systems used on smart phones or vehicles are usually required to process speech queries from very different domains. In such situations, a vanilla ASR system usually fails to perform well on every domain. This paper proposes a multi-domain ASR framework for Tencent Map, a navigation app used on smart phones and in-vehicle infotainment systems. The proposed framework consists of three core parts: a basic ASR module to generate n-best lists of a speech query, a text classification module to determine which domain the speech query belongs to, and a reranking module to rescore n-best lists using domain-specific language models. In addition, an instance sampling based method to training neural network language models (NNLMs) is proposed to address the data imbalance problem in multi-domain ASR. In experiments, the proposed framework was evaluated on navigation domain and music domain, since navigating and playing music are two main features of Tencent Map. Compared to a general ASR system, the proposed framework achieves a relative 13% $\sim$ 22% character error rate reduction on several test sets collected from Tencent Map and our in-car voice assistant.

preprint2022arXiv

An efficient multi-modes Monte Carlo homogenization method for random materials

In this paper, we propose and analyze a new stochastic homogenization method for diffusion equations with random and fast oscillatory coefficients. In the proposed method, the homogenized solutions are sought through a two-stage procedure. In the first stage, the original oscillatory diffusion equation is approximated, for each fixed random sample w, by a spatially homogenized diffusion equation with piecewise constant coefficients, resulting a random diffusion equation. In the second stage, the resulted random diffusion equation is approximated and computed by using an efficient multi-modes Monte Carlo method which only requires to solve a diffusion equation with a constant diffusion coefficient and a random right-hand side. The main advantage of the proposed method is that it separates the computational difficulty caused by the spatial fast oscillation of the solution and that caused by the randomness of the solution, so they can be overcome separately using different strategies. The convergence of the solution of the spatially homogenized equation (from the first stage) to the solution of the original random diffusion equation is established and the optimal rate of convergence is also obtained for the proposed multi-modes Monte Carlo method. Numerical experiments on some benchmark test problems for random composite materials are also presented to gauge the efficiency and accuracy of the proposed two-stage stochastic homogenization method.

preprint2022arXiv

Higher order time discretization for the stochastic semilinear wave equation with multiplicative noise

In this paper, a higher-order time-discretization scheme is proposed, where the iterates approximate the solution of the stochastic semilinear wave equation driven by multiplicative noise with general drift and diffusion. We employ a variational method for its error analysis and prove an improved convergence order of 3/2 for the approximates of the solution. The core of the analysis is Holder continuity in time and moment bounds for the solutions of the continuous and the discrete problem. Computational experiments are also presented.

preprint2021arXiv

High-order mass- and energy-conserving SAV-Gauss collocation finite element methods for the nonlinear Schrödinger equation

A family of arbitrarily high-order fully discrete space-time finite element methods are proposed for the nonlinear Schrödinger equation based on the scalar auxiliary variable formulation, which consists of a Gauss collocation temporal discretization and the finite element spatial discretization. The proposed methods are proved to be well-posed and conserving both mass and energy at the discrete level. An error bound of the form $O(h^p+τ^{k+1})$ in the $L^\infty(0,T;H^1)$-norm is established, where $h$ and $τ$ denote the spatial and temporal mesh sizes, respectively, and $(p,k)$ is the degree of the space-time finite elements. Numerical experiments are provided to validate the theoretical results on the convergence rates and conservation properties. The effectiveness of the proposed methods in preserving the shape of a soliton wave is also demonstrated by numerical results.

preprint2020arXiv

A new theory of fractional differential calculus

This paper presents a self-contained new theory of weak fractional differential calculus in one-dimension. The crux of this new theory is the introduction of a weak fractional derivative notion which is a natural generalization of integer order weak derivatives; it also helps to unify multiple existing fractional derivative definitions and characterize what functions are fractionally differentiable. Various calculus rules including a fundamental theorem calculus, product and chain rules, and integration by parts formulas are established for weak fractional derivatives. Additionally, relationships with classical fractional derivatives and detailed characterizations of weakly fractional differentiable functions are also established. Furthermore, the notion of weak fractional derivatives is also systematically extended to general distributions instead of only to some special distributions. This new theory lays down a solid theoretical foundation for systematically and rigorously developing new theories of fractional Sobolev spaces, fractional calculus of variations, and fractional PDEs as well as their numerical solutions in subsequent works. This paper is a concise presentation of the materials of Sections 1-4 and 6 of reference [9].

preprint2020arXiv

A new theory of fractional differential calculus and fractional Sobolev spaces: One-dimensional case

This paper presents a self-contained new theory of weak fractional differential calculus and fractional Sobolev spaces in one-dimension. The crux of this new theory is the introduction of a weak fractional derivative notion which is a natural generalization of integer order weak derivatives; it also helps to unify multiple existing fractional derivative definitions and characterize what functions are fractionally differentiable. Various calculus rules including a fundamental theorem of calculus, product and chain rules, and integration by parts formulas are established for weak fractional derivatives and relationships with classical derivatives are also obtained. Based on the weak fractional derivative notion, new fractional order Sobolev spaces are introduced and many important theorems and properties, such as density/approximation theorem, extension theorems, trace theorem, and various embedding theorems in these Sobolev spaces are established. Moreover, a few relationships with existing fractional Sobolev spaces are also established. Furthermore, the notion of weak fractional derivatives is also systematically extended to general distributions instead of only to some special distributions. The new theory lays down a solid theoretical foundation for systematically and rigorously developing a fractional calculus of variations theory and a fractional PDE theory as well as their numerical solutions in subsequent works.

preprint2020arXiv

Accelerating Deep Learning Inference with Cross-Layer Data Reuse on GPUs

Accelerating the deep learning inference is very important for real-time applications. In this paper, we propose a novel method to fuse the layers of convolutional neural networks (CNNs) on Graphics Processing Units (GPUs), which applies data reuse analysis and access optimization in different levels of the memory hierarchy. To achieve the balance between computation and memory access, we explore the fusion opportunities in the CNN computation graph and propose three fusion modes of convolutional neural networks: straight, merge and split. Then, an approach for generating efficient fused code is designed, which goes deeper in multi-level memory usage for cross-layer data reuse. The effectiveness of our method is evaluated with the network layers from state-of-the-art CNNs on two different GPU platforms, NVIDIA TITAN Xp and Tesla P4. The experiments show that the average speedup is 2.02x on representative structures of CNNs, and 1.57x on end-to-end inference of SqueezeNet.

preprint2020arXiv

Analysis of Fully Discrete Mixed Finite Element Methods for Time-dependent Stochastic Stokes Equations with Multiplicative Noise

This paper is concerned with fully discrete mixed finite element approximations of the time-dependent stochastic Stokes equations with multiplicative noise. A prototypical method, which comprises of the Euler-Maruyama scheme for time discretization and the Taylor-Hood mixed element for spatial discretization is studied in detail. Strong convergence with rates is established not only for the velocity approximation but also for the pressure approximation (in a time-averaged fashion). A stochastic inf-sup condition is established and used in a nonstandard way to obtain the error estimate for the pressure approximation in the time-averaged fashion. Numerical results are also provided to validate the theoretical results and to gauge the performance of the proposed fully discrete mixed finite methods.

preprint2020arXiv

LANCE: Efficient Low-Precision Quantized Winograd Convolution for Neural Networks Based on Graphics Processing Units

Accelerating deep convolutional neural networks has become an active topic and sparked an interest in academia and industry. In this paper, we propose an efficient low-precision quantized Winograd convolution algorithm, called LANCE, which combines the advantages of fast convolution and quantization techniques. By embedding linear quantization operations into the Winograd-domain, the fast convolution can be performed efficiently under low-precision computation on graphics processing units. We test neural network models with LANCE on representative image classification datasets, including SVHN, CIFAR, and ImageNet. The experimental results show that our 8-bit quantized Winograd convolution improves the performance by up to 2.40x over the full-precision convolution with trivial accuracy loss.

preprint2020arXiv

On New Families of Fractional Sobolev Spaces

This paper presents three new families of fractional Sobolev spaces and their accompanying theory in one-dimension. The new construction and theory are based on a newly developed notion of weak fractional derivatives, which are natural generalizations of the well-established integer order Sobolev spaces and theory. In particular, two new families of one-sided fractional Sobolev spaces are introduced and analyzed, they reveal more insights about another family of so-called symmetric fractional Sobolev spaces. Many key theorems/properties, such as density/approximation theorem, extension theorems, one-sided trace theorem, and various embedding theorems and Sobolev inequalities in those Sobolev spaces are established. Moreover, a few relationships with existing fractional Sobolev spaces are also discovered. The results of this paper lay down a solid theoretical foundation for systematically developing a fractional calculus of variations theory and a fractional PDE theory as well as their numerical solutions in subsequent works. This paper is a concise presentation of the materials of Sections 1, 4 and 5 of reference [7].

preprint2020arXiv

Optimally Convergent Mixed Finite Element Methods for the Stochastic Stokes Equations

We propose some new mixed finite element methods for the time dependent stochastic Stokes equations with multiplicative noise, which use the Helmholtz decomposition of the driving multiplicative noise. It is known [16] that the pressure solution has a low regularity, which manifests in sub-optimal convergence rates for well-known inf-sup stable mixed finite element methods in numerical simulations, see [10]. We show that eliminating this gradient part from the noise in the numerical scheme leads to optimally convergent mixed finite element methods, and that this conceptual idea may be used to retool numerical methods that are well-known in the deterministic setting, including pressure stabilization methods, so that their optimal convergence properties can still be maintained in the stochastic setting. Computational experiments are also provided to validate the theoretical results and to illustrate the conceptional usefulness of the proposed numerical approach.

preprint2016arXiv

An efficient Monte Carlo interior penalty discontinuous Galerkin method for elastic wave scattering in random media

This paper develops and analyzes an efficient Monte Carlo interior penalty discontinuous Galerkin (MCIP-DG) method for elastic wave scattering in random media. The method is constructed based on a multi-modes expansion of the solution of the governing random partial differential equations. It is proved that the mode functions satisfy a three-term recurrence system of partial differential equations (PDEs) which are nearly deterministic in the sense that the randomness only appears in the right-hand side source terms, not in the coefficients of the PDEs. Moreover, the same differential operator applies to all mode functions. A proven unconditionally stable and optimally convergent IP-DG method is used to discretize the deterministic PDE operator, an efficient numerical algorithm is proposed based on combining the Monte Carlo method and the IP-DG method with the $LU$ direct linear solver. It is shown that the algorithm converges optimally with respect to both the mesh size $h$ and the sampling number $M$, and practically its total computational complexity is only amount to solving very few deterministic elastic Helmholtz equations using the $LU$ direct linear solver. Numerically experiments are also presented to demonstrate the performance and key features of the proposed MCIP-DG method.

preprint2016arXiv

An enhanced finite element method for a class of variational problems exhibiting the Lavrentiev gap phenomenon

This paper develops an enhanced finite element method for approximating a class of variational problems which exhibit the \textit{Lavrentiev gap phenomenon} in the sense that the minimum values of the energy functional have a nontrivial gap when the functional is minimized on spaces $W^{1,1}$ and $W^{1,\infty}$. To remedy the standard finite element method, which fails to converge for such variational problems, a simple and effective cut-off procedure is utilized to design the (enhanced finite element) discrete energy functional. In essence the proposed discrete energy functional curbs the gap phenomenon by capping the derivatives of its input on a scale of $O(h^{-α})$ (where $h$ denotes the mesh size) for some positive constant $α$. A sufficient condition is proposed for determining the problem-dependent parameter $\a$. Extensive 1-D and 2-D numerical experiment results are provided to show the convergence behavior and the performance of the proposed enhanced finite element method.

preprint2016arXiv

Convergent semi-Lagrangian methods for the Monge-Ampère equation on unstructured grids

This paper is concerned with developing and analyzing convergent semi-Lagrangian methods for the fully nonlinear elliptic Monge-Ampère equation on general triangular grids. This is done by establishing an equivalent (in the viscosity sense) Hamilton-Jacobi-Bellman formulation of the Monge-Ampère equation. A significant benefit of the reformulation is the removal of the convexity constraint from the admissible space as convexity becomes a built-in property of the new formulation. Moreover, this new approach allows one to tap the wealthy numerical methods, such as semi-Lagrangian schemes, for Hamilton-Jacobi-Bellman equations to solve Monge-Ampère type equations. It is proved that the considered numerical methods are monotone, pointwise consistent and uniformly stable. Consequently, its solutions converge uniformly to the unique convex viscosity solution of the Monge-Ampère Dirichlet problem. A super-linearly convergent Howard's algorithm, which is a Newton type method, is utilized as the nonlinear solver to take advantage of the monotonicity of the scheme. Numerical experiments are also presented to gauge the performance of the proposed numerical method and the nonlinear solver.

preprint2016arXiv

Interior Penalty Discontinuous Galerkin Methods for Second Order Linear Non-Divergence Form Elliptic PDEs

This paper develops interior penalty discontinuous Galerkin (IP-DG) methods to approximate $W^{2,p}$ strong solutions of second order linear elliptic partial differential equations (PDEs) in non-divergence form with continuous coefficients. The proposed IP-DG methods are closely related to the IP-DG methods for advection-diffusion equations, and they are easy to implement on existing standard IP-DG software platforms. It is proved that the proposed IP-DG methods have unique solutions and converge with optimal rate to the $W^{2,p}$ strong solution in a discrete $W^{2,p}$-norm. The crux of the analysis is to establish a DG discrete counterpart of the Calderon-Zygmund estimate and to adapt a freezing coefficient technique used for the PDE analysis at the discrete level. As a byproduct of our analysis, we also establish broken $W^{1,p}$-norm error estimates for IP-DG approximations of constant coefficient elliptic PDEs. Numerical experiments are provided to gauge the performance of the proposed IP-DG methods and to validate the theoretical convergence results.

preprint2015arXiv

$C^0$ discontinuous Galerkin finite element methods for second order linear elliptic partial differential equations in non-divergence form

This paper is concerned with finite element approximations of $W^{2,p}$ strong solutions of second-order linear elliptic partial differential equations (PDEs) in non-divergence form with continuous coefficients. A nonstandard (primal) finite element method, which uses finite-dimensional subspaces consisting globally continuous piecewise polynomial functions, is proposed and analyzed. The main novelty of the finite element method is to introduce an interior penalty term, which penalizes the jump of the flux across the interior element edges/faces, to augment a nonsymmetric piecewise defined and PDE-induced bilinear form. Existence, uniqueness and error estimate in a discrete $W^{2,p}$ energy norm are proved for the proposed finite element method. This is achieved by establishing a discrete Calderon-Zygmund-type estimate and mimicking strong solution PDE techniques at the discrete level. Numerical experiments are provided to test the performance of proposed finite element method and to validate the convergence theory.

preprint2015arXiv

An unconditionally stable discontinuous Galerkin method for the elastic Helmholtz equations with large frequency

In this paper we propose and analyze an interior penalty discontinuous Galerkin (IP-DG) method using piecewise linear polynomials for the elastic Helmholtz equations with the first order absorbing boundary condition. It is proved that the sesquilinear form for the problem satisfies a generalized weak coercivity property, which immediately infers a stability estimate for the solution of the differential problem in all frequency regimes. It is also proved that the proposed IP-DG method is unconditionally stable with respect to both frequency $ω$ and mesh size $h$. Sub-optimal order (with respect to $h$) error estimates in the broken $H^1$-norm and in the $L^2$-norm are obtained in all mesh regimes. These estimate improve to optimal order when the mesh size $h$ is restricted to the pre-asymptotic regime (i.e., $ω^βh =O(1)$ for some $1\leq β<2$). The novelties of the proposed IP-DG method include: first, the method penalizes not only the jumps of the function values across the element edges but also the jumps of the normal derivatives; second, the penalty parameters are taken as complex numbers with positive imaginary parts. In order to establish the desired unconditional stability estimate for the numerical solution, the main idea is to exploit a (simple) property of linear functions to overcome the main difficulty caused by non-Hermitian nature and strong indefiniteness of the Helmholtz-type problem. The error estimate is then derived using a nonstandard technique adapted from \cite{Feng_Wu09}. Numerical experiments are also presented to validate the theoretical results and to numerically examine the pollution effect (with respect to $ω$) in the error bounds.

preprint2015arXiv

Analysis of mixed interior penalty discontinuous Galerkin methods for the Cahn-Hilliard equation and the Hele-Shaw flow

This paper proposes and analyzes two fully discrete mixed interior penalty discontinuous Galerkin (DG) methods for the fourth order nonlinear Cahn-Hilliard equation. Both methods use the backward Euler method for time discretization and interior penalty discontinuous Galerkin methods for spatial discretization. They differ from each other on how the nonlinear term is treated, one of them is based on fully implicit time-stepping and the other uses the energy-splitting time-stepping. The primary goal of the paper is to prove the convergence of the numerical interfaces of the DG methods to the interface of the Hele-Shaw flow. This is achieved by establishing error estimates that depend on $ε^{-1}$ only in some low polynomial orders, instead of exponential orders. Similar to [14], the crux is to prove a discrete spectrum estimate in the discontinuous Galerkin finite element space. However, the validity of such a result is not obvious because the DG space is not a subspace of the (energy) space $H^1$ and it is larger than the finite element space. This difficult is overcome by a delicate perturbation argument which relies on the discrete spectrum estimate in the finite element space proved in \cite{Feng_Prohl04}. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete mixed DG methods.

preprint2015arXiv

Finite Element Methods for the Stochastic Allen-Cahn Equation with Gradient-type Multiplicative Noises

This paper studies finite element approximations of the stochastic Allen-Cahn equation with gradient-type multiplicative noises that are white in time and correlated in space. The sharp interface limit as the parameter $ε\rightarrow 0$ of the stochastic equation formally approximates a stochastic mean curvature flow which is described by a stochastically perturbed geometric law of the deterministic mean curvature flow. Both the stochastic Allen-Cahn equation and the stochastic mean curvature flow arise from materials science, fluid mechanics and cell biology applications. Two fully discrete finite element methods which are based on different time-stepping strategies for the nonlinear term are proposed. Strong convergence with sharp rates for both fully discrete finite element methods is proved. This is done with a crucial help of the Hölder continuity in time with respect to the spatial $L^2$-norm and $H^1$-seminorm for the strong solution of the stochastic Allen-Cahn equation, which are key technical lemmas proved in paper. It also relies on the fact that high moments of the strong solution are bounded in various spatial and temporal norms. Numerical experiments are provided to gauge the performance of the proposed fully discrete finite element methods and to study the interplay of the geometric evolution and gradient-type noises.

preprint2014arXiv

An efficient numerical method for acoustic wave scattering in random media

This paper is concerned with developing efficient numerical methods for acoustic wave scattering in random media which can be expressed as random perturbations of homogeneous media. We first analyze the random Helmholtz problem by deriving some wave-number-explicit solution estimates. We then establish a multi-modes representation of the solution as a power series of the perturbation parameter and analyze its finite modes approximations. Based on this multi-modes representation, we develop a Monte Carlo interior penalty discontinuous Galerkin (MCIP-DG) method for approximating the mode functions, which are governed by recursively defined nearly deterministic Helmholtz equations. Optimal order error estimates are derived for the method and an efficient algorithm, which is based on the LU direct solver, is also designed for efficiently implementing the proposed multi-modes MCIP-DG method. It is proved that the computational complexity of the whole algorithm is comparable to that of solving one deterministic Helmholtz problem using the LU director solver. Numerical experiments are provided to validate the theoretical results and to gauge the performance of the proposed numerical method and algorithm.

preprint2014arXiv

Analysis of interior penalty discontinuous Galerkin methods for the Allen-Cahn equation and the mean curvature flow

This paper develops and analyzes two fully discrete interior penalty discontinuous Galerkin (IP-DG) methods for the Allen-Cahn equation, which is a nonlinear singular perturbation of the heat equation and originally arises from phase transition of binary alloys in materials science, and its sharp interface limit (the mean curvature flow) as the perturbation parameter tends to zero. Both fully implicit and energy-splitting time-stepping schemes are proposed. The primary goal of the paper is to derive sharp error bounds which depend on the reciprocal of the perturbation parameter $ε$ (also called "interaction length") only in some lower polynomial order, instead of exponential order, for the proposed IP-DG methods. The derivation is based on a refinement of the nonstandard error analysis technique first introduced in [12]. The centerpiece of this new technique is to establish a spectrum estimate result in totally discontinuous DG finite element spaces with a help of a similar spectrum estimate result in the conforming finite element spaces which was established in [12]. As a nontrivial application of the sharp error estimates, they are used to establish convergence and the rates of convergence of the zero level sets of the fully discrete IP-DG solutions to the classical and generalized mean curvature flow. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete IP-DG methods.

preprint2014arXiv

Group Orbit Optimization: A Unified Approach to Data Normalization

In this paper we propose and study an optimization problem over a matrix group orbit that we call \emph{Group Orbit Optimization} (GOO). We prove that GOO can be used to induce matrix decomposition techniques such as singular value decomposition (SVD), LU decomposition, QR decomposition, Schur decomposition and Cholesky decomposition, etc. This gives rise to a unified framework for matrix decomposition and allows us to bridge these matrix decomposition methods. Moreover, we generalize GOO for tensor decomposition. As a concrete application of GOO, we devise a new data decomposition method over a special linear group to normalize point cloud data. Experiment results show that our normalization method is able to obtain recovery well from distortions like shearing, rotation and squeezing.

preprint2014arXiv

International Transmission of Shocks and Fragility of a Bank Network

The weighted and directed network of countries based on the number of overseas banks is analyzed in terms of its fragility to the banking crisis of one country. We use two different models to describe transmission of shocks, one local and the other global. Depending on the original source of the crisis, the overall size of crisis impacts is found to differ country by country. For the two-step local spreading model, it is revealed that the scale of the first impact is determined by the out-strength, the total number of overseas branches of the country at the origin of the crisis, while the second impact becomes more serious if the in-strength at the origin is increased. For the global spreading model, some countries named "triggers" are found to play important roles in shock transmission, and the importance of the feed-forward-loop mechanism is pointed out. We also discuss practical policy implications of the present work.

preprint2014arXiv

Measurement and Internalization of Systemic Risk in a Global Banking Network

The negative externalities from an individual bank failure to the whole system can be huge. One of the key purposes of bank regulation is to internalize the social costs of potential bank failures via capital charges. This study proposes a method to evaluate and allocate the systemic risk to different countries/regions using a SIR type of epidemic spreading model and the Shapley value in game theory. The paper also explores features of a constructed bank network using real globe-wide banking data.

preprint2014arXiv

Multiphysics Finite Element Methods for a Poroelasticity Model

This paper concerns with finite element approximations of a quasi-static poroelasticity model in displacement-pressure formulation which describes the dynamics of poro-elastic materials under an applied mechanical force on the boundary. To better describe the multiphysics process of deformation and diffusion for poro-elastic materials, we first present a reformulation of the original model by introducing two pseudo-pressures, one of them is shown to satisfy a diffusion equation, we then propose a time-stepping algorithm which decouples (or couples) the reformulated PDE problem at each time step into two sub-problems, one of which is a generalized Stokes problem for the displacement vector field (of the solid network of the poro-elastic material) along with one pseudo-pressure field and the other is a diffusion problem for the other pseudo-pressure field (of the solvent of the material). In the paper, the Taylor-Hood mixed finite element method combined with the $P_1$-conforming finite element method is used as an example to demonstrate the viability of the proposed multiphysics approach. It is proved that the solutions of the fully discrete finite element methods fulfill a discrete energy law which mimics the differential energy law satisfied by the PDE solution and converges optimally in the energy norm. Moreover, it is showed that the proposed formulation also has a built-in mechanism to overcome so-called "locking phenomenon" associated with the numerical approximations of the poroelasticity model. Numerical experiments are presented to show the performance of the proposed approach and methods and to demonstrate the absence of "locking phenomenon" in our numerical experiments.

preprint2013arXiv

Analysis of a Darcy-Cahn-Hilliard Diffuse Interface Model for the Hele-Shaw Flow and its Fully Discrete Finite Element Approximation

In this paper we present PDE and finite element analyses for a system of partial differential equations (PDEs) consisting of the Darcy equation and the Cahn-Hilliard equation, which arises as a diffuse interface model for the two phase Hele-Shaw flow. We propose a fully discrete implicit finite element method for approximating the PDE system, which consists of the implicit Euler method combined with a convex splitting energy strategy for the temporal discretization, the standard finite element discretization for the pressure and a split (or mixed) finite element discretization for the fourth order Cahn-Hilliard equation. It is shown that the proposed numerical method satisfies a mass conservation law in addition to a discrete energy law that mimics the basic energy law for the Darcy-Cahn-Hilliard phase field model and holds uniformly in the phase field parameter $ε$. With help of the discrete energy law, we first prove that the fully discrete finite method is unconditionally energy stable and uniquely solvable at each time step. We then show that, using the compactness method, the finite element solution has an accumulation point that is a weak solution of the PDE system. As a result, the convergence result also provides a constructive proof of the existence of global-in-time weak solutions to the Darcy-Cahn-Hilliard phase field model in both two and three dimensions. Finally, we propose a nonlinear multigrid iterative algorithm to solve the finite element equations at each time step. Numerical experiments based on the overall solution method of combining the proposed finite element discretization and the nonlinear multigrid solver are presented to validate the theoretical results and to show the effectiveness of the proposed fully discrete finite element method for approximating the Darcy-Cahn-Hilliard phase field model.

preprint2013arXiv

Convergent finite difference methods for one-dimensional fully nonlinear second order partial differential equations

This paper develops a new framework for designing and analyzing convergent finite difference methods for approximating both classical and viscosity solutions of second order fully nonlinear partial differential equations (PDEs) in 1-D. The goal of the paper is to extend the successful framework of monotone, consistent, and stable finite difference methods for first order fully nonlinear Hamilton-Jacobi equations to second order fully nonlinear PDEs such as Monge-Ampère and Bellman type equations. New concepts of consistency, generalized monotonicity, and stability are introduced; among them, the generalized monotonicity and consistency, which are easier to verify in practice, are natural extensions of the corresponding notions of finite difference methods for first order fully nonlinear Hamilton-Jacobi equations. The main component of the proposed framework is the concept of "numerical operator", and the main idea used to design consistent, monotone and stable finite difference methods is the concept of "numerical moment". These two new concepts play the same roles as the "numerical Hamiltonian" and the "numerical viscosity" play in the finite difference framework for first order fully nonlinear Hamilton-Jacobi equations. In the paper, two classes of consistent and monotone finite difference methods are proposed for second order fully nonlinear PDEs. The first class contains Lax-Friedrichs-like methods which also are proved to be stable and the second class contains Godunov-like methods. Numerical results are also presented to gauge the performance of the proposed finite difference methods and to validate the theoretical results of the paper.

preprint2013arXiv

Discontinuous Galerkin finite element differential calculus and applications to numerical solutions of linear and nonlinear partial differential equations

This paper develops a discontinuous Galerkin (DG) finite element differential calculus theory for approximating weak derivatives of Sobolev functions and piecewise Sobolev functions. By introducing numerical one-sided derivatives as building blocks, various first and second order numericaloperators such as the gradient, divergence, Hessian, and Laplacian operator are defined, and their corresponding calculus rules are established. Among the calculus rules are product and chain rules, integration by parts formulas and the divergence theorem. Approximation properties and the relationship between the proposed DG finite element numerical derivatives and some well-known finite difference numerical derivative formulas on Cartesian grids are also established. Efficient implementation of the DG finite element numerical differential operators is also proposed. Besides independent interest in numerical differentiation, the primary motivation and goal of developing the DG finite element differential calculus is to solve partial differential equations. It is shown that several existing finite element, finite difference and DG methods can be rewritten compactly using the proposed DG finite element differential calculus framework. Moreover, new DG methods for linear and nonlinear PDEs are also obtained from the framework.

preprint2013arXiv

Finite element approximations of the stochastic mean curvature flow of planar curves of graphs

This paper develops and analyzes a semi-discrete and a fully discrete finite element method for a one-dimensional quasilinear parabolic stochastic partial differential equation (SPDE) which describes the stochastic mean curvature flow for planar curves of graphs. To circumvent the difficulty caused by the low spatial regularity of the SPDE solution, a regularization procedure is first proposed to approximate the SPDE, and an error estimate for the regularized problem is derived. A semi-discrete finite element method, and a space-time fully discrete method are then proposed to approximate the solution of the regularized SPDE problem. Strong convergence with rates are established for both, semi- and fully discrete methods. Computational experiments are provided to study the interplay of the geometric evolution and gradient type-noises.

preprint2013arXiv

On Schwarz Methods for Nonsymmetric and Indefinite Problems

In this paper we introduce a new Schwarz framework and theory, based on the well-known idea of space decomposition, for nonsymmetric and indefinite linear systems arising from continuous and discontinuous Galerkin approximations of general nonsymmetric and indefinite elliptic partial differential equations. The proposed Schwarz framework and theory are presented in a variational setting in Banach spaces instead of Hilbert spaces which is the case for the well-known symmetric and positive definite (SPD) Schwarz framework and theory. Condition number estimates for the additive and hybrid Schwarz preconditioners are established. The main idea of our nonsymmetric and indefinite Schwarz framework and theory is to use weak coercivity (satisfied by the nonsymmetric and indefinite bilinear form) induced norms to replace the standard bilinear form induced norm in the SPD Schwarz framework and theory. Applications of the proposed nonsymmetric and indefinite Schwarz framework and theory. Applications of the proposed nonsymmetric and indefinite Schwarz framework to solutions of discontinuous Galerkin approximations of convection-diffusion problems are also discussed. Extensive 1-D numerical experiments are also provided to gauge the performance of the proposed Schwarz methods.

preprint2012arXiv

An absolutely stable discontinuous Galerkin method for the indefinite time-harmonic Maxwell equations with large wave number

This paper develops and analyzes an interior penalty discontinuous Galerkin (IPDG) method using piecewise linear polynomials for the indefinite time harmonic Maxwell equations with the impedance boundary condition in the three dimensional space. The main novelties of the proposed IPDG method include the following: first, the method penalizes not only the jumps of the tangential component of the electric field across the element faces but also the jumps of the tangential component of its vorticity field; second, the penalty parameters are taken as complex numbers of negative imaginary parts. For the differential problem, we prove that the sesquilinear form associated with the Maxwell problem satisfies a generalized weak stability (i.e., inf-sup condition) for star-shaped domains.Such a generalized weak stability readily infers wave-number explicit a priori estimates for the solution of the Maxwell problem, which plays an important role in the error analysis for the IPDG method. For the proposed IPDG method, we show that the discrete sesquilinear form satisfies a coercivity for all positive mesh size $h$ and wave number $k$ and for general domains including non-star-shaped ones. In turn, the coercivity easily yields the well-posedness and stability estimates (i.e., a priori estimates) for the discrete problem without imposing any mesh constraint. Based on these discrete stability estimates, by adapting a nonstandard error estimate technique of Fung and Wu (2009), we derive both the energy-norm and the $L^2$-norm error estimates for the IPDG method in all mesh parameter regimes including pre-asymptotic regime (i.e., $k^2 h\gtrsim 1$). Numerical experiments are also presented to gauge the theoretical results and to numerically examine the pollution effect (with respect to $k$) in the error bounds.

preprint2012arXiv

Convergence of a fourth order singular perturbation of the $n$-dimensional radially symmetric Monge-Ampère equation

This paper concerns with the convergence analysis of a fourth order singular perturbation of the Dirichlet Monge-Ampère problem in the $n$-dimensional radial symmetric case. A detailed study of the fourth order problem is presented. In particular, various {\em a priori} estimates with explicit dependence on the perturbation parameter $\vepsi$ are derived, and a crucial convexity property is also proved for the solution of the fourth order problem. Using these estimates and the convexity property, we prove that the solution of the perturbed problem converges uniformly and compactly to the unique convex viscosity solution of the Dirichlet Monge-Ampère problem. Rates of convergence in the $H^k$-norm for $k=0,1,2$ are established, and illustrating numerical experiment results are also presented in the paper.

preprint2012arXiv

Local discontinuous Galerkin methods for one-dimensional second order fully nonlinear elliptic and parabolic equations

This paper is concerned with developing accurate and efficient discontinuous Galerkin methods for fully nonlinear second order elliptic and parabolic partial differential equations (PDEs) in the case of one spatial dimension. The primary goal of the paper to develop a general framework for constructing high order local discontinuous Galerkin (LDG) methods for approximating viscosity solutions of these fully nonlinear PDEs which are merely continuous functions by definition. In order to capture discontinuities of the first order derivative $u_x$ of the solution $u$, two independent functions $q_1$ and $q_2$ are introduced to approximate one-sided derivatives of $u$. Similarly, to capture the discontinuities of the second order derivative $u_{xx}$, four independent functions $p_{1}$, $p_{2}$, $p_{3}$, and $p_{4}$ are used to approximate one-sided derivatives of $q_1$ and $q_2$. The proposed LDG framework, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a given fully nonlinear problem into a mostly linear system of equations where the given nonlinear differential operator must be replaced by a numerical operator which allows multiple value inputs of the first and second order derivatives $u_x$ and $u_{xx}$. An easy to verify criterion for constructing "good" numerical operators is also proposed. It consists of a consistency and a generalized monotonicity. To ensure such a generalized monotonicity, the crux of the construction is to introduce the numerical moment in the numerical operator. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes.

preprint2012arXiv

Mixed Interior Penalty Discontinuous Galerkin Methods for One-dimensional Fully Nonlinear Second Order Elliptic and Parabolic Equations

This paper is concerned with developing accurate and efficient numerical methods for one-dimensional fully nonlinear second order elliptic and parabolic partial differential equations (PDEs). In the paper we present a general framework for constructing high order interior penalty discontinuous Galerkin (IP-DG) methods for approximating viscosity solutions of these fully nonlinear PDEs. In order to capture discontinuities of the second order derivative $u_{xx}$ of the solution $u$, three independent functions $p_1, p_2$ and $p_3$ are introduced to represent numerical derivatives using various one-sided limits. The proposed DG framework, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a nonlinear problem into a mostly linear system of equations where the nonlinearity has been modified to include multiple values of the second order derivative $u_{xx}$. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes.In addition to the nonstandard mixed formulation setting, another main idea is to replace the fully nonlinear differential operator by a numerical operator, which is consistent with the differential operator and satisfies certain monotonicity (called g-monotonicity) properties. To ensure such a g-monotonicity, the crux of the construction is to introduce the numerical moment, which plays a critical role in the proposed DG framework. This paper also presents and analyzes numerical results for several numerical test problems which are used to guage the accuracy and efficiency of the proposed DG methods.

preprint2011arXiv

The Vanishing Moment Method for Fully Nonlinear Second Order Partial Differential Equations: Formulation, Theory, and Numerical Analysis

The vanishing moment method was introduced by the authors in [37] as a reliable methodology for computing viscosity solutions of fully nonlinear second order partial differential equations (PDEs), in particular, using Galerkin-type numerical methods such as finite element methods, spectral methods, and discontinuous Galerkin methods, a task which has not been practicable in the past. The crux of the vanishing moment method is the simple idea of approximating a fully nonlinear second order PDE by a family (parametrized by a small parameter $\vepsi$) of quasilinear higher order (in particular, fourth order) PDEs. The primary objectives of this book are to present a detailed convergent analysis for the method in the radial symmetric case and to carry out a comprehensive finite element numerical analysis for the vanishing moment equations (i.e., the regularized fourth order PDEs). Abstract methodological and convergence analysis frameworks of conforming finite element methods and mixed finite element methods are first developed for fully nonlinear second order PDEs in general settings. The abstract frameworks are then applied to three prototypical nonlinear equations, namely, the Monge-Ampère equation, the equation of prescribed Gauss curvature, and the infinity-Laplacian equation. Numerical experiments are also presented for each problem to validate the theoretical error estimate results and to gauge the efficiency of the proposed numerical methods and the vanishing moment methodology.

preprint2010arXiv

Absolutely stable local discontinuous Galerkin methods for the Helmholtz equation with large wave number

Two local discontinuous Galerkin (LDG) methods using some non-standard numerical fluxes are developed for the Helmholtz equation with the first order absorbing boundary condition in the high frequency regime. It is shown that the proposed LDG methods are absolutely stable (hence well-posed) with respect to both the wave number and the mesh size. Optimal order (with respect to the mesh size) error estimates are proved for all wave numbers in the preasymptotic regime. To analyze the proposed LDG methods, they are recasted and treated as (non-conforming) mixed finite element methods. The crux of the analysis is to establish a generalized {\em inf-sup} condition, which holds without any mesh constraint, for each LDG method. The generalized {\em inf-sup} conditions then easily infer the desired absolute stability of the proposed LDG methods. In return, the stability results not only guarantee the well-posedness of the LDG methods but also play a crucial role in the derivation of the error estimates. Numerical experiments, which confirm the theoretical results and compare the proposed two LDG methods, are also presented in the paper.

preprint2008arXiv

Structural relaxation effects on interface and transport properties of Fe/MgO(001) tunnel junctions

The interface structure of Fe/MgO(100) magnetic tunnel junctions predicted by density functional theory (DFT) depends significantly on the choice of exchange and correlation functional. Bader analysis reveals that structures obtained by relaxing the cell with the local spin-density approximation (LSDA) display a different charge transfer than those relaxed with the generalized gradient approximation (GGA). As a consequence, the electronic transport is found to be extremely sensitive to the interface structure. In particular, the conductance for the LSDA-relaxed geometry is about one order of magnitude smaller than that of the GGA-relaxed one. The high sensitivity of the electronic current to the details of the interface might explain the discrepancy between the experimental and calculated values of magnetoresistance.