DLPrimitives Blog
Development Blog
Posts in category ‘Internals’.
Attempt to integrate with OneDNN
Intel's OneDNN is great project that provides cudnn/inference/training like tools for Intel's GPU.
Also it is called OneDNN... it should be called IntelDNN since it supports only Intel gpus and cpus.
Bottom line I tried to add OneDNN based convolutions for Intel GPU just to discover that my simple GEMM based convolution works better. Why? Apparently Intel's implementation seems to be optimized for Channel Last format only.
https://github.com/oneapi-src/oneDNN/issues/1194
A simple convolution with 3x3 kernel with 64 input and output channels with image dimension of 56 on Intel HD 530 with 400 GFlops capacity gives:
- 295.6 GFlops for OneDNN's channels last format
- 144.7 GFlops for dlprimitive's channel first format
- 33.4(!) GFlops for OneDNN's channels first format.
The problem is that channels first is the most common format used by pytorch, mxnet, caffe and many other tools (including dlprimitives)
Ok... I'll check it later when one of two happens:
- They fix channel first performance
- I'll support channel last format internally
Pytorch Updates
In order to improve the progress I started validating all pretrained torchvision models one by one. I found several features I needed to implement but what is more important I found several critical bugs I could fix.
https://pytorch.org/vision/stable/models.html#classification
At this point following networks are validated against CPU version in both forward and backward propagation:
alexnetresnet18resnet50vgg16densenet161googlenetsqueezenet1_0inception_v3(fwd only - backward fails on cuda/cpu)shufflenet_v2_x1_0mobilenet_v2mobilenet_v3_largemobilenet_v3_small(fwd only - same failure on bwd on cuda)resnext50_32x4dwide_resnet50_2mnasnet1_0efficientnet_b0efficientnet_b4regnet_y_400mf
To be continued...
Update Nov 17, 2021: I implemneted ceil rounding pooling mode, thus googlenet and squeezenet1_0 now pass validation
Pointwise Broadcast Reduce
Lots of deep learning operations can be implemented as simple element-by-element operations over different tensors with numpy broadcasting and reduction afterwards. For example:
Adding Bias [C] to [B,C,H,W] image is can be seen in numpy as:
x + bias.reshape((C,1,1))
Gradient of bias can be calculated as:
np.sum(dy,dims=(0,2,3))
That is simple reduction operations. Calculation of mean and variance in batch normalisation requires calculation of x and x*x over all dims but C.
Observing this I implemented a broadcast/reduce templates API to simplify development. http://dlprimitives.org/docs/pointwise_8hpp_source.html
The idea is following:
- You provide input tensors and scalar parameters
- You define the operation need to performed on each operand
- You provide reduction operation
The OpenCL kernel code is auto-generated for you. For example calculations of x and x*x sums over all dims but channels would look like:
auto op = dlprim::core::PointwiseOperationBroadcastReduce::create(
ctx,
{X.specs()},{Xsum.specs(),X2sum.specs()},
0,dlprim::float_data,
"y0=x0; y1=x0*x0;", // operations
"reduce_y0 = 0; reduce_y1 = 0", // reduce init
"reduce_y0 += y0; reduce_y1 += y1"
);
op->enqueue({X},{Xsum,X2sum},s,{},{1,1},{0,0},q);
So - 1st output is just x - sum and second is x*x - sum. So if you provide X in shape of [B,C,H,W] and Xsum, X2sum in shape [C,1,1] that is broadcast-able to X you'll get the sums you need without writing custom reduction code of manually writing kernels.
This vastly simplified writing multiple operators especially ones that are expected to support numpy style broadcasting in pytorch.
Hello Pytorch OpenCL
TL;DR: I managed to run an inference of alexnet using OpenCL/DLPrimitives based pytorch backend!
Details
I started from this tutorial to implement out of source backend for pytorch. It wasn't that simple and I had to do small changes in original pytorch source code but finally something is working:
Now I implemented only handful of ops and mostly for forward computations: github backend code. However, I managed to do forward computations and get correct result on pretrained alexnet.
$ python validate_network.py --model alexnet --device cuda *.ppm
cat.ppm,281,tabby
dog.ppm,207,golden retriever
parrot.ppm,87,African grey
$ python validate_network.py --model alexnet --device opencl:1 *.ppm
Accessing device #1:GeForce GTX 960 on NVIDIA CUDA
cat.ppm,281,tabby
dog.ppm,207,golden retriever
parrot.ppm,87,African grey
Performance for this tiny task isn't not brilliant, but not horrible either, GTX 960, alexnet batch size of 16 images 224x224:
- Pytorch Cuda/CUDNN: 15.317 ms - updated 2021-10-10
- Pytorch OpenCL/DLPrimitives: 22.932 ms - updated 2021-10-10
- DLPrim - microframework: 22.401 ms
- Caffe/CuDNN: 16.1812 ms
- Caffe/OpenCL: 41.072 ms
- Caffe/OpenCL+DLPrimitives: 28.618 ms
- Keras/CuDNN: 23.341 ms
- Keras/PlaidML: 44.041 ms
Now, one of the issues that I currently have is synchronous execution that gives significant penalty for every operation. I need to understand an asynchronous execution and memory management stuff before I continue. The penalty for NVidia OpenCL backend isn't horrible but it is devastating for AMD OpenCL driver. Need to dive in.
Keep updated.
Edit, Oct 10, 2021
I found the way to implement asynchronous execution + initial GPU memory caching. That allowed to be bring the performance of pytorch OpenCL to same level as vanilla dlprimitives. This also solved the performance issues I had with AMD GPU.
Additionally I found that I didn't take in an account host-to-device transfer in pytorch benchmarks - so original CUDA run time for pytorch increased.
Why do we need OpenCL based deep learning?
Why do we need OpenCL based solution for deep learning?
NVidia provides high performance tools like cuDNN and TensorRT that power AI industry running CUDA API. AMD does the same for their compute cards with ROCm/MIOpen using their own CUDA clone called "hip".
Why should we care? We have a working solution why do we need to reimplement something that already exists with something that likely going to have lower performance?
The problem is somewhat deeper than that. I'll talk about 3 points.
- We need open source high performance low level algorithms especially for research purposes.
- We need to use unified GPU API that is standard, open and vendor independent.
- We need a worthy competition that would lead to better and more affordable products.
Open source algorithms are mandatory
A good example is a convolution algorithm. In the beginning of DL boom the most common approach was to run im2col with GEMM combination that gave decent performance thanks to very efficient cublas library.
However it was found that it is much better to merge two GPU kernels - one that converts image to a matrix and another that computes matrix multiplication to a single kernel. However, if you want to implement this technique and check its efficiency you need to have an efficient implementation of core matrix multiplication algorithm in first place. Also matrix multiplication seems like a simple task, in fact it is one of the most complicated alorithms to implement in GPU
If a researcher wants to add a new method of convolution, lets say dilation with variable steps or any other kind of modification it wouldn't be possible to do this efficiently without easily available source code. On the other hand if a researcher wants to improve some performance aspects of the Conv+GEMM algorithm he need to implement what nVidia engineers did in assembly to be able to compare his new implementation to what is done in cuDNN.
I would say a good example of such a case would be Winograd Convolution. There is a widely cited paper that describes the algorithm. In reality, there are virtually 0 details about specific implementation. How do you load/store data to shared memory without bank conflicts? how do you perform computations? How do you store and convert results between shared memory to main one? And so on. I personally attempted to implement it several times and couldn't even reach a GEMM based convolution performance.
Only this paper from 2020 had cleared many low level details an allowed me writing a relatively efficient kernel. It still does not reach the performance of cuDNN implementation but it is way better than GEMM based convolution.
Unfortunately, this is a typical case, if tomorrow somebody wants to implement an additional variant of convolution or any other algorithm and prove its efficiency, it is first needed to overcome the performance of unknown highly optimised code to prove that you did something efficient.
Thus having open source low level algorithms is very critical for research and the field itself. It is very clear to me why NVidia do not want to release the source code - it contradicts their interest as a monopoly in the field.
However the community that works in the field is hurt by this policy. We need faster and more efficient software to solve real world problems, without having the source code of core components the progress in the field may be much slower.
Open and Standardised API
I recently discovered that I can't run a code that was compiled with cuda10 and cudnn7 on RTX 30 series of GPU. I didn't want any of the new GPU features - just wanted to run the program. I couldn't. I need to rebuild the software from the scratch with a new API.
What do you think would happen in the gaming industry if NVidia would told the game developers who released their games a year ago to rebuild them from the scratch to be able to run on RTX 3060?
It isn't only user unfriendly, it is also sometimes impossible because the team isn't actively developing a product any more and moved to a new one.
Of course this does not happen in the gaming world, there are standard open or proprietary APIs: OpenGL, Vulkan and Direct3D. To get a new card you need to install a driver. The rest is the same. You can run today a classic old game like Jane's IAF: Israeli Air Force on Windows 10 on modern hardware.
There is a standard GPU compute API: OpenCL. It works very well, it works on virtually every modern GPU from 1st release date: NVidia, AMD, Intel. It even runs on your Android phone's GPU.
Kernels that are written for OpenCL or CUDA can be mechanically converted between the platforms as is. Virtually all concepts of OpenCL and CUDA are identical, code is similar.
This is highly critical for software maintainability and long-term support.
Of course, there are many cases you can't have same kernel running equally efficient on AMD, NVidia, Intel and Mali GPUs. It is true for GEMM and Convolution kernels, but it true to much lesser degree for vast majority of other kernels that do basic but important stuff. Kernels for activation, normalization, element-wise operators and many many others can run same code on different platforms with only minute tweaks.
So as long as you have optimised kernels for critical parts, rest can be shared across platforms. It can be compared to writing a code in C or C++ for all platforms like Intel, ARM and MIPS and only small computationally intensive parts are implemented using platform specific assembly on intrinsic code.
Finally, if you look into gaming industry dealing with GPUs on daily basis you would clearly understand how the requirement of rebuilding your entire code base upon arrival of new GPU generation is ridiculous. But today it is the sad reality.
Open standards are always better than best proprietary ones. Even if you loose something now you are winning in long-term run.
Monopoly Issue
It isn't a secret that NVidia controls virtually all deep-learning and GPU-compute market. It goes to that level that AMD creates "hip" - virtually implemented cuda API with fast replacement of s/cuda/hip/ just to be able to run NVidia's code.
As long as vendor specific API is used, you will remain under the control of this vendor.
It isn't the case in Gaming. AMD, Intel and NVidia always try to challenge each other to this degree or other while we as customers enjoy better and more affordable products, enjoy new features like Ray Tracing etc.
It isn't the case for deep-learning. NVidia can charge premium for their "professional" grade devices that provide virtually same silicone. Why? Because they really can. If I was NVidia's CEO I would do the same.
Just to make it clear, currently NVidia is superior DL platform over AMD and Intel, mostly due to absolute ignorance of software by AMD and lack of any kind of powerful GPUs from Intel.
But even if you rightfully use NVidia, relying on open API and open source libraries will allow you to both challenge NVidia that does not take competition seriously and have much lesser headache when NVidia releases a new GPU.
And what is more important you are empowered to choose whatever product you or your customer want or has and not be limited to a single vendor.
Additional Points
Why not SYCL?
The answer is actually very simple. There are several points that make OpenCL superior:
- The fact that the kernel source code is separate from C++ code actually makes it much more portable. Mixing
C++and GPU code the way it is done in CUDA it is nice to start with but brings lots of issues in long term where you need to compile for numerous platforms and combinations. - Many SYCL Implementation target one platform only instead of being cross platform.
- There are still no serious working open source SYCL implementation, while OpenCL is working today on virtually every device
Why wouldn't AMD or Intel create an alternative?
Good question. And this is should be addressed to them.
AMD created hip and MiOpen that aren't even compatible with their own RDNA hardware, started dropping support of older GCN (like rx580) and limited to Linux only. Their investment into deep learning is minimalist to none.
Intel played along and actually created OpenCL implementations. However their product oneDNN is optimised for Intel and wrapped cuDNN for NVidia.
Bottom line, each vendor cares about his own interests, probably they decided that this is a lost battle.
However the good thing that both AMD and Intel published almost everything as open source so it makes easier in future to use their work in critical training paths.