Mercurial > repos > public > sbplib_julia
diff src/SbpOperators/volumeops/volume_operator.jl @ 995:1ba8a398af9c refactor/lazy_tensors
Rename types
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Fri, 18 Mar 2022 21:14:47 +0100 |
parents | 469ed954b493 |
children | 52f07c77299d 3bb94ce74697 2b6298905692 |
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/volume_operator.jl Fri Mar 18 20:44:17 2022 +0100 +++ b/src/SbpOperators/volumeops/volume_operator.jl Fri Mar 18 21:14:47 2022 +0100 @@ -6,7 +6,7 @@ the stencils `inner_stencil` and `closure_stencils`. When `Dim=1`, the corresponding `VolumeOperator` tensor mapping is returned. When `Dim>1`, the returned operator is the appropriate outer product of a one-dimensional -operators and `IdentityMapping`s, e.g for `Dim=3` the volume operator in the +operators and `IdentityTensor`s, e.g for `Dim=3` the volume operator in the y-direction is `I⊗op⊗I`. """ function volume_operator(grid::EquidistantGrid, inner_stencil, closure_stencils, parity, direction) @@ -14,9 +14,9 @@ # Create 1D volume operator in along coordinate direction op = VolumeOperator(restrict(grid, direction), inner_stencil, closure_stencils, parity) - # Create 1D IdentityMappings for each coordinate direction + # Create 1D IdentityTensors for each coordinate direction one_d_grids = restrict.(Ref(grid), Tuple(1:dimension(grid))) - Is = IdentityMapping{eltype(grid)}.(size.(one_d_grids)) + Is = IdentityTensor{eltype(grid)}.(size.(one_d_grids)) # Formulate the correct outer product sequence of the identity mappings and # the volume operator parts = Base.setindex(Is, op, direction) @@ -27,7 +27,7 @@ VolumeOperator{T,N,M,K} <: TensorOperator{T,1} Implements a one-dimensional constant coefficients volume operator """ -struct VolumeOperator{T,N,M,K} <: TensorMapping{T,1,1} +struct VolumeOperator{T,N,M,K} <: LazyTensor{T,1,1} inner_stencil::Stencil{T,N} closure_stencils::NTuple{M,Stencil{T,K}} size::NTuple{1,Int}