diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/README.md
new file mode 100644
index 000000000000..1bc8c67dde33
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/README.md
@@ -0,0 +1,146 @@
+
+
+# dgemm
+
+> Perform the matrix-matrix operation `C = alpha*op(A)*op(B) + beta*C`.
+
+
+
+
+
+
+
+## Usage
+
+```javascript
+var dgemm = require( '@stdlib/blas/base/ndarray/dgemm' );
+```
+
+#### dgemm( arrays )
+
+Performs the matrix-matrix operation `C = alpha*op(A)*op(B) + beta*C`, where `op(X)` is either `op(X) = X` or `op(X) = X^T`, `alpha` and `beta` are scalars, `A`, `B`, and `C` are matrices, with `op(A)` an `M` by `K` matrix, `op(B)` a `K` by `N` matrix, and `C` an `M` by `N` matrix.
+
+```javascript
+var Float64Matrix = require( '@stdlib/ndarray/matrix/float64' );
+var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
+var resolveEnum = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
+
+var A = new Float64Matrix( [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ] );
+var B = new Float64Matrix( [ [ 1.0, 1.0 ], [ 0.0, 1.0 ] ] );
+var C = new Float64Matrix( [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ] );
+
+var transA = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+});
+var transB = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+});
+var alpha = scalar2ndarray( 1.0, {
+ 'dtype': 'float64'
+});
+var beta = scalar2ndarray( 1.0, {
+ 'dtype': 'float64'
+});
+
+var out = dgemm( [ A, B, C, transA, transB, alpha, beta ] );
+// returns [ [ 2.0, 5.0 ], [ 6.0, 11.0 ] ]
+
+var bool = ( out === C );
+// returns true
+```
+
+The function has the following parameters:
+
+- **arrays**: array-like object containing the following ndarrays:
+
+ - a two-dimensional input ndarray corresponding to `A`.
+ - a two-dimensional input ndarray corresponding to `B`.
+ - a two-dimensional input/output ndarray corresponding to `C`.
+ - a zero-dimensional ndarray specifying whether `A` should be transposed, conjugate-transposed, or not transposed.
+ - a zero-dimensional ndarray specifying whether `B` should be transposed, conjugate-transposed, or not transposed.
+ - a zero-dimensional ndarray containing a scalar constant corresponding to `alpha`.
+ - a zero-dimensional ndarray containing a scalar constant corresponding to `beta`.
+
+
+
+
+
+
+
+
+
+
+
+## Examples
+
+
+
+```javascript
+var discreteUniform = require( '@stdlib/random/discrete-uniform' );
+var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
+var resolveEnum = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
+var ndarray2array = require( '@stdlib/ndarray/to-array' );
+var dgemm = require( '@stdlib/blas/base/ndarray/dgemm' );
+
+var opts = {
+ 'dtype': 'float64'
+};
+
+var A = discreteUniform( [ 3, 4 ], 0, 10, opts );
+var B = discreteUniform( [ 4, 2 ], 0, 10, opts );
+var C = discreteUniform( [ 3, 2 ], 0, 10, opts );
+
+var transA = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+});
+var transB = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+});
+var alpha = scalar2ndarray( 1.0, opts );
+var beta = scalar2ndarray( 1.0, opts );
+
+var out = dgemm( [ A, B, C, transA, transB, alpha, beta ] );
+console.log( ndarray2array( out ) );
+```
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/benchmark/benchmark.js
new file mode 100644
index 000000000000..2ad514bf7f6e
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/benchmark/benchmark.js
@@ -0,0 +1,125 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2026 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+* http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var bench = require( '@stdlib/bench' );
+var uniform = require( '@stdlib/random/uniform' );
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var pow = require( '@stdlib/math/base/special/pow' );
+var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
+var resolveEnum = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
+var format = require( '@stdlib/string/format' );
+var pkg = require( './../package.json' ).name;
+var dgemm = require( './../lib' );
+
+
+// VARIABLES //
+
+var options = {
+ 'dtype': 'float64'
+};
+
+
+// FUNCTIONS //
+
+/**
+* Creates a benchmark function.
+*
+* @private
+* @param {PositiveInteger} len - array length
+* @returns {Function} benchmark function
+*/
+function createBenchmark( len ) {
+ var transA;
+ var transB;
+ var alpha;
+ var beta;
+ var A;
+ var B;
+ var C;
+
+ A = uniform( [ len, len ], -100.0, 100.0, options );
+ B = uniform( [ len, len ], -100.0, 100.0, options );
+ C = uniform( [ len, len ], -100.0, 100.0, options );
+
+ alpha = scalar2ndarray( 1.0, options );
+ beta = scalar2ndarray( 1.0, options );
+ transA = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+ });
+ transB = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+ });
+
+ return benchmark;
+
+ /**
+ * Benchmark function.
+ *
+ * @private
+ * @param {Benchmark} b - benchmark instance
+ */
+ function benchmark( b ) {
+ var z;
+ var i;
+
+ b.tic();
+ for ( i = 0; i < b.iterations; i++ ) {
+ z = dgemm( [ A, B, C, transA, transB, alpha, beta ] );
+ if ( typeof z !== 'object' ) {
+ b.fail( 'should return an ndarray' );
+ }
+ }
+ b.toc();
+ if ( isnan( z.get( i%len, i%len ) ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ b.pass( 'benchmark finished' );
+ b.end();
+ }
+}
+
+
+// MAIN //
+
+/**
+* Main execution sequence.
+*
+* @private
+*/
+function main() {
+ var len;
+ var min;
+ var max;
+ var f;
+ var i;
+
+ min = 1; // 10^min
+ max = 3; // 10^max
+
+ for ( i = min; i <= max; i++ ) {
+ len = pow( 10, i );
+ f = createBenchmark( len );
+ bench( format( '%s:len=%d', pkg, len ), f );
+ }
+}
+
+main();
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/repl.txt b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/repl.txt
new file mode 100644
index 000000000000..62f1588fec66
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/repl.txt
@@ -0,0 +1,54 @@
+
+{{alias}}( arrays )
+ Performs the matrix-matrix operation `C = alpha*op(A)*op(B) + beta*C`,
+ where `op(X)` is either `op(X) = X` or `op(X) = X^T`, `alpha` and `beta`
+ are scalars, `A`, `B`, and `C` are matrices, with `op(A)` an `M` by `K`
+ matrix, `op(B)` a `K` by `N` matrix, and `C` an `M` by `N` matrix.
+
+ Parameters
+ ----------
+ arrays: ArrayLikeObject
+ Array-like object containing the following ndarrays:
+
+ - a two-dimensional input ndarray corresponding to `A`.
+ - a two-dimensional input ndarray corresponding to `B`.
+ - a two-dimensional input/output ndarray corresponding to `C`.
+ - a zero-dimensional ndarray specifying whether `A` should be
+ transposed, conjugate-transposed, or not transposed.
+ - a zero-dimensional ndarray specifying whether `B` should be
+ transposed, conjugate-transposed, or not transposed.
+ - a zero-dimensional ndarray containing a scalar constant corresponding
+ to `alpha`.
+ - a zero-dimensional ndarray containing a scalar constant corresponding
+ to `beta`.
+
+ Returns
+ -------
+ out: ndarray
+ Output ndarray.
+
+ Examples
+ --------
+ > var abuf = [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ];
+ > var A = new {{alias:@stdlib/ndarray/matrix/float64}}( abuf );
+
+ > var bbuf = [ [ 1.0, 1.0 ], [ 0.0, 1.0 ] ];
+ > var B = new {{alias:@stdlib/ndarray/matrix/float64}}( bbuf );
+
+ > var cbuf = [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ];
+ > var C = new {{alias:@stdlib/ndarray/matrix/float64}}( cbuf );
+
+ > var opts = { 'dtype': 'float64' };
+ > var alpha = {{alias:@stdlib/ndarray/from-scalar}}( 1.0, opts );
+ > var beta = {{alias:@stdlib/ndarray/from-scalar}}( 1.0, opts );
+
+ > var transA = {{alias:@stdlib/ndarray/from-scalar}}( 'no-transpose' );
+ > var transB = {{alias:@stdlib/ndarray/from-scalar}}( 'no-transpose' );
+
+ > {{alias}}( [ A, B, C, transA, transB, alpha, beta ] );
+ > C
+ [ [ 2.0, 5.0 ], [ 6.0, 11.0 ] ]
+
+ See Also
+ --------
+
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/types/index.d.ts
new file mode 100644
index 000000000000..57d7cd9dbb62
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/types/index.d.ts
@@ -0,0 +1,76 @@
+/*
+* @license Apache-2.0
+*
+* Copyright (c) 2026 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+* http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+// TypeScript Version: 4.1
+
+///
+
+import { float64ndarray, ndarray } from '@stdlib/types/ndarray';
+
+/**
+* Performs the matrix-matrix operation `C = alpha*op(A)*op(B) + beta*C`, where `op(X)` is either `op(X) = X` or `op(X) = X^T`, `alpha` and `beta` are scalars, `A`, `B`, and `C` are matrices, with `op(A)` an `M` by `K` matrix, `op(B)` a `K` by `N` matrix, and `C` an `M` by `N` matrix.
+*
+* ## Notes
+*
+* - The function expects the following ndarrays:
+*
+* - a two-dimensional input ndarray corresponding to `A`.
+* - a two-dimensional input ndarray corresponding to `B`.
+* - a two-dimensional input/output ndarray corresponding to `C`.
+* - a zero-dimensional ndarray specifying whether `A` should be transposed, conjugate-transposed, or not transposed.
+* - a zero-dimensional ndarray specifying whether `B` should be transposed, conjugate-transposed, or not transposed.
+* - a zero-dimensional ndarray containing a scalar constant corresponding to `alpha`.
+* - a zero-dimensional ndarray containing a scalar constant corresponding to `beta`.
+*
+* @param arrays - array-like object containing ndarrays
+* @returns output ndarray
+*
+* @example
+* var Float64Matrix = require( '@stdlib/ndarray/matrix/float64' );
+* var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
+* var resolveEnum = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
+*
+* var A = new Float64Matrix( [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ] );
+* var B = new Float64Matrix( [ [ 1.0, 1.0 ], [ 0.0, 1.0 ] ] );
+* var C = new Float64Matrix( [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ] );
+*
+* var transA = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+* 'dtype': 'int32'
+* });
+* var transB = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+* 'dtype': 'int32'
+* });
+* var alpha = scalar2ndarray( 1.0, {
+* 'dtype': 'float64'
+* });
+* var beta = scalar2ndarray( 1.0, {
+* 'dtype': 'float64'
+* });
+*
+* var z = dgemm( [ A, B, C, transA, transB, alpha, beta ] );
+* // returns [ [ 2.0, 5.0 ], [ 6.0, 11.0 ] ]
+*
+* var bool = ( z === C );
+* // returns true
+*/
+declare function dgemm( arrays: [ float64ndarray, float64ndarray, float64ndarray, ndarray, ndarray, float64ndarray, float64ndarray ] ): float64ndarray;
+
+
+// EXPORTS //
+
+export = dgemm;
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/types/test.ts b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/types/test.ts
new file mode 100644
index 000000000000..7e5afa0d84bd
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/docs/types/test.ts
@@ -0,0 +1,93 @@
+/*
+* @license Apache-2.0
+*
+* Copyright (c) 2026 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+* http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+/* eslint-disable space-in-parens */
+
+import zeros = require( '@stdlib/ndarray/zeros' );
+import dgemm = require( '@stdlib/blas/base/ndarray/dgemm' );
+
+
+// TESTS //
+
+// The function returns an ndarray...
+{
+ const A = zeros( [ 2, 2 ], {
+ 'dtype': 'float64'
+ });
+ const B = zeros( [ 2, 2 ], {
+ 'dtype': 'float64'
+ });
+ const C = zeros( [ 2, 2 ], {
+ 'dtype': 'float64'
+ });
+ const transA = zeros( [], {
+ 'dtype': 'int32'
+ });
+ const transB = zeros( [], {
+ 'dtype': 'int32'
+ });
+ const alpha = zeros( [], {
+ 'dtype': 'float64'
+ });
+ const beta = zeros( [], {
+ 'dtype': 'float64'
+ });
+
+ dgemm( [ A, B, C, transA, transB, alpha, beta ] ); // $ExpectType float64ndarray
+}
+
+// The compiler throws an error if the function is provided a first argument which is not an array of ndarrays...
+{
+ dgemm( '10' ); // $ExpectError
+ dgemm( 10 ); // $ExpectError
+ dgemm( true ); // $ExpectError
+ dgemm( false ); // $ExpectError
+ dgemm( null ); // $ExpectError
+ dgemm( undefined ); // $ExpectError
+ dgemm( [] ); // $ExpectError
+ dgemm( {} ); // $ExpectError
+ dgemm( ( x: number ): number => x ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided an unsupported number of arguments...
+{
+ const A = zeros( [ 2, 2 ], {
+ 'dtype': 'float64'
+ });
+ const B = zeros( [ 2, 2 ], {
+ 'dtype': 'float64'
+ });
+ const C = zeros( [ 2, 2 ], {
+ 'dtype': 'float64'
+ });
+ const transA = zeros( [], {
+ 'dtype': 'int32'
+ });
+ const transB = zeros( [], {
+ 'dtype': 'int32'
+ });
+ const alpha = zeros( [], {
+ 'dtype': 'float64'
+ });
+ const beta = zeros( [], {
+ 'dtype': 'float64'
+ });
+
+ dgemm(); // $ExpectError
+ dgemm( [ A, B, C, transA, transB, alpha, beta ], {} ); // $ExpectError
+}
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/examples/index.js b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/examples/index.js
new file mode 100644
index 000000000000..21377c39b161
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/examples/index.js
@@ -0,0 +1,45 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2026 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+* http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+var discreteUniform = require( '@stdlib/random/discrete-uniform' );
+var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
+var resolveEnum = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
+var ndarray2array = require( '@stdlib/ndarray/to-array' );
+var dgemm = require( './../lib' );
+
+var opts = {
+ 'dtype': 'float64'
+};
+
+var A = discreteUniform( [ 3, 4 ], 0, 10, opts );
+var B = discreteUniform( [ 4, 2 ], 0, 10, opts );
+var C = discreteUniform( [ 3, 2 ], 0, 10, opts );
+
+var transA = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+});
+var transB = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+ 'dtype': 'int32'
+});
+var alpha = scalar2ndarray( 1.0, opts );
+var beta = scalar2ndarray( 1.0, opts );
+
+var out = dgemm( [ A, B, C, transA, transB, alpha, beta ] );
+console.log( ndarray2array( out ) );
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/lib/index.js b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/lib/index.js
new file mode 100644
index 000000000000..ca4ead3d7541
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/lib/index.js
@@ -0,0 +1,63 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2026 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+* http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+/**
+* BLAS level 3 routine to perform the matrix-matrix operation `C = alpha*op(A)*op(B) + beta*C`.
+*
+* @module @stdlib/blas/base/ndarray/dgemm
+*
+* @example
+* var Float64Matrix = require( '@stdlib/ndarray/matrix/float64' );
+* var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
+* var resolveEnum = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
+* var dgemm = require( '@stdlib/blas/base/ndarray/dgemm' );
+*
+* var A = new Float64Matrix( [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ] );
+* var B = new Float64Matrix( [ [ 1.0, 1.0 ], [ 0.0, 1.0 ] ] );
+* var C = new Float64Matrix( [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ] );
+*
+* var transA = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+* 'dtype': 'int32'
+* });
+* var transB = scalar2ndarray( resolveEnum( 'no-transpose' ), {
+* 'dtype': 'int32'
+* });
+* var alpha = scalar2ndarray( 1.0, {
+* 'dtype': 'float64'
+* });
+* var beta = scalar2ndarray( 1.0, {
+* 'dtype': 'float64'
+* });
+*
+* var out = dgemm( [ A, B, C, transA, transB, alpha, beta ] );
+* // returns [ [ 2.0, 5.0 ], [ 6.0, 11.0 ] ]
+*
+* var bool = ( out === C );
+* // returns true
+*/
+
+// MODULES //
+
+var main = require( './main.js' );
+
+
+// EXPORTS //
+
+module.exports = main;
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/lib/main.js b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/lib/main.js
new file mode 100644
index 000000000000..00f94f4b5c88
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dgemm/lib/main.js
@@ -0,0 +1,134 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2026 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+* http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var getShape = require( '@stdlib/ndarray/base/shape' );
+var getStrides = require( '@stdlib/ndarray/base/strides' );
+var getOffset = require( '@stdlib/ndarray/base/offset' );
+var getData = require( '@stdlib/ndarray/base/data-buffer' );
+var ndarraylike2scalar = require( '@stdlib/ndarray/base/ndarraylike2scalar' );
+var resolveEnum = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
+var strided = require( '@stdlib/blas/base/dgemm' ).ndarray;
+
+
+// VARIABLES //
+
+var NO_TRANSPOSE = resolveEnum( 'no-transpose' );
+
+
+// MAIN //
+
+/**
+* Performs the matrix-matrix operation `C = alpha*op(A)*op(B) + beta*C`, where `op(X)` is either `op(X) = X` or `op(X) = X^T`, `alpha` and `beta` are scalars, `A`, `B`, and `C` are matrices, with `op(A)` an `M` by `K` matrix, `op(B)` a `K` by `N` matrix, and `C` an `M` by `N` matrix.
+*
+* ## Notes
+*
+* - The function expects the following ndarrays:
+*
+* - a two-dimensional input ndarray corresponding to `A`.
+* - a two-dimensional input ndarray corresponding to `B`.
+* - a two-dimensional input/output ndarray corresponding to `C`.
+* - a zero-dimensional ndarray specifying whether `A` should be transposed, conjugate-transposed, or not transposed.
+* - a zero-dimensional ndarray specifying whether `B` should be transposed, conjugate-transposed, or not transposed.
+* - a zero-dimensional ndarray containing a scalar constant corresponding to `alpha`.
+* - a zero-dimensional ndarray containing a scalar constant corresponding to `beta`.
+*
+* @param {ArrayLikeObject