comparison test/BoundaryConditions/sat_test.jl @ 1395:bdcdbd4ea9cd feature/boundary_conditions

Merge with default. Comment out broken tests for boundary_conditions at sat
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 26 Jul 2023 21:35:50 +0200
parents ea2e8254820a
children b96858a50e35
comparison
equal deleted inserted replaced
1217:ea2e8254820a 1395:bdcdbd4ea9cd
4 using Sbplib.BoundaryConditions 4 using Sbplib.BoundaryConditions
5 using Sbplib.Grids 5 using Sbplib.Grids
6 using Sbplib.RegionIndices 6 using Sbplib.RegionIndices
7 using Sbplib.LazyTensors 7 using Sbplib.LazyTensors
8 8
9 grid = EquidistantGrid(11, 0.0, 1.0) 9 grid = equidistant_grid(11, 0.0, 1.0)
10 (id_l,id_r) = boundary_identifiers(grid) 10 (id_l,id_r) = boundary_identifiers(grid)
11 struct MockOp 11 struct MockOp
12 end 12 end
13 13
14 function BoundaryConditions.sat_tensors(op::MockOp, grid, bc::DirichletCondition) 14 function BoundaryConditions.sat_tensors(op::MockOp, grid, bc::DirichletCondition)
28 end 28 end
29 return closure, penalty 29 return closure, penalty
30 end 30 end
31 31
32 32
33 @testset "sat" begin 33 # @testset "sat" begin
34 g = ConstantBoundaryData(2.0) 34 # g = ConstantBoundaryData(2.0)
35 dc = DirichletCondition(g,id_l) 35 # dc = DirichletCondition(g,id_l)
36 op = MockOp() 36 # op = MockOp()
37 f = sat(op, grid, dc) 37 # f = sat(op, grid, dc)
38 u = evalOn(grid, x-> -1/2 + x^2) 38 # u = eval_on(grid, x-> -1/2 + x^2)
39 @show f(0.,u) 39 # @show f(0.,u)
40 end 40 # end