This is an archive of the discontinued LLVM Phabricator instance.

Add support for numpy arrays to memref conversions.
ClosedPublic

Authored by pashu123 on Apr 7 2021, 4:21 PM.

Details

Summary

This offers the ability to pass numpy arrays to the corresponding
memref argument.

Diff Detail

Event Timeline

pashu123 created this revision.Apr 7 2021, 4:21 PM
pashu123 requested review of this revision.Apr 7 2021, 4:21 PM
pashu123 updated this revision to Diff 335948.Apr 7 2021, 4:29 PM

Minor Changes.

pashu123 updated this revision to Diff 335950.Apr 7 2021, 4:40 PM

Python formating.

mehdi_amini added inline comments.Apr 7 2021, 4:42 PM
mlir/lib/Bindings/Python/mlir/runtime/memref/np_to_memref.py
27 ↗(On Diff #335948)
mlir/test/Bindings/Python/execution_engine.py
144

This will not work with a memref that has strides I believe. You'll need to use something like np.lib.stride_tricks.as_strided.
I think this should be a utility in the runtime as well memref_to_numpy_view or something like that.

(If you can add a test with an array and a view with strides that could show this)

183

Debug left-over? (same below)

pashu123 updated this revision to Diff 335952.Apr 7 2021, 4:45 PM

Clean left over.

Harbormaster completed remote builds in B97615: Diff 335950.
bondhugula requested changes to this revision.Apr 7 2021, 8:35 PM

I'd really like to see this support! Thanks for implementing this.

mlir/lib/Bindings/Python/mlir/runtime/memref/np_to_memref.py
4 ↗(On Diff #335952)

Please add a line on what this file is about.

41 ↗(On Diff #335952)

Doc comment.

mlir/test/Bindings/Python/execution_engine.py
136

Terminate all comments with a full stop.

144

Is there a way to restrict this patch to identity layout map? This would enable and unlock *many* things even if it just handles the default/row major identity layout. Can we check and bail out / assert on non-default strides?

191

Unrelated to this PR: If the wrong number of arguments are provided, will ExecutionEngine catch it? I think from a developer standpoint, this is quite useful since folks might sometimes lower the MLIR module to a form where there is a mismatch in the number of function arguments.

This revision now requires changes to proceed.Apr 7 2021, 8:35 PM
bondhugula added inline comments.Apr 7 2021, 8:36 PM
mlir/lib/Bindings/Python/mlir/runtime/memref/np_to_memref.py
47 ↗(On Diff #335952)

get_ranked_memref_descriptor?

Here and elsewhere, can we name these things descriptively?

mehdi_amini added inline comments.Apr 7 2021, 8:50 PM
mlir/test/Bindings/Python/execution_engine.py
144

The stride is really easy to add though, maybe easier than bailing out (and adding the tests for all the bail out).

But if you prefer to add these tests, I can implement the strides in a follow up.

191

This is a layer of consistency that is important, but we'll have to build above this, these are all quite low-level APIs: the execution engine does not have any knowledge of the APIs present in the Module.

The only comment I have (not covered elsewhere) is a weak preference to unnest by one directory level (i.e. remove the memref directory and put this in the parent). All of the exported APIs already have memref in the name and fewer imports/less repeating for these kind of things is better, imo. Also fine to see how things evolve and rework the namespace once/if there is more there.

pashu123 updated this revision to Diff 336556.Apr 9 2021, 1:29 PM

Added support for strided memrefs.

mlir/lib/Bindings/Python/mlir/runtime/memref/np_to_memref.py
68 ↗(On Diff #336556)

Note that numpy uses byte quantities to express strides.

MLIR OTOH uses the torch abstraction which specifies strides in terms of elements; the ConversionToLLVM takes care of generating the right addresses (which will also requirer furrther hooking the data layout better).

Bottom line, I experimented with memref<...xf32> by hacking the following on top of this:

+    # x.strides = nparray.ctypes.strides
+   strides_ctype_t = ctypes.c_longlong * nparray.ndim
+    x.strides = strides_ctype_t(*[t // 4 for t in nparray.strides])

It seemed to help a bit but I still saw issues that I have not yet debugged.
I am not sure whether the ctype lifetime is reasonable the code I wrote.

mlir/test/Bindings/Python/execution_engine.py
144

Also note the discrepancy between np.strides and memref / torch strides I signal above.

pashu123 updated this revision to Diff 336557.Apr 9 2021, 1:41 PM

Remove the memref directory and put files in the parent.

mlir/lib/Bindings/Python/mlir/runtime/memref/np_to_memref.py
110 ↗(On Diff #336556)

nice, you also want to pointwise-multiply by the element size in bytes.

pashu123 updated this revision to Diff 336686.Apr 11 2021, 1:14 PM

Cleaning up.

pashu123 marked 3 inline comments as done.Apr 11 2021, 1:15 PM
pashu123 marked an inline comment as done.Apr 11 2021, 1:23 PM
nicolasvasilache requested changes to this revision.Apr 12 2021, 12:08 AM
nicolasvasilache added inline comments.
mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
75

The 4 is only valid for 4-byte type like f32, you need to get the size in bytes of the type (here and other places you updated).

This revision now requires changes to proceed.Apr 12 2021, 12:08 AM
pashu123 added inline comments.Apr 12 2021, 12:15 AM
mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
75

I see, thanks for pointing this out.

pashu123 updated this revision to Diff 336761.Apr 12 2021, 12:24 AM

Fix itemsize.

pashu123 updated this revision to Diff 336763.Apr 12 2021, 12:26 AM

Formatting.

pashu123 marked an inline comment as done.Apr 12 2021, 12:27 AM
bondhugula requested changes to this revision.Apr 12 2021, 5:21 AM
bondhugula added inline comments.
mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
75

This also means that test cases were missing - you don't seem to have f64 test cases.

mlir/test/Bindings/Python/execution_engine.py
159

You are missing test cases on float64 and so this isn't exercising all of the code well.

This revision now requires changes to proceed.Apr 12 2021, 5:21 AM
mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
60

This is not the element size but the displacement from the "base pointer address" at which the relevant data lives.
You probably picked the value here: https://numpy.org/doc/stable/reference/generated/numpy.ndarray.html ?
This shows that similary to strides, numpy represents this "offset in bytes" whereas MLIR is in number of elements.

I verified that setting this to 0 in my local experiment (here and below), worked in simple cases.

In the general case you need to translate between MLIR and NP offsets similarly as you do for strides.

pashu123 added inline comments.Apr 12 2021, 11:30 AM
mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
60

I tried setting this to multiple different values and it worked. If offset is represented as a number of elements in MLIR, then the offset should be 1. Similarly, it should be 1 * nparray.itemsize while converting back to NumPy from memrefs world.

mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
60

I am unclear what "it worked" means in the absence of more context on what you tried?

Please note that lowering will do the right thing in the static case here: https://github.com/llvm/llvm-project/blob/main/mlir/lib/Conversion/StandardToLLVM/StandardToLLVM.cpp#L3399.
The problem I mention is only exhibited by dynamic offsets which I do not see your tests implementing.

For instance, this commit will fail with an offset different than 0 (at runtime, in this case I hardcoded 0 for expediency while awaiting for this PR to land) : https://github.com/google/iree-llvm-sandbox/commit/87015b29c5f7cbc445bd85e1ce4a5d7597e80361.

The underlying reason is that after tiling, the offset becomes dependence on the loop IV (i.e. the memref has a static ? in its type) and then the dynamic branch of the StandardToLLVM.cpp will kick in.
Address computation will compute an offset that is 1 + actual value and will segfault at runtime.

pashu123 added inline comments.Apr 13 2021, 4:45 AM
mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
60

Thanks for pointing this out. I didn't try on Dynamic memrefs.

pashu123 updated this revision to Diff 337108.Apr 13 2021, 5:21 AM

Redirecting prints to stderr.

mehdi_amini accepted this revision.Apr 13 2021, 3:50 PM

LGTM, seems like a good basis to iterate on!

@bondhugula @nicolasvasilache : can you confirm the recent changes address your previous comments?

@mehdi_amini not yet, the author should at the very least set the 2 offsets to 0 to as described in the commit message of https://github.com/google/iree-llvm-sandbox/commit/87015b29c5f7cbc445bd85e1ce4a5d7597e80361, so that it passes without crashing.

I'm fine with iterating in-tree once that is fixed.

mehdi_amini added inline comments.Apr 13 2021, 4:04 PM
mlir/lib/Bindings/Python/mlir/runtime/np_to_memref.py
60

I don't quite indeed this x.offset = ctypes.c_longlong(4)? Shouldn't it always be zero? I suspect Numpy is baking the offset into the pointer we get with nparray.ctypes.data_as.

@mehdi_amini indeed, in https://numpy.org/doc/stable/reference/generated/numpy.ndarray.html we can see:

>>> np.ndarray((2,), buffer=np.array([1,2,3]),
...            offset=np.int_().itemsize,
...            dtype=int) # offset = 1*itemsize, i.e. skip first element
array([2, 3])

as well as the default = 0 for offset:

class numpy.ndarray(shape, dtype=float, buffer=None, offset=0, strides=None, order=None)[source]¶

0 is the right value to use for now.
When we do morer fancy numpy / linalg level subviews etc it will be a different story; but for now anything else than 0 will crash non-trivial cases.

@mehdi_amini indeed, in https://numpy.org/doc/stable/reference/generated/numpy.ndarray.html we can see:

>>> np.ndarray((2,), buffer=np.array([1,2,3]),
...            offset=np.int_().itemsize,
...            dtype=int) # offset = 1*itemsize, i.e. skip first element
array([2, 3])

Right, but the ctypes interface does not expose the offset, so I wonder if in such cases the offset that you provide on construction will not just be computed in the "data" field internally (I could check the source code for ndarray...)

pashu123 updated this revision to Diff 337336.Apr 13 2021, 10:35 PM

Changing the default value of offset to zero.

pashu123 added a comment.EditedApr 13 2021, 10:36 PM

@mehdi_amini not yet, the author should at the very least set the 2 offsets to 0 to as described in the commit message of https://github.com/google/iree-llvm-sandbox/commit/87015b29c5f7cbc445bd85e1ce4a5d7597e80361, so that it passes without crashing.

I'm fine with iterating in-tree once that is fixed.

Done. I have changed it to zero. I had pushed the wrong changes.

@bondhugula @nicolasvasilache : can you confirm the recent changes address your previous comments?

Is this revision also meant to handle dynamic memrefs? If yes, a test case is missing. If not, we should check and assert (and add a TODO)?

@bondhugula @nicolasvasilache : can you confirm the recent changes address your previous comments?

Is this revision also meant to handle dynamic memrefs? If yes, a test case is missing. If not, we should check and assert (and add a TODO)?

It supports a subset of cases (e.g. anything that takes subviews and crosses the NP -> MLIR or MLIR -> NP boundary should be considered broken atm).
+1 on adding a simple mixed static / dynamic 2-D add and compare against NP's.

pashu123 updated this revision to Diff 337860.Apr 15 2021, 12:00 PM

Adding test case of element wise addition of dynamic and static memrefs.

Great, thanks for your contribution, let's land this ! :)

This revision was not accepted when it landed; it landed in state Needs Review.Apr 15 2021, 4:41 PM
This revision was automatically updated to reflect the committed changes.