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

Yuto Miyatake

Yuto Miyatake contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical AnalysisDistributed, Parallel, and Cluster ComputingMachine Learningphysics.data-an

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
5topics
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

7 published item(s)

preprint2022arXiv

Composing a surrogate observation operator for sequential data assimilation

In data assimilation, state estimation is not straightforward when the observation operator is unknown. This study proposes a method for composing a surrogate operator when the true operator is unknown. A neural network is used to improve the surrogate model iteratively to decrease the difference between the observations and the results of the surrogate model. A twin experiment suggests that the proposed method outperforms approaches that tentatively use a specific operator throughout the data assimilation process.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High-order linearly implicit schemes conserving quadratic invariants

In this paper, we propose linearly implicit and arbitrary high-order conservative numerical schemes for ordinary differential equations with a quadratic invariant. Many differential equations have invariants, and numerical schemes for preserving them have been extensively studied. Since linear invariants can be easily kept after discretisation, quadratic invariants are essentially the simplest ones. Quadratic invariants are important objects that appear not only in many physical examples but also in the computationally efficient conservative schemes for general invariants such as scalar auxiliary variable approach, which have been studied in recent years. It is known that quadratic invariants can be kept relatively easily compared to general invariants, and indeed can be preserved by canonical Runge--Kutta methods. However, there is no unified method for constructing linearly implicit and high order conservative schemes. In this paper, we construct such schemes based on canonical Runge--Kutta methods and prove some properties involving accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Structure-preserving numerical methods for constrained gradient flows of planar closed curves with explicit tangential velocities

In this paper, we consider numerical approximation of constrained gradient flows of planar closed curves, including the Willmore and the Helfrich flows. These equations have energy dissipation and the latter has conservation properties due to the constraints. We will develop structure-preserving methods for these equations that preserve both the dissipation and the constraints. To preserve the energy structures, we introduce the discrete version of gradients according to the discrete gradient method and determine the Lagrange multipliers appropriately. We directly address higher order derivatives by using the Galerkin method with B-spline curves to discretize curves. Moreover, we will consider stabilization of the schemes by adding tangential velocities. We introduce a new Lagrange multiplier to obtain both the energy structures and the stability. Several numerical examples are presented to verify that the proposed schemes preserve the energy structures with good distribution of control points.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A fully discrete curve-shortening polygonal evolution law for moving boundary problems

We consider the numerical integration of moving boundary problems with the curve-shortening property, such as the mean curvature flow and Hele-Shaw flow. We propose a fully discrete curve-shortening polygonal evolution law. The proposed evolution law is fully implicit, and the key to the derivation is to devise the definitions of tangent and normal vectors and tangential and normal velocities at each vertex in an implicit manner. Numerical experiments show that the proposed method allows the use of relatively large time step sizes and also captures the area-preserving or dissipative property in good accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Parallelizable Energy-Preserving Integrator MB4 and Its Application to Quantum-Mechanical Wavepacket Dynamics

In simulating physical systems, conservation of the total energy is often essential, especially when energy conversion between different forms of energy occurs frequently. Recently, a new fourth order energy-preserving integrator named MB4 was proposed based on the so-called continuous stage Runge--Kutta methods (Y.~Miyatake and J.~C.~Butcher, SIAM J.~Numer.~Anal., 54(3), 1993-2013). A salient feature of this method is that it is parallelizable, which makes its computational time for one time step comparable to that of second order methods. In this paper, we illustrate how to apply the MB4 method to a concrete ordinary differential equation using the nonlinear Schrödinger-type equation on a two-dimensional grid as an example. This system is a prototypical model of two-dimensional disordered organic material and is difficult to solve with standard methods like the classical Runge--Kutta methods due to the nonlinearity and the $δ$-function like potential coming from defects. Numerical tests show that the method can solve the equation stably and preserves the total energy to 16-digit accuracy throughout the simulation. It is also shown that parallelization of the method yields up to 2.8 times speedup using 3 computational nodes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adjoint-based exact Hessian computation

We consider a scalar function depending on a numerical solution of an initial value problem, and its second-derivative (Hessian) matrix for the initial value. The need to extract the information of the Hessian or to solve a linear system having the Hessian as a coefficient matrix arises in many research fields such as optimization, Bayesian estimation, and uncertainty quantification. From the perspective of memory efficiency, these tasks often employ a Krylov subspace method that does not need to hold the Hessian matrix explicitly and only requires computing the multiplication of the Hessian and a given vector. One of the ways to obtain an approximation of such Hessian-vector multiplication is to integrate the so-called second-order adjoint system numerically. However, the error in the approximation could be significant even if the numerical integration to the second-order adjoint system is sufficiently accurate. This paper presents a novel algorithm that computes the intended Hessian-vector multiplication exactly and efficiently. For this aim, we give a new concise derivation of the second-order adjoint system and show that the intended multiplication can be computed exactly by applying a particular numerical method to the second-order adjoint system. In the discussion, symplectic partitioned Runge--Kutta methods play an essential role.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generalization of partitioned Runge--Kutta methods for adjoint systems

This study computes the gradient of a function of numerical solutions of ordinary differential equations (ODEs) with respect to the initial condition. The adjoint method computes the gradient approximately by solving the corresponding adjoint system numerically. In this context, Sanz-Serna [SIAM Rev., 58 (2016), pp. 3--33] showed that when the initial value problem is solved by a Runge--Kutta (RK) method, the gradient can be exactly computed by applying an appropriate RK method to the adjoint system. Focusing on the case where the initial value problem is solved by a partitioned RK (PRK) method, this paper presents a numerical method, which can be seen as a generalization of PRK methods, for the adjoint system that gives the exact gradient.

0 citations0 reviewsNo code linkedOpen access