Mercurial > repos > public > sbplib_julia

comparison src/Grids/mapped_grid.jl @ 1705:4870fc3faa25 feature/grids/curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add tests for equallity of mapped grids
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 04 Sep 2024 15:38:10 +0200
parents e5e76c8e52c5
children 11640aa3e348
comparison
equal deleted inserted replaced
1704:e5e76c8e52c5 1705:4870fc3faa25
1 struct MappedGrid{T,D, GT<:Grid{<:Any,D}, CT<:AbstractArray{T,D}, JT<:AbstractArray{<:AbstractArray{<:Any, 2}, D}} <: Grid{T,D} 1 struct MappedGrid{T,D, GT<:Grid{<:Any,D}, CT<:AbstractArray{T,D}, JT<:AbstractArray{<:AbstractArray{<:Any, 2}, D}} <: Grid{T,D}
2 logicalgrid::GT 2 logicalgrid::GT
3 physicalcoordinates::CT 3 physicalcoordinates::CT
4 jacobian::JT 4 jacobian::JT
5 end
6
7 function Base.:(==)(a::MappedGrid, b::MappedGrid)
8 same_logicalgrid = logicalgrid(a) == logicalgrid(b)
9 same_coordinates = collect(a) == collect(b)
10 same_jacobian = jacobian(a) == jacobian(b)
11
12 return same_logicalgrid && same_coordinates && same_jacobian
5 end 13 end
6 14
7 jacobian(g::MappedGrid) = g.jacobian 15 jacobian(g::MappedGrid) = g.jacobian
8 logicalgrid(g::MappedGrid) = g.logicalgrid 16 logicalgrid(g::MappedGrid) = g.logicalgrid
9 17