fftsvd/SVector.c at master · jbardhan/fftsvd · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
#include "FFTSVD.h"

/* Constructors and Destructors */

SVector SVector_allocate(unsigned int length) {
   return (SVector)calloc(length, sizeof(sreal));
}

void SVector_free(SVector vector) {
   free(vector);
}

/* Initialization and Copying */

void SVector_copy(SVector vectordest, SVector vectorsrc, unsigned int length) {
   unsigned int i;

#pragma ivdep
   for (i = 0; i < length; i++)
      vectordest[i] = vectorsrc[i];
}

void SVector_copypiece(SVector vectordest, unsigned int destStart, SVector vectorsrc, unsigned int srcStart,unsigned int length) {
   unsigned int i;
#pragma ivdep
   for (i = 0; i < length; i++)
      vectordest[destStart + i] = vectorsrc[srcStart + i];
}

void SVector_copyscaledpiece(SVector vectordest, unsigned int destStart, SVector vectorsrc, unsigned int srcStart,unsigned int length, sreal scalefactor) {
   unsigned int i;
#pragma ivdep
   for (i = 0; i < length; i++)
      vectordest[destStart + i] = scalefactor * vectorsrc[srcStart + i];
}

void SVector_zero(SVector vector, unsigned int length) {
   unsigned int i;

#pragma ivdep
   for (i = 0; i < length; i++)
      vector[i] = (sreal)0;
}

/* Arithmetic Operations */

void SVector_scale(SVector vector, sreal scale, unsigned int length) {
   unsigned int i;

#pragma ivdep
   for (i = 0; i < length; i++)
      vector[i] *= scale;
}

void SVector_subtractscaledvector(SVector vector1, sreal scale, SVector vector2, unsigned int length) {
   unsigned int i;

#pragma ivdep
   for (i = 0; i < length; i++)
      vector1[i] -= scale * vector2[i];
}

void SVector_addscaledvector(SVector vector1, sreal scale, SVector vector2, unsigned int length) {
   unsigned int i;

#pragma ivdep
   for (i = 0; i < length; i++)
      vector1[i] += scale * vector2[i];
}

void SVector_addvector(SVector vector1, SVector vector2, unsigned int length) {
   unsigned int i;

#pragma ivdep
   for (i = 0; i < length; i++)
      vector1[i] += vector2[i];
}

void SVector_subtractvector(SVector vector1, SVector vector2, unsigned int length) {
   unsigned int i;

#pragma ivdep
   for (i = 0; i < length; i++)
      vector1[i] -= vector2[i];
}

/* Linear Algebra Operations */

sreal SVector_dot(SVector vector1, SVector vector2, unsigned int length) {
   unsigned int i;
   sreal dot = (sreal)0;

#pragma ivdep
   for (i = 0; i < length; i++)
      dot += vector1[i] * vector2[i];

   return dot;
}

sreal SVector_norm(SVector vector, unsigned int length) {
   unsigned int i;
   sreal norm = (sreal)0;

#pragma ivdep
   for (i = 0; i < length; i++)
      norm += vector[i] * vector[i];

   norm = sqrt(norm);

   return norm;
}

void SVector_normalize(SVector vector, unsigned int length) {
   SVector_scale(vector, (sreal)1 / SVector_norm(vector, length), length);
}