From 01d6f065bd1601d0a4a1365aac2a91c06defb22d Mon Sep 17 00:00:00 2001
From: Shabareesh Shetty <shabareeshshetty@Shabareeshs-MacBook-Air.local>
Date: Fri, 25 Jul 2025 13:48:49 +0530
Subject: [PATCH 1/4] feat: add blas/base/gsyr

---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: passed
  - task: lint_package_json
    status: passed
  - task: lint_repl_help
    status: passed
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: passed
  - task: lint_javascript_tests
    status: passed
  - task: lint_javascript_benchmarks
    status: passed
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: na
  - task: lint_c_examples
    status: na
  - task: lint_c_benchmarks
    status: na
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: passed
  - task: lint_typescript_tests
    status: passed
  - task: lint_license_headers
    status: passed
---
---
 .../@stdlib/blas/base/gsyr/README.md          |  196 +++
 .../blas/base/gsyr/benchmark/benchmark.js     |   93 ++
 .../base/gsyr/benchmark/benchmark.ndarray.js  |   93 ++
 .../@stdlib/blas/base/gsyr/docs/repl.txt      |  123 ++
 .../blas/base/gsyr/docs/types/index.d.ts      |  114 ++
 .../@stdlib/blas/base/gsyr/docs/types/test.ts |  342 ++++
 .../@stdlib/blas/base/gsyr/examples/index.js  |   42 +
 .../@stdlib/blas/base/gsyr/lib/accessors.js   |  141 ++
 .../@stdlib/blas/base/gsyr/lib/base.js        |  125 ++
 .../@stdlib/blas/base/gsyr/lib/index.js       |   61 +
 .../@stdlib/blas/base/gsyr/lib/main.js        |   97 ++
 .../@stdlib/blas/base/gsyr/lib/ndarray.js     |   83 +
 .../@stdlib/blas/base/gsyr/package.json       |   65 +
 .../column_major_complex_access_pattern.json  |   20 +
 .../gsyr/test/fixtures/column_major_l.json    |   20 +
 .../gsyr/test/fixtures/column_major_oa.json   |   20 +
 .../gsyr/test/fixtures/column_major_ox.json   |   20 +
 .../test/fixtures/column_major_sa1_sa2.json   |   20 +
 .../test/fixtures/column_major_sa1_sa2n.json  |   20 +
 .../test/fixtures/column_major_sa1n_sa2.json  |   20 +
 .../test/fixtures/column_major_sa1n_sa2n.json |   20 +
 .../gsyr/test/fixtures/column_major_u.json    |   20 +
 .../gsyr/test/fixtures/column_major_xn.json   |   20 +
 .../gsyr/test/fixtures/column_major_xp.json   |   20 +
 .../row_major_complex_access_pattern.json     |   20 +
 .../base/gsyr/test/fixtures/row_major_l.json  |   20 +
 .../base/gsyr/test/fixtures/row_major_oa.json |   20 +
 .../base/gsyr/test/fixtures/row_major_ox.json |   20 +
 .../gsyr/test/fixtures/row_major_sa1_sa2.json |   20 +
 .../test/fixtures/row_major_sa1_sa2n.json     |   20 +
 .../test/fixtures/row_major_sa1n_sa2.json     |   20 +
 .../test/fixtures/row_major_sa1n_sa2n.json    |   20 +
 .../base/gsyr/test/fixtures/row_major_u.json  |   20 +
 .../base/gsyr/test/fixtures/row_major_xn.json |   20 +
 .../base/gsyr/test/fixtures/row_major_xp.json |   20 +
 .../@stdlib/blas/base/gsyr/test/test.js       |   38 +
 .../@stdlib/blas/base/gsyr/test/test.main.js  |  798 ++++++++++
 .../blas/base/gsyr/test/test.ndarray.js       | 1417 +++++++++++++++++
 38 files changed, 4268 insertions(+)
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/README.md
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.ndarray.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/docs/repl.txt
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/docs/types/index.d.ts
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/lib/accessors.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/lib/index.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/package.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_complex_access_pattern.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_l.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_oa.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_ox.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2n.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2n.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_u.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xn.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xp.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_complex_access_pattern.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_l.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_oa.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_ox.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2n.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2n.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_u.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xn.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xp.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/test.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js
 create mode 100644 lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js

diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/README.md b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
new file mode 100644
index 000000000000..6061f6ab252e
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
@@ -0,0 +1,196 @@
+<!--
+
+@license Apache-2.0
+
+Copyright (c) 2024 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.
+
+-->
+
+# gsyr
+
+> Perform the symmetric rank 1 operation `A = α*x*x^T + A`.
+
+<section class="usage">
+
+## Usage
+
+```javascript
+var gsyr = require( '@stdlib/blas/base/gsyr' );
+```
+
+#### gsyr( order, uplo, N, α, x, sx, A, LDA )
+
+Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
+var x = [ 1.0, 2.0, 3.0 ];
+
+gsyr( 'row-major', 'upper', 3, 1.0, x, 1, A, 3 );
+// A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+```
+
+The function has the following parameters:
+
+-   **order**: storage layout.
+-   **uplo**: specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced.
+-   **N**: number of elements along each dimension of `A`.
+-   **α**: scalar constant.
+-   **x**: input array.
+-   **sx**: stride length for `x`.
+-   **A**: input matrix stored in linear memory.
+-   **LDA**: stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`).
+
+The stride parameters determine how elements in the input arrays are accessed at runtime. For example, to iterate over the elements of `x` in reverse order,
+
+```javascript
+var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
+var x = [ 3.0, 2.0, 1.0 ];
+
+gsyr( 'row-major', 'upper', 3, 1.0, x, -1, A, 3 );
+// A => <Float64Array>[ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+```
+
+Note that indexing is relative to the first index. To introduce an offset, use [`typed array`][mdn-typed-array] views.
+
+<!-- eslint-disable stdlib/capitalized-comments -->
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+// Initial arrays...
+var x0 = new Float64Array( [ 0.0, 3.0, 2.0, 1.0 ] );
+var A = new Float64Array( [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ] );
+
+// Create offset views...
+var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element
+
+gsyr( 'row-major', 'upper', 3, 1.0, x1, -1, A, 3 );
+// A => <Float64Array>[ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+```
+
+#### gsyr.ndarray( uplo, N, α, x, sx, ox, A, sa1, sa2, oa )
+
+Performs the symmetric rank 1 operation `A = α*x*x^T + A`, using alternative indexing semantics and where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
+var x = [ 1.0, 2.0, 3.0 ];
+
+gsyr.ndarray( 'upper', 3, 1.0, x, 1, 0, A, 3, 1, 0 );
+// A => <Float64Array>[ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+```
+
+The function has the following additional parameters:
+
+-   **ox**: starting index for `x`.
+-   **sa1**: stride of the first dimension of `A`.
+-   **sa2**: stride of the second dimension of `A`.
+-   **oa**: starting index for `A`.
+
+While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example,
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
+var x = [ 1.0, 2.0, 3.0, 4.0, 5.0 ];
+
+gsyr.ndarray( 'upper', 3, 1.0, x, -2, 4, A, 3, 1, 0 );
+// A => [ 26.0, 17.0, 8.0, 2.0, 10.0, 5.0, 3.0, 2.0, 2.0 ]
+```
+
+</section>
+
+<!-- /.usage -->
+
+<section class="notes">
+
+## Notes
+
+-   `gsyr()` corresponds to the [BLAS][blas] level 2 function [`dsyr`][dsyr] with the exception that this implementation works with any array type, not just Float64Arrays. Depending on the environment, the typed versions ([`dsyr`][@stdlib/blas/base/dsyr], [`ssyr`][@stdlib/blas/base/ssyr], etc.) are likely to be significantly more performant.
+-   Both functions support array-like objects having getter and setter accessors for array element access (e.g., [`@stdlib/array/base/accessor`][@stdlib/array/base/accessor]).
+
+</section>
+
+<!-- /.notes -->
+
+<section class="examples">
+
+## Examples
+
+<!-- eslint no-undef: "error" -->
+
+```javascript
+var discreteUniform = require( '@stdlib/random/array/discrete-uniform' );
+var ones = require( '@stdlib/array/ones' );
+var gsyr = require( '@stdlib/blas/base/gsyr' );
+
+var opts = {
+    'dtype': 'gsyr'
+};
+
+var N = 3;
+
+// Create N-by-N symmetric matrices:
+var A1 = ones( N*N, opts.dtype );
+var A2 = ones( N*N, opts.dtype );
+
+// Create a random vector:
+var x = discreteUniform( N, -10.0, 10.0, opts );
+
+gsyr( 'row-major', 'upper', 3, 1.0, x, 1, A1, 3 );
+console.log( A1 );
+
+gsyr.ndarray( 'upper', 3, 1.0, x, 1, 0, A2, 3, 1, 0 );
+console.log( A2 );
+```
+
+</section>
+
+<!-- /.examples -->
+
+<!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
+
+<section class="related">
+
+</section>
+
+<!-- /.related -->
+
+<!-- Section for all links. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="links">
+
+[blas]: http://www.netlib.org/blas
+
+[dsyr]: https://www.netlib.org/lapack/explore-html/dc/d82/group__her_ga07f0e3f8592107877f12a554a41c7413.html#ga07f0e3f8592107877f12a554a41c7413
+
+[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray
+
+[@stdlib/blas/base/dsyr]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/base/dsyr
+
+[@stdlib/blas/base/ssyr]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/base/ssyr
+
+[@stdlib/array/base/accessor]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/base/accessor
+
+
+</section>
+
+<!-- /.links -->
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.js
new file mode 100644
index 000000000000..1bb6cc804ff9
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.js
@@ -0,0 +1,93 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 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/array/uniform' );
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var pow = require( '@stdlib/math/base/special/pow' );
+var floor = require( '@stdlib/math/base/special/floor' );
+var pkg = require( './../package.json' ).name;
+var gsyr = require( './../lib' );
+
+
+// VARIABLES //
+
+var options = {
+	'dtype': 'float64'
+};
+
+
+// FUNCTIONS //
+
+/**
+* Create a benchmark function.
+*
+* @private
+* @param {PositiveInteger} N - array dimension size
+* @returns {Function} benchmark function
+*/
+function createBenchmark( N ) {
+	var x = uniform( N, -10.0, 10.0, options );
+	var A = uniform( N*N, -10.0, 10.0, options );
+	return benchmark;
+
+	function benchmark( b ) {
+		var z;
+		var i;
+
+		b.tic();
+		for ( i = 0; i < b.iterations; i++ ) {
+			z = gsyr( 'row-major', 'upper', N, 1.0, x, 1, A, N );
+			if ( isnan( z[ i%z.length ] ) ) {
+				b.fail( 'should not return NaN' );
+			}
+		}
+		b.toc();
+		if ( isnan( z[ i%z.length ] ) ) {
+			b.fail( 'should not return NaN' );
+		}
+		b.pass( 'benchmark finished' );
+		b.end();
+	}
+}
+
+
+// MAIN //
+
+function main() {
+	var min;
+	var max;
+	var N;
+	var f;
+	var i;
+
+	min = 1; // 10^min
+	max = 6; // 10^max
+
+	for ( i = min; i <= max; i++ ) {
+		N = floor( pow( pow( 10, i ), 1.0/2.0 ) );
+		f = createBenchmark( N );
+		bench( pkg+':size='+(N*N), f );
+	}
+}
+
+main();
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.ndarray.js
new file mode 100644
index 000000000000..483f2f7bf35c
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/benchmark/benchmark.ndarray.js
@@ -0,0 +1,93 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 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/array/uniform' );
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var pow = require( '@stdlib/math/base/special/pow' );
+var floor = require( '@stdlib/math/base/special/floor' );
+var pkg = require( './../package.json' ).name;
+var gsyr = require( './../lib' ).ndarray;
+
+
+// VARIABLES //
+
+var options = {
+	'dtype': 'float64'
+};
+
+
+// FUNCTIONS //
+
+/**
+* Create a benchmark function.
+*
+* @private
+* @param {PositiveInteger} N - array dimension size
+* @returns {Function} benchmark function
+*/
+function createBenchmark( N ) {
+	var x = uniform( N, -10.0, 10.0, options );
+	var A = uniform( N*N, -10.0, 10.0, options );
+	return benchmark;
+
+	function benchmark( b ) {
+		var z;
+		var i;
+
+		b.tic();
+		for ( i = 0; i < b.iterations; i++ ) {
+			z = gsyr( 'upper', N, 1.0, x, 1, 0, A, N, 1, 0 );
+			if ( isnan( z[ i%z.length ] ) ) {
+				b.fail( 'should not return NaN' );
+			}
+		}
+		b.toc();
+		if ( isnan( z[ i%z.length ] ) ) {
+			b.fail( 'should not return NaN' );
+		}
+		b.pass( 'benchmark finished' );
+		b.end();
+	}
+}
+
+
+// MAIN //
+
+function main() {
+	var min;
+	var max;
+	var N;
+	var f;
+	var i;
+
+	min = 1; // 10^min
+	max = 6; // 10^max
+
+	for ( i = min; i <= max; i++ ) {
+		N = floor( pow( pow( 10, i ), 1.0/2.0 ) );
+		f = createBenchmark( N );
+		bench( pkg+':ndarray:size='+(N*N), f );
+	}
+}
+
+main();
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/docs/repl.txt b/lib/node_modules/@stdlib/blas/base/gsyr/docs/repl.txt
new file mode 100644
index 000000000000..4f78e2d51796
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/docs/repl.txt
@@ -0,0 +1,123 @@
+
+{{alias}}( order, uplo, N, α, x, sx, A, lda )
+    Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a
+    scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric
+    matrix.
+
+    Indexing is relative to the first index. To introduce an offset, use typed
+    array views.
+
+    If `N` is equal to `0` or `α` is equal to `0`, the function returns `A`
+    unchanged.
+
+    Parameters
+    ----------
+    order: string
+        Row-major (C-style) or column-major (Fortran-style) order.
+
+    uplo: string
+        Specifies whether to reference the upper or lower triangular part of
+        `A`. Must be either 'upper' or 'lower'.
+
+    N: integer
+        Number of elements along each dimension of `A`.
+
+    α: number
+        Scalar constant.
+
+    x: Array<number>|TypedArray
+        Input vector.
+
+    sx: integer
+        Stride length for `x`.
+
+    A: Array<number>|TypedArray
+        Input matrix.
+
+    lda: integer
+        Stride of the first dimension of `A` (a.k.a., leading dimension of the
+        matrix `A`).
+
+    Returns
+    -------
+    A: Array<number>|TypedArray
+        Input matrix.
+
+    Examples
+    --------
+    // Standard usage:
+    > var x = [ 1.0, 1.0 ];
+    > var A = [ 1.0, 2.0, 2.0, 1.0 ];
+    > {{alias}}( 'row-major', 'upper', 2, 1.0, x, 1, A, 2 )
+    [ 2.0, 3.0, 2.0, 2.0 ]
+
+    // Advanced indexing:
+    > x = [ 1.0, 1.0 ];
+    > A = [ 1.0, 2.0, 2.0, 1.0 ];
+    > {{alias}}( 'row-major', 'upper', 2, 1.0, x, -1, A, 2 )
+    [ 2.0, 3.0, 2.0, 2.0 ]
+
+    // Using typed array views:
+    > var x0 = new {{alias:@stdlib/array/float64}}( [ 0.0, 1.0, 1.0 ] );
+    > A = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0, 2.0, 1.0 ] );
+    > var x1 = new {{alias:@stdlib/array/float64}}( x0.buffer, x0.BYTES_PER_ELEMENT*1 );
+    > {{alias}}( 'row-major', 'upper', 2, 1.0, x, 1, A, 2 )
+    <Float64Array>[ 2.0, 3.0, 2.0, 2.0 ]
+
+
+{{alias}}.ndarray( uplo, N, α, x, sx, ox, A, sa1, sa2, oa )
+    Performs the symmetric rank 1 operation `A = α*x*x^T + A`, using alternative
+    indexing semantics and where `α` is a scalar, `x` is an `N` element vector,
+    and `A` is an `N` by `N` symmetric matrix.
+
+    While typed array views mandate a view offset based on the underlying
+    buffer, the offset parameters support indexing semantics based on starting
+    indices.
+
+    Parameters
+    ----------
+    uplo: string
+        Specifies whether to reference the upper or lower triangular part of
+        `A`. Must be either 'upper' or 'lower'.
+
+    N: integer
+        Number of elements along each dimension of `A`.
+
+    α: number
+        Scalar constant.
+
+    x: Array<number>|TypedArray
+        Input vector.
+
+    sx: integer
+        Stride length for `x`.
+
+    ox: integer
+        Starting index for `x`.
+
+    A: Array<number>|TypedArray
+        Input matrix.
+
+    sa1: integer
+        Stride of the first dimension of `A`.
+
+    sa2: integer
+        Stride of the second dimension of `A`.
+
+    oa: integer
+        Starting index for `A`.
+
+    Returns
+    -------
+    A: Array<number>|TypedArray
+        Input matrix.
+
+    Examples
+    --------
+    > var x = [ 1.0, 1.0 ];
+    > var A = [ 1.0, 2.0, 2.0, 1.0 ];
+    > {{alias}}.ndarray( 'upper', 2, 1.0, x, 1, 0, A, 2, 1, 0 )
+    [ 2.0, 3.0, 2.0, 2.0 ]
+
+    See Also
+    --------
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/index.d.ts
new file mode 100644
index 000000000000..4ec0f6d61a15
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/index.d.ts
@@ -0,0 +1,114 @@
+/*
+* @license Apache-2.0
+*
+* Copyright (c) 2025 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
+
+/// <reference types="@stdlib/types"/>
+
+import { NumericArray, Collection, AccessorArrayLike } from '@stdlib/types/array';
+import { Layout, MatrixTriangle } from '@stdlib/types/blas';
+
+/**
+* Input array.
+*/
+type InputArray = NumericArray | Collection<number> | AccessorArrayLike<number>;
+
+/**
+* Interface describing `gsyr`.
+*/
+interface Routine {
+	/**
+	* Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+	*
+	* @param order - storage layout
+	* @param uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced
+	* @param N - number of elements along each dimension in the matrix `A`
+	* @param alpha - scalar constant
+	* @param x - input vector
+	* @param strideX - stride length for `x`
+	* @param A - input matrix
+	* @param LDA - stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`)
+	* @returns `A`
+	*
+	* @example
+	* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+	* var x = [ 1.0, 2.0, 3.0 ];
+	*
+	* gsyr( 'row-major', 'upper', 3, 1.0, x, 1, A, 3 );
+	* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+	*/
+	<T extends InputArray = InputArray>( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, x: InputArray, strideX: number, A: T, LDA: number ): T;
+
+	/**
+	* Performs the symmetric rank 1 operation `A = α*x*x^T + A` using alternative indexing semantics and where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+	*
+	* @param uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced
+	* @param N - number of elements along each dimension in the matrix `A`
+	* @param alpha - scalar constant
+	* @param x - input vector
+	* @param strideX - stride length for `x`
+	* @param offsetX - starting index for `x`
+	* @param A - input matrix
+	* @param strideA1 - stride of the first dimension of `A`
+	* @param strideA2 - stride of the second dimension of `A`
+	* @param offsetA - starting index for `A`
+	* @returns `A`
+	*
+	* @example
+	* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+	* var x = [ 1.0, 2.0, 3.0 ];
+	*
+	* gsyr.ndarray( 'upper', 3, 1.0, x, 1, 0, A, 3, 1, 0 );
+	* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+	*/
+	ndarray<T extends InputArray = InputArray>( uplo: MatrixTriangle, N: number, alpha: number, x: InputArray, strideX: number, offsetX: number, A: T, strideA1: number, strideA2: number, offsetA: number ): T;
+}
+
+/**
+* Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+*
+* @param order - storage layout
+* @param uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced
+* @param N - number of elements along each dimension in the matrix `A`
+* @param alpha - scalar constant
+* @param x - input vector
+* @param strideX - stride length for `x`
+* @param A - input matrix
+* @param LDA - stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`)
+* @returns `A`
+*
+* @example
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr( 'row-major', 'upper', 3, 1.0, x, 1, A, 3 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*
+* @example
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr.ndarray( 'upper', 3, 1.0, x, 1, 0, A, 3, 1, 0 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*/
+declare var gsyr: Routine;
+
+
+// EXPORTS //
+
+export = gsyr;
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts b/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts
new file mode 100644
index 000000000000..290aea1ab7f2
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts
@@ -0,0 +1,342 @@
+/*
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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.
+*/
+
+import gsyr = require( './index' );
+
+
+// TESTS //
+
+// The function returns an array...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectType Float64Array
+}
+
+// The compiler throws an error if the function is provided a first argument which is not a string...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr( 10, 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( true, 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( false, 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( null, 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( undefined, 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( [], 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( {}, 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( ( x: number ): number => x, 'upper', 10, 1.0, x, 1, A, 10 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a second argument which is not a string...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr( 'row-major', 10, 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', true, 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', false, 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', null, 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', undefined, 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', [ '1' ], 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', {}, 10, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', ( x: number ): number => x, 10, 1.0, x, 1, A, 10 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a third argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr( 'row-major', 'upper', '10', 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', true, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', false, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', null, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', undefined, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', [], 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', {}, 1.0, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', ( x: number ): number => x, 1.0, x, 1, A, 10 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fourth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr( 'row-major', 'upper', 10, '10', x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, true, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, false, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, null, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, undefined, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, [], x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, {}, x, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, ( x: number ): number => x, x, 1, A, 10 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fifth argument which is not a numeric array...
+{
+	const A = new Float64Array( 20 );
+
+	gsyr( 'row-major', 'upper', 10, 1.0, 10, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, '10', 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, true, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, false, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, null, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, undefined, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, [ '1' ], 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, {}, 1, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, ( x: number ): number => x, 1, A, 10 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a sixth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr( 'row-major', 'upper', 10, 1.0, x, '10', A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, true, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, false, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, null, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, undefined, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, [], A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, {}, A, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, ( x: number ): number => x, A, 10 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a seventh argument which is not a numeric array...
+{
+	const x = new Float64Array( 10 );
+
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, 10, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, '10', 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, true, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, false, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, null, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, undefined, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, [ '1' ], 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, {}, 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, ( x: number ): number => x, 10 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided an eighth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, '10' ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, true ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, false ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, null ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, undefined ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, [] ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, {} ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, ( x: number ): number => x ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided an unsupported number of arguments...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr(); // $ExpectError
+	gsyr( 'row-major' ); // $ExpectError
+	gsyr( 'row-major', 'upper' ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1 ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A ); // $ExpectError
+	gsyr( 'row-major', 'upper', 10, 1.0, x, 1, A, 10, 1 ); // $ExpectError
+}
+
+// Attached to main export is an `ndarray` method which returns an array...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectType Float64Array
+}
+
+// The compiler throws an error if the function is provided a first argument which is not a string...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 10, 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( true, 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( false, 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( null, 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( undefined, 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( [ '1' ], 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( {}, 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( ( x: number ): number => x, 10, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a second argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', '10', 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', true, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', false, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', null, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', undefined, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', [], 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', {}, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', ( x: number ): number => x, 1.0, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a third argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, '10', x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, true, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, false, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, null, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, undefined, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, [], x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, {}, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, ( x: number ): number => x, x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fourth argument which is not a numeric array...
+{
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, 10, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, '10', 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, true, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, false, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, null, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, undefined, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, [ '1' ], 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, {}, 1, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, ( x: number ): number => x, 1, 0, A, 10, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fifth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, x, '10', 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, true, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, false, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, null, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, undefined, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, [], 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, {}, 0, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, ( x: number ): number => x, 0, A, 10, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a sixth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, '10', A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, true, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, false, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, null, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, undefined, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, [], A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, {}, A, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, ( x: number ): number => x, A, 10, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a seventh argument which is not a numeric array...
+{
+	const x = new Float64Array( 10 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, 10, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, '10', 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, true, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, false, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, null, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, undefined, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, [ '1' ], 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, {}, 10, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, ( x: number ): number => x, 10, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided an eighth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, '10', 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, true, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, false, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, null, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, undefined, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, [], 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, {}, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, ( x: number ): number => x, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a ninth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, '10', 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, true, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, false, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, null, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, undefined, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, [], 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, {}, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, ( x: number ): number => x, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a tenth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, '10' ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, true ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, false ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, null ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, undefined ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, [] ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, {} ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, ( x: number ): number => x ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments...
+{
+	const x = new Float64Array( 10 );
+	const A = new Float64Array( 20 );
+
+	gsyr.ndarray(); // $ExpectError
+	gsyr.ndarray( 'upper' ); // $ExpectError
+	gsyr.ndarray( 'upper', 10 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1 ); // $ExpectError
+	gsyr.ndarray( 'upper', 10, 1.0, x, 1, 0, A, 10, 1, 0, 10 ); // $ExpectError
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js b/lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js
new file mode 100644
index 000000000000..6c072a6e9b7e
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js
@@ -0,0 +1,42 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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/array/discrete-uniform' );
+var ones = require( '@stdlib/array/ones' );
+var gsyr = require( './../lib' );
+
+var opts = {
+	'dtype': 'generic'
+};
+
+var N = 3;
+
+// Create N-by-N symmetric matrices:
+var A1 = ones( N*N, opts.dtype );
+var A2 = ones( N*N, opts.dtype );
+
+// Create a random vector:
+var x = discreteUniform( N, -10.0, 10.0, opts );
+
+gsyr( 'row-major', 'upper', 3, 1.0, x, 1, A1, 3 );
+console.log( A1 );
+
+gsyr.ndarray( 'upper', 3, 1.0, x, 1, 0, A2, 3, 1, 0 );
+console.log( A2 );
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/accessors.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/accessors.js
new file mode 100644
index 000000000000..d66849a0be6d
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/accessors.js
@@ -0,0 +1,141 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 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 isRowMajor = require( '@stdlib/ndarray/base/assert/is-row-major' );
+
+
+// MAIN //
+
+/**
+* Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+*
+* @private
+* @param {string} uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced
+* @param {NonNegativeInteger} N - number of columns in the matrix `A`
+* @param {number} alpha - scalar constant
+* @param {Object} x - input vector object
+* @param {Collection} x.data - input vector data
+* @param {Array<Function>} x.accessors - array element accessors
+* @param {integer} strideX - stride length for `x`
+* @param {NonNegativeInteger} offsetX - starting index for `x`
+* @param {Object} A - output matrix object
+* @param {Collection} A.data - output matrix data
+* @param {Array<Function>} A.accessors - array element accessors
+* @param {integer} strideA1 - stride of the first dimension of `A`
+* @param {integer} strideA2 - stride of the second dimension of `A`
+* @param {NonNegativeInteger} offsetA - starting index for `A`
+* @returns {Object} output array object
+*
+* @example
+* var toAccessorArray = require( '@stdlib/array/base/to-accessor-array' );
+* var arraylike2object = require( '@stdlib/array/base/arraylike2object' );
+*
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr( 'upper', 3, 1.0, arraylike2object( toAccessorArray( x ) ), 1, 0, arraylike2object( toAccessorArray( A ) ), 3, 1, 0 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*/
+function gsyr( uplo, N, alpha, x, strideX, offsetX, A, strideA1, strideA2, offsetA ) { // eslint-disable-line max-len
+	var getX;
+	var getA;
+	var setA;
+	var xbuf;
+	var Abuf;
+	var isrm;
+	var tmp;
+	var ix0;
+	var ix1;
+	var sa0;
+	var sa1;
+	var i0;
+	var i1;
+	var ia;
+	var ox;
+	var v;
+
+	// Cache references to array data:
+	xbuf = x.data;
+	Abuf = A.data;
+
+	// Cache references to element accessors:
+	getX = x.accessors[ 0 ];
+	getA = A.accessors[ 0 ];
+	setA = A.accessors[ 1 ];
+
+	// Note on variable naming convention: S#, da#, ia#, i# where # corresponds to the loop number, with `0` being the innermost loop...
+
+	isrm = isRowMajor( [ strideA1, strideA2 ] );
+	if ( isrm ) {
+		// For row-major matrices, the last dimension has the fastest changing index...
+		sa0 = strideA2; // stride for innermost loop
+		sa1 = strideA1; // stride for outermost loop
+	} else { // isColMajor
+		// For column-major matrices, the first dimension has the fastest changing index...
+		sa0 = strideA1; // stride for innermost loop
+		sa1 = strideA2; // stride for outermost loop
+	}
+	ox = offsetX;
+	if (
+		( !isrm && uplo === 'upper' ) ||
+		( isrm && uplo === 'lower' )
+	) {
+		ix1 = ox;
+		for ( i1 = 0; i1 < N; i1++ ) {
+			if ( getX( xbuf, ix1 ) !== 0.0 ) {
+				tmp = alpha * getX( xbuf, ix1 );
+				ia = offsetA + (sa1*i1);
+				ix0 = ox;
+				for ( i0 = 0; i0 <= i1; i0++ ) {
+					v = getX( xbuf, ix0 ) * tmp;
+					setA( Abuf, ia, getA( Abuf, ia ) + v );
+					ix0 += strideX;
+					ia += sa0;
+				}
+			}
+			ix1 += strideX;
+		}
+		return A;
+	}
+	// ( isrm && uplo === 'upper' ) || ( !isrm && uplo === 'lower' )
+	ix1 = ox;
+	for ( i1 = 0; i1 < N; i1++ ) {
+		if ( getX( xbuf, ix1 ) !== 0.0 ) {
+			tmp = alpha * getX( xbuf, ix1 );
+			ia = offsetA + (sa1*i1) + (sa0*i1);
+			ix0 = ix1;
+			for ( i0 = i1; i0 < N; i0++ ) {
+				v = getX( xbuf, ix0 ) * tmp;
+				setA( Abuf, ia, getA( Abuf, ia ) + v );
+				ix0 += strideX;
+				ia += sa0;
+			}
+		}
+		ix1 += strideX;
+	}
+	return A;
+}
+
+
+// EXPORTS //
+
+module.exports = gsyr;
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js
new file mode 100644
index 000000000000..e84a1e9f8bdb
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js
@@ -0,0 +1,125 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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 arraylike2object = require( '@stdlib/array/base/arraylike2object' );
+var isRowMajor = require( '@stdlib/ndarray/base/assert/is-row-major' );
+var accessors = require( './accessors.js' );
+
+
+// MAIN //
+
+/**
+* Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+*
+* @private
+* @param {string} uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced
+* @param {NonNegativeInteger} N - number of elements along each dimension of `A`
+* @param {number} alpha - scalar constant
+* @param {NumericArray} x - input vector
+* @param {integer} strideX - stride length for `x`
+* @param {NonNegativeInteger} offsetX - starting index for `x`
+* @param {NumericArray} A - input matrix
+* @param {integer} strideA1 - stride of the first dimension of `A`
+* @param {integer} strideA2 - stride of the second dimension of `A`
+* @param {NonNegativeInteger} offsetA - starting index for `A`
+* @returns {NumericArray} `A`
+*
+* @example
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr( 'upper', 3, 1.0, x, 1, 0, A, 3, 1, 0 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*/
+function gsyr( uplo, N, alpha, x, strideX, offsetX, A, strideA1, strideA2, offsetA ) { // eslint-disable-line max-len
+	var isrm;
+	var tmp;
+	var ix0;
+	var ix1;
+	var sa0;
+	var sa1;
+	var i0;
+	var i1;
+	var ia;
+	var ix;
+	var oa;
+	var ox;
+
+	ox = arraylike2object( x );
+	oa = arraylike2object( A );
+	if ( ox.accessorProtocol || oa.accessorProtocol ) {
+		accessors( uplo, N, alpha, ox, strideX, offsetX, oa, strideA1, strideA2, offsetA );
+		return A;
+	}
+	isrm = isRowMajor( [ strideA1, strideA2 ] );
+	if ( isrm ) {
+		// For row-major matrices, the last dimension has the fastest changing index...
+		sa0 = strideA2; // stride for innermost loop
+		sa1 = strideA1; // stride for outermost loop
+	} else { // isColMajor
+		// For column-major matrices, the first dimension has the fastest changing index...
+		sa0 = strideA1; // stride for innermost loop
+		sa1 = strideA2; // stride for outermost loop
+	}
+	ix = offsetX;
+	if (
+		( !isrm && uplo === 'upper' ) ||
+		( isrm && uplo === 'lower' )
+	) {
+		ix1 = ix;
+		for ( i1 = 0; i1 < N; i1++ ) {
+			if ( x[ ix1 ] !== 0.0 ) {
+				tmp = alpha * x[ ix1 ];
+				ia = offsetA + (sa1*i1);
+				ix0 = ix;
+				for ( i0 = 0; i0 <= i1; i0++ ) {
+					A[ ia ] += x[ ix0 ] * tmp;
+					ix0 += strideX;
+					ia += sa0;
+				}
+			}
+			ix1 += strideX;
+		}
+		return A;
+	}
+	// ( isrm && uplo === 'upper' ) || ( !isrm && uplo === 'lower' )
+	ix1 = ix;
+	for ( i1 = 0; i1 < N; i1++ ) {
+		if ( x[ ix1 ] !== 0.0 ) {
+			tmp = alpha * x[ ix1 ];
+			ia = offsetA + (sa1*i1) + (sa0*i1);
+			ix0 = ix1;
+			for ( i0 = i1; i0 < N; i0++ ) {
+				A[ ia ] += x[ ix0 ] * tmp;
+				ix0 += strideX;
+				ia += sa0;
+			}
+		}
+		ix1 += strideX;
+	}
+	return A;
+}
+
+
+// EXPORTS //
+
+module.exports = gsyr;
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/index.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/index.js
new file mode 100644
index 000000000000..0e95bec0738d
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/index.js
@@ -0,0 +1,61 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 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 2 routine to perform the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+*
+* @module @stdlib/blas/base/gsyr
+*
+* @example
+* var gsyr = require( '@stdlib/blas/base/gsyr' );
+*
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr( 'row-major', 'upper', 3, 1.0, x, 1, A, 3 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*
+* @example
+* var gsyr = require( '@stdlib/blas/base/gsyr' );
+*
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr.ndarray( 'upper', 3, 1.0, x, 1, 0, A, 3, 1, 0 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*/
+
+// MODULES //
+
+var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' );
+var main = require( './main.js' );
+var ndarray = require( './ndarray.js' );
+
+
+// MAIN //
+
+setReadOnly( main, 'ndarray', ndarray );
+
+
+// EXPORTS //
+
+module.exports = main;
+
+// exports: { "ndarray": "main.ndarray" }
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js
new file mode 100644
index 000000000000..b5620466d6ba
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js
@@ -0,0 +1,97 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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 isLayout = require( '@stdlib/blas/base/assert/is-layout' );
+var isMatrixTriangle = require( '@stdlib/blas/base/assert/is-matrix-triangle' );
+var isColumnMajor = require( '@stdlib/ndarray/base/assert/is-column-major-string' );
+var stride2offset = require( '@stdlib/strided/base/stride2offset' );
+var max = require( '@stdlib/math/base/special/fast/max' );
+var format = require( '@stdlib/string/format' );
+var base = require( './base.js' );
+
+
+// MAIN //
+
+/**
+* Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+*
+* @param {string} order - storage layout
+* @param {string} uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced
+* @param {NonNegativeInteger} N - number of elements along each dimension of `A`
+* @param {number} alpha - scalar constant
+* @param {NumericArray} x - input vector
+* @param {integer} strideX - stride length for `x`
+* @param {NumericArray} A - input matrix
+* @param {integer} LDA - stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`)
+* @throws {TypeError} first argument must be a valid order
+* @throws {TypeError} second argument must specify whether to reference the lower or upper triangular matrix
+* @throws {RangeError} third argument must be a nonnegative integer
+* @throws {RangeError} sixth argument must be non-zero
+* @throws {RangeError} eighth argument must be a valid stride
+* @returns {NumericArray} `A`
+*
+* @example
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr( 'row-major', 'upper', 3, 1.0, x, 1, A, 3 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*/
+function gsyr( order, uplo, N, alpha, x, strideX, A, LDA ) {
+	var sa1;
+	var sa2;
+	var ox;
+
+	if ( !isLayout( order ) ) {
+		throw new TypeError( format( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ) );
+	}
+	if ( !isMatrixTriangle( uplo ) ) {
+		throw new TypeError( format( 'invalid argument. Second argument must specify whether to reference the lower or upper triangular matrix. Value: `%s`.', uplo ) );
+	}
+	if ( N < 0 ) {
+		throw new RangeError( format( 'invalid argument. Third argument must be a nonnegative integer. Value: `%d`.', N ) );
+	}
+	if ( strideX === 0 ) {
+		throw new RangeError( format( 'invalid argument. Sixth argument must be non-zero. Value: `%d`.', strideX ) );
+	}
+	if ( LDA < max( 1, N ) ) {
+		throw new RangeError( format( 'invalid argument. Eighth argument must be greater than or equal to max(1,%d). Value: `%d`.', N, LDA ) );
+	}
+	// Check if we can early return...
+	if ( N === 0 || alpha === 0.0 ) {
+		return A;
+	}
+	if ( isColumnMajor( order ) ) {
+		sa1 = 1;
+		sa2 = LDA;
+	} else { // order === 'row-major'
+		sa1 = LDA;
+		sa2 = 1;
+	}
+	ox = stride2offset( N, strideX );
+	return base( uplo, N, alpha, x, strideX, ox, A, sa1, sa2, 0 );
+}
+
+
+// EXPORTS //
+
+module.exports = gsyr;
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js
new file mode 100644
index 000000000000..e42a667f638c
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js
@@ -0,0 +1,83 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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 isMatrixTriangle = require( '@stdlib/blas/base/assert/is-matrix-triangle' );
+var format = require( '@stdlib/string/format' );
+var base = require( './base.js' );
+
+
+// MAIN //
+
+/**
+* Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
+*
+* @param {string} uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced
+* @param {NonNegativeInteger} N - number of elements along each dimension of `A`
+* @param {number} alpha - scalar constant
+* @param {NumericArray} x - input vector
+* @param {integer} strideX - stride length for `x`
+* @param {NonNegativeInteger} offsetX - starting index for `x`
+* @param {NumericArray} A - input matrix
+* @param {integer} strideA1 - stride of the first dimension of `A`
+* @param {integer} strideA2 - stride of the second dimension of `A`
+* @param {NonNegativeInteger} offsetA - starting index for `A`
+* @throws {TypeError} first argument must specify whether to reference the lower or upper triangular matrix
+* @throws {RangeError} second argument must be a nonnegative integer
+* @throws {RangeError} fifth argument must be non-zero
+* @throws {RangeError} eighth argument must be non-zero
+* @throws {RangeError} ninth argument must be non-zero
+* @returns {NumericArray} `A`
+*
+* @example
+* var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ]; // => [ [ 1.0, 2.0, 3.0 ], [ 2.0, 1.0, 2.0 ], [ 3.0, 2.0, 1.0 ] ]
+* var x = [ 1.0, 2.0, 3.0 ];
+*
+* gsyr( 'upper', 3, 1.0, x, 1, 0, A, 3, 1, 0 );
+* // A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+*/
+function gsyr( uplo, N, alpha, x, strideX, offsetX, A, strideA1, strideA2, offsetA ) { // eslint-disable-line max-len
+	if ( !isMatrixTriangle( uplo ) ) {
+		throw new TypeError( format( 'invalid argument. First argument must specify whether to reference the lower or upper triangular matrix. Value: `%s`.', uplo ) );
+	}
+	if ( N < 0 ) {
+		throw new RangeError( format( 'invalid argument. Second argument must be a nonnegative integer. Value: `%d`.', N ) );
+	}
+	if ( strideX === 0 ) {
+		throw new RangeError( format( 'invalid argument. Fifth argument must be non-zero. Value: `%d`.', strideX ) );
+	}
+	if ( strideA1 === 0 ) {
+		throw new RangeError( format( 'invalid argument. Eighth argument must be non-zero. Value: `%d`.', strideA1 ) );
+	}
+	if ( strideA2 === 0 ) {
+		throw new RangeError( format( 'invalid argument. Ninth argument must be non-zero. Value: `%d`.', strideA2 ) );
+	}
+	// Check if we can early return...
+	if ( N === 0 || alpha === 0.0 ) {
+		return A;
+	}
+	return base( uplo, N, alpha, x, strideX, offsetX, A, strideA1, strideA2, offsetA ); // eslint-disable-line max-len
+}
+
+
+// EXPORTS //
+
+module.exports = gsyr;
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/package.json b/lib/node_modules/@stdlib/blas/base/gsyr/package.json
new file mode 100644
index 000000000000..60e639b9bf75
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/package.json
@@ -0,0 +1,65 @@
+{
+  "name": "@stdlib/blas/base/gsyr",
+  "version": "0.0.0",
+  "description": "Perform the symmetric rank 1 operation `A = α*x*x^T + A`.",
+  "license": "Apache-2.0",
+  "author": {
+    "name": "The Stdlib Authors",
+    "url": "https://github.com/stdlib-js/stdlib/graphs/contributors"
+  },
+  "contributors": [
+    {
+      "name": "The Stdlib Authors",
+      "url": "https://github.com/stdlib-js/stdlib/graphs/contributors"
+    }
+  ],
+  "main": "./lib",
+  "directories": {
+    "benchmark": "./benchmark",
+    "doc": "./docs",
+    "example": "./examples",
+    "lib": "./lib",
+    "test": "./test"
+  },
+  "types": "./docs/types",
+  "scripts": {},
+  "homepage": "https://github.com/stdlib-js/stdlib",
+  "repository": {
+    "type": "git",
+    "url": "git://github.com/stdlib-js/stdlib.git"
+  },
+  "bugs": {
+    "url": "https://github.com/stdlib-js/stdlib/issues"
+  },
+  "dependencies": {},
+  "devDependencies": {},
+  "engines": {
+    "node": ">=0.10.0",
+    "npm": ">2.7.0"
+  },
+  "os": [
+    "aix",
+    "darwin",
+    "freebsd",
+    "linux",
+    "macos",
+    "openbsd",
+    "sunos",
+    "win32",
+    "windows"
+  ],
+  "keywords": [
+    "stdlib",
+    "stdmath",
+    "mathematics",
+    "math",
+    "blas",
+    "level 2",
+    "gsyr",
+    "linear",
+    "algebra",
+    "subroutines",
+    "array",
+    "ndarray"
+  ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_complex_access_pattern.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_complex_access_pattern.json
new file mode 100644
index 000000000000..7498487815a1
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_complex_access_pattern.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 3.0, 2.0, 1.0 ],
+  "strideX": -1,
+  "offsetX": 2,
+  "A": [ 1.0, 999.0, 0.0, 999.0, 0.0, 999.0, 2.0, 999.0, 1.0, 999.0, 0.0, 999.0, 3.0, 999.0, 2.0, 999.0, 1.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": 2,
+  "strideA2": 6,
+  "offsetA": 0,
+  "A_out": [ 2.0, 999.0, 0.0, 999.0, 0.0, 999.0, 4.0, 999.0, 5.0, 999.0, 0.0, 999.0, 6.0, 999.0, 8.0, 999.0, 10.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_l.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_l.json
new file mode 100644
index 000000000000..088fa13b1590
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_l.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "lower",
+  "N": 3,
+  "alpha": 2.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 1.0, 1.0, 1.0, 0.0, 2.0, 2.0, 0.0, 0.0, 3.0 ],
+  "A_mat": [
+      [ 1.0, 0.0, 0.0 ],
+      [ 1.0, 2.0, 0.0 ],
+      [ 1.0, 2.0, 3.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 1,
+  "strideA2": 3,
+  "offsetA": 0,
+  "A_out": [ 3.0, 5.0, 7.0, 0.0, 10.0, 14.0, 0.0, 0.0, 21.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_oa.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_oa.json
new file mode 100644
index 000000000000..60925d070e13
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_oa.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "lower",
+  "N": 3,
+  "alpha": 2.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 999.0, 1.0, 1.0, 1.0, 0.0, 2.0, 2.0, 0.0, 0.0, 3.0 ],
+  "A_mat": [
+      [ 1.0, 0.0, 0.0 ],
+      [ 1.0, 2.0, 0.0 ],
+      [ 1.0, 2.0, 3.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 1,
+  "strideA2": 3,
+  "offsetA": 1,
+  "A_out": [ 999.0, 3.0, 5.0, 7.0, 0.0, 10.0, 14.0, 0.0, 0.0, 21.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_ox.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_ox.json
new file mode 100644
index 000000000000..89045e3412da
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_ox.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "lower",
+  "N": 3,
+  "alpha": 2.0,
+  "x": [ 0.0, 0.0, 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 2,
+  "A": [ 1.0, 1.0, 1.0, 0.0, 2.0, 2.0, 0.0, 0.0, 3.0 ],
+  "A_mat": [
+      [ 1.0, 0.0, 0.0 ],
+      [ 1.0, 2.0, 0.0 ],
+      [ 1.0, 2.0, 3.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 1,
+  "strideA2": 3,
+  "offsetA": 0,
+  "A_out": [ 3.0, 5.0, 7.0, 0.0, 10.0, 14.0, 0.0, 0.0, 21.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2.json
new file mode 100644
index 000000000000..7ff0a090a92a
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 1.0, 999.0, 0.0, 999.0, 0.0, 999.0, 2.0, 999.0, 1.0, 999.0, 0.0, 999.0, 3.0, 999.0, 2.0, 999.0, 1.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": 2,
+  "strideA2": 6,
+  "offsetA": 0,
+  "A_out": [ 2.0, 999.0, 0.0, 999.0, 0.0, 999.0, 4.0, 999.0, 5.0, 999.0, 0.0, 999.0, 6.0, 999.0, 8.0, 999.0, 10.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2n.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2n.json
new file mode 100644
index 000000000000..5ba9c7a29796
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1_sa2n.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 3.0, 2.0, 1.0, 999.0, 999.0, 999.0, 2.0, 1.0, 0.0, 999.0, 999.0, 999.0, 1.0, 0.0, 0.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": 1,
+  "strideA2": -6,
+  "offsetA": 12,
+  "A_out": [ 6.0, 8.0, 10.0, 999.0, 999.0, 999.0, 4.0, 5.0, 0.0, 999.0, 999.0, 999.0, 2.0, 0.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2.json
new file mode 100644
index 000000000000..a9b50f734cc5
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 0.0, 0.0, 1.0, 999.0, 999.0, 999.0, 0.0, 1.0, 2.0, 999.0, 999.0, 999.0, 1.0, 2.0, 3.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": -1,
+  "strideA2": 6,
+  "offsetA": 2,
+  "A_out": [ 0.0, 0.0, 2.0, 999.0, 999.0, 999.0, 0.0, 5.0, 4.0, 999.0, 999.0, 999.0, 10.0, 8.0, 6.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2n.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2n.json
new file mode 100644
index 000000000000..fddbd49751a8
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_sa1n_sa2n.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 1.0, 2.0, 3.0, 0.0, 1.0, 2.0, 0.0, 0.0, 1.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": -1,
+  "strideA2": -3,
+  "offsetA": 8,
+  "A_out": [ 10.0, 8.0, 6.0, 0.0, 5.0, 4.0, 0.0, 0.0, 2.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_u.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_u.json
new file mode 100644
index 000000000000..c65f4f26fbfe
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_u.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 1.0, 0.0, 0.0, 2.0, 1.0, 0.0, 3.0, 2.0, 1.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 1,
+  "strideA2": 3,
+  "offsetA": 0,
+  "A_out": [ 2.0, 0.0, 0.0, 4.0, 5.0, 0.0, 6.0, 8.0, 10.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xn.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xn.json
new file mode 100644
index 000000000000..ba0f363e2fa4
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xn.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0, 4.0, 5.0 ],
+  "strideX": -2,
+  "offsetX": 4,
+  "A": [ 1.0, 0.0, 0.0, 2.0, 1.0, 0.0, 3.0, 2.0, 1.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 1,
+  "strideA2": 3,
+  "offsetA": 0,
+  "A_out": [ 26.0, 0.0, 0.0, 17.0, 10.0, 0.0, 8.0, 5.0, 2.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xp.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xp.json
new file mode 100644
index 000000000000..41ffcc1e291f
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/column_major_xp.json
@@ -0,0 +1,20 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0, 4.0, 5.0 ],
+  "strideX": 2,
+  "offsetX": 0,
+  "A": [ 1.0, 0.0, 0.0, 2.0, 1.0, 0.0, 3.0, 2.0, 1.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 1,
+  "strideA2": 3,
+  "offsetA": 0,
+  "A_out": [ 2.0, 0.0, 0.0, 5.0, 10.0, 0.0, 8.0, 17.0, 26.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_complex_access_pattern.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_complex_access_pattern.json
new file mode 100644
index 000000000000..43b7f6158961
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_complex_access_pattern.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 3.0, 2.0, 1.0 ],
+  "strideX": -1,
+  "offsetX": 2,
+  "A": [ 999.0, 1.0, 999.0, 2.0, 999.0, 3.0, 999.0, 0.0, 999.0, 1.0, 999.0, 2.0, 999.0, 0.0, 999.0, 0.0, 999.0, 1.0 ],
+  "A_mat": [
+    [ 1.0, 2.0, 3.0 ],
+    [ 0.0, 1.0, 2.0 ],
+    [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": 6,
+  "strideA2": 2,
+  "offsetA": 1,
+  "A_out": [ 999.0, 2.0, 999.0, 4.0, 999.0, 6.0, 999.0, 0.0, 999.0, 5.0, 999.0, 8.0, 999.0, 0.0, 999.0, 0.0, 999.0, 10.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_l.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_l.json
new file mode 100644
index 000000000000..1b51b23b043b
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_l.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "lower",
+  "N": 3,
+  "alpha": 2.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 1.0, 0.0, 0.0, 1.0, 2.0, 0.0, 1.0, 2.0, 3.0 ],
+  "A_mat": [
+    [ 1.0, 0.0, 0.0 ],
+    [ 1.0, 2.0, 0.0 ],
+    [ 1.0, 2.0, 3.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 3,
+  "strideA2": 1,
+  "offsetA": 0,
+  "A_out": [ 3.0, 0.0, 0.0, 5.0, 10.0, 0.0, 7.0, 14.0, 21.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_oa.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_oa.json
new file mode 100644
index 000000000000..f65891a426d4
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_oa.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "lower",
+  "N": 3,
+  "alpha": 2.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 999.0, 999.0, 999.0, 1.0, 0.0, 0.0, 1.0, 2.0, 0.0, 1.0, 2.0, 3.0 ],
+  "A_mat": [
+    [ 1.0, 0.0, 0.0 ],
+    [ 1.0, 2.0, 0.0 ],
+    [ 1.0, 2.0, 3.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 3,
+  "strideA2": 1,
+  "offsetA": 3,
+  "A_out": [ 999.0, 999.0, 999.0, 3.0, 0.0, 0.0, 5.0, 10.0, 0.0, 7.0, 14.0, 21.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_ox.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_ox.json
new file mode 100644
index 000000000000..94a3538145c5
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_ox.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "lower",
+  "N": 3,
+  "alpha": 2.0,
+  "x": [ 0.0, 0.0, 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 2,
+  "A": [ 1.0, 0.0, 0.0, 1.0, 2.0, 0.0, 1.0, 2.0, 3.0 ],
+  "A_mat": [
+    [ 1.0, 0.0, 0.0 ],
+    [ 1.0, 2.0, 0.0 ],
+    [ 1.0, 2.0, 3.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 3,
+  "strideA2": 1,
+  "offsetA": 0,
+  "A_out": [ 3.0, 0.0, 0.0, 5.0, 10.0, 0.0, 7.0, 14.0, 21.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2.json
new file mode 100644
index 000000000000..82ed634406f4
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 999.0, 1.0, 999.0, 2.0, 999.0, 3.0, 999.0, 0.0, 999.0, 1.0, 999.0, 2.0, 999.0, 0.0, 999.0, 0.0, 999.0, 1.0 ],
+  "A_mat": [
+    [ 1.0, 2.0, 3.0 ],
+    [ 0.0, 1.0, 2.0 ],
+    [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": 6,
+  "strideA2": 2,
+  "offsetA": 1,
+  "A_out": [ 999.0, 2.0, 999.0, 4.0, 999.0, 6.0, 999.0, 0.0, 999.0, 5.0, 999.0, 8.0, 999.0, 0.0, 999.0, 0.0, 999.0, 10.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2n.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2n.json
new file mode 100644
index 000000000000..7ef5794e0fc8
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1_sa2n.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 3.0, 2.0, 1.0, 999.0, 999.0, 999.0, 2.0, 1.0, 0.0, 999.0, 999.0, 999.0, 1.0, 0.0, 0.0 ],
+  "A_mat": [
+      [ 1.0, 2.0, 3.0 ],
+      [ 0.0, 1.0, 2.0 ],
+      [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": 6,
+  "strideA2": -1,
+  "offsetA": 2,
+  "A_out": [ 6.0, 4.0, 2.0, 999.0, 999.0, 999.0, 8.0, 5.0, 0.0, 999.0, 999.0, 999.0, 10.0, 0.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2.json
new file mode 100644
index 000000000000..d9f1f314215a
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 0.0, 0.0, 1.0, 999.0, 999.0, 999.0, 0.0, 1.0, 2.0, 999.0, 999.0, 999.0, 1.0, 2.0, 3.0 ],
+  "A_mat": [
+    [ 1.0, 2.0, 3.0 ],
+    [ 0.0, 1.0, 2.0 ],
+    [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 6,
+  "strideA1": -6,
+  "strideA2": 1,
+  "offsetA": 12,
+  "A_out": [  0.0, 0.0, 10.0, 999.0, 999.0, 999.0, 0.0, 5.0, 8.0, 999.0, 999.0, 999.0, 2.0, 4.0, 6.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2n.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2n.json
new file mode 100644
index 000000000000..48869bbcb9ff
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_sa1n_sa2n.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 1.0, 0.0, 0.0, 2.0, 1.0, 0.0, 3.0, 2.0, 1.0 ],
+  "A_mat": [
+    [ 1.0, 2.0, 3.0 ],
+    [ 0.0, 1.0, 2.0 ],
+    [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": -3,
+  "strideA2": -1,
+  "offsetA": 8,
+  "A_out": [ 10.0, 0.0, 0.0, 8.0, 5.0, 0.0, 6.0, 4.0, 2.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_u.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_u.json
new file mode 100644
index 000000000000..7162797cd9a8
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_u.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "strideX": 1,
+  "offsetX": 0,
+  "A": [ 1.0, 2.0, 3.0, 0.0, 1.0, 2.0, 0.0, 0.0, 1.0 ],
+  "A_mat": [
+    [ 1.0, 2.0, 3.0 ],
+    [ 0.0, 1.0, 2.0 ],
+    [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 3,
+  "strideA2": 1,
+  "offsetA": 0,
+  "A_out": [ 2.0, 4.0, 6.0, 0.0, 5.0, 8.0, 0.0, 0.0, 10.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xn.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xn.json
new file mode 100644
index 000000000000..ecd88396bc38
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xn.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0, 4.0, 5.0 ],
+  "strideX": -2,
+  "offsetX": 4,
+  "A": [ 1.0, 2.0, 3.0, 0.0, 1.0, 2.0, 0.0, 0.0, 1.0 ],
+  "A_mat": [
+    [ 1.0, 2.0, 3.0 ],
+    [ 0.0, 1.0, 2.0 ],
+    [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 3,
+  "strideA2": 1,
+  "offsetA": 0,
+  "A_out": [ 26.0, 17.0, 8.0, 0.0, 10.0, 5.0, 0.0, 0.0, 2.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xp.json b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xp.json
new file mode 100644
index 000000000000..67873efe2d97
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/fixtures/row_major_xp.json
@@ -0,0 +1,20 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "x": [ 1.0, 2.0, 3.0, 4.0, 5.0 ],
+  "strideX": 2,
+  "offsetX": 0,
+  "A": [ 1.0, 2.0, 3.0, 0.0, 1.0, 2.0, 0.0, 0.0, 1.0 ],
+  "A_mat": [
+    [ 1.0, 2.0, 3.0 ],
+    [ 0.0, 1.0, 2.0 ],
+    [ 0.0, 0.0, 1.0 ]
+  ],
+  "lda": 3,
+  "strideA1": 3,
+  "strideA2": 1,
+  "offsetA": 0,
+  "A_out": [ 2.0, 5.0, 8.0, 0.0, 10.0, 17.0, 0.0, 0.0, 26.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/test.js b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.js
new file mode 100644
index 000000000000..a1699b586a50
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.js
@@ -0,0 +1,38 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 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 tape = require( 'tape' );
+var gsyr = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof gsyr, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'attached to the main export is a method providing an ndarray interface', function test( t ) {
+	t.strictEqual( typeof gsyr.ndarray, 'function', 'method is a function' );
+	t.end();
+});
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js
new file mode 100644
index 000000000000..7afa6a7fa4a1
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js
@@ -0,0 +1,798 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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 max-len */
+
+'use strict';
+
+// MODULES //
+
+var tape = require( 'tape' );
+var toAccessorArray = require( '@stdlib/array/base/to-accessor-array' );
+var copy = require( '@stdlib/array/base/copy' );
+var gsyr = require( './../lib/main.js' );
+var ru = require( './fixtures/row_major_u.json' );
+var rl = require( './fixtures/row_major_l.json' );
+var rxp = require( './fixtures/row_major_xp.json' );
+var rxn = require( './fixtures/row_major_xn.json' );
+var cu = require( './fixtures/column_major_u.json' );
+var cl = require( './fixtures/column_major_l.json' );
+var cxp = require( './fixtures/column_major_xp.json' );
+var cxn = require( './fixtures/column_major_xn.json' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof gsyr, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function has an arity of 8', function test( t ) {
+	t.strictEqual( gsyr.length, 8, 'returns expected value' );
+	t.end();
+});
+
+tape( 'the function throws an error if provided an invalid first argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		'foo',
+		'bar',
+		'beep',
+		'boop'
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( value, data.uplo, data.N, data.alpha, data.x, data.strideX, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid first argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		'foo',
+		'bar',
+		'beep',
+		'boop'
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( value, data.uplo, data.N, data.alpha, data.x, data.strideX, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid second argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		'foo',
+		'bar',
+		'beep',
+		'boop'
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, value, data.N, data.alpha, data.x, data.strideX, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid second argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		'foo',
+		'bar',
+		'beep',
+		'boop'
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, value, data.N, data.alpha, data.x, data.strideX, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid third argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		-1,
+		-2,
+		-3
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, data.uplo, value, data.alpha, data.x, data.strideX, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid third argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		-1,
+		-2,
+		-3
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, data.uplo, value, data.alpha, data.x, data.strideX, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid sixth argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, data.uplo, data.N, data.alpha, data.x, value, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid sixth argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, data.uplo, data.N, data.alpha, data.x, value, data.A, data.lda );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid eighth argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		2,
+		1,
+		0,
+		-1,
+		-2,
+		-3
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, data.uplo, data.N, data.alpha, data.x, data.strideX, data.A, value );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid eighth argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		2,
+		1,
+		0,
+		-1,
+		-2,
+		-3
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			gsyr( data.order, data.uplo, data.N, data.alpha, data.x, data.strideX, data.A, value );
+		};
+	}
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, upper)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, upper) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, upper)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cu;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, upper) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cu;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, lower)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, lower) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, lower)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, lower) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function returns a reference to the input matrix `A`', function test( t ) {
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function returns a reference to the input matrix `A` (accessors)', function test( t ) {
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A;
+
+	out = gsyr( data.order, data.uplo, 0, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	out = gsyr( data.order, data.uplo, data.N, 0.0, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A;
+
+	out = gsyr( data.order, data.uplo, 0, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	out = gsyr( data.order, data.uplo, data.N, 0.0, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A;
+
+	out = gsyr( data.order, data.uplo, 0, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	out = gsyr( data.order, data.uplo, data.N, 0.0, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A;
+
+	out = gsyr( data.order, data.uplo, 0, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	out = gsyr( data.order, data.uplo, data.N, 0.0, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `x` stride (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxp;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `x` stride (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxp;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `x` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxp;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `x` stride (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxp;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxn;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxn;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxn;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxn;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = gsyr( data.order, data.uplo, data.N, data.alpha, x, data.strideX, a, data.lda );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js
new file mode 100644
index 000000000000..3c2f3478f813
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js
@@ -0,0 +1,1417 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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 max-len */
+
+'use strict';
+
+// MODULES //
+
+var tape = require( 'tape' );
+var toAccessorArray = require( '@stdlib/array/base/to-accessor-array' );
+var copy = require( '@stdlib/array/base/copy' );
+var dsyr = require( './../lib/ndarray.js' );
+
+
+// FIXTURES //
+
+var ru = require( './fixtures/row_major_u.json' );
+var rl = require( './fixtures/row_major_l.json' );
+var rxp = require( './fixtures/row_major_xp.json' );
+var rxn = require( './fixtures/row_major_xn.json' );
+var roa = require( './fixtures/row_major_oa.json' );
+var rox = require( './fixtures/row_major_ox.json' );
+var rsa1sa2 = require( './fixtures/row_major_sa1_sa2.json' );
+var rsa1nsa2 = require( './fixtures/row_major_sa1n_sa2.json' );
+var rsa1sa2n = require( './fixtures/row_major_sa1_sa2n.json' );
+var rsa1nsa2n = require( './fixtures/row_major_sa1n_sa2n.json' );
+var rcap = require( './fixtures/row_major_complex_access_pattern.json' );
+
+var cu = require( './fixtures/column_major_u.json' );
+var cl = require( './fixtures/column_major_l.json' );
+var cxp = require( './fixtures/column_major_xp.json' );
+var cxn = require( './fixtures/column_major_xn.json' );
+var coa = require( './fixtures/column_major_oa.json' );
+var cox = require( './fixtures/column_major_ox.json' );
+var csa1sa2 = require( './fixtures/column_major_sa1_sa2.json' );
+var csa1nsa2 = require( './fixtures/column_major_sa1n_sa2.json' );
+var csa1sa2n = require( './fixtures/column_major_sa1_sa2n.json' );
+var csa1nsa2n = require( './fixtures/column_major_sa1n_sa2n.json' );
+var ccap = require( './fixtures/column_major_complex_access_pattern.json' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof dsyr, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function has an arity of 10', function test( t ) {
+	t.strictEqual( dsyr.length, 10, 'returns expected value' );
+	t.end();
+});
+
+tape( 'the function throws an error if provided an invalid first argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		'foo',
+		'bar',
+		'beep',
+		'boop'
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( value, data.N, data.alpha, data.x, data.strideX, data.offsetX, data.A, data.strideA1, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid first argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		'foo',
+		'bar',
+		'beep',
+		'boop'
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( value, data.N, data.alpha, data.x, data.strideX, data.offsetX, data.A, data.strideA1, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid second argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		-1,
+		-2,
+		-3
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, value, data.alpha, data.x, data.strideX, data.offsetX, data.A, data.strideA1, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid second argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		-1,
+		-2,
+		-3
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, value, data.alpha, data.x, data.strideX, data.offsetX, data.A, data.strideA1, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid fifth argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, data.N, data.alpha, data.x, value, data.offsetX, data.A, data.strideA1, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid fifth argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, data.N, data.alpha, data.x, value, data.offsetX, data.A, data.strideA1, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid eighth argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, data.N, data.alpha, data.x, data.strideX, data.offsetX, data.A, value, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid eighth argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, data.N, data.alpha, data.x, data.strideX, data.offsetX, data.A, value, data.strideA2, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid ninth argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, data.N, data.alpha, data.x, data.strideX, data.offsetX, data.A, data.strideA1, value, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function throws an error if provided an invalid ninth argument (accessors)', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = ru;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dsyr( data.uplo, data.N, data.alpha, data.x, data.strideX, data.offsetX, data.A, data.strideA1, value, data.offsetA );
+		};
+	}
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, upper)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, upper) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, upper)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cu;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, upper) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cu;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, lower)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (row-major, lower) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, lower)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function performs the symmetric rank 1 operation `A = α*x*x^T + A` (column-major, lower) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function returns a reference to the input matrix `A`', function test( t ) {
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function returns a reference to the input matrix `A` (accessors)', function test( t ) {
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ru;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A;
+
+	out = dsyr( data.uplo, 0, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	out = dsyr( data.uplo, data.N, 0.0, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A;
+
+	out = dsyr( data.uplo, 0, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	out = dsyr( data.uplo, data.N, 0.0, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A;
+
+	out = dsyr( data.uplo, 0, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	out = dsyr( data.uplo, data.N, 0.0, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( a, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'if `N` is zero or the scalar constant is zero, the function returns the input matrix `A` unchanged (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cl;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A;
+
+	out = dsyr( data.uplo, 0, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	out = dsyr( data.uplo, data.N, 0.0, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying the strides of the first and second dimensions of `A` (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1sa2;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying the strides of the first and second dimensions of `A` (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1sa2;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying the strides of the first and second dimensions of `A` (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1sa2;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying the strides of the first and second dimensions of `A` (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1sa2;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the first dimension of `A` (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1nsa2;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the first dimension of `A` (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1nsa2;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the first dimension of `A` (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1nsa2;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the first dimension of `A` (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1nsa2;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the second dimension of `A` (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1sa2n;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the second dimension of `A` (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1sa2n;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the second dimension of `A` (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1sa2n;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports a negative stride for the second dimension of `A` (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1sa2n;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports negative strides for both dimensions of `A` (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1nsa2n;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports negative strides for both dimensions of `A` (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rsa1nsa2n;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports negative strides for both dimensions of `A` (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1nsa2n;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports negative strides for both dimensions of `A` (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = csa1nsa2n;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `A` offset (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = roa;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `A` offset (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = roa;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `A` offset (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = coa;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `A` offset (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = coa;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` stride (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxp;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` stride (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxp;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxp;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` stride (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxp;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxn;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rxn;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxn;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cxn;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` offset (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rox;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` offset (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rox;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` offset (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cox;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `x` offset (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = cox;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports complex access patterns (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rcap;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports complex access patterns (row-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = rcap;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports complex access patterns (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ccap;
+
+	a = copy( data.A );
+	x = copy( data.x );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports complex access patterns (column-major) (accessors)', function test( t ) {
+	var expected;
+	var abuf;
+	var data;
+	var out;
+	var a;
+	var x;
+
+	data = ccap;
+
+	abuf = copy( data.A );
+	a = toAccessorArray( abuf );
+	x = toAccessorArray( copy( data.x ) );
+
+	expected = data.A_out;
+
+	out = dsyr( data.uplo, data.N, data.alpha, x, data.strideX, data.offsetX, a, data.strideA1, data.strideA2, data.offsetA );
+	t.strictEqual( out, a, 'returns expected value' );
+	t.deepEqual( abuf, expected, 'returns expected value' );
+
+	t.end();
+});

From 44797299ca7ab1362dd142131c1babf9f5397ac0 Mon Sep 17 00:00:00 2001
From: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com>
Date: Mon, 28 Jul 2025 11:33:05 +0530
Subject: [PATCH 2/4] chore: remove Float64Array

Signed-off-by: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com>
---
 lib/node_modules/@stdlib/blas/base/gsyr/README.md | 10 ++--------
 1 file changed, 2 insertions(+), 8 deletions(-)

diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/README.md b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
index 6061f6ab252e..1a85f698fed2 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/README.md
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
@@ -35,8 +35,6 @@ var gsyr = require( '@stdlib/blas/base/gsyr' );
 Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
 
 ```javascript
-var Float64Array = require( '@stdlib/array/float64' );
-
 var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
 var x = [ 1.0, 2.0, 3.0 ];
 
@@ -62,7 +60,7 @@ var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
 var x = [ 3.0, 2.0, 1.0 ];
 
 gsyr( 'row-major', 'upper', 3, 1.0, x, -1, A, 3 );
-// A => <Float64Array>[ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+// A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
 ```
 
 Note that indexing is relative to the first index. To introduce an offset, use [`typed array`][mdn-typed-array] views.
@@ -88,13 +86,11 @@ gsyr( 'row-major', 'upper', 3, 1.0, x1, -1, A, 3 );
 Performs the symmetric rank 1 operation `A = α*x*x^T + A`, using alternative indexing semantics and where `α` is a scalar, `x` is an `N` element vector, and `A` is an `N` by `N` symmetric matrix.
 
 ```javascript
-var Float64Array = require( '@stdlib/array/float64' );
-
 var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
 var x = [ 1.0, 2.0, 3.0 ];
 
 gsyr.ndarray( 'upper', 3, 1.0, x, 1, 0, A, 3, 1, 0 );
-// A => <Float64Array>[ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
+// A => [ 2.0, 4.0, 6.0, 2.0, 5.0, 8.0, 3.0, 2.0, 10.0 ]
 ```
 
 The function has the following additional parameters:
@@ -107,8 +103,6 @@ The function has the following additional parameters:
 While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example,
 
 ```javascript
-var Float64Array = require( '@stdlib/array/float64' );
-
 var A = [ 1.0, 2.0, 3.0, 2.0, 1.0, 2.0, 3.0, 2.0, 1.0 ];
 var x = [ 1.0, 2.0, 3.0, 4.0, 5.0 ];
 

From 0147902f1bb691b7eece5915ea41ef42ae2b4a16 Mon Sep 17 00:00:00 2001
From: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com>
Date: Mon, 28 Jul 2025 11:38:20 +0530
Subject: [PATCH 3/4] chore: update dtype

Signed-off-by: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com>
---
 lib/node_modules/@stdlib/blas/base/gsyr/README.md | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/README.md b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
index 1a85f698fed2..d68e0a536eea 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/README.md
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
@@ -137,7 +137,7 @@ var ones = require( '@stdlib/array/ones' );
 var gsyr = require( '@stdlib/blas/base/gsyr' );
 
 var opts = {
-    'dtype': 'gsyr'
+    'dtype': 'generic'
 };
 
 var N = 3;

From 83509e3b0ce38c25c61248722c4ff4b51e434d36 Mon Sep 17 00:00:00 2001
From: stdlib-bot <82920195+stdlib-bot@users.noreply.github.com>
Date: Mon, 28 Jul 2025 06:11:45 +0000
Subject: [PATCH 4/4] chore: update copyright years

---
 lib/node_modules/@stdlib/blas/base/gsyr/README.md            | 2 +-
 lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts   | 2 +-
 lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js    | 2 +-
 lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js          | 2 +-
 lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js          | 2 +-
 lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js       | 2 +-
 lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js    | 2 +-
 lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js | 2 +-
 8 files changed, 8 insertions(+), 8 deletions(-)

diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/README.md b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
index d68e0a536eea..78300385b0d6 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/README.md
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/README.md
@@ -2,7 +2,7 @@
 
 @license Apache-2.0
 
-Copyright (c) 2024 The Stdlib Authors.
+Copyright (c) 2025 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.
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts b/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts
index 290aea1ab7f2..2b27831b15d1 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/docs/types/test.ts
@@ -1,7 +1,7 @@
 /*
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 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.
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js b/lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js
index 6c072a6e9b7e..277adc828c8d 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/examples/index.js
@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 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.
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js
index e84a1e9f8bdb..7c3a413eb4ee 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/base.js
@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 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.
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js
index b5620466d6ba..be34a72b1df1 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/main.js
@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 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.
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js
index e42a667f638c..fcbe90edddc7 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/lib/ndarray.js
@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 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.
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js
index 7afa6a7fa4a1..43246bffd174 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.main.js
@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 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.
diff --git a/lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js
index 3c2f3478f813..d0a36b482552 100644
--- a/lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js
+++ b/lib/node_modules/@stdlib/blas/base/gsyr/test/test.ndarray.js
@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 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.