{"title":"Time domain boundary integral equations and convolution quadrature for scattering by composite media","authors":"A. Rieder, F. Sayas, J. Melenk","doi":"10.1090/mcom/3730","DOIUrl":"https://doi.org/10.1090/mcom/3730","url":null,"abstract":"We consider acoustic scattering in heterogeneous media with piecewise constant wave number. The discretization is carried out using a Galerkin boundary element method in space and Runge-Kutta convolution quadrature in time. We prove well-posedness of the scheme and provide a priori estimates for the convergence in space and time.","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"1 1","pages":"2165-2195"},"PeriodicalIF":0.0,"publicationDate":"2020-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86494174","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A generalized Faulhaber inequality, improved bracketing covers, and applications to discrepancy","authors":"M. Gnewuch, Hendrik Pasing, Christian Weiss","doi":"10.1090/mcom/3666","DOIUrl":"https://doi.org/10.1090/mcom/3666","url":null,"abstract":"<p>We prove a generalized Faulhaber inequality to bound the sums of the <inline-formula content-type=\"math/mathml\">\u0000<mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" alttext=\"j\">\u0000 <mml:semantics>\u0000 <mml:mi>j</mml:mi>\u0000 <mml:annotation encoding=\"application/x-tex\">j</mml:annotation>\u0000 </mml:semantics>\u0000</mml:math>\u0000</inline-formula>-th powers of the first <inline-formula content-type=\"math/mathml\">\u0000<mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" alttext=\"n\">\u0000 <mml:semantics>\u0000 <mml:mi>n</mml:mi>\u0000 <mml:annotation encoding=\"application/x-tex\">n</mml:annotation>\u0000 </mml:semantics>\u0000</mml:math>\u0000</inline-formula> (possibly shifted) natural numbers. With the help of this inequality we are able to improve the known bounds for bracketing numbers of <inline-formula content-type=\"math/mathml\">\u0000<mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" alttext=\"d\">\u0000 <mml:semantics>\u0000 <mml:mi>d</mml:mi>\u0000 <mml:annotation encoding=\"application/x-tex\">d</mml:annotation>\u0000 </mml:semantics>\u0000</mml:math>\u0000</inline-formula>-dimensional axis-parallel boxes anchored in <inline-formula content-type=\"math/mathml\">\u0000<mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" alttext=\"0\">\u0000 <mml:semantics>\u0000 <mml:mn>0</mml:mn>\u0000 <mml:annotation encoding=\"application/x-tex\">0</mml:annotation>\u0000 </mml:semantics>\u0000</mml:math>\u0000</inline-formula> (or, put differently, of lower left orthants intersected with the <inline-formula content-type=\"math/mathml\">\u0000<mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" alttext=\"d\">\u0000 <mml:semantics>\u0000 <mml:mi>d</mml:mi>\u0000 <mml:annotation encoding=\"application/x-tex\">d</mml:annotation>\u0000 </mml:semantics>\u0000</mml:math>\u0000</inline-formula>-dimensional unit cube <inline-formula content-type=\"math/mathml\">\u0000<mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" alttext=\"left-bracket 0 comma 1 right-bracket Superscript d\">\u0000 <mml:semantics>\u0000 <mml:mrow>\u0000 <mml:mo stretchy=\"false\">[</mml:mo>\u0000 <mml:mn>0</mml:mn>\u0000 <mml:mo>,</mml:mo>\u0000 <mml:mn>1</mml:mn>\u0000 <mml:msup>\u0000 <mml:mo stretchy=\"false\">]</mml:mo>\u0000 <mml:mi>d</mml:mi>\u0000 </mml:msup>\u0000 </mml:mrow>\u0000 <mml:annotation encoding=\"application/x-tex\">[0,1]^d</mml:annotation>\u0000 </mml:semantics>\u0000</mml:math>\u0000</inline-formula>). We use these improved bracketing numbers to establish new bounds for the star-discrepancy of negatively dependent random point sets and its expectation. We apply our findings also to the weighted star-discrepancy.</p>","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"437 1","pages":"2873-2898"},"PeriodicalIF":0.0,"publicationDate":"2020-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78783459","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Approximation properties of sum-up rounding in the presence of vanishing constraints","authors":"Paul Manns, C. Kirches, F. Lenders","doi":"10.1090/mcom/3606","DOIUrl":"https://doi.org/10.1090/mcom/3606","url":null,"abstract":"Approximation algorithms like sum-up rounding that allow to compute integer-valued approximations of the continuous controls in a weak∗ sense have attracted interest recently. They allow to approximate (optimal) feasible solutions of continuous relaxations of mixed-integer control problems (MIOCPs) with integer controls arbitrarily close. To this end, they use compactness properties of the underlying state equation, a feature that is tied to the infinite-dimensional vantage point. In this work, we consider a class of MIOCPs that are constrained by pointwise mixed state-control constraints. We show that a continuous relaxation that involves so-called vanishing constraints has beneficial properties for the described approximation methodology. Moreover, we complete recent work on a variant of the sum-up rounding algorithm for this problem class. In particular, we prove that the observed infeasibility of the produced integer-valued controls vanishes in an L∞-sense with respect to the considered relaxation. Moreover, we improve the bound on the control approximation error to a value that is asymptotically tight.","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"4058 2 1","pages":"1263-1296"},"PeriodicalIF":0.0,"publicationDate":"2020-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86743200","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Imaginary multiquadratic number fields with class group of exponent 3 and 5","authors":"Jürgen Klüners, T. Komatsu","doi":"10.1090/mcom/3609","DOIUrl":"https://doi.org/10.1090/mcom/3609","url":null,"abstract":"","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"53 1","pages":"1483-1497"},"PeriodicalIF":0.0,"publicationDate":"2020-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83103469","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Optimal error estimates of discontinuous Galerkin methods with generalized fluxes for wave equations on unstructured meshes","authors":"Zheng Sun, Y. Xing","doi":"10.1090/mcom/3605","DOIUrl":"https://doi.org/10.1090/mcom/3605","url":null,"abstract":"L2 stable discontinuous Galerkin method with a family of numerical fluxes was studied for the one-dimensional wave equation by Cheng, Chou, Li, and Xing in [Math. Comp. 86 (2017), pp. 121–155]. Although optimal convergence rates were numerically observed with wide choices of parameters in the numerical fluxes, their error estimates were only proved for a sub-family with the construction of a local projection. In this paper, we first complete the one-dimensional analysis by providing optimal error estimates that match all numerical observations in that paper. The key ingredient is to construct an optimal global projection with the characteristic decomposition. We then extend the analysis on optimal error estimate to multidimensions by constructing a global projection on unstructured meshes, which can be considered as a perturbation away from the local projection studied by Cockburn, Gopalakrishnan, and Sayas in [Math. Comp. 79 (2010), pp. 1351–1367] for hybridizable discontinuous Galerkin methods. As a main contribution, we use a novel energy argument to prove the optimal approximation property of the global projection. This technique does not require explicit assembly of the matrix for the perturbed terms and hence can be easily used for unstructured meshes in multidimensions. Finally, numerical tests in two dimensions are provided to validate our analysis is sharp and at least one of the unknowns will degenerate to suboptimal rates if the assumptions are not satisfied.","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"28 1","pages":"1741-1772"},"PeriodicalIF":0.0,"publicationDate":"2020-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80436843","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Sampling and homology via bottlenecks","authors":"S. Rocco, David Eklund, Oliver Gäfvert","doi":"10.1090/mcom/3757","DOIUrl":"https://doi.org/10.1090/mcom/3757","url":null,"abstract":"In this paper we present an efficient algorithm to produce a provably dense sample of a smooth compact variety. The procedure is partly based on computing $textit{bottlenecks}$ of the variety. Using geometric information such as the bottlenecks and the $textit{local reach}$ we also provide bounds on the density of the sample needed in order to guarantee that the homology of the variety can be recovered from the sample. An implementation of the algorithm is provided together with numerical experiments and a computational comparison to the algorithm by Dufresne et. al.","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"89 ","pages":"2969-2995"},"PeriodicalIF":0.0,"publicationDate":"2020-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91444765","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sualeh Asif, Francesc Fit'e, Dylan Pentland, Andrew V. Sutherland
{"title":"Computing L-polynomials of Picard curves from Cartier-Manin matrices","authors":"Sualeh Asif, Francesc Fit'e, Dylan Pentland, Andrew V. Sutherland","doi":"10.1090/mcom/3675","DOIUrl":"https://doi.org/10.1090/mcom/3675","url":null,"abstract":"We study the sequence of zeta functions $Z(C_p,T)$ of a generic Picard curve $C:y^3=f(x)$ defined over $mathbb{Q}$ at primes $p$ of good reduction for $C$. By determining the density of the set of primes of ordinary reduction, we prove that, for all but a density zero subset of primes, the Zeta function $Z(C_p,T)$ is uniquely determined by the Cartier--Manin matrix $A_p$ of $C$ modulo $p$, the irreducibility of $f$ modulo $p$ (or the failure thereof), and the exponent of the Jacobian of $C$ modulo $p$; we also show that for primes $equiv 1 pmod{3}$ the matrix $A_p$ suffices and that for primes $equiv 2 pmod{3}$ the genericity assumption on $C$ is unnecessary. By combining this with recent work of Sutherland, we obtain a practical probabilistic algorithm of Las Vegas type that computes $Z(C_p,T)$ for almost all primes $p le N$ using $Nlog(N)^{3+o(1)}$ expected bit operations. This is the first practical result of this type for curves of genus greater than 2.","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"10 1","pages":"943-971"},"PeriodicalIF":0.0,"publicationDate":"2020-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81918599","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An exponent one-fifth algorithm for deterministic integer factorisation","authors":"David Harvey","doi":"10.1090/MCOM/3658","DOIUrl":"https://doi.org/10.1090/MCOM/3658","url":null,"abstract":"Hittmeir recently presented a deterministic algorithm that provably computes the prime factorisation of a positive integer $N$ in $N^{2/9+o(1)}$ bit operations. Prior to this breakthrough, the best known complexity bound for this problem was $N^{1/4+o(1)}$, a result going back to the 1970s. In this paper we push Hittmeir's techniques further, obtaining a rigorous, deterministic factoring algorithm with complexity $N^{1/5+o(1)}$.","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":"28 1","pages":"2937-2950"},"PeriodicalIF":0.0,"publicationDate":"2020-10-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77977199","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}