Skip to content

Spgemm #72

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 35 commits into from
Apr 10, 2023
Merged

Spgemm #72

merged 35 commits into from
Apr 10, 2023

Conversation

IgorErin
Copy link
Contributor

@IgorErin IgorErin commented Apr 7, 2023

Proposed Changes

  • spGeMM (naive ESC)
  • Scatter and Gather refactor
  • ClArray.choose, *.choose2

Types of changes

What types of changes does your code introduce to GraphBLAS-sharp?
Put an x in the boxes that apply

  • Bugfix (non-breaking change which fixes an issue)
  • New feature (non-breaking change which adds functionality)
  • Breaking change (fix or feature that would cause existing functionality to not work as expected)

Checklist

Put an x in the boxes that apply. You can also fill these out after creating the PR. If you're unsure about any of them, don't hesitate to ask. We're here to help! This is simply a reminder of what we are going to look for before merging your code.

  • Build and tests pass locally
  • I have added tests that prove my fix is effective or that my feature works (if appropriate)
  • I have added necessary documentation (if appropriate)

Further comments

If this is a relatively large or complex change, kick off the discussion by explaining why you chose the solution you did and what alternatives you considered, etc...


let getUniqueBitmapLastOccurrence clContext =
getUniqueBitmapGeneral
<| <@ fun (gid: int) (length: int) (inputArray: ClArray<'a>) ->
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Code duplication with this

@IgorErin IgorErin requested a review from gsvgit April 8, 2023 13:23
else
isUniqueBitmap.[i] <- 1 @>
if gid < inputLength then
let isUnique = (%predicate) gid inputLength inputArray // brahma error
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is it reported?

// add and filter
|> Array.map Some
|> Array.fold add None
|> function
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

if Utils.isFloat64Available context.ClDevice then
createTest ((=) 0.0) Expand.getSegmentPointers

createTest ((=) 0f) Expand.getSegmentPointers
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is it planned to use createTest with other testFun-ctions? If no, what is the reason to pass it as parameter?

let (total: ClCell<_>) = scan q clArray zero

let actual = Array.zeroCreate<'a> clArray.Length
let actual = clArray.ToHostAndFree q
let actualSum = total.ToHostAndFree(q)
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should q be in brackets or not?

@@ -24,7 +36,7 @@ let checkResult isEqual actualKeys actualValues keys values reduceOp =
"Values must the same"
|> Utils.compareArrays isEqual actualValues expectedValues

let makeTest isEqual reduce reduceOp (arrayAndKeys: (int * 'a) []) =
let private makeTest isEqual reduce reduceOp (arrayAndKeys: (int * 'a) []) =
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Why `makeTest' sometimes private but sometimes -- not


fun (queue: MailboxProcessor<_>) (matrix1: ClMatrix<'a>) (matrix2: ClMatrix<'b>) (mask: ClMatrix<_>) ->
match matrix1, matrix2, mask with
| ClMatrix.CSR m1, ClMatrix.CSC m2, ClMatrix.COO mask -> runCSRnCSC queue m1 m2 mask |> ClMatrix.COO
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

How can one to multiply two CSR matrix and get CSR matrix? Can one use CSR matrix as a mask?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Or it is an old version form Artiom?

fun (processor: MailboxProcessor<_>) lengths (segmentsPointers: Indices) (leftMatrix: ClMatrix.CSR<'a>) (rightMatrix: ClMatrix.CSR<'b>) ->

// Compute A positions
let APositions = zeroCreate processor DeviceOnly lengths
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

leftMatrixPositions? Same for other names below.

match sum with
| Some value ->
let result =
((%reduceOp) value values.[currentPosition]) // brahma error
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is it reported?

@IgorErin IgorErin requested a review from gsvgit April 9, 2023 08:09
@gsvgit gsvgit merged commit e0485c5 into SparseLinearAlgebra:dev Apr 10, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants