Mercurial > repos > public > sbplib_julia

view ext/SbplibSparseArraysExt.jl @ 1731:3684db043add feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add Base.@constprop :aggressive to src/SbpOperators/volumeops/derivatives/second_derivative_variable.jl
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 09 Sep 2024 09:01:57 +0200
parents 3e9c3986930d
children 8787b2f9f5f8
line wrap: on
line source

module SbplibSparseArraysExt

using Sbplib
using Sbplib.LazyTensors

using SparseArrays
using Tokens

function SparseArrays.sparse(t::LazyTensor)
    v = ArrayToken(:v, prod(domain_size(t)))

    v̄ = reshape(v,domain_size(t)...)
    tv = reshape(t*v̄, :)
    return Tokens._to_matrix(tv, prod(range_size(t)), prod(domain_size(t)))
end

end