Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/laplace/laplace.jl @ 756:1970ebceabe4 feature/laplace_opset

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add suggestion for pretty printing of Laplace
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 02 Jul 2021 11:13:14 +0200
parents 36adc15d3935
children 1784b1c0af3e 9929c99754fb
comparison
equal deleted inserted replaced
755:36adc15d3935 756:1970ebceabe4
49 Hᵧ_pairs = ntuple(i -> ids[i] => inner_product(boundary_grid(grid,ids[i]),H_closure_stencils),n_ids) 49 Hᵧ_pairs = ntuple(i -> ids[i] => inner_product(boundary_grid(grid,ids[i]),H_closure_stencils),n_ids)
50 50
51 return Laplace(Δ, H, H⁻¹, StaticDict(e_pairs), StaticDict(d_pairs), StaticDict(Hᵧ_pairs)) 51 return Laplace(Δ, H, H⁻¹, StaticDict(e_pairs), StaticDict(d_pairs), StaticDict(Hᵧ_pairs))
52 end 52 end
53 53
54 54 # TODO: Consider pretty printing of the following form
55 Base.show(io::IO, L::Laplace{T, Dim}) where {T,Dim} = print(io, "Laplace{$T, $Dim}(", L.D, L.H, L.H_inv, L.e, L.d, L.H_boundary, ")") 55 # Base.show(io::IO, L::Laplace{T, Dim}) where {T,Dim,TM} = print(io, "Laplace{$T, $Dim, $TM}(", L.D, L.H, L.H_inv, L.e, L.d, L.H_boundary, ")")
56 56
57 LazyTensors.range_size(L::Laplace) = LazyTensors.range_size(L.D) 57 LazyTensors.range_size(L::Laplace) = LazyTensors.range_size(L.D)
58 LazyTensors.domain_size(L::Laplace) = LazyTensors.domain_size(L.D) 58 LazyTensors.domain_size(L::Laplace) = LazyTensors.domain_size(L.D)
59 LazyTensors.apply(L::Laplace, v::AbstractArray, I...) = LazyTensors.apply(L.D,v,I...) 59 LazyTensors.apply(L::Laplace, v::AbstractArray, I...) = LazyTensors.apply(L.D,v,I...)
60 60