Special pointsets


  • For special pointsets over intervals, triangles (simplex), disk, square, and polygons, see the following url.

Recent codes (last update: June 7, 2022)





» Object: CQMC: a software package for low-dimensional Compressed Quasi-MonteCarlo cubature.

The main folder contains the following compressed folders:

  • ROUTINES_2D.zip, that include many routines for computation of integrals via QMC/CQMC methods on domains obtained by set operation from two bivariate domains whose boundary is defined by piecewise NURBS;
  • ROUTINES_3D.zip, that include many routines for computation of integrals via QMC/CQMC methods on domains obtained by set operation from a trivariate polyhedron and a sphere;
In both cases, the main routines are:

  • make_rule, that defines a QMC rule over the domain and then extracts the CQMC rule via cqmc_01;
  • cqmc_01, routine that from a sufficiently dense QMC rule computes a low cardinality CQMC rule that matches the QMC moments up to a certain total polynomial degree.

» Sources:
» First update: June 07, 2022. Last update: June 07, 2022.





» Object: IN_RS: Matlab routine for points indomain on NURBS-shaped planar domains.

The main folder IN_RS contains:

  • inRS.m, a new routine for indomain over NURBS-shaped domains;
  • inRS_21.m, an old routine for indomain over NURBS-shaped domains;
  • inRSpoly2.m, a routine for indomain over NURBS-shaped domains (based on approximating the domain with a polygon, it requires inpoly2 by D. Engwirda);
  • demo_nurbs_indomain_01.m, demo_nurbs_indomain_02.m, demo_nurbs_indomain_03.m are demos, useful to define the domains and to understand how to use the indomain routine.
  • demo_nurbs_indomain_04.m, evaluates the performance of our two indomain routines for NURBS presented in [1] and [2], comparing them to one in which the NURBS boundary is approximated by a polygon, finally using inpoly2 by D. Engwirda. To make that comparison, please download the routine inpoly2 and save it in the IN_RS folder.
  • the routine demo_nurbs_indomain_05.m shows how to define the problem when only some vertices of the curved polygon are available;
  • the routine demo_nurbs_indomain_06.m considers the case of a domain whose boundary is defined only by a point cloud of possibly inexact data; this routine requires spcsp.m (cubic smoothing periodic spline), by M. Zanetti, available at the Matlab Central File Exchange;
  • the routine demo_nurbs_indomain_07.m shows how to define the problem when the whole domain is defined via point cloud and NURBS use the vertices, provided by alphashape polygonal approximation of the boundary, as control points;
  • the routine demo_nurbs_indomain_08.m shows how to define the problem when the whole domain is defined via point cloud and parametric cubic smoothing periodic spline are used on the data provided by alphashape polygonal approximation of the boundary; this routine requires spcsp.m, by M. Zanetti, available at the Matlab Central File Exchange;

» Sources:

    Papers:
  • [1] A. Sommariva and M. Vianello, Low-cardinality Positive Interior cubature on NURBS-shaped domains, [PDF] (submitted),
  • [2] A. Sommariva and M. Vianello, inRS: implementing the indicator function for NURBS-shaped planar domains (submitted).

    • Codes:
  • Last version: [MATLAB CODES (zip file)],

    • Presentation:
  • Presentation at SA 2022: [PDF (Beamer)]
» First update: January 07, 2022. Last update: February 22, 2022.




» Object: CUB_RS: Matlab toolbox for algebraic cubature on NURBS-shaped domains.

The main folder CUB_RS contains:

  • PAPER_EXPERIMENTS that includes all the Matlab routines useful for reproducing the experiments described in the work [1] and [2] (see references below);
  • BEZIER_ROUTINES that includes all the Matlab routines useful for integration over domains whose boundary is defined by composite Bezier arcs;
  • NURBS_ROUTINES that includes all the Matlab routines useful for integration over domains whose boundary is defined by NURBS;
  • SPLINE_ROUTINES that includes all the Matlab routines useful for integration over domains whose boundary is defined by general parametric splines.
All these folders contain several demos, useful to define the domains, to understand how to use the indomain routine as well as the cubature routine. Some examples on hyperinterpolation are also available.

The main routines are:

  • inRS that implements a fast indomain algorithm for these instances; the last update presents a faster indomain algorithm inRS, presented in [2]; the older version introduced in [1] is present in the file inRS_21;
  • cubRS that determines an algebraic cubature formula of the desired degree of exactness, with low cardinality, internal nodes and positive weights.

» Sources:

    Paper:
  • [1] A. Sommariva and M. Vianello, Low-cardinality Positive Interior cubature on NURBS-shaped domains, [PDF] (submitted),
  • [2] A. Sommariva and M. Vianello, inRS: implementing the indicator function for NURBS-shaped planar domains (submitted).

    • Codes:
  • Last version: [MATLAB CODES (zip file)],
» Last update: January 05, 2022.
Presentation:
  • Presentation at SA 2022: [PDF (Beamer)]




    • » Object: Cubature on a domain defined by the union/difference of a disk with a polygon..

    The main routine is
  • polygcirc_2021.m that computes an algebraic cubature rule on a domain defined by the union/difference of a disk with a polygon.

    The software includes the demo file:
  • demo_plot that tests the cubature results on some domains.

    » Sources:

      Paper:
    • E. Artioli, A. Sommariva and M. Vianello, Algebraic cubature on polygonal elements with a circular edge, [PDF],

      • Codes:
    • Last version: [MATLAB CODES (zip file)],
    » Last update: December 02, 2021.



    • » Object: Cubature on a domain defined by the union/difference of a disk with a polygon..

    The main routine is
  • polygcirc_2021.m that computes an algebraic cubature rule on a domain defined by the union/difference of a disk with a polygon.

    The software includes the demo file:
  • demo_plot that tests the cubature results on some domains.

    » Sources:
    • Papers: E. Artioli, A. Sommariva and M. Vianello, Algebraic cubature on polygonal elements with a circular edge, [PDF],
    • Last version: [MATLAB CODES (zip file)],
    » Last update: December 02, 2021.



    • » Object: Basic adaptive code for quadrature of an integrand over polygonal regions.

    The main routine is
  • cub_polygon_adaptive.m that computes adaptive quadrature of an integrand over polygonal regions.

    The software includes the demos file:
  • demo_adaptive_polygons.m that tests the cubature results on some domains and integrands.

    » Sources:
    » Last update: July 13, 2021.



    • » Object: RBF cubature on polygonal regions (with almost optimal shape parameter).

    The main routines are:
  • RBF_cub_polygon_OPT.m that determines the RBF based cubature rule on the polygonal region;
  • RBF_moms_polar.m and RBF_moms_tri.m that determines the RBF moments on the polygonal region;
  • RBF_optimal_scale.m that chooses the optimal RBF scale parameter;

    The software includes several demos:
  • demo_RBF_cub_OPT.m that performs a battery of tests, allowing different integrands, polygonal regions, RBF, number of centers and choices of the RBF shape parameter.
  • demo_RBF_cub_OPT_basic that performs a simple example, on a certain integrand, polygonal region, RBF, number of centers and with almost optimal RBF shape parameter (useful to adapt the code to user experiments).

    For the demos of the previous version published in RBF moment computation and meshless cubature on general polygonal regions, see:
  • RBF_cub_polygon.m that determines the RBF based cubature rule on the polygonal region (not optimized shape parameter);

    » Source: » Last update: July 9, 2021.



    • » Object: Near-algebraic quadrature rules on spherical polygons..

    The main routines are
  • cub_sphpgon, based on cub_sphtri.m ,that determines full and compressed quadrature rules with fixed (numerical) algebraic degree of precision on spherical polygons; it requires dCATCHsph.m, dCHEBVAND.m, gqcircsect.m, LHDM.m, trigauss.m.
  • adaptive_cub_sphpgon.m, that implements a general purpose adaptive code for quadrature over spherical polygons.

    The software includes the demos files:
  • demo_cub_sphpgon.m that tests the cubature rules on a battery of tests functions,
  • demo_cub_continents.m that tests the cubature of several functions on approximation of continents.

    » Object: Polynomial hyperinterpolation on spherical polygons.

    The software includes the demos files:
  • demo_hyp_sphpgon that tests the hyperinterpolants of several functions on some spherical polygons, also approximation of continents.

    » Sources:
    » Last update: June 17, 2021.



    • » Object: Approximation/Interpolation on spherical triangles.

    The zip file contains three folders:
    • BASIC_USAGE: here we have introduced a demo demo_am_00.m on the norming mesh and its properties as well as the main routine am_sphtri.m for the determination of the norming mesh, depending on a spherical triangle contained in a certain hemisphere.
    • ADVANCED_USAGE: here there are two demos for the computation of norming meshes on spherical triangle T, as well as a set of Approximate Fekete Points, useful for polynomial interpolation.
      1. demo_am_01.m demo on using "dPOLYFIT", "dPOLYVAL", "dLEB" for algebraic polynomial interpolation/least squares approximation on spherical triangle. The approximation points are taken from AMs on the spherical triangle; the new routine introduced here is "am_sphtri" that determines an Admissible Mesh over a spherical triangle.
      2. demo_am_02.m this routine computes a AM on a spherical triangle and extracts AFEK pointset; next it determines its main features and plots AM and AFEK.
    • FIGURES: it contains all the Matlab files useful for determining the figures and their experiments, as cited in the reference paper Chebyshev-Dubiner norming grids and Fekete-like interpolation on spherical triangles.

    » Last update: February 18, 2021.



    » Object: Near-algebraic quadrature rules on spherical triangles..

    The main routine is
  • cub_sphtri.m that determines full and compressed quadrature rules with fixed (numerical) algebraic degree of precision on spherical triangles; it requires dCATCHsph.m, dCHEBVAND.m, gqcircsect.m, LHDM.m, trigauss.m.

    The software includes the demos files:
  • demo_01.m that tests the cubature rules on spherical harmonics (it requires sphharmVAND.m and plot_s2.m),
  • demo_02.m that tests the cubature rules on a battery of tests functions.

    » Object: Polynomial hyperinterpolation on spherical triangles.

    The software includes the demos files:
  • demo_hyperinterpolation_01.m that tests the hyperinterpolants of several functions on some spherical triangles,
  • demo_hyperinterpolation_02.m that shows how to determine the hyperinterpolant on a single example.
  • demo_hyperinterpolation_03.m that shows estimates of the Lesbegue constant based on Christoffel function.

    » Sources:
    » Last update: March 22, 2021.



    • » Object: RBF moments and cubature rules on polygonal regions.

    The main routine are
  • RBF_moms_polar.m and RBF_moms_tri.m that determines the RBF moments on the polygonal region;
  • RBF_cub_polygon.m that determines the RBF based cubature rule on the polygonal region;

    The software includes several demos that test moments as well as the cubature rules. In particular the tests in the paper can be obtained by means of demo_RBF_moms.m and demo_RBF_cub.m.

    » Source: » Last update: December 28, 2020.



    • » Object: cubature rules on spherical rectangles of the 2-sphere.

    The main routine is cub_sphrect.m that computes rules with a fixed algebraic degree of precision on a spherical rectangle.

    Another relevant function is cub_sphrect_adaptive_2020.m that is a basic adaptive code for cubature on a spherical rectangle.

    The software includes two demos, on integration of a general function as well as test on the algebraic degree of precision of the rule.

    » Source: » Last update: December 20, 2020.



    Incurvpolygon_cheb: indomain routine on bivariate Jordan domains whose boundary is defined by chebfuns. The files are stored in the compressed folder, that contains:

    » incurvpolygon_cheb.m: it is an indomain routine on domains defined whose boundary is defined by chebfuns;

    » in the folder there is a demo demo_incurvpolygon_cheb.m, that tests the provided software.

    » Source:
    » Last update: October 10, 2020.



    ChebfunGauss: a Matlab code for cubature on Jordan bivariate domains defined analytically and approximated by Chebfuns. The files are stored in the compressed folder, that contains:

    » chebcub_2019.m: it is a more stable version of the previous software and quickly produces a rule with many points, of the choosen degree of precision;

    » chebfuncub_compression.m: by this routine one can compress the previous rule to one with the same degree of precision but fewer nodes; first it tries to compute one with internal nodes and positive weights, then one with internal nodes and possible negative weights, then one with possible external nodes and positive weights and finally one with possible external nodes and possible negative weights;

    » chebcub_2020.m: it is a version that produces a rule with few points, of the choosen degree of precision (the rule has internal nodes and usually positive weights);

    » incurvpolygon_cheb.m: it is an indomain routine on bivariate Jordan domains defined whose boundary is defined by chebfuns;

    » in the folder there are several demo,
    • demo_2019.m,
    • demo_2020.m,
    • demo_random_polynomials.m,
    • demo_19vs20.m,
    • demo_incurvpolygon_cheb.m,
    that test the provided software.

    » Source:
    » Last update: October 10, 2020.



    In this software, we define a Matlab code for establishing numerically if a point is inside or on the boundary of a spline-curvilinear polygon (the so called "in domain" algorithms). The main routine is incurvpolygon.m. Several demos are available, that show how to use incurvpolygon.m.

    » Paper: Tchakaloff-like cubature rules on spline curvilinear polygons with M. Vianello.
    » Source: [.zip],
    » First Version: April 2019, Last update: September 29, 2020.



    Fast Matlab codes for computing algebraic cubature rules on spline-curvilinear polygons. In general the cubature formula may have points outside the domain as well as negative weights.
    The main routine is the m. file splinegauss_2020.m that performs cubature over domains with spline curvilinear boundaries. In case one wants to recover rules with less points (possibly with points in the domain and positive weights) one can try to achieve such rule by the routine splinegauss_compression.m.
    Several demos are available, that show how to use splinegauss_2020.m and splinegauss_compression.m to produce cubature rules over spline-curvilinear polygons.

    » Paper: Gauss-Green cubature and moment computation over arbitrary geometries, with M. Vianello.
    » Source: [.zip],
    » First Version: March 2009, Last update: September 20, 2020.



    Object: this software computes algebraic cubature rules with fixed degree of precision over spline-curvilinear polygons. It includes an in-domain routine incurvpolygon.m that determines if a point is inside one of such regions.
    The cubature formula has all points inside the domain as well as positive weights.

    » Paper: Tchakaloff-like cubature rules on spline curvilinear polygons with M. Vianello.
    » Source: [.zip],
    » First Version: May 2020, Last update: May 2020.



    Object: this software computes algebraic cubature rules with fixed degree of precision over general polygons (convex or not convex, connected or disconnected, simply or not simply connected). It needs Matlab-built in polyshape environment.

    Applications to optics are given.
    » Paper: Compressed cubature over polygons with applications to optical design
    » Source: [.zip]
    » Last update: May 12, 2019.



    Object: this software computes cubature rules over the sphere, via Radial Basis Functions and Spherical Radial Basis Functions. We point out that as described in demoRBFweights.m and demoSRBFweights.m , cubature weights are obtained, without the need of the RBF interpolant (though via moments and transpose of Vandermonde matrix).
    » Paper: Integration by RBF over the sphere
    » Source: [.zip],
    » First Version: May 2005, Last update: September 26, 2019.



    Object: this software computes cubature rules over disks and annular sectors, via Radial Basis Functions. We point out that as described in demoRBFweights.m, cubature weights are obtained, without the need of the RBF interpolant (though via moments and transpose of Vandermonde matrix).
    » Paper: Meshless cubature over the disk by Thin-Plates Splines
    » Source: [.zip],
    » Old sources: [GreenDisk_v1_0.zip], [GreenDisk_v1_0a.zip]
    » First Version: 2008, Last update: September 24, 2019.



    Object: this software computes cubature rules over rectangles, via Radial Basis Functions. The first .zip file considers RBF with the same scale for all the centers, the second zip file manages the variable scale case. We point out that as described in demo_2019_weights.m, cubature weights can be obtained, without the need of the RBF interpolant.
    » Paper: Numerical cubature on scattered datas by Radial Basis Functions
    » Source: [.zip], [.zip]
    » First Version: May 2004, Last update: September 23, 2019.



    Object: this software computes cubature rules over polygons (convex or not convex), via Radial Basis Functions (Thin Plate Splines)
    » Paper: Meshless cubature by Green's formula.
    » Source: [.zip]
    » First Version: May 2006, Last update: September 11, 2019.



    » Source: [ZIP]
    » Last update: July 27, 2012.



    » Source: [ZIP]
    » Last update: July 27, 2012 (Tested August 22, 2016)]



    » Object: The MULTIVOP_220816.zip suite allows to compute hyperinterpolants over general regions. Here we show examples on some circular regions or polygons. The more important routine is multivop.m (see the paper http://www.math.unipd.it/~marcov/pdf/multivop.pdf) that allows the numerical computation of orthogonal polynomials over general regions.
    » Source: Matlab code for Gauss-like cubature over polygons [MULTIVOP_220816]
    » Last update: August 25, 2016.



    » Object: The routine trigauss.m computes a quadrature rule over an interval, exact for trigonometric polynomial degree N. Several versions, including the Gauss-Kronrod version (it is not known if they always exist!), are available.
    » Source: » Last update: September 13, 2016.



    » Object: Software for 1D trigonometric hyperinterpolation and hyperinterpolation over spherical rectangles. HYPERTRIG is based on trigonometric tensorial hyperinterpolation, while HYPERPOL consists in polynomial hyperinterpolation over spherical rectangles (of total degree). For details, see the paper G. Da Fies, A. Sommariva and M. Vianello, Subperiodic trigonometric hyperinterpolation.
    » Source: » Last update: October 24, 2016.



    » Object: Demo of least squares on a 2D mesh and least squares on a compressed mesh. Examples on union of disks and polygons. For details, see the paper F. Piazzon, A. Sommariva and M. Vianello, Caratheodory-Tchakaloff Subsampling.
    » Source: » Last update: November 20, 2016.



    » Object: Computation of optimal meshes on sphere, ball and torus. For details, see the paper A. Sommariva and M. Vianello, Discrete norming inequalities on sections of sphere, ball and torus.
    » Source: » Last update: April 24, 2017.



    » Object: Computation of nested cubature rules on polygons, from bivariate data. For details, see the paper A. Sommariva and M. Vianello, Nearly optimal nested sensors location for polynomial regression on complex geometriess.
    » Source: » Last update: May 23, 2017.



    Matlab code for cubature over a square [Padua_CC.zip] (using non tensorial Padua and Morrow-Patterson Xu points). The main routine is the m. file main.m that performs cubature over a square using different formulas. All numerical experiments have been tested on Intel Centrino Duo T2400, with Windows XP, Matlab 6.1.0.450 Release 12.1 and GNU Octave 2.1.73. [Last update: October 09, 2007]



    Software for Padua points interpolation and cubature. [Padua2DM.tar.gz].



    Some routines that can be used to make experiments with WAMs can be found in the following compressed files:
    • [WAM_2D_280210.zip]: a complete version that works on several bivariate domains, not easily adaptable to add own WAMs but right for experiments on general 2D domains.
    • [WAM_2De_010310.zip]: a simple version that does not work on many bivariate domains, but that is easily adaptable to add 2D WAMs.
    • [WAM_3D_010310.zip]: a simple version that does not work on many trivariate domains, but that is easily adaptable to add 3D WAMs.
    [Last update: March 01, 2010]




    Old versions of our software



    ChebfunGauss: a Matlab code for Gauss-Green cubature by the Chebfun package: Matlab code that uses a suitable version of Splinegauss for bivariate domains defined analytically and approximated by Chebfuns. The files are stored in the compressed folder [CHEBFUNGAUSS_2017.zip]. For a demo use the file demo.m, for plotting cubature points plot_pts.m.
    » Source:
    » Last update: November 26, 2017.




    » Object: The polygauss.m routine allows the computation of cubature rules with a certain degree of exactness on a polygon (without self-intersections). The new routine works via quadrangles panelling or triangularization. Some new demos show its usage. We point out that these rules can be compressed, providing a rule with less nodes/weights but the same algebraic degree of exactness (see comprexcub.m).
    » Source: » Last update: August 25, 2016.



    Matlab code for Gauss-like cubature over polygons [POLYGAUSS_2013_300413.zip]

    Previous version: Matlab code for Gauss-like cubature over polygons [GAUSSCUB_2013.zip]. The main routine is the m. file polygauss.m. As demo, please use demo_gausscub.m. A different version, that uses less points but is more time expensive is polygaussj.m. We point out that the routines work in the GNU Octave version 2.1.73 (i686-pc-cygwin) environment. For some numerical experiments in Matlab/Octave see [readme_polygauss.pdf].
  • April 30, 2013: Modification of polygauss.m outputs.
    • For a numerical comparison, see the routine that calls "dblquad" [dblquad_polygons.zip]. [Last update: March 23, 2013]



    Matlab code for Gauss-like cubature over spline curvilinear polygons [SPLINEGAUSS_2009.zip]. The main routine is the m. file splinegauss_2009b that performs cubature over domains with spline curvilinear boundaries.. As demo, please use demo.m for computing cubature over spline curvilinear polygons. In addition, try demo_splinegauss_moments.m for tests on computing moments w.r.t. some specific bases on such domains. All numerical experiments have been done on MAC Book, with MAC OS X, Matlab 7.6.0.324 Release 2008a. [Latest update: November 19, 2009].
    If you are interested in this software, see the more recent routine in SPLINEGAUSS_2020.zip. [Last update: March 27, 2009]



    [POLYGINT_010411.zip] contains several Matlab files for cubature/interpolation over polygons based on AFP (Approximate Fekete Points) and DLP (Discrete Leja Points). The zip file contains also Matlab routines for a minimal triangulation of a polygon and a quadrangulation of a polygon using few quadrangles and triangles. [Last Update: April 01, 2011]



    LEBFEK: Matlab codes and data sets for bivariate Lebesgue and Fekete interpolation.

    In M. Briani, A. Sommariva, M. Vianello, Computing Fekete and Lebesgue points: simplex, square, disk we have described some algorithms for computing (quasi-)Lebesgue and (quasi-)Fekete points on the unit simplex, a square and on the unit disk. The respective codes and sets for these domains can be found in the links below. [Last Update: February 14, 2011]



    Computation of (almost-)minimal rules on the square.

    Matlab codes for computing (almost-)minimal cubature rules on the square.
    1. Asymmetric rule: [ZIP]
    2. Pi symmetric rule: [ZIP]
    3. Pi/2 symmetric rule: [ZIP]
    [Last Update: December 11, 2011]



    Gaussian Rules for Symmetric weight functions.

    Matlab codes for computing Gaussian Rules for Symmetric weight functions.
    1. [ZIP]
    [Last Update: July 27, 2012]



    Matlab codes for computing Fejer and Clenshaw-Curtis rules for general weight functions.
    1. [ZIP]
    [Last Update: July 27, 2012]



    Matlab codes for computing trigonometric quadrature rules.

    Matlab codes for computing trigonometric quadrature rules. A comparison.
    1. [ZIP]
    [Last Update: February 27, 2013]

    Matlab codes for computing good interpolation points or quadrature rules on geographic subregions of the sphere.
      [MATLAB CODES (zip file), interpolation] [MATLAB CODES (zip file), cubature]
    [Last Update: November 20, 2013]