BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jongho Park

Jongho Park contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEmath.NANumerical Analysismath.OCastro-ph.GAMachine Learningastro-ph.IMcond-mat.mtrl-sciData Structures and Algorithmsphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2026arXiv

A high-order augmented Lagrangian method with arbitrarily fast convergence

We propose a high-order version of the augmented Lagrangian method for solving convex optimization problems with linear constraints, which achieves arbitrarily fast -- and even superlinear -- convergence rates. First, we analyze the convergence rates of the high-order proximal point method under certain uniform convexity assumptions on the energy functional. We then introduce the high-order augmented Lagrangian method and analyze its convergence by leveraging the convergence results of the high-order proximal point method. Finally, we present applications of the high-order augmented Lagrangian method to various problems arising in the sciences, including data fitting, flow in porous media, and scientific machine learning.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast Gradient Methods for Uniformly Convex and Weakly Smooth Problems

In this paper, acceleration of gradient methods for convex optimization problems with weak levels of convexity and smoothness is considered. Starting from the universal fast gradient method which was designed to be an optimal method for weakly smooth problems whose gradients are Hölder continuous, its momentum is modified appropriately so that it can also accommodate uniformly convex and weakly smooth problems. Different from the existing works, fast gradient methods proposed in this paper do not use the restarting technique but use momentums that are suitably designed to reflect both the uniform convexity and the weak smoothness information of the target energy function. Both theoretical and numerical results that support the superiority of proposed methods are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

ReLU Regression with Massart Noise

We study the fundamental problem of ReLU regression, where the goal is to fit Rectified Linear Units (ReLUs) to data. This supervised learning task is efficiently solvable in the realizable setting, but is known to be computationally hard with adversarial label noise. In this work, we focus on ReLU regression in the Massart noise model, a natural and well-studied semi-random noise model. In this model, the label of every point is generated according to a function in the class, but an adversary is allowed to change this value arbitrarily with some probability, which is {\em at most} $η< 1/2$. We develop an efficient algorithm that achieves exact parameter recovery in this model under mild anti-concentration assumptions on the underlying distribution. Such assumptions are necessary for exact recovery to be information-theoretically possible. We demonstrate that our algorithm significantly outperforms naive applications of $\ell_1$ and $\ell_2$ regression on both synthetic and real data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Photon Ring in M87*

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called &#34;photon rings&#34; -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Unravelling the Innermost Jet Structure of OJ 287 with the First GMVA+ALMA Observations

We present the first very-long-baseline interferometric (VLBI) observations of the blazar OJ287 carried out jointly with the Global Millimeter VLBI Array (GMVA) and the phased Atacama Large Millimeter/submillimeter Array (ALMA) at 3.5 mm on April 2, 2017. Participation of phased-ALMA not only has improved the GMVA north-south resolution by a factor of ~3, but also has enabled fringe detection with signal-to-noise ratios up to 300 at baselines longer than 2 Gλ. The high sensitivity has motivated us to image the data with the newly developed regularized maximum likelihood imaging methods, revealing the innermost jet structure with unprecedentedly high angular resolution. Our images reveal a compact and twisted jet extending along the northwest direction with two bends within the inner 200 μas that resembles a precessing jet in projection. The component at the southeastern end shows a compact morphology and high brightness temperature, and is identified as the VLBI core. An extended jet feature that lies at ~200 μas northwest of the core shows a conical shape in both total and linearly polarized intensity, and a bimodal distribution of the linear polarization electric vector position angle. We discuss the nature of this feature by comparing our observations with models and simulations of oblique and recollimation shocks with various magnetic field configurations. Our high-fidelity images also enabled us to search for possible jet features from the secondary supermassive black hole (SMBH) and test the SMBH binary hypothesis proposed for this source.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Accelerated Additive Schwarz Methods for Convex Optimization with Adaptive Restart

Based on an observation that additive Schwarz methods for general convex optimization can be interpreted as gradient methods, we propose an acceleration scheme for additive Schwarz methods. Adopting acceleration techniques developed for gradient methods such as momentum and adaptive restarting, the convergence rate of additive Schwarz methods is greatly improved. The proposed acceleration scheme does not require any a priori information on the levels of smoothness and sharpness of a target energy functional, so that it can be applied to various convex optimization problems. Numerical results for linear elliptic problems, nonlinear elliptic problems, nonsmooth problems, and nonsharp problems are provided to highlight the superiority and the broad applicability of the proposed scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Additive Schwarz Methods for Convex Optimization with Backtracking

This paper presents a novel backtracking strategy for additive Schwarz methods for general convex optimization problems as an acceleration scheme. The proposed backtracking strategy is independent of local solvers, so that it can be applied to any algorithms that can be represented in an abstract framework of additive Schwarz methods. Allowing for adaptive increasing and decreasing of the step size along the iterations, the convergence rate of an algorithm is greatly improved. Improved convergence rate of the algorithm is proven rigorously. In addition, combining the proposed backtracking strategy with a momentum acceleration technique, we propose a further accelerated additive Schwarz method. Numerical results for various convex optimization problems that support our theory are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The intrinsic structure of Sagittarius A* at 1.3 cm and 7 mm

Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets to resolve the innermost region of SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multi-wavelength campaign of the Event Horizon Telescope (EHT) in 2017 April. To mitigate scattering effects, the physically motivated scattering kernel model from Psaltis et al. (2018) and the scattering parameters from Johnson et al. (2018) have been applied. As a result, a single, symmetric Gaussian model well describes the intrinsic structure of Sgr A* at both wavelengths. From closure amplitudes, the major-axis sizes are ~704$\pm$102 $μ$as (axial ratio $\sim$1.19$^{+0.24}_{-0.19}$) and $\sim$300$\pm$25 $μ$as (axial ratio $\sim$1.28$\pm$0.2) at 1.349 cm and 6.95 mm respectively. Together with a quasi-simultaneous observation at 3.5 mm (86 GHz) by Issaoun et al. (2019), we show that the intrinsic size scales with observing wavelength as a power-law, with an index $\sim$1.2$\pm$0.2. Our results also provide estimates of the size and compact flux density at 1.3 mm, which can be incorporated into the analysis of the EHT observations. In terms of the origin of radio emission, we have compared the intrinsic structures with the accretion flow scenario, especially the radiatively inefficient accretion flow based on the Keplerian shell model. With this, we show that a nonthermal electron population is necessary to reproduce the source sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Two-level Group Convolution

Group convolution has been widely used in order to reduce the computation time of convolution, which takes most of the training time of convolutional neural networks. However, it is well known that a large number of groups significantly reduce the performance of group convolution. In this paper, we propose a new convolution methodology called ``two-level&#39;&#39; group convolution that is robust with respect to the increase of the number of groups and suitable for multi-GPU parallel computation. We first observe that the group convolution can be interpreted as a one-level block Jacobi approximation of the standard convolution, which is a popular notion in the field of numerical analysis. In numerical analysis, there have been numerous studies on the two-level method that introduces an intergroup structure that resolves the performance degradation issue without disturbing parallel computation. Motivated by these, we introduce a coarse-level structure which promotes intergroup communication without being a bottleneck in the group convolution. We show that all the additional work induced by the coarse-level structure can be efficiently processed in a distributed memory system. Numerical results that verify the robustness of the proposed method with respect to the number of groups are presented. Moreover, we compare the proposed method to various approaches for group convolution in order to highlight the superiority of the proposed method in terms of execution time, memory efficiency, and performance.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Additive Schwarz Methods for Convex Optimization as Gradient Methods

This paper gives a unified convergence analysis of additive Schwarz methods for general convex optimization problems. Resembling to the fact that additive Schwarz methods for linear problems are preconditioned Richardson methods, we prove that additive Schwarz methods for general convex optimization are in fact gradient methods. Then an abstract framework for convergence analysis of additive Schwarz methods is proposed. The proposed framework applied to linear elliptic problems agrees with the classical theory. We present applications of the proposed framework to various interesting convex optimization problems such as nonlinear elliptic problems, nonsmooth problems, and nonsharp problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Overlapping Domain Decomposition Framework without Dual Formulation for Variational Imaging Problems

In this paper, we propose a novel overlapping domain decomposition method that can be applied to various problems in variational imaging such as total variation minimization. Most of recent domain decomposition methods for total variation minimization adopt the Fenchel--Rockafellar duality, whereas the proposed method is based on the primal formulation. Thus, the proposed method can be applied not only to total variation minimization but also to those with complex dual problems such as higher order models. In the proposed method, an equivalent formulation of the model problem with parallel structure is constructed using a custom overlapping domain decomposition scheme with the notion of essential domains. As a solver for the constructed formulation, we propose a decoupled augmented Lagrangian method for untying the coupling of adjacent subdomains. Convergence analysis of the decoupled augmented Lagrangian method is provided. We present implementation details and numerical examples for various model problems including total variation minimizations and higher order models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GPCAL: a generalized calibration pipeline for instrumental polarization in VLBI data

We present the Generalized Polarization CALibration pipeline (GPCAL), an automated pipeline for instrumental polarization calibration of very long baseline interferometry (VLBI) data. The pipeline is designed to achieve a high calibration accuracy by means of fitting the instrumental polarization model, including the second-order terms, to multiple calibrators data simultaneously. It also allows using more accurate linear polarization models of calibrators for D-term estimation compared to the conventional way that assumes similar linear polarization and total intensity structures. This assumption has widely been used in the existing packages for instrumental polarization calibration but could be a source of significant uncertainties when there is no suitable calibrator satisfying the assumption. We demonstrate the capabilities of GPCAL by using simulated data, archival Very Long Baseline Array (VLBA) data of many active galactic nuclei (AGN) jets at 15 and 43 GHz, and our Korean VLBI Network (KVN) observations of many AGN jets at 86, 95, 130, and 142 GHz. The pipeline could reproduce the complex linear polarization structures of several sources shown in the previous studies using the same VLBA data. GPCAL also reveals a complex linear polarization structure in the flat-spectrum radio quasar 3C 273 from the KVN data at all four frequencies. These results demonstrate that GPCAL can achieve a high calibration accuracy for various VLBI arrays.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Off-centered Pb interstitials in PbTe

In this work, we calculate the defect properties of low-symmetry Pb interstitials in PbTe using first-principles density-functional theory calculations. We break the symmetry imposed on on-centered interstitial defects and show that the lowest ground state of Pb interstitial defects is off-centered along the [111] directions. Due to the four multi-stable structures with low defect formation energies, the defect density of Pb interstitials is expected to be ~5.6 times larger than previous predictions when PbTe is synthesized at 900 K. In contrast to the on-centered Pbinterstitial, the off-centered Pb interstitials in PbTe can exhibit long-range lattice relaxation toward [111] direction beyond distance of 1 nm, indicating the potential formation of weak local dipoles. This result provides an alternative explanation for the emphanitic anharmonicity of PbTe.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pseudo-linear Convergence of an Additive Schwarz Method for Dual Total Variation Minimization

In this paper, we propose an overlapping additive Schwarz method for total variation minimization based on a dual formulation. The $O(1/n)$-energy convergence of the proposed method is proven, where $n$ is the number of iterations. In addition, we introduce an interesting convergence property called pseudo-linear convergence of the proposed method; the energy of the proposed method decreases as fast as linearly convergent algorithms until it reaches a particular value. It is shown that such the particular value depends on the overlapping width $δ$, and the proposed method becomes as efficient as linearly convergent algorithms if $δ$ is large. As the latest domain decomposition methods for total variation minimization are sublinearly convergent, the proposed method outperforms them in the sense of the energy decay. Numerical experiments which support our theoretical results are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Exploring the Morphology and Origins of the 4C 38.41 Jet

We study the properties of the innermost jet of the flat spectrum radio quasar 1633+382 (4C~38.41) based on VLBI data from the radio monitoring observations of the Boston University VLBI program at 43~GHz. Analysis of the components suggests a semi-parabolic jet geometry with jet radius $R$ following the relation $R\propto r^{0.7}$ with distance $r$, with indications of a jet geometry break towards a conical geometry. Brightness temperature falls with distance following $T_B\propto r^{-2.1}$. Combining this information, magnetic field and electron densities are found to fall along the jet as $B\propto r^{-1.5}$ and $n\propto r^{-1.1}$ respectively, suggesting that the magnetic configuration in the jet may be dominated by the poloidal component. Our analysis of the jet structure suggests that the innermost jet regions do not follow a ballistic trajectory and, instead, match a sinusoidal morphology which could be due to jet precession from a helical pattern or Kelvin-Helmholtz instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Kinematics of the M87 jet in the collimation zone: gradual acceleration and velocity stratification

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($Γ\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

0 citations0 reviewsNo code linkedOpen access