Mercurial > repos > public > sbplib_julia
view ext/DiffinitiveSparseArraysExt.jl @ 1831:27e74f524c4f refactor/grids/iterable_boundary_indices
Some cleanup in tensor_grid_test.jl
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Sun, 20 Oct 2024 21:50:18 +0200 |
parents | 471a948cd2b2 |
children |
line wrap: on
line source
module DiffinitiveSparseArraysExt using Diffinitive using Diffinitive.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