Source author record

Jingzhi Li

Jingzhi Li 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

23works
9topics
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

23 published item(s)

preprint2026arXiv

Post-hoc Provider Fairness Adaptation via Hierarchical Exposure Alignment

Provider exposure fairness is crucial for sustaining a healthy content ecosystem and preventing monopolization in recommender systems. Yet, most existing methods either incorporate fairness constraints during model training, requiring expensive retraining when fairness objectives change, or rely on post-hoc reranking with fixed criteria, which lacks adaptability to diverse fairness requirements. To overcome these limitations, we propose Post-hoc Fairness Adaptation (PFA), a lightweight framework that equips a frozen recommender with a fairness adapter, enabling flexible fairness control without retraining the backbone model. Specifically, the fairness adapter learns personalized additive score adjustments from user-item embeddings, which are injected into the original ranking scores to steer provider exposure toward fairness. To train the adapter, we minimize the KL divergence between the actual and the target fair exposure distributions. However, this global objective implicitly treats all providers equally, ignoring structural disparities such as imbalanced provider group sizes and heterogeneous exposure within groups. Consequently, fairness may appear satisfied at an aggregate level while severe inter-group and intra-group exposure imbalances persist, undermining practical fairness. To address this, we design Hierarchical Exposure Fairness Alignment (HEFA), which explicitly balances inter- and intra-group provider exposure disparities, enabling flexible adaptation to diverse fairness requirements. To mitigate potential accuracy degradation, PFA jointly optimizes HEFA with a differentiable NDCG loss, enabling end-to-end fairness optimization while preserving ranking quality. Extensive experiments on three public datasets demonstrate that PFA achieves substantial fairness gains with negligible accuracy loss, consistently outperforming strong baselines.

preprint2024arXiv

Spatiotemporal Monitoring of Epidemics via Solution of a Coefficient Inverse Problem

Let S,I and R be susceptible, infected and recovered populations in a city affected by an epidemic. The SIR model of Lee, Liu, Tembine, Li and Osher, \emph{SIAM J. Appl. Math.},~81, 190--207, 2021 of the spatiotemoral spread of epidemics is considered. This model consists of a system of three nonlinear coupled parabolic Partial Differential Equations with respect to the space and time dependent functions S,I and R. For the first time, a Coefficient Inverse Problem (CIP) for this system is posed. The so-called \textquotedblleft convexification" numerical method for this inverse problem is constructed. The presence of the Carleman Weight Function (CWF) in the resulting regularization functional ensures the global convergence of the gradient descent method of the minimization of this functional to the true solution of the CIP, as long as the noise level tends to zero. The CWF is the function, which is used as the weight in the Carleman estimate for the corresponding Partial Differential Operator. Numerical studies demonstrate an accurate reconstruction of unknown coefficients as well as S,I,R functions inside of that city. As a by-product, uniqueness theorem for this CIP is proven. Since the minimal measured input data are required, then the proposed methodology has a potential of a significant decrease of the cost of monitoring of epidemics.

preprint2022arXiv

A $C^{0}$ interior penalty method for $m$th-Laplace equation

In this paper, we propose a $C^{0}$ interior penalty method for $m$th-Laplace equation on bounded Lipschitz polyhedral domain in $\mathbb{R}^{d}$, where $m$ and $d$ can be any positive integers. The standard $H^{1}$-conforming piecewise $r$-th order polynomial space is used to approximate the exact solution $u$, where $r$ can be any integer greater than or equal to $m$. Unlike the interior penalty method in [T.~Gudi and M.~Neilan, {\em An interior penalty method for a sixth-order elliptic equation}, IMA J. Numer. Anal., \textbf{31(4)} (2011), pp. 1734--1753], we avoid computing $D^{m}$ of numerical solution on each element and high order normal derivatives of numerical solution along mesh interfaces. Therefore our method can be easily implemented. After proving discrete $H^{m}$-norm bounded by the natural energy semi-norm associated with our method, we manage to obtain stability and optimal convergence with respect to discrete $H^{m}$-norm. Numerical experiments validate our theoretical estimate.

preprint2022arXiv

An Iterative Decoupled Algorithm with Unconditional Stability for Biot Model

This paper is concerned with numerical algorithms for Biot model. By introducing an intermediate variable, the classical 2-field Biot model is written into a 3-field formulation. Based on such a 3-field formulation, we propose a coupled algorithm, some time-extrapolation based decoupled algorithms, and an iterative decoupled algorithm. Our focus is the analysis of the iterative decoupled algorithm. It is shown that the convergence of the iterative decoupled algorithm requires no extra assumptions on physical parameters or stabilization parameters. Numerical experiments are provided to demonstrate the accuracy and efficiency of the proposed method.

preprint2022arXiv

Convexification for a CIP for the RTE]{Convexification Numerical Method for a Coefficient Inverse Problem for the Radiative Transport Equation

An $\left( n+1\right) -$D coefficient inverse problem for the radiative stationary transport equation is considered for the first time. A globally convergent so-called convexification numerical \ method is developed and its convergence analysis is provided. The analysis is based on a Carleman estimate. In particular, convergence analysis implies a certain uniqueness theorem. Extensive numerical studies in the 2-D case are presented.

preprint2021arXiv

Computation of transmission eigenvalues by the regularized Schur complement for the boundary integral operators

This paper is devoted to the computation of transmission eigenvalues in the inverse acoustic scattering theory. This problem is first reformulated as a two by two boundary system of boundary integral equations. Next, utilizing the Schur complement technique, we develop a Schur complement operator with regularization to obtain a reduced system of boundary integral equations. The Nyström discretization is then used to obtain an eigenvalue problem for a matrix. We employ the recursive integral method for the numerical computation of the matrix eigenvalue. Numerical results show that the proposed method is efficient and reduces computational costs.

preprint2020arXiv

Convexification for an Inverse Parabolic Problem

A convexification-based numerical method for a Coefficient Inverse Problem for a parabolic PDE is presented. The key element of this method is the presence of the so-called Carleman Weight Function in the numerical scheme. Convergence analysis ensures the global convergence of this method, as opposed to the local convergence of the conventional least squares minimization techniques. Numerical results demonstrate a good performance.

preprint2020arXiv

Inverse source problems in an inhomogeneous medium with a single far-field pattern

This paper concerns time-harmonic inverse source problems with a single far-field pattern in two dimensions, where the source term is compactly supported in an a priori given inhomogeneous background medium. For convex-polygonal source terms, we prove that the source support together with the zeroth and first order derivatives of the source function at corner points can be uniquely determined. Further, we prove that an admissible set of source functions (including harmonic functions) having a convex-polygonal support can be uniquely identified by a single far-field pattern. A class of radiating sources is characterized and the extension of the radiated field across a corner point is proven impossible. The corner scattering theory leads to a data-driven inversion scheme for imaging an arbitrarily convex-polygonal source support.

preprint2020arXiv

Modeling the Control of COVID-19: Impact of Policy Interventions and Meteorological Factors

In this paper, we propose a dynamical model to describe the transmission of COVID-19, which is spreading in China and many other countries. To avoid a larger outbreak in the worldwide, Chinese government carried out a series of strong strategies to prevent the situation from deteriorating. Home quarantine is the most important one to prevent the spread of COVID-19. In order to estimate the effect of population quarantine, we divide the population into seven categories for simulation. Based on a Least-Squares procedure and officially published data, the estimation of parameters for the proposed model is given. Numerical simulations show that the proposed model can describe the transmission of COVID-19 accurately, the corresponding prediction of the trend of the disease is given. The home quarantine strategy plays an important role in controlling the disease spread and speeding up the decline of COVID-19. The control reproduction number of most provinces in China are analyzed and discussed adequately. We should pay attention to that, though the epidemic is in decline in China, the disease still has high risk of human-to-human transmission continuously. Once the control strategy is removed, COVID-19 may become a normal epidemic disease just like flu. Further control for the disease is still necessary, we focus on the relationship between the spread rate of the virus and the meteorological conditions. A comprehensive meteorological index is introduced to represent the impact of meteorological factors on both high and low migration groups. As the progress on the new vaccine, we design detail vaccination strategies for COVID-19 in different control phases and show the effectiveness of efficient vaccination. Once the vaccine comes into use, the numerical simulation provide a promptly prospective research.

preprint2019arXiv

Determining a random Schrödinger equation with unknown source and potential

We are concerned with the direct and inverse scattering problems associated with a time-harmonic random Schrödinger equation with unknown source and potential terms. The well-posedness of the direct scattering problem is first established. Three uniqueness results are then obtained for the corresponding inverse problems in determining the variance of the source, the potential and the expectation of the source, respectively, by the associated far-field measurements. First, a single realization of the passive scattering measurement can uniquely recover the variance of the source without the a priori knowledge of the other unknowns. Second, if active scattering measurement can be further obtained, a single realization can uniquely recover the potential function without knowing the source. Finally, both the potential and the first two statistic moments of the random source can be uniquely recovered with full measurement data. The major novelty of our study is that on the one hand, both the random source and the potential are unknown, and on the other hand, both passive and active scattering measurements are used for the recovery in different scenarios.

preprint2019arXiv

Uniqueness in inverse cavity scattering problems with phaseless near-field data

This paper is concerned with the uniqueness of inverse acoustic scattering problem for cavities with the modulus of the near-fields. With the aid of the reference ball technique and the superpositions of two point sources as the incident waves, we rigorously prove that the location and shape of the cavity as well as its boundary condition can be uniquely determined by the modulus of near-fields at an admissible surface. To our knowledge, this is the first uniqueness result in inverse cavity scattering problems with phaseless near-field data. In this paper, we make use of the phaseless near-field data incurred by the cavity and the point sources, and thus the configuration is more feasible in practice.

preprint2016arXiv

Mathematical design of a novel input/instruction device using a moving emitter

This paper is concerned with the mathematical design of a novel input/instruction device using a moving emitter. The emitter generates a point source and can be installed on a digit pen or worn on the finger of the human being who wants to interact/communicate with the computer. The input/instruction can be recognized by identifying the motion trajectory of the emitter performed by the human being from the collected wave field data. The identification process is modelled as an in- verse source problem where one intends to identify the trajectory of a moving point source. There are several salient features of our study which distinguish our result from the existing ones in the literature. First, the point source is moving in an inhomogeneous background medium, which models the human body. Second, the dynamical wave field data are collected in a limited aperture. Third, the recon- struction method is independent of the background medium, and it is totally direct without any matrix inversion. Hence, it is efficient and robust with respect to the measurement noise. Both theoretical justifications and computational experiments are presented to verify our novel findings.

preprint2016arXiv

Numerical solution of an ill-posed Cauchy problem for a quasilinear parabolic equation using a Carleman weight function

This is the first publication in which an ill-posed Cauchy problem for a quasi- linear PDE is solved numerically by a rigorous method. More precisely, we solve the side Cauchy problem for a 1-d quasilinear parabolc equation. The key idea is to minimize a strictly convex cost functional with the Carleman Weight Function in it. Previous publications about numerical solutions of ill-posed Cauchy problems were considering only linear equations.

preprint2015arXiv

On quasi-static cloaking due to anomalous localized resonance in $\mathbb{R}^3$

This work concerns the cloaking due to anomalous localized resonance (CALR) in the quasi-static regime. We extend the related two-dimensional studies in [2,10] to the three-dimensional setting. CALR is shown not to take place for the plasmonic configuration considered in [2,10] in the three-dimensional case. We give two different constructions which ensure the occurrence of CALR. There may be no core or an arbitrary shape core for the cloaking device. If there is a core, then the dielectric distribution inside it could be arbitrary.

preprint2015arXiv

Recovering a polyhedral obstacle by a few backscattering measurements

We propose an inverse scattering scheme of recovering a polyhedral obstacle in $\mathbb{R}^n$, $n=2,3$, by only a few high-frequency acoustic backscattering measurements. The obstacle could be sound-soft or sound-hard. It is shown that the modulus of the far-field pattern in the backscattering aperture possesses a certain local maximum behavior, from which one can determine the exterior normal directions of the front sides/faces. Then by using the phaseless backscattering data corresponding to a few incident plane waves with suitably chosen incident directions, one can determine the exterior unit normal vector of each side/face of the obstacle. After the determination of the exterior unit normals, the recovery is reduced to a finite-dimensional problem of determining a location point of the obstacle and the distance of each side/face away from the location point. For the latter reconstruction, we need make use of the far-field data with phases. Numerical experiments are also presented to illustrate the effectiveness of the proposed scheme.

preprint2014arXiv

First order least squares method with weakly imposed boundary condition for convection dominated diffusion problems

We present and analyze a first order least squares method for convection dominated diffusion problems, which provides robust L2 a priori error estimate for the scalar variable even if the given data f in L2 space. The novel theoretical approach is to rewrite the method in the framework of discontinuous Petrov - Galerkin (DPG) method, and then show numerical stability by using a key equation discovered by J. Gopalakrishnan and W. Qiu [Math. Comp. 83(2014), pp. 537-552]. This new approach gives an alternative way to do numerical analysis for least squares methods for a large class of differential equations. We also show that the condition number of the global matrix is independent of the diffusion coefficient. A key feature of the method is that there is no stabilization parameter chosen empirically. In addition, Dirichlet boundary condition is weakly imposed. Numerical experiments verify our theoretical results and, in particular, show our way of weakly imposing Dirichlet boundary condition is essential to the design of least squares methods - numerical solutions on subdomains away from interior layers or boundary layers have remarkable accuracy even on coarse meshes, which are unstructured quasi-uniform.

preprint2014arXiv

Robust a Posteriori Error Estimates for HDG method for Convection-Diffusion Equations

We propose a robust a posteriori error estimator for the hybridizable discontinuous Galerkin (HDG) method for convection-diffusion equations with dominant convection. The reliability and efficiency of the estimator are established for the error measured in an energy norm. The energy norm is uniformly bounded even when the diffusion coefficient tends to zero. The estimators are robust in the sense that the upper and lower bounds of error are uniformly bounded with respect to the diffusion coefficient. A weighted test function technique and the Oswald interpolation are key ingredients in the analysis. Numerical results verify the robustness of the proposed a posteriori error estimator. In numerical experiments, optimal convergence is observed.

preprint2013arXiv

Inverse Elastic Scattering for Multiscale Rigid Bodies with A Single Far-field Pattern

We develop three inverse elastic scattering schemes for locating multiple small, extended and multiscale rigid bodies, respectively. There are some salient and promising features of the proposed methods. The cores of those schemes are certain indicator functions, which are obtained by using only a single far-field pattern of the pressure (longitudinal) wave, or the shear (transversal) wave, or the total wave field. Though the inverse scattering problem is known to be nonlinear and ill-posed, the proposed reconstruction methods are totally "direct" and there are no inversions involved. Hence, the methods are very efficient and robust against noisy data. Both rigorous mathematical justifications and numerical simulations are presented in our study.

preprint2013arXiv

Locating Multiple Multi-scale Electromagnetic Scatterers by A Single Far-field Measurement

Two inverse scattering schemes were recently developed in \cite{LiLiuShangSun} for locating multiple electromagnetic (EM) scatterers, respectively, of small size and regular size compared to the detecting EM wavelength. Both schemes make use of a single far-field measurement. The scheme of locating regular-size scatterers requires the {\it a priori} knowledge of the possible shapes, orientations and sizes of the underlying scatterer components. In this paper, we extend that imaging scheme to a much more practical setting by relaxing the requirement on the orientations and sizes. We also develop an imaging scheme of locating multiple multi-scale EM scatterers, which may include at the same time, both components of regular size and small size. For the second scheme, a novel local re-sampling technique is developed. Furthermore, more robust and accurate reconstruction can be achieved for the second scheme if an additional far-field measurement is used. Rigorous mathematical justifications are provided and numerical results are presented to demonstrate the effectiveness and the promising features of the proposed imaging schemes.

preprint2013arXiv

Two Single-shot Methods for Locating Multiple Electromagnetic Scatterers

We develop two inverse scattering schemes for locating multiple electromagnetic (EM) scatterers by the electric far-field measurement corresponding to a single incident/detecting plane wave. The first scheme is for locating scatterers of small size compared to the wavelength of the detecting plane wave. The multiple scatterers could be extremely general with an unknown number of components, and each scatterer component could be either an impenetrable perfectly conducting obstacle or a penetrable inhomogeneous medium with an unknown content. The second scheme is for locating multiple perfectly conducting obstacles of regular size compared to the detecting EM wavelength. The number of the obstacle components is not required to be known in advance, but the shape of each component must be from a certain known admissible class. The admissible class may consist of multiple different reference obstacles. The second scheme could also be extended to include the medium components if a certain generic condition is satisfied. Both schemes are based on some novel indicator functions whose indicating behaviors could be used to locate the scatterers. No inversion will be involved in calculating the indicator functions, and the proposed methods are every efficient and robust to noise. Rigorous mathematical justifications are provided and extensive numerical experiments are conducted to illustrate the effectiveness of the imaging schemes.

preprint2012arXiv

A Direct Sampling Method for Inverse Scattering Using Far-Field Data

This work is concerned with a direct sampling method (DSM) for inverse acoustic scattering problems using far-field data. The method characterizes some unknown obstacles, inhomogeneous media or cracks, directly through an indicator function computed from the measured data. Using one or very few incident waves, the DSM provides quite reasonable profiles of scatterers in time-harmonic inverse acoustic scattering without a priori knowledge of either the physical properties or the number of disconnected components of the scatterer. We shall first derive the DSM using far-field data, then carry out a systematic evaluation of the performances and distinctions of the DSM using both near-field and far-field data. The numerical simulations are shown to demonstrate interesting and promising potentials of the DSM: a) ability to identify not only medium scatterers, but also obstacles, and even cracks, using measurement data from one or few incident directions, b) robustness with respect to large noise, and c) computational efficiency with only inner products involved.

preprint2011arXiv

Enhanced Approximate Cloaking by SH and FSH Lining

We consider approximate cloaking from a regularization viewpoint introduced in [13] for EIT and further investigated in [12] [17] for the Helmholtz equation. The cloaking schemes in [12] and [17] are shown to be (optimally) within $|\lnρ|^{-1}$ in 2D and $ρ$ in 3D of perfect cloaking, where $ρ$ denotes the regularization parameter. In this paper, we show that by employing a sound-hard layer right outside the cloaked region, one could (optimally) achieve $ρ^N$ in $\mathbb{R}^N,\ N\geq 2$, which significantly enhances the near-cloak. We then develop a cloaking scheme by making use of a lossy layer with well-chosen parameters. The lossy-layer cloaking scheme is shown to possess the same cloaking performance as the one with a sound-hard layer. Moreover, it is shown that the lossy layer could be taken as a finite realization of the sound-hard layer. Numerical experiments are also presented to assess the cloaking performances of all the cloaking schemes for comparisons.