Mercurial > repos > public > sbplib_julia

diff src/Grids/mapped_grid.jl @ 1743:49bd573ab07f feature/grids/curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Implement checking of sizes for inputs to MappedGrid
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 11 Sep 2024 13:46:33 +0200
parents 248304f88e49
children c38eead8be17
line wrap: on
line diff
--- a/src/Grids/mapped_grid.jl	Wed Sep 11 12:09:59 2024 +0200
+++ b/src/Grids/mapped_grid.jl	Wed Sep 11 13:46:33 2024 +0200
@@ -2,6 +2,14 @@
     logicalgrid::GT
     physicalcoordinates::CT
     jacobian::JT
+
+    function MappedGrid(logicalgrid::GT, physicalcoordinates::CT, jacobian::JT) where {T,D, GT<:Grid{<:Any,D}, CT<:AbstractArray{T,D}, JT<:AbstractArray{<:AbstractArray{<:Any, 2}, D}}
+        if !(size(logicalgrid) == size(physicalcoordinates) == size(jacobian))
+            throw(ArgumentError("Sizes must match"))
+        end
+
+        return new{T,D,GT,CT,JT}(logicalgrid, physicalcoordinates, jacobian)
+    end
 end
 
 function Base.:(==)(a::MappedGrid, b::MappedGrid)