Mercurial > repos > public > sbplib_julia

diff src/SbpOperators/volumeops/quadratures/inverse_quadrature.jl @ 675:1ce3a104afc8 feature/boundary_quads

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge in default
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 07 Feb 2021 21:28:53 +0100
parents e14627e79a54
children 0bec3c4e78c0
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/quadratures/inverse_quadrature.jl	Sat Feb 06 12:04:06 2021 +0100
+++ b/src/SbpOperators/volumeops/quadratures/inverse_quadrature.jl	Sun Feb 07 21:28:53 2021 +0100
@@ -32,7 +32,7 @@
 the closure stencils are those of the quadrature operator (and not the inverse).
 """
 function InverseDiagonalQuadrature(grid::EquidistantGrid, closure_stencils::NTuple{M,Stencil{T,1}}) where {T,M}
-    inv_inner_stencil = Stencil(Tuple{T}(1),center=1)
+    inv_inner_stencil = Stencil(one(T), center=1)
     inv_closure_stencils = reciprocal_stencil.(closure_stencils)
     return InverseQuadrature(grid, inv_inner_stencil, inv_closure_stencils)
 end