Matrix multiplication and triangle counting by Danil-Zaripov · Pull Request #7 · Lamagraph/QTreeFSharp

Danil-Zaripov · 2026-02-04T15:53:35Z

No description provided.

gsvgit · 2026-02-04T18:13:28Z

QuadTree/LinearAlgebra.fs

+            (Matrix.mkNode tree (Matrix.qtree.Leaf Dummy) (Matrix.qtree.Leaf Dummy) (Matrix.qtree.Leaf Dummy))
+            (expandRatio / 2UL)
+
+let mxm op_add op_mult (m1: Matrix.SparseMatrix<'a>) (m2: Matrix.SparseMatrix<'a>) =


Типы матриц могут быть разными.

gsvgit · 2026-02-04T18:14:14Z

QuadTree/LinearAlgebra.fs

    else
        (Error.InconsistentSizeOfArguments(vector, matrix)) |> Result.Failure
+
+let rec shrink tree (size: uint64<storageSize>) =


Правда ли эти функции должны быть на верхнем уровне?

gsvgit · 2026-02-05T05:48:19Z

QuadTree.Tests/Tests.LinearAlgebra.fs

+    // DDDD
+    let tree =
+        qtree.Node(
+            qtree.Leaf << UserValue <| Some 2,


А может сделать для таких вещей вспомогательные функции?
Типа

leaf_v v = qtree.Leaf << UserValue <| Some v leaf_n = qtree.Leaf << UserValue <| None

и т.д.

gsvgit · 2026-02-05T05:50:32Z

QuadTree/LinearAlgebra.fs

+            let halfSize = size / 2UL
+
+            // Double check this code
+            let nw1xnw2, ne1xsw2, nw1xne2, ne1xse2, sw1xnw2, se1xsw2, sw1xne2, se1xse2 =


Может nw1_x_nw2?

gsvgit · 2026-02-05T05:51:12Z

QuadTree/LinearAlgebra.fs

+                (multiply halfSize se1 se2)
+
+            match nw1xnw2, ne1xsw2, nw1xne2, ne1xse2, sw1xnw2, se1xsw2, sw1xne2, se1xse2 with
+            | Result.Success(tnw1xnw2, nvals_nw1xnw2),


И тут тоже t_nw1_x_nw2

gsvgit · 2026-02-05T05:55:19Z

QuadTree/LinearAlgebra.fs

+
+        match m1, m2 with
+        | Matrix.qtree.Leaf(UserValue v1), Matrix.qtree.Leaf(UserValue v2) ->
+            let res = multScalar op_add (uint64 size) (op_mult v1 v2)


Чот меня тип size смущает. Если я правильно понимаю, то тут он уже uint64. Тогда каст не нужен. Хотя я и не уверен, что это хорошая затея. Я бы предпочёл до последнего тянуть типизированный uint, и только в самом конце конвертировать к другим единицам. Больше кода, но и больше контроля со стороны компилятора.

gsvgit · 2026-02-05T05:57:05Z

QuadTree/LinearAlgebra.fs

+        let m1_tree, m2_tree =
+            if m1.storage.size < m2.storage.size then
+                expand (m1.storage.data) (m2.storage.size / (uint64 m1.storage.size)), m2.storage.data
+            else if m1.storage.size > m2.storage.size then


gsvgit · 2026-02-05T12:30:33Z

QuadTree.Tests/Tests.LinearAlgebra.fs

-                    )
+                    Matrix.qtree.Node(leaf_n (), leaf_v 1, leaf_v 3, leaf_v 2),
+                    Matrix.qtree.Node(leaf_v 1, leaf_n (), leaf_v 2, leaf_v 3),
+                    Matrix.qtree.Node(leaf_n (), leaf_n (), Matrix.qtree.Leaf(Dummy), Matrix.qtree.Leaf(Dummy)),


Dummy тоже вынести через хелпер?

gsvgit · 2026-02-05T12:32:44Z

QuadTree.Tests/Tests.TriangleCount.fs

+
+[<Fact>]
+let ``7V Triangle count`` () =
+    // Lower triangle


То есть Вы считаете, что в функцию подсчёта треугольиков сразу подаётся треугольная матрица?

В коде самой функции, вроде, есть получение треугольника. Может тогда тут всё же полноценная матрица смежности графа пусть будет?

Писать долго руками, и не понятно что это дает

Даёт отсутствие ложных ожиданий от функции. Тесты - тоже документация. Можно, например, граф попроще взять.

gsvgit · 2026-02-05T12:39:06Z

QuadTree/Matrix.fs

-        Leaf (s, op v1)
-    | Node (s,x1,x2,x3,x4) ->
-        mkNode s (map op x1) (map op x2) (map op x3) (map op x4)
+            for _ in 1UL .. area do


Не очень нравится такой код. В сети хорошо транслируется и параллелится код в духе рекурсивного "разделяй и властвуй". А с такимим циклами не понятно, что делать.

gsvgit · 2026-02-05T12:40:51Z

QuadTree/Matrix.fs

    else
        (Error.InconsistentSizeOfArguments(matrix1, matrix2)) |> Result.Failure

+let fold (folder: 'State option -> 'T option -> 'State option) (state: 'State option) (matrix: SparseMatrix<'T>) =


Надо не забыть потом написать что-то в духе fold_values, чтобы пустые куски не обходить. Вообще, надо потом подумать, как такое глобально дизайнить. Может и для map должна быть пара, которая только значеня обходит.

gsvgit · 2026-02-05T12:43:37Z

QuadTree.Tests/Tests.Matrix.fs

+    Assert.Equal(expected, actual)
+
+[<Fact>]
+let ``4x4 lower triangle`` () =


А для нечётного размера? А для каких-нибудь "пустых" блоков со странной структурой?

gsvgit · 2026-02-05T12:43:53Z

QuadTree.Tests/Tests.Matrix.fs


    Assert.Equal(expected.storage.data, actual.storage.data)
+
+[<Fact>]


На транспонирование ещё надо бы.

gsvgit · 2026-02-05T12:44:39Z

QuadTree/TriangleCount.fs

+// Assume non-oriented graph adjacency matrix
+// Some () -> edge, None -> no edge
+// Computes triangle count
+let triangle_count (graph: Matrix.SparseMatrix<unit>) =


Почему unit? Вам же всё равно, что там на рёбрах.

gsvgit · 2026-02-05T12:45:37Z

QuadTree/TriangleCount.fs

+
+    let op_mult o1 o2 =
+        match o1, o2 with
+        | Some(), Some() -> Some 1UL


Во. Вам точно всё равно, что на рёбрах.

gsvgit · 2026-02-06T06:16:14Z

QuadTree/Matrix.fs

    else
        (Error.InconsistentSizeOfArguments(matrix1, matrix2)) |> Result.Failure

+let foldAssociative (folder: 'T option -> 'T option -> 'T option) (state: 'T option) (matrix: SparseMatrix<'T>) =


Получился в тоге reduce. Давайте так и назовём. А над тем, как fold нормальный делать, я ещё подумаю.

Danil-Zaripov added 2 commits February 4, 2026 18:02

Add mxm

b166100

Add expanding and shrinking mxm tests

0d141d5

Danil-Zaripov requested a review from gsvgit February 4, 2026 15:53

gsvgit reviewed Feb 4, 2026

View reviewed changes

gsvgit reviewed Feb 5, 2026

View reviewed changes

Danil-Zaripov added 6 commits February 5, 2026 10:35

Move helper functions into inner block

7b443f3

Make mxm more generic

c431f1c

Add leaf instantiation helper functions

3a039e6

else if -> elif

cd11138

Use storageSize inside mxm

ecca8cb

nw1xnw2 -> nw1_x_nw2

cad5fac

Danil-Zaripov requested a review from gsvgit February 5, 2026 08:06

Add triangle counting

0ca94b4

Danil-Zaripov changed the title ~~Matrix multiplication~~ Matrix multiplication and triangle counting Feb 5, 2026

gsvgit reviewed Feb 5, 2026

View reviewed changes

Danil-Zaripov added 2 commits February 5, 2026 16:03

Make triangle_count generic

167cb29

Make more use of helpers

419195b

Danil-Zaripov force-pushed the mxm branch from 5d70f17 to 419195b Compare February 5, 2026 13:52

Danil-Zaripov added 4 commits February 5, 2026 16:57

Add upper triangle part in count test

544900e

Remove unhelpful comment

7f9c74d

Make fold use recursion

9d5c478

Assume folder is associative

8b8d36a

Danil-Zaripov requested a review from gsvgit February 5, 2026 15:06

gsvgit approved these changes Feb 6, 2026

View reviewed changes

gsvgit merged commit ecd5ba0 into Lamagraph:main Feb 6, 2026
2 checks passed


		Assert.Equal(expected.storage.data, actual.storage.data)

		[<Fact>]

Conversation

Danil-Zaripov commented Feb 4, 2026

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

2 participants