call dpnp.kron and dpnp.dot for special cases in dpnp_matmul (#1815)

* reorganize the structure of dpnp_matmul function

reorganizing for easier future code development
calling dpnp.kron and dpnp.dot for special cases in dpnp.matmul

* address comments

* improve coverage

* fix pre-commit

---------

Co-authored-by: Anton <100830759+antonwolfy@users.noreply.github.com>

Author: vtavana

dpnp/dpnp_iface_histograms.py

@@ -173,6 +173,7 @@ def _get_bin_edges(a, bins, range, usm_type):
         # numpy's gh-10322 means that type resolution rules are dependent on
         # array shapes. To avoid this causing problems, we pick a type now and
         # stick with it throughout.
+        # pylint: disable=possibly-used-before-assignment
         bin_type = dpnp.result_type(first_edge, last_edge, a)
         if dpnp.issubdtype(bin_type, dpnp.integer):
             bin_type = dpnp.result_type(
@@ -332,6 +333,7 @@ def histogram(a, bins=10, range=None, density=None, weights=None):
 
     if density:
         db = dpnp.diff(bin_edges).astype(dpnp.default_float_type())
+        # pylint: disable=possibly-used-before-assignment
         return n / db / n.sum(), bin_edges
 
     return n, bin_edges