Mercurial > repos > public > sbplib_julia

comparison DiffOps/src/laplace.jl @ 250:863a98d9e798 boundary_conditions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove obsolete comment
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 27 Jun 2019 14:24:09 +0200
parents 05e7bbe0af97
children 89a101a63e7a
comparison
equal deleted inserted replaced
249:7cb4492ccd60 250:863a98d9e798
81 # TODO: This is obviouly strange. Is domain_size just discarded? Is there a way to avoid storing grid in BoundaryValue? 81 # TODO: This is obviouly strange. Is domain_size just discarded? Is there a way to avoid storing grid in BoundaryValue?
82 # Can we give special treatment to TensorMappings that go to a higher dim? 82 # Can we give special treatment to TensorMappings that go to a higher dim?
83 LazyTensors.range_size(e::NormalDerivative{T}, domain_size::NTuple{1,Integer}) where T = size(e.grid) 83 LazyTensors.range_size(e::NormalDerivative{T}, domain_size::NTuple{1,Integer}) where T = size(e.grid)
84 LazyTensors.domain_size(e::NormalDerivative{T}, range_size::NTuple{2,Integer}) where T = (range_size[3-dim(e.bId)],) 84 LazyTensors.domain_size(e::NormalDerivative{T}, range_size::NTuple{2,Integer}) where T = (range_size[3-dim(e.bId)],)
85 85
86 # Not correct abstraction level
87 # TODO: Not type stable D:< 86 # TODO: Not type stable D:<
88 function LazyTensors.apply(d::NormalDerivative, v::AbstractArray, I::NTuple{2,Int}) 87 function LazyTensors.apply(d::NormalDerivative, v::AbstractArray, I::NTuple{2,Int})
89 i = I[dim(d.bId)] 88 i = I[dim(d.bId)]
90 j = I[3-dim(d.bId)] 89 j = I[3-dim(d.bId)]
91 N_i = size(d.grid)[dim(d.bId)] 90 N_i = size(d.grid)[dim(d.bId)]