Mercurial > repos > public > sbplib_julia

view test/BoundaryConditions/boundary_condition_test.jl @ 1417:5f79549f60ae feature/boundary_conditions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Test constant dirichlet data
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 23 Aug 2023 09:27:32 +0200
parents b4ec84190e6b
children fefc81a9c155
line wrap: on
line source

using Test

using Sbplib.BoundaryConditions
using Sbplib.Grids

@testset "BoundaryCondition" begin
    grid_2d = equidistant_grid((11,15), (0.0, 0.0), (1.0,1.0))
    grid_3d = equidistant_grid((11,15,13), (0.0, 0.0, 0.0), (1.0,1.0, 1.0))
    (_,_,_,id_n) = boundary_identifiers(grid_2d)
    (_,_,_,_,id_b,_) = boundary_identifiers(grid_3d)

    g = 3.14
    f(x,y) = x^2+y^2
    @test DirichletCondition(g,id_n) isa BoundaryCondition{Float64}
    @test NeumannCondition(f,id_n) isa BoundaryCondition{<:Function}

   @test g*ones(11,1) ≈ discretize_data(grid_2d,DirichletCondition(g,id_n))
end