view ext/SbplibSparseArraysExt.jl @ 1721:7379f492c4b3 feature/lazy_tensors/sparse_conversions

Fix typos in docstring
author Vidar Stiernström <vidar.stiernstrom@gmail.com>
date Thu, 05 Sep 2024 08:31:03 -0700
parents 8787b2f9f5f8
children
line wrap: on
line source

module SbplibSparseArraysExt

using Sbplib
using Sbplib.LazyTensors

using SparseArrays
using Tokens

"""
    sparse(t::LazyTensor)

The sparse matrix representation of `t`.

If `L` is a `LazyTensor` and `v` a tensor, then `A = sparse(L)` is constructed
so that `A*reshape(v,:) == reshape(L*v,:)`.
"""
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