Support 3D input for MKL-DNN softmax operator

[TaoLv]

Description

1. to support 3D softmax layers in GluonNLP BERT (need Update MKL-DNN submodule to v0.19 #14783 for better performance)
2. fix in-place softmax
3. remove ctx.is_train check so the cpp test for softmax can work
4. enhance the checks in SupportMKLDNNSoftmax

Checklist

Essentials

Please feel free to remove inapplicable items for your PR.

• The PR title starts with [MXNET-$JIRA_ID], where $JIRA_ID refers to the relevant JIRA issue created (except PRs with tiny changes)
• Changes are complete (i.e. I finished coding on this PR)
• All changes have test coverage:
• Unit tests are added for small changes to verify correctness (e.g. adding a new operator)
• Nightly tests are added for complicated/long-running ones (e.g. changing distributed kvstore)
• Build tests will be added for build configuration changes (e.g. adding a new build option with NCCL)
• Code is well-documented:
• For user-facing API changes, API doc string has been updated.
• For new C++ functions in header files, their functionalities and arguments are documented.
• For new examples, README.md is added to explain the what the example does, the source of the dataset, expected performance on test set and reference to the original paper if applicable
• Check the API doc at http://mxnet-ci-doc.s3-accelerate.dualstack.amazonaws.com/PR-$PR_ID/$BUILD_ID/index.html
• To the my best knowledge, examples are either not affected by this change, or have been fixed to be compatible with this change

Changes

• Feature1, tests, (and when applicable, API doc)
• Feature2, tests, (and when applicable, API doc)

Comments

• If this change is a backward incompatible change, why must this change be made.
• Interesting edge cases to note here

[TaoLv]

@pengzhao-intel @juliusshufan

[pengzhao-intel]

@TaoLv thanks for the PR.

Is there a test for the 1D softmax and could you show the performance of MKL-DNN primitive against original implementation?

[pengzhao-intel]

LGTM :)

[TaoLv]

Tests should be covered by
https://github.com/apache/incubator-mxnet/blob/master/tests/python/unittest/test_operator.py#L4697
and
https://github.com/apache/incubator-mxnet/blob/master/tests/cpp/operator/mkldnn_operator_test.cc#L1288.
I used below code snippet for performance benchmarking:

def test_performance():
    shapes = [(1024,), (96, 512), (96, 128, 128), (96, 256, 256), (1, 8, 1024, 1024)]
    for sh in shapes:
        a = mx.nd.random.uniform(shape=sh)
        # warm up
        b = mx.nd.softmax(a, axis=-1)
        b.wait_to_read()

        tic = time.time()
        for i in range(1000):
            b = mx.nd.softmax(a, axis=-1)
            b.wait_to_read()

        toc = time.time()
        print("softmax %s, take %f ms" % (sh, (toc - tic)/1000*1000.0))

Some performance numbers as following:
mxnet==1.5.0b20190426

softmax (1024,), take 0.103340 ms
softmax (96, 512), take 0.127465 ms
softmax (96, 128, 128), take 1.655400 ms
softmax (96, 256, 256), take 6.369653 ms
softmax (1, 8, 1024, 1024), take 11.450656 ms

This PR with MKL-DNN backend:

softmax (1024,), take 0.062743 ms
softmax (96, 512), take 0.104104 ms
softmax (96, 128, 128), take 0.385350 ms
softmax (96, 256, 256), take 0.463220 ms
softmax (1, 8, 1024, 1024), take 1.704757 ms

[TaoLv]

Pending on MKL-DNN update for better performance~

[TaoLv]

Fallback all softmax operations when axis != last dimension because they are not optimized in MKL-DNN.

[pengzhao-intel]

@TaoLv I have merged the MKL-DNN 0.19 and please rebase the code and see if everything is OK :)

[pengzhao-intel]

@TaoLv please rebase and retrigger again the CI issue is fixed now.

[pengzhao-intel]

Merging now :) Thanks for your contribution.