example.c

/**
* @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.
*/

#include "stdlib/blas/base/drotm.h"
#include <stdio.h>

int main( void ) {
	// Create strided arrays:
	double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0 };
	double y[] = { 6.0, 7.0, 8.0, 9.0, 10.0 };

	// Specify the number of elements:
	const int N = 5;

	// Specify stride lengths:
	const int strideX = 1;
	const int strideY = 1;

	// Specify parameters for the modified Givens transformation:
	const double param[5] = { 0.0, 0.0, 2.0, -3.0, 0.0 };

	// Apply plane rotation:
	c_drotm( N, x, strideX, y, strideY, param );

	// Print the result:
	for ( int i = 0; i < 5; i++ ) {
		printf( "x[ %i ] = %lf, y[ %i ] = %lf\n", i, x[ i ], i, y[ i ] );
	}

	// Apply plane rotation:
	c_drotm_ndarray( N, x, -strideX, N-1, y, -strideY, N-1, param );

	// Print the result:
	for ( int i = 0; i < 5; i++ ) {
		printf( "x[ %i ] = %lf, y[ %i ] = %lf\n", i, x[ i ], i, y[ i ] );
	}
}