feat: add C implementation for blas/base/ssyr2

lib/node_modules/@stdlib/blas/base/ssyr2/README.md

@@ -193,21 +193,36 @@ console.log( A );
 ### Usage
 
 ```c
-TODO
+#include "stdlib/blas/base/ssyr2.h"
 ```
 
-#### TODO
+#### c_ssyr2( order, uplo, N, alpha, \*X, strideX, \*Y, strideY, \*A, LDA )
 
-TODO.
+Performs the symmetric rank 2 operation `A = α*x*y^T + α*y*x^T + A`.
 
 ```c
-TODO
+float A[] = { 1.0f, 0.0f, 0.0f, 2.0f, 1.0f, 0.0f, 3.0f, 2.0f, 1.0f };
+float x[] = { 1.0f, 2.0f, 3.0f };
+float y[] = { 1.0f, 2.0f, 3.0f };
+
+c_ssyr2( 102, 121, 3, 1.0f, x, 1, y, 1, A, 3 );
 ```
 
-TODO
+The function accepts the following arguments:
+
+- **order**: `[in] CBLAS_LAYOUT` storage layout.
+- **uplo**: `[in] CBLAS_UPLO` specifies whether the upper or lower triangular part of the symmetric matrix `A` should be referenced.
+- **N**: `[in] CBLAS_INT` number of elements along each dimension of `A`.
+- **alpha**: `[in] float` scalar.
+- **X**: `[in] float*` first input array.
+- **strideX**: `[in] CBLAS_INT` index increment for `X`.
+- **Y**: `[in] float*` second input array.
+- **strideY**: `[in] CBLAS_INT` index increment for `Y`.
+- **A**: `[inout] float*` input matrix.
+- **LDA**: `[in] CBLAS_INT` stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`).
 
 ```c
-TODO
+void c_ssyr2( const CBLAS_LAYOUT order, const CBLAS_UPLO uplo, const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY, float *A, const CBLAS_INT LDA )
 ```
 
 </section>
@@ -229,7 +244,26 @@ TODO
 ### Examples
 
 ```c
-TODO
+#include "stdlib/blas/base/ssyr2.h"
+#include <stdio.h>
+
+int main( void ) {
+ // Create strided arrays:
+ float A[] = { 1.0f, 0.0f, 0.0f, 2.0f, 1.0f, 0.0f, 3.0f, 2.0f, 1.0f };
+ float x[] = { 1.0f, 2.0f, 3.0f };
+ float y[] = { 1.0f, 2.0f, 3.0f };
+
+ // Specify the number of elements:
+ const int N = 3;
+
+ // Perform the symmetric rank 2 operation `A = α*x*y^T + α*y*x^T + A`:
+ c_ssyr2( 102, 121, N, 1.0f, x, 1, y, 1, A, N );
+
+ // Print the result:
+ for ( int i = 0; i < N*N; i++ ) {
+ printf( "A[ %i ] = %f\n", i, A[ i ] );
+ }
+}
 ```
 
 </section>

lib/node_modules/@stdlib/blas/base/ssyr2/benchmark/benchmark.native.js

@@ -0,0 +1,110 @@
+/**
+* @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 resolve = require( 'path' ).resolve;
+var bench = require( '@stdlib/bench' );
+var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
+var ones = require( '@stdlib/array/ones' );
+var pow = require( '@stdlib/math/base/special/pow' );
+var floor = require( '@stdlib/math/base/special/floor' );
+var tryRequire = require( '@stdlib/utils/try-require' );
+var pkg = require( './../package.json' ).name;
+
+
+// VARIABLES //
+
+var ssyr2 = tryRequire( resolve( __dirname, './../lib/ssyr2.native.js' ) );
+var opts = {
+'skip': ( ssyr2 instanceof Error )
+};
+var options = {
+'dtype': 'float32'
+};
+
+
+// FUNCTIONS //
+
+/**
+* Creates a benchmark function.
+*
+* @private
+* @param {PositiveInteger} N - number of elements along each dimension
+* @returns {Function} benchmark function
+*/
+function createBenchmark( N ) {
+var x = ones( N, options.dtype );
+var y = ones( N, options.dtype );
+var A = ones( N*N, options.dtype );
+return benchmark;
+
+/**
+* Benchmark function.
+*
+* @private
+* @param {Benchmark} b - benchmark instance
+*/
+function benchmark( b ) {
+var z;
+var i;
+
+b.tic();
+for ( i = 0; i < b.iterations; i++ ) {
+z = ssyr2( 'row-major', 'upper', N, 1.0, x, 1, y, 1, A, N );
+if ( isnanf( z[ i%z.length ] ) ) {
+b.fail( 'should not return NaN' );
+}
+}
+b.toc();
+if ( isnanf( z[ i%z.length ] ) ) {
+b.fail( 'should not return NaN' );
+}
+b.pass( 'benchmark finished' );
+b.end();
+}
+}
+
+
+// MAIN //
+
+/**
+* Main execution sequence.
+*
+* @private
+*/
+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), opts, f );
+}
+}
+
+main();

lib/node_modules/@stdlib/blas/base/ssyr2/benchmark/c/Makefile

@@ -0,0 +1,146 @@
+#/
+# @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.
+#/
+
+# VARIABLES #
+
+ifndef VERBOSE
+QUIET := @
+else
+QUIET :=
+endif
+
+# Determine the OS ([1][1], [2][2]).
+#
+# [1]: https://en.wikipedia.org/wiki/Uname#Examples
+# [2]: http://stackoverflow.com/a/27776822/2225624
+OS ?= $(shell uname)
+ifneq (, $(findstring MINGW,$(OS)))
+OS := WINNT
+else
+ifneq (, $(findstring MSYS,$(OS)))
+OS := WINNT
+else
+ifneq (, $(findstring CYGWIN,$(OS)))
+OS := WINNT
+else
+ifneq (, $(findstring Windows_NT,$(OS)))
+OS := WINNT
+endif
+endif
+endif
+endif
+
+# Define the program used for compiling C source files:
+ifdef C_COMPILER
+CC := $(C_COMPILER)
+else
+CC := gcc
+endif
+
+# Define the command-line options when compiling C files:
+CFLAGS ?= \
+-std=c99 \
+-O3 \
+-Wall \
+-pedantic
+
+# Determine whether to generate position independent code ([1][1], [2][2]).
+#
+# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options
+# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option
+ifeq ($(OS), WINNT)
+fPIC ?=
+else
+fPIC ?= -fPIC
+endif
+
+# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`):
+INCLUDE ?=
+
+# List of source files:
+SOURCE_FILES ?=
+
+# List of libraries (e.g., `-lopenblas -lpthread`):
+LIBRARIES ?=
+
+# List of library paths (e.g., `-L /foo/bar -L /beep/boop`):
+LIBPATH ?=
+
+# List of C targets:
+c_targets := benchmark.length.out
+
+
+# RULES #
+
+#/
+# Compiles source files.
+#
+# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`)
+# @param {string} [CFLAGS] - C compiler options
+# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`)
+# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`)
+# @param {string} [SOURCE_FILES] - list of source files
+# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`)
+# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`)
+#
+# @example
+# make
+#
+# @example
+# make all
+#/
+all: $(c_targets)
+
+.PHONY: all
+
+#/
+# Compiles C source files.
+#
+# @private
+# @param {string} CC - C compiler (e.g., `gcc`)
+# @param {string} CFLAGS - C compiler options
+# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`)
+# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`)
+# @param {string} SOURCE_FILES - list of source files
+# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`)
+# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`)
+#/
+$(c_targets): %.out: %.c
+$(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES)
+
+#/
+# Runs compiled benchmarks.
+#
+# @example
+# make run
+#/
+run: $(c_targets)
+$(QUIET) ./$<
+
+.PHONY: run
+
+#/
+# Removes generated files.
+#
+# @example
+# make clean
+#/
+clean:
+$(QUIET) -rm -f *.o *.out
+
+.PHONY: clean