Source author record

Min Tang

Min Tang 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

27works
17topics
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

27 published item(s)

preprint2026arXiv

LFD: Layer Fused Decoding to Exploit External Knowledge in Retrieval-Augmented Generation

Retrieval-augmented generation (RAG) incorporates external knowledge into large language models (LLMs), improving their adaptability to downstream tasks and enabling information updates. Surprisingly, recent empirical evidence demonstrates that injecting noise into retrieved relevant documents paradoxically facilitates exploitation of external knowledge and improves generation quality. Although counterintuitive and challenging to apply in practice, this phenomenon enables granular control and rigorous analysis of how LLMs integrate external knowledge. Therefore, in this paper, we intervene on noise injection and establish a layer-specific functional demarcation within the LLM: shallow layers specialize in local context modeling, intermediate layers focus on integrating long-range external factual knowledge, and deeper layers primarily rely on parametric internal knowledge. Building on this insight, we propose Layer Fused Decoding (LFD), a simple decoding strategy that directly combines representations from an intermediate layer with final-layer decoding outputs to fully exploit the external factual knowledge. To identify the optimal intermediate layer, we introduce an internal knowledge score (IKS) criterion that selects the layer with the lowest IKS value in the latter half of layers. Experimental results across multiple benchmarks demonstrate that LFD helps RAG systems more effectively surface retrieved context knowledge with minimal cost.

preprint2024arXiv

A fast offline/online forward solver for stationary transport equation with multiple inflow boundary conditions and varying coefficients

It is of great interest to solve the inverse problem of stationary radiative transport equation (RTE) in optical tomography. The standard way is to formulate the inverse problem into an optimization problem, but the bottleneck is that one has to solve the forward problem repeatedly, which is time-consuming. Due to the optical property of biological tissue, in real applications, optical thin and thick regions coexist and are adjacent to each other, and the geometry can be complex. To use coarse meshes and save the computational cost, the forward solver has to be asymptotic preserving across the interface (APAL). In this paper, we propose an offline/online solver for RTE. The cost at the offline stage is comparable to classical methods, while the cost at the online stage is much lower. Two cases are considered. One is to solve the RTE with fixed scattering and absorption cross sections while the boundary conditions vary; the other is when cross sections vary in a small domain and the boundary conditions change many times. The solver can be decomposed into offline/online stages in these two cases. One only needs to calculate the offline stage once and update the online stage when the parameters vary. Our proposed solver is much cheaper when one needs to solve RTE with multiple right-hand sides or when the cross sections vary in a small domain, thus can accelerate the speed of solving inverse RTE problems. We illustrate the online/offline decomposition based on the Tailored Finite Point Method (TFPM), which is APAL on general quadrilateral meshes.

preprint2024arXiv

Confined run-and-tumble model with boundary aggregation: long time behavior and convergence to the confined Fokker-Planck model

The motile micro-organisms such as E. coli, sperm, or some seaweed are usually modelled by self-propelled particles that move with the run-and-tumble process. Individual-based stochastic models are usually employed to model the aggregation phenomenon at the boundary, which is an active research field that has attracted a lot of biologists and biophysicists. Self-propelled particles at the microscale have complex behaviors, while characteristics at the population level are more important for practical applications but rely on individual behaviors. Kinetic PDE models that describe the time evolution of the probability density distribution of the motile micro-organisms are widely used. However, how to impose the appropriate boundary conditions that take into account the boundary aggregation phenomena is rarely studied. In this paper, we propose the boundary conditions for a 2D confined run-and-tumble model (CRTM) for self-propelled particle populations moving between two parallel plates with a run-and-tumble process. The proposed model satisfies the relative entropy inequality and thus long-time convergence. We establish the relation between CRTM and the confined Fokker-Planck model (CFPM) studied in [22]. We prove theoretically that when the tumble is highly forward peaked and frequent enough, CRTM converges asymptotically to the CFPM. A numerical comparison of the CRTM with aggregation and CFPM is given. The time evolution of both the deterministic PDE model and individual-based stochastic simulations are displayed, which match each other well.

preprint2024arXiv

Reconstructing the kinetic chemotaxis kernel using macroscopic data: well-posedness and ill-posedness

Bacterial motion is steered by external stimuli (chemotaxis), and the motion described on the mesoscopic scale is uniquely determined by a parameter $K$ that models velocity change response from the bacteria. This parameter is called chemotaxis kernel. In a practical setting, it is inferred by experimental data. We deploy a PDE-constrained optimization framework to perform this reconstruction using velocity-averaged, localized data taken in the interior of the domain. The problem can be well-posed or ill-posed depending on the data preparation and the experimental setup. In particular, we propose one specific design that guarantees numerical reconstructability and local convergence. This design is adapted to the discretization of $K$ in space and decouples the reconstruction of local values of $K$ into smaller cell problems, opening up parallelization opportunities. Numerical evidences support the theoretical findings.

preprint2022arXiv

N-Cloth: Predicting 3D Cloth Deformation with Mesh-Based Networks

We present a novel mesh-based learning approach (N-Cloth) for plausible 3D cloth deformation prediction. Our approach is general and can handle cloth or obstacles represented by triangle meshes with arbitrary topologies. We use graph convolution to transform the cloth and object meshes into a latent space to reduce the non-linearity in the mesh space. Our network can predict the target 3D cloth mesh deformation based on the initial state of the cloth mesh template and the target obstacle mesh. Our approach can handle complex cloth meshes with up to 100K triangles and scenes with various objects corresponding to SMPL humans, non-SMPL humans or rigid bodies. In practice, our approach can be used to generate plausible cloth simulation at 30-45 fps on an NVIDIA GeForce RTX 3090 GPU. We highlight its benefits over prior learning-based methods and physically-based cloth simulators.

preprint2021arXiv

A Spatial-Temporal asymptotic preserving scheme for radiation magnetohydrodynamics in the equilibrium and non-equilibrium diffusion limit

The radiation magnetohydrodynamics (RMHD) system couples the ideal magnetohydrodynamics equations with a gray radiation transfer equation. The main challenge is that the radiation travels at the speed of light while the magnetohydrodynamics changes with the time scale of the fluid. The time scales of these two processes can vary dramatically. In order to use mesh sizes and time steps that are independent of the speed of light, asymptotic preserving (AP) schemes in both space and time are desired. In this paper, we develop an AP scheme in both space and time for the RMHD system. Two different scalings are considered. One results in an equilibrium diffusion limit system, while the other results in a non-equilibrium system. The main idea is to decompose the radiative intensity into three parts, each part is treated differently with suitable combinations of explicit and implicit discretizations guaranteeing the favorable stability conditionand computational efficiency. The performance of the AP method is presented, for both optically thin and thick regions, as well as for the radiative shock problem.

preprint2021arXiv

Multi-scale GCN-assisted two-stage network for joint segmentation of retinal layers and disc in peripapillary OCT images

An accurate and automated tissue segmentation algorithm for retinal optical coherence tomography (OCT) images is crucial for the diagnosis of glaucoma. However, due to the presence of the optic disc, the anatomical structure of the peripapillary region of the retina is complicated and is challenging for segmentation. To address this issue, we developed a novel graph convolutional network (GCN)-assisted two-stage framework to simultaneously label the nine retinal layers and the optic disc. Specifically, a multi-scale global reasoning module is inserted between the encoder and decoder of a U-shape neural network to exploit anatomical prior knowledge and perform spatial reasoning. We conducted experiments on human peripapillary retinal OCT images. The Dice score of the proposed segmentation network is 0.820$\pm$0.001 and the pixel accuracy is 0.830$\pm$0.002, both of which outperform those from other state-of-the-art techniques.

preprint2020arXiv

Hierarchical Optimization Time Integration for CFL-rate MPM Stepping

We propose Hierarchical Optimization Time Integration (HOT) for efficient implicit time-stepping of the Material Point Method (MPM) irrespective of simulated materials and conditions. HOT is an MPM-specialized hierarchical optimization algorithm that solves nonlinear time step problems for large-scale MPM systems near the CFL-limit. HOT provides convergent simulations "out-of-the-box" across widely varying materials and computational resolutions without parameter tuning. As an implicit MPM time stepper accelerated by a custom-designed Galerkin multigrid wrapped in a quasi-Newton solver, HOT is both highly parallelizable and robustly convergent. As we show in our analysis, HOT maintains consistent and efficient performance even as we grow stiffness, increase deformation, and vary materials over a wide range of finite strain, elastodynamic and plastic examples. Through careful benchmark ablation studies, we compare the effectiveness of HOT against seemingly plausible alternative combinations of MPM with standard multigrid and other Newton-Krylov models. We show how these alternative designs result in severe issues and poor performance. In contrast, HOT outperforms the existing state-of-the-art, heavily optimized implicit MPM codes with an up to 10x performance speedup across a wide range of challenging benchmark test simulations.

preprint2020arXiv

On a subset sums problem of Chen and Wu

For a set $A$, let $P(A)$ be the set of all finite subset sums of $A$. We prove that if a sequence $B=\{11\leq b_1<b_2<\cdots\}$ satisfies $b_2=3b_1+5$, $b_3=3b_2+2$ and $b_{n+1}=3b_n+4b_{n-1}$ for all $n\geq 3$, then there is a sequence of positive integers $A=\{a_1<a_2<\cdots\}$ such that $P(A)=\mathbb{N}\setminus B$. This result shows that the answer to the problem of Chen and Wu [`The inverse problem on subset sums', European. J. Combin. 34(2013), 841-845] is negative.

preprint2020arXiv

On an inverse problem in additive number theory

For a set $A$, let $P(A)$ be the set of all finite subset sums of $A$. In this paper, for a sequence of integers $B=\{1<b_1<b_2<\cdots\}$ and $3b_1+5\leq b_2\leq 6b_1+10$, we determine the critical value for $b_3$ such that there exists an infinite sequence $A$ of positive integers for which $P(A)=\mathbb{N}\setminus B$. This result shows that we partially solve the problem of Fang and Fang [`On an inverse problem in additive number theory', Acta Math. Hungar. 158(2019), 36-39].

preprint2020arXiv

P-Cloth: Interactive Complex Cloth Simulation on Multi-GPU Systems using Dynamic Matrix Assembly and Pipelined Implicit Integrators

We present a novel parallel algorithm for cloth simulation that exploits multiple GPUs for fast computation and the handling of very high resolution meshes. To accelerate implicit integration, we describe new parallel algorithms for sparse matrix-vector multiplication (SpMV) and for dynamic matrix assembly on a multi-GPU workstation. Our algorithms use a novel work queue generation scheme for a fat-tree GPU interconnect topology. Furthermore, we present a novel collision handling scheme that uses spatial hashing for discrete and continuous collision detection along with a non-linear impact zone solver. Our parallel schemes can distribute the computation and storage overhead among multiple GPUs and enable us to perform almost interactive simulation on complex cloth meshes, which can hardly be handled on a single GPU due to memory limitations. We have evaluated the performance with two multi-GPU workstations (with 4 and 8 GPUs, respectively) on cloth meshes with 0.5-1.65M triangles. Our approach can reliably handle the collisions and generate vivid wrinkles and folds at 2-5 fps, which is significantly faster than prior cloth simulation systems. We observe almost linear speedups with respect to the number of GPUs.

preprint2016arXiv

An Asymptotic Preserving method for strongly anisotropic diffusion equations based on field line integration

In magnetized plasma, the magnetic field confines the particles around the field lines. The anisotropy intensity in the viscosity and heat conduction may reach the order of $10^{12}$. When the boundary conditions are periodic or Neumann, the strong diffusion leads to an ill-posed limiting problem. To remove the ill-conditionedness in the highly anisotropic diffusion equations, we introduce a simple but very efficient asymptotic preserving reformulation in this paper. The key idea is that, instead of discretizing the Neumann boundary conditions locally, we replace one of the Neumann boundary condition by the integration of the original problem along the field line, the singular $1/ε$ terms can be replaced by $O(1)$ terms after the integration, so that yields a well-posed problem. Small modifications to the original code are required and no change of coordinates nor mesh adaptation are needed. Uniform convergence with respect to the anisotropy strength $1/ε$ can be observed numerically and the condition number does not scale with the anisotropy.

preprint2016arXiv

Macroscopic limits of pathway-based kinetic models for E.coli chemotaxis in large gradient environments

It is of great biological interest to understand the molecular origins of chemotactic behavior of E. coli by developing population-level models based on the underlying signaling pathway dynamics. We derive macroscopic models for E.coli chemotaxis that match quantitatively with the agent-based model (SPECS) for all ranges of the spacial gradient, in particular when the chemical gradient is large such that the standard Keller-Segel model is no longer valid. These equations are derived both formally and rigorously as asymptotic limits for pathway-based kinetic equations. We also present numerical results that show good agreement between the macroscopic models and SPECS. Our work provides an answer to the question of how to determine the population-level diffusion coefficient and drift velocity from the molecular mechanisms of chemotaxis, for both shallow gradients and large gradients environments.

preprint2016arXiv

Uniform convergent scheme for strongly anisotropic diffusion equations with closed field lines

In magnetized plasma, the magnetic field confines particles around field lines. The ratio between the intensity of the parallel and perpendicular viscosity or heat conduction may reach the order of $10^{12}$. When the magnetic fields have closed field lines and form a "magnetic island", the convergence order of most known schemes depends on the anisotropy strength. In this paper, by integration of the original differential equation along each closed field line, we introduce a simple but very efficient asymptotic preserving reformulation, which yields uniform convergence with respect to the anisotropy strength. Only slight modification to the original code is required and neither change of coordinates nor mesh adaptation is needed. Numerical examples demonstrating the performance of the new scheme are presented.

preprint2016arXiv

Well-balanced and asymptotic preserving schemes for kinetic models

In this paper, we propose a general framework for designing numerical schemes that have both well-balanced (WB) and asymptotic preserving (AP) properties, for various kinds of kinetic models. We are interested in two different parameter regimes, 1) When the ratio between the mean free path and the characteristic macroscopic length $ε$ tends to zero, the density can be described by (advection) diffusion type (linear or nonlinear) macroscopic models; 2) When $ε$ = O(1), the models behave like hyperbolic equations with source terms and we are interested in their steady states. We apply the framework to three different kinetic models: neutron transport equation and its diffusion limit, the transport equation for chemotaxis and its Keller-Segel limit, and grey radiative transfer equation and its nonlinear diffusion limit. Numerical examples are given to demonstrate the properties of the schemes.

preprint2015arXiv

Derivation of the bacterial run-and-tumble kinetic equation from a model with biochemical pathway

Kinetic-transport equations are, by now, standard models to describe the dynamics of populations of bacteria moving by run-and-tumble. Experimental observations show that bacteria increase their run duration when encountering an increasing gradient of chemotactic molecules. This led to a first class of models which heuristically include tumbling frequencies depending on the path-wise gradient of chemotactic signal. More recently, the biochemical pathways regulating the flagellar motors were uncovered. This knowledge gave rise to a second class of kinetic-transport equations, that takes into account an intra-cellular molecular content and which relates the tumbling frequency to this information. It turns out that the tumbling frequency depends on the chemotactic signal, and not on its gradient. For these two classes of models, macroscopic equations of Keller-Segel type, have been derived using diffusion or hyperbolic rescaling. We complete this program by showing how the first class of equations can be derived from the second class with molecular content after appropriate rescaling. The main difficulty is to explain why the path-wise gradient of chemotactic signal can arise in this asymptotic process. Randomness of receptor methylation events can be included, and our approach can be used to compute the tumbling frequency in presence of such a noise.

preprint2015arXiv

PolyDepth: Real-time Penetration Depth Computation using Iterative Contact-Space Projection

We present a real-time algorithm that finds the Penetration Depth (PD) between general polygonal models based on iterative and local optimization techniques. Given an in-collision configuration of an object in configuration space, we find an initial collision-free configuration using several methods such as centroid difference, maximally clear configuration, motion coherence, random configuration, and sampling-based search. We project this configuration on to a local contact space using a variant of continuous collision detection algorithm and construct a linear convex cone around the projected configuration. We then formulate a new projection of the in-collision configuration onto the convex cone as a Linear Complementarity Problem (LCP), which we solve using a type of Gauss-Seidel iterative algorithm. We repeat this procedure until a locally optimal PD is obtained. Our algorithm can process complicated models consisting of tens of thousands triangles at interactive rates.

preprint2014arXiv

Derivation of a Hele-Shaw type system from a cell model with active motion

We formulate a Hele-Shaw type free boundary problem for a tumor growing under the combined effects of pressure forces, cell multiplication and active motion, the latter being the novelty of the present paper. This new ingredient is considered here as a standard diffusion process. The free boundary model is derived from a description at the cell level using the asymptotic of a stiff pressure limit. Compared to the case when active motion is neglected, the pressure satisfies the same complementarity Hele-Shaw type formula. However, the cell density is smoother (Lipschitz continuous), while there is a deep change in the free boundary velocity, which is no longer given by the gradient of the pressure, because some kind of 'mushy region' prepares the tumor invasion.

preprint2014arXiv

Traveling wave solution of the Hele-Shaw model of tumor growth with nutrient

Several mathematical models of tumor growth are now commonly used to explain medical observations and predict cancer evolution based on images. These models incorporate mechanical laws for tissue compression combined with rules for nutrients availability which can differ depending on the situation under consideration, in vivo or in vitro. Numerical solutions exhibit, as expected from medical observations, a proliferative rim and a necrotic core. However, their precise profiles are rather complex, both in one and two dimensions.

preprint2013arXiv

A pathway-based mean-field model for E. coli chemotaxis: Mathematical derivation and Keller-Segel limit

A pathway-based mean-field theory (PBMFT) was recently proposed for E. coli chemotaxis in [G. Si, T. Wu, Q. Quyang and Y. Tu, Phys. Rev. Lett., 109 (2012), 048101]. In this paper, we derived a new moment system of PBMFT by using the moment closure technique in kinetic theory under the assumption that the methylation level is locally concentrated. The new system is hyperbolic with linear convection terms. Under certain assumptions, the new system can recover the original model. Especially the assumption on the methylation difference made there can be understood explicitly in this new moment system. We obtain the Keller-Segel limit by taking into account the different physical time scales of tumbling, adaptation and the experimental observations. We also present numerical evidence to show the quantitative agreement of the moment system with the individual based E. coli chemotaxis simulator.

preprint2013arXiv

Jeśmanowicz' conjecture and Fermat numbers

Let $a,b,c$ be relatively prime positive integers such that $a^{2}+b^{2}=c^{2}.$ In 1956, Jeśmanowicz conjectured that for any positive integer $n$, the only solution of $(an)^{x}+(bn)^{y}=(cn)^{z}$ in positive integers is $(x,y,z)=(2,2,2)$. Let $k\geq 1$ be an integer and $F_k=2^{2^k}+1$ be a Fermat number. In this paper, we show that Jeśmanowicz' conjecture is true for Pythagorean triples $(a,b,c)=(F_k-2,2^{2^{k-1}+1},F_k)$.

preprint2010arXiv

Can a traveling wave connect two unstable states? The case of the nonlocal Fisher equation

This note investigates the properties of the traveling waves solutions of the nonlocal Fisher equation. The existence of such solutions has been proved recently in \cite{BNPR} but their asymptotic behavior was still unclear. We use here a new numerical approximation of these traveling waves which shows that some traveling waves connect the two homogeneous steady states $0$ and $1$, which is a striking fact since $0$ is dynamically unstable and $1$ is unstable in the sense of Turing.

preprint2010arXiv

Traveling plateaus for a hyperbolic Keller-Segel system with attraction and repulsion: existence and branching instabilities

How can repulsive and attractive forces, acting on a conservative system, create stable traveling patterns or branching instabilities? We have proposed to study this question in the framework of the hyperbolic Keller-Segel system with logistic sensitivity. This is a model system motivated by experiments on cell communities auto-organization, a field which is also called socio-biology. We continue earlier modeling work, where we have shown numerically that branching patterns arise for this system and we have analyzed this instability by formal asymptotics for small diffusivity of the chemo-repellent. Here we are interested in the more general situation, where the diffusivities of both the chemo-attractant and the chemo-repellent are positive. To do so, we develop an appropriate functional analysis framework. We apply our method to two cases. Firstly we analyze steady states. Secondly we analyze traveling waves when neglecting the degradation coefficient of the chemo-repellent; the unique wave speed appears through a singularity cancelation which is the main theoretical difficulty. This shows that in different situations the cell density takes the shape of a plateau. The existence of steady states and traveling plateaus are a symptom of how rich the system is and why branching instabilities can occur. Numerical tests show that large plateaus may split into smaller ones, which remain stable.

preprint2009arXiv

All speed scheme for the low mach number limit of the Isentropic Euler equation

An all speed scheme for the Isentropic Euler equation is presented in this paper. When the Mach number tends to zero, the compressible Euler equation converges to its incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficulty in the low Mach regime. The key idea of our all speed scheme is the special semi-implicit time discretization, in which the low Mach number stiff term is divided into two parts, one being treated explicitly and the other one implicitly. Moreover, the flux of the density equation is also treated implicitly and an elliptic type equation is derived to obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared with previous semi-implicit methods, nonphysical oscillations can be suppressed. We develop this semi-implicit time discretization in the framework of a first order local Lax-Friedrich (LLF) scheme and numerical tests are displayed to demonstrate its performances.