MTL vector rotation stretching function

// Solving by Robinkin

#include "mtl / mtl.h" #include "mtl / utils.h" #include "mtl / tinalg_vec.h"

/ * Eliminates the last element of vector y.

Sample Output

[1, 2, 3, 4, 5,] [2, 4, 8, 16, 32,] [2.1304, 14.2608, 8.36723, 16.4257, 32.3883,] [-0.679258, -1.35852, -1.72902, -1.48202, 0 ,]

* /

Using namespace MTL;

TYPEDEF EXTERNAL_VEC VEC;

INTMAIN () {// begin const INT n = 5; double dx [] = {1, 2, 3, 4, 5}; double dy [] = {2, 4, 8, 16, 32};

VEC X (DX, N), Y (DY, N); // end print_vector (x); print_vector (y); // begin double a = x [n-1]; double b = y [n-1] ;

Givens_Rotation Rot (A, B);

/ / For the 5 pairs (x, y) described above, the rotation is plugged, and the rotated reading is ArctG (A / B) stretched SQRT (A * a, b * b)

Rot.Apply (x, y); // end

Print_vector (x); print_vector (y); return 0;}