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

Makoto Matsumoto

Makoto Matsumoto contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.COmath.NAmath.CTMathematical Softwaremath.AGmath.RANumerical Analysis

Trust snapshot

Quick read

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

11 published item(s)

preprint2026arXiv

Some Patterns of Duplications in the outputs of Mersenne Twister Pseudorandom Number Generator MT19937

The Mersenne Twister MT19937 pseudorandom number generator, introduced by the last two authors in 1998, is still widely used. It passes all existing statistical tests, except for the linear complexity test, which measures the ratio of the even-odd of the number of 1's among specific bits (and hence should not be important for most applications). Harase reported that MT19937 is rejected by some birthday-spacing tests, which are rather artificially designed. In this paper, we report that MT19937 fails in a natural test based on the distribution of run-lengths on which we found an identical value in the output 32-bit integers. The number of observations of the run-length 623 is some 40 times larger than the expectation (and than the numbers of the observations of 622 and 624, etc.), which implies that the corresponding p-value is almost 0. We mathematically analyze the phenomena, and obtain a theorem which explains these failures. It seems not to be a serious defect of MT19937, because finding the defect requires astronomical efforts. Still, the phenomena should be reported to the academic society relating to pseudorandom number generation.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

A categorical proof of the nonexistence of (120, 35, 10)-difference sets

A difference set with parameters $(v, k, λ)$ is a subset $D$ of cardinality $k$ in a finite group $G$ of order $v$, such that the number $λ$ of occurrences of $g \in G$ as the ratio $d^{-1}d'$ in distinct pairs $(d, d')\in D\times D$ is independent of $g$. We prove the nonexistence of $(120, 35, 10)$-difference sets, which has been an open problem for 70 years since Bruck introduced the notion of nonabelian difference sets. Our main tools are 1. a generalization of the category of finite groups to that of association schemes (actually, to that of relation partitions), 2. a generalization of difference sets to equi-distributed functions and its preservation by pushouts along quotients, 3. reduction to a linear programming in the nonnegative integer lattice with quadratic constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Functoriality of Bose-Mesner algebras and profinite association schemes

We show that taking the set of primitive idempotents of commutative association schemes is a functor from the category of commutative association schemes with surjective morphisms to the category of finite sets with surjective partial functions. We then consider projective systems of commutative association schemes consisting of surjections (which we call profinite association schemes), for which Bose-Mesner algebra is defined, and describe a Delsarte theory on such schemes. This is another method for generalizing association schemes to those on infinite sets, related with the approach by Barg and Skriganov. Relation with $(t,m,s)$-nets and $(t,s)$-sequences is studied. We reprove some of the results of Martin-Stinson from this viewpoint.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Locally-finite extensive categories, their semi-rings, and decomposition to connected objects

Let $\mathcal C$ be the category of finite graphs. Lovàsz shows that the semi-ring of isomorphism classes of $\mathcal C$ (with coproduct as sum, and product as multiplication) is embedded into the direct product of the semi-ring of natural numbers. Our aim is to generalize this result to other categories. For this, one crucial property is that every object decomposes to a finite coproduct of connected objects. We show that a locally-finite extensive category satisfies this condition. Conversely, a category where any object is decomposed into a finite coproduct of connected objects is shown to be extensive. The decomposition turns out to be unique. Using these results, we give some sufficient conditions that the semi-ring (the ring) of isomorphism classes of a locally finite category embeds to the direct product of natural numbers (integers, respectively). Such a construction of rings from a category is a most primitive form of Burnside rings and Grothendieck rings.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lovàsz's hom-counting theorem by inclusion-exclusion principle

Let ${\mathcal C}$ be the category of finite graphs. Lovàsz (1967) shows that if $|\mathrm{Hom}(X,A)|=|\mathrm{Hom}(X,B)|$ holds for any $X$, then $A$ is isomorphic to $B$. Pultr (1973) gives a categorical generalization using a similar argument. Both proofs assume that each object has a finite number of isomorphism classes of subobjects. Generalizations without this assumption are given by Dawar, Jakl, and Reggio (2021) and Regio (2021). Here another generalization without this assumption is given, with a shorter proof. Examples of categories are given, for which our theorem is applicable, but the existing theorems are not.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Approximation of integration over finite groups, difference sets and association schemes

Let $G$ be a finite group and $f:G \to {\mathbb C}$ be a function. For a non-empty finite subset $Y\subset G$, let $I_Y(f)$ denote the average of $f$ over $Y$. Then, $I_G(f)$ is the average of $f$ over $G$. Using the decomposition of $f$ into irreducible components of ${\mathbb C}^G$ as a representation of $G\times G$, we define non-negative real numbers $V(f)$ and $D(Y)$, each depending only on $f$, $Y$, respectively, such that an inequality of the form $|I_G(f)-I_Y(f)|\leq V(f)\cdot D(Y)$ holds. We give a lower bound of $D(Y)$ depending only on $\#Y$ and $\#G$. We show that the lower bound is achieved if and only if $\#\{(x,y)\in Y^2 \mid x^{-1}y \in [a]\}/\#[a]$ is independent of the choice of the conjugacy class $[a]\subset G$ for $a \neq 1$. We call such a $Y\subset G$ as a pre-difference set in $G$, since the condition is satisfied if $Y$ is a difference set. If $G$ is abelian, the condition is equivalent to that $Y$ is a difference set. We found a non-trivial pre-difference set in the dihedral group of order 16, where no non-trivial difference set exists. The pre-difference sets in non-abelian groups of order 16 are classified. A generalization to commutative association schemes is also given.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Universal Mixed Elliptic Motives

In this paper we construct a Q-linear tannakian category MEM_1 of universal mixed elliptic motives over the moduli space M_{1,1} of elliptic curves. It contains MTM, the category of mixed Tate motives unramified over the integers. Each object of MEM_1 is an object of MTM endowed with an action of SL_2(Z) that is compatible with its structure. Universal mixed elliptic motives can be thought of as motivic local systems over M_{1,1} whose fiber over the tangential base point d/dq at the cusp is a mixed Tate motive. The basic structure of the tannakian fundamental group of MEM is determined and the lowest order terms of all relations are found (using computations of Francis Brown), including the arithmetic relations, which describe the "infinitesimal Galois action". We use the presentation to give a new and more conceptual proof of the Ihara-Takao congruences.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Walsh Figure of Merit for Digital Nets: An Easy Measure for Higher Order Convergent QMC

Fix an integer $s$. Let $f:[0,1)^s \to \mathbb R$ be an integrable function. Let $P\subset [0,1]^s$ be a finite point set. Quasi-Monte Carlo integration of $f$ by $P$ is the average value of $f$ over $P$ that approximates the integration of $f$ over the $s$-dimensional cube. Koksma-Hlawka inequality tells that, by a smart choice of $P$, one may expect that the error decreases roughly $O(N^{-1}(\log N)^s)$. For any $α\geq 1$, J.\ Dick gave a construction of point sets such that for $α$-smooth $f$, convergence rate $O(N^{-α}(\log N)^{sα})$ is assured. As a coarse version of his theory, M-Saito-Matoba introduced Walsh figure of Merit (WAFOM), which gives the convergence rate $O(N^{-C\log N/s})$. WAFOM is efficiently computable. By a brute-force search of low WAFOM point sets, we observe a convergence rate of order $N^{-α}$ with $α>1$, for several test integrands for $s=4$ and $8$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Computable Figure of Merit for Quasi-Monte Carlo Point Sets

Let $\mathcal{P} \subset [0,1)^S$ be a finite point set of cardinality $N$ in an $S$-dimensional cube, and let $f:[0,1)^S \to \mathbb{R}$ be an integrable function. A QMC integration of $f$ by $\mathcal{P}$ is the average of values of $f$ at each point in $\mathcal{P}$, which approximates the integration of $f$ over the cube. Assume that $\mathcal{P}$ is constructed from an $\mathbb{F}2$-vector space $P\subset (\F2^n)^S$ by means of a digital net with $n$-digit precision. As an $n$-digit discretized version of Josef Dick's method, we introduce Walsh figure of merit (WAFOM) $\textnormal{WF}(P)$ of $P$, which satisfies a Koksma-Hlawka type inequality, namely, QMC integration error is bounded by $C_{S,n}||f||_n \textnormal{WF}(P)$ under $n$-smoothness of $f$, where $C_{S,n}$ is a constant depending only on $S,n$. We show a Fourier inversion formula for $\textnormal{WF}(P)$ which is computable in $O(n SN)$ steps. This effectiveness enables us a random search for $P$ with small value of $\textnormal{WF}(P)$, which would be difficult for other figures of merit such as discrepancy. From an analogy to coding theory, we expect that random search may find better point sets than mathematical constructions. In fact, a naïve search finds point sets $P$ with small $\textnormal{WF}(P)$. In experiments, we show better performance of these point sets in QMC integration than widely used QMC rules. We show some experimental evidence on the effectiveness of our point sets to even non-smooth integrands appearing in finance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the fast computation of the weight enumerator polynomial and the $t$ value of digital nets over finite abelian groups

In this paper we introduce digital nets over finite abelian groups which contain digital nets over finite fields and certain rings as a special case. We prove a MacWilliams type identity for such digital nets. This identity can be used to compute the strict $t$-value of a digital net over finite abelian groups. If the digital net has $N$ points in the $s$ dimensional unit cube $[0,1]^s$, then the $t$-value can be computed in $\mathcal{O}(N s \log N)$ operations and the weight enumerator polynomial can be computed in $\mathcal{O}(N s (\log N)^2)$ operations, where operations mean arithmetic of integers. By precomputing some values the number of operations of computing the weight enumerator polynomial can be reduced further.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Variants of Mersenne Twister Suitable for Graphic Processors

This paper proposes a type of pseudorandom number generator, Mersenne Twister for Graphic Processor (MTGP), for efficient generation on graphic processessing units (GPUs). MTGP supports large state sizes such as 11213 bits, and uses the high parallelism of GPUs in computing many steps of the recursion in parallel. The second proposal is a parameter-set generator for MTGP, named MTGP Dynamic Creator (MTGPDC). MT- GPDC creates up to 2^32 distinct parameter sets which generate sequences with high-dimensional uniformity. This facility is suitable for a large grid of GPUs where each GPU requires separate random number streams. MTGP is based on linear recursion over the two-element field, and has better high-dimensional equidistribution than the Mersenne Twister pseudorandom number generator.

0 citations0 reviewsNo code linkedOpen access