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
|
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#define asize1 1000
double rnd(double width);
int main(int argc, char *argv[])
{
int i,j,k;
double sum;
double **a;
double **b;
double **c;
a = (double**)malloc(sizeof(double*)*asize1);
for (i=0; i<asize1; i++){
a[i] = (double*)malloc(sizeof(double)*asize1);
}
b = (double**)malloc(sizeof(double*)*asize1);
for (i=0; i<asize1; i++){
b[i] = (double*)malloc(sizeof(double)*asize1);
}
c = (double**)malloc(sizeof(double*)*asize1);
for (i=0; i<asize1; i++){
c[i] = (double*)malloc(sizeof(double)*asize1);
}
for (i=0; i<asize1; i++){
for (j=0; j<asize1; j++){
a[i][j] = 1.0e-10*(double)(i + j);
b[i][j] = 1.0e-10*(double)(i - j);
}
}
for (i=0; i<asize1; i++){
for (j=0; j<asize1; j++){
sum = 0.0;
for (k=0; k<asize1; k++){
sum += a[i][k]*b[j][k];
}
c[i][j] = sum;
}
}
}
double rnd(double width)
{
/****************************************************
This rnd() function generates random number
-width/2 to width/2
****************************************************/
double result;
result = rand();
while (width<result){
result = result/2.0;
}
result = result - width*0.75;
return result;
}
|
