Mercurial > repos > public > sbplib_julia

changeset 1662:bdb4becac704 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update inner_product and inverse_inner_product to work on boundary grids
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 28 Jun 2024 21:35:33 +0200
parents 6d196fb85133
children de6300bd36cc
files src/SbpOperators/volumeops/inner_products/inner_product.jl src/SbpOperators/volumeops/inner_products/inverse_inner_product.jl
diffstat 2 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/inner_products/inner_product.jl	Fri Jun 28 17:04:05 2024 +0200
+++ b/src/SbpOperators/volumeops/inner_products/inner_product.jl	Fri Jun 28 21:35:33 2024 +0200
@@ -52,6 +52,6 @@
 
 
 function inner_product(g::MappedGrid, stencil_set)
-    J = jacobian_determinant(g)
+    J = map(sqrt∘det, geometric_tensor(g))
     DiagonalTensor(J)∘inner_product(logicalgrid(g), stencil_set)
 end
--- a/src/SbpOperators/volumeops/inner_products/inverse_inner_product.jl	Fri Jun 28 17:04:05 2024 +0200
+++ b/src/SbpOperators/volumeops/inner_products/inverse_inner_product.jl	Fri Jun 28 21:35:33 2024 +0200
@@ -51,6 +51,6 @@
 inverse_inner_product(g::ZeroDimGrid, stencil_set::StencilSet) = IdentityTensor{component_type(g)}()
 
 function inverse_inner_product(g::MappedGrid, stencil_set)
-    J⁻¹ = map(inv, jacobian_determinant(g))
+    J⁻¹ = map(inv∘sqrt∘det, geometric_tensor(g))
     DiagonalTensor(J⁻¹)∘inner_product(logicalgrid(g), stencil_set)
 end