QR鍒嗚В瀹為獙鎶ュ憡.doc - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

骞惰�岃�＄畻璇剧▼鑰冩牳

 

瀹為獙鎶ュ憡

 

鑰冩牳棰樼洰锛?/p>

QR

鍒嗚В绠楁硶鐨勫苟琛屽疄鐜?/p>

 

骞惰�屽疄鐜版柟寮忥�?/p>

OpenMP 

 

1. 

闂�棰樻�傝堪

 

QR

鍒嗚В娉曟槸鐩�鍓嶆眰涓�鑸�鐭╅樀鍏ㄩ儴鐗瑰緛鍊肩殑鏈�鏈夋晥骞跺箍娉涘簲鐢ㄧ殑鏂规硶锛屼竴鑸�鐭╅樀鍏堢粡杩囨�ｄ氦

鐩镐技鍙樺寲鎴愪负

householder

鐭╅樀锛岀劧鍚庡啀搴旂敤

QR

鏂规硶姹傜壒寰佸�煎拰鐗瑰緛鍚戦噺銆傚畠鏄�灏嗙煩闃靛�?/p>

瑙ｆ垚涓�涓�姝ｈ�勬�ｄ氦鐭╅�?/p>

Q

涓庝笂涓夎�掑舰鐭╅�?/p>

R

锛屾墍浠ョО涓?/p>

QR

鍒嗚В娉曘�?/p>

 

2. 

涓茶�屼唬鐮佹弿杩�

 

涓茶�屼富瑕佷唬鐮佸�備笅锛?/p>

 

#include

 <stdio.h> 

#include

 <math.h> 

#include

 <stdlib.h> 

#include

 <time.h> 

#define

 

n

 100 

//maxN

 

void

 Matrix_print(

double

 

A

[

n

][

n

]) 

{ 

 

for

 (

int

 i=0;i<

n

;i++) 

 

{ 

 

 

for

 (

int

 j=0;j<

n

;j++) 

 

 

 

printf("%8.4lf\t",

A

[i][j]); 

 

 

printf("\n"); 

 

} 

} 

double

 Matrix_norm(

double

 

a

[

n

]) 

 

{ 

 

double

 d=0; 

 

for

 (

int

 i=0;i<

n

;i++) 

 

 

d+=

a

[i]*

a

[i]; 

 

return

 sqrt(d); 

} 

 

void

 Matrix_multiply(

double

 

A

[

n

][

n

],

double

 

B

[

n

][

n

],

double

 

C

[

n

][

n

]) 

{ 

 

 

for

 (

int

 i=0;i<

n

;i++) 

 

 

for

 (

int

 j=0;j<

n

;j++) 

 

 

{ 

 

 

 

C

[i][j]=0; 

 

 

 

for

 (

int

 t=0;t<

n

;t++) 

 

 

 

 

C

[i][j]+=

A

[i][t]*

B

[t][j]; 

 

 

}