Mercurial > repos > public > sbplib_julia

diff ext/SbplibSparseArraysExt.jl @ 1700:87d603499fc3 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge feature/lazy_tensors/sparse_conversions
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 03 Sep 2024 08:38:27 +0200
parents 3e9c3986930d
children 8787b2f9f5f8
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/ext/SbplibSparseArraysExt.jl	Tue Sep 03 08:38:27 2024 +0200
@@ -0,0 +1,17 @@
+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