blob: 4e005674f5315f262c33634348d07b1a00989932 (
plain)
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
|
#include "lib9.h"
#include "mathi.h"
double
dot(int n, double *x, double *y)
{
double sum = 0;
if (n <= 0)
return 0;
while (n--) {
sum += *x++ * *y++;
}
return sum;
}
int
iamax(int n, double *x)
{
int i, m;
double xm, a;
if (n <= 0)
return 0;
m = 0;
xm = fabs(*x);
for (i = 1; i < n; i++) {
a = fabs(*++x);
if (xm < a) {
m = i;
xm = a;
}
}
return m;
}
double
norm1(int n, double *x)
{
double sum = 0;
if (n <= 0)
return 0;
while (n--) {
sum += fabs(*x);
x++;
}
return sum;
}
double
norm2(int n, double *x)
{
double sum = 0;
if (n <= 0)
return 0;
while (n--) {
sum += *x * *x;
x++;
}
return sqrt(sum);
}
|
